Diﬀusion of the carbon dioxide–ethanol mixture in the extended critical region †

The eﬀect of traces of ethanol in supercritical carbon dioxide on the mixture’s thermodynamic properties is studied by molecular simulations and Taylor dispersion measurements. This mixture is investigated along the isobar p = 10 MPa in the temperature range between T = 304 and 343 K. Along this path, the mixture undergoes two transitions: First, the Widom line is crossed, marking the transition from liquid-like to gas-like conditions. A second transition occurs from the supercritical gas-like domain to a subcritical gas. The Widom line crossover entails inflection points for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell–Stefan and intradiffusion coefficients. On the other hand, the transition between the super- and subcritical regions is found to be generally smooth, an observation that is qualitatively confirmed by experimental Taylor dispersion measurements. Dedicated atomistic simulations show the presence of microheterogeneities due to ethanol self-association along the investigated path, which lead to the mixture’s anomalous behavior in its extended critical region. schematic of the employed Taylor dispersion apparatus, graphical results for the enthalpy h , specific volume v , c p , c v , a v and b T , as well as the center-of-mass radial distribution function and average coordination number of the CO 2 –ethanol pair. Predictions of the Fick diffusion coefficient of the CO 2 + ethanol mixture with the regarded equations are in comparison with experimental and simulation results. equilibrium molecular dynamics simulation also


Introduction
Investigating the thermophysical properties of supercritical carbon dioxide (scCO 2 ) has grown to a topic of great interest. 1-4 While the majority of work had focused on pure fluids in sufficient distance from the critical point, 5,6 only some special case scenarios with respect to a second component have been examined under near-critical conditions. [7][8][9] CO 2 is commonly used in various technological applications, ranging from environmental, mechanical, chemical, geothermal to pharmaceutical industries. 10 It constitutes an attractive alternative to organic or aqueous solvents by having a remarkably mild critical point, being non-toxic, non-flammable, widely available and largely inert. However, due to its low dielectric constant in combination with a zero dipole moment, 11 it is a poor solvent for polar substances. Yet, the solubility of high molecular weight solutes is readily enhanced by adding small amounts of strongly polar entrainers. 12,13 Consequently, ethanol emerges as a widespread commodity in supercritical extraction processes. Further, ethanol forms a non-ideal mixture with CO 2 due to the presence of strong solute-solute interactions and hydrogen bonding making this mixture particularly interesting.
In contrast to the classical perception of the supercritical region as a featureless domain, experimental evidence [14][15][16][17] and molecular dynamics simulations 3,4,9,[18][19][20] indicate a partition into several subdomains. The continuous dynamic crossover between different domains occurs across transition lines, i.e. Fisher-Widom, 21 Nishikawa,14,16 Frenkel 22,23 and Widom. 24 For moderate pressures, i.e. up to three times the critical pressure 25 p r 3p c , the crossover between gas-like and liquidlike regimes is accomplished over a delta shaped area 20,26,27 that is delimited by the loci of extrema of particular thermodynamic response functions, i.e. isobaric and isochoric heat capacities c p , c v , thermal expansion a v , isothermal compressibility b T , density r and speed of sound c. These response functions eventually confluence to a single line in close vicinity to the critical point, the so-called Widom line, which can be considered as an extension of the vapor pressure curve. Crossing the Widom line is additionally associated with a minimum of the thermodynamic factor and with large density fluctuations. 28 While the critical point exhibits critical opalescence, 29 its extended vicinity is characterized by emerging microscopic clusters that have been observed experimentally 28 and with equilibrium molecular dynamics simulations. [30][31][32][33][34] The present work contributes to the understanding of how ethanol diluted into CO 2 affects the mixture's microscopic structure and a Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany. E-mail: vrabec@tu-berlin.de; Tel: +49 30 314 22755 b Microgravity Research Center, Université Libre de Bruxelles, 1050 Bruxelles, Belgium † Electronic supplementary information (ESI) available: It includes a schematic of the employed Taylor dispersion apparatus, graphical results for the enthalpy h, specific volume v, c p , c v , a v and b T , as well as the center-of-mass radial distribution function and average coordination number of the CO 2 -ethanol pair. Predictions of the Fick diffusion coefficient of the CO 2 + ethanol mixture with the regarded equations are shown in comparison with experimental and simulation results. The numerical data from equilibrium molecular dynamics simulation are also listed. See DOI: 10.1039/d0cp04985a consequently its transport properties in the extended critical region. The CO 2 + ethanol mixture was examined twofold. While the microscopic structure was investigated exclusively by equilibrium molecular dynamics (EMD) simulations, mutual diffusion was sampled atomistically and was measured in the lab with the Taylor dispersion technique. The present mixture has already been investigated experimentally [35][36][37] and by molecular simulations. In particular, the vapor-liquid equilibrium (VLE), 38 solubility parameters, 39 microscopic structure 31,32,[40][41][42] and diffusion coefficients 43 have been studied by molecular simulation. However, the effects of structural changes on mutual diffusion across the Widom line in the extended critical region of this mixture have not yet been considered neither experimentally, nor by molecular simulation. This is of special interest for the understanding and description of supercritical extraction processes which are performed in this region.

Molecular simulation
Classical molecular simulations are based on molecular force fields models that are typically optimized to describe the pure component phase behavior. The transferability to a given mixture is validated by comparing simulation predictions to experimental data along the VLE curves. While the employed force fields were based on the Lennard-Jones potential for both components, a predictive mode was used to describe the interactions between unlike Lennard-Jones sites, i.e. the Lorentz-Berthelot combining rules were assumed. The force field for ethanol consisted of three Lennard-Jones sites with three superimposed point charges. 44 In contrast, the two Lennard-Jones sites of carbon dioxide's force field included a superimposed point quadrupole. 45 States along the VLE curves were sampled by Monte Carlo (MC) simulations with the grand equilibrium method containing N = 1372 molecules for the liquid phase and N = 460 molecules for the gas phase simulations.
The MC simulation results for the VLE curves are in excellent agreement with experimental literature data 46-52 over the entire investigated temperature and pressure ranges, cf. Fig. 1. In order to properly appreciate the simulation results, the coexistence curves were additionally computed by three equations of state, allowing for a temperature-dependent binary parameter k ij , cf. The intra-D i and Maxwell-Stefan diffusion coefficients as well as the shear viscosity were sampled directly with EMD simulation and the Green-Kubo formalism. 57,58 The working equations for the determination of these transport properties have been published e.g. in ref. 59 and are not repeated here. EMD simulations for transport properties were made in two steps. In the first step, a simulation in the isobaric-isothermal (NpT) ensemble was performed to calculate the density at the desired temperature and pressure. In the second step, EMD simulations were performed at this temperature and density in the canonical (NVT) ensemble containing N = 3000 molecules in cubic volumes with periodic boundary conditions and specifying an integrator time step of Dt = 1 fs. The associated finite size effects were corrected with a modified Yeh-Hummer approach 60,61 employing the sampled shear viscosity values. To improve statistics, a total of 2.5 Â 10 5 correlation functions was averaged. The thermodynamic factor G was sampled directly with Kirkwood-Buff integration based on the methodology proposed by Ganguly and van der Vegt, 62 which was found to be the most adequate in previous work. 63 Extrapolation to the thermodynamic limit was not necessary. All molecular simulations were performed with the fully open source software ms2. 64 Statistical uncertainties were calculated with the error propagation law and a coverage factor of k = 1. Throughout, the cut-off radius was set to r c = 17.5 Å. Electrostatic long-range corrections were made using the

Taylor dispersion technique
The employed Taylor apparatus has been optimized for the use of scCO 2 as described in preceding works. 35,65,66 The experimental determination of the mixture's mutual diffusion coefficient D constitutes a four stage procedure utilizing the Taylor dispersion technique. The experimental schematic is disclosed in the ESI. † During the first step, the pure CO 2 carrier fluid which was stored under VLE conditions, i.e. at T = 288.15 K, p B 5 MPa with a purity 0.99998 mol mol À1 (Air Liquide), was initially liquefied through a cryostat reducing its temperature to T = 269.15 K. Liquid CO 2 was subsequently pumped above the critical pressure and eventually heated to its target temperature before it reached the injection valve (Knauer model D-14163). In the second step, scCO 2 was delivered to the dispersion tube with a constant flow rate. The carrier stream was thermostated by means of a heat reservoir ensuring a constant target temperature ranging from T exp = 304 to 343 K with an accuracy of AE0.1 K and barostated with a high pressure pump at p exp = 10 MPa with AE0.05 MPa accuracy. The ethanol sample (purity 99.9% (GC) in volume fraction, CAS 64-17-5; purchased from VWR) adopts to the respective target temperature T exp within the valve's loop prior to its injection into the scCO 2 stream. In order to ensure temperature homogenization in this section of the experiment, the dispersion tube and injection valve were placed inside a polyurethane foam insulated housing with an additional air fan. In the third step, the sample was injected into the scCO 2 stream. The resulting strongly diluted mixture was fed to a l = 30.916 AE 0.001 m long dispersion capillary with a circular cross section of radius r = 0.375 mm. The capillary was coiled around a grooved, hollow aluminum cylinder with radius R c = 0.175 m providing stability and fixation. The cylinder was additionally thermostated with an internal circular flow ensuring a temperature stabilization of AE0.1 K. In order to minimize pressure and density disturbances during injection, an ethanol sample volume of V 0 = 2 Â 10 À6 dm 3 was selected, with smaller volumes having a negative effect on the signal-to-noise ratio. 35,65 The pressure of the system was controlled by a back pressure regulator (Jasco BP-2080) and measured by means of a pressure sensor (JUMO dTrans p30) with an accuracy of AE0.05 MPa. In the experiment's final step, the Taylor peak was monitored at the outlet of the dispersion tube by means of a FT-IR spectrometer (Jasco FT-IR 4100) with AE 0.01 cm À1 accuracy and 4 cm À1 resolution. In contrast to its nominal operation, the employed FT-IR was equipped with a custom-built high pressure demountable cell (Harrick) that was optimized for the best possible signal-to-noise ratio. 35 The ZnSe cell had a thickness 150 mm 35 and allowed for a maximum working pressure of 25 MPa. The data generated by the FTIR were digitally read out by a specific software (Spectra Manager by Jasco) and the variation of the solute concentration over time was monitored through the absorbance spectra at wavenumbers corresponding to different vibration modes. The procedure to select the working wavenumbers, the experimental protocol and the fitting procedure have been reported in ref. 65 and 66.

Critical line and Widom line
The CO 2 + ethanol mixture was investigated by EMD simulation along the isobar p = 10 MPa in the temperature range between T = 305 and 340 K with a composition x CO2 = 0.97 mol mol À1 , cf. Fig. 2. This composition was chosen as the closest to the infinite dilution limit that allowed for adequate statistics. Along this path, the mixture undergoes two transitions, indicated as I and II in Fig. 2. First, the Widom line is crossed at T B 323 K (point I) marking the transition from liquid-like to gas-like conditions, as indicated by the maxima of the response functions c p , a v , b T , cf. Fig. S2-S4 in the ESI. † Each function's inflection point corresponds to a maximum of a thermodynamic response function and consequently determines the mixture's Widom line. The inflection points of the enthalpy and density, which correspond to maxima of the isobaric heat capacity c p and the thermal expansion a v , respectively, occur according to the employed force field at T = 317 K for pure CO 2 and T = 323 K for the mixture, cf. Fig. 3. Thus, the addition of a small amount of ethanol shifts the Widom line up by B 6 K. To assess the capability of the employed force field to predict the studied properties, the temperature dependence of density and enthalpy of pure CO 2 was compared to the Span-Wagner equation of state, 67 which is of reference quality. In general, the employed CO 2 force field is able to predict the density of CO 2 in the studied temperature and pressure ranges with a good accuracy. The average deviation between simulation results and the Span-Wagner equation of state 67 is 5.4%.  Table 2 (black crosses). The red bullet symbolizes the molar composition x CO2 = 0.97 mol mol À1 and pressure p = 10 MPa at which all molecular dynamics simulations were carried out.
However, the predicted inflection points, which define the Widom line, occur approximately 2 K below those according to the Span-Wagner equation of state. 67 Therefore, the related uncertainty of the present simulation results is expected to have a similar magnitude, which is marked as a shaded region around the calculated transition points I and II. After crossing the Widom line, the mixture undergoes a transition from the supercritical gas-like domain to a gas at subcritical conditions at T B 328 K (point II), as indicated by the critical line of the mixture, cf. Table 2. More specifically, while for T = 318 K all states above p = 8.64 MPa and x CO2 = 0.938 mol mol À1 are located in the supercritical region, for T = 328 K only states above p = 10.09 MPa and x CO2 = 0.863 mol mol À1 lie in the supercritical region. This can be clearly seen in the inset of Fig. 2, while for temperatures below T B 328 K the simulated pressure-composition pair is clearly located above the critical line, for higher temperatures, the simulated pressure is below that of the critical line. The pressure-temperature projection of the vapor pressure curves for the CO 2 + ethanol mixture is shown in Fig. 2. The critical line of the mixture was determined by joining the mixture's critical points from experimental data. [67][68][69][70] A distinctive feature of the curve is that it goes through a pronounced pressure maximum at p B 15.5 MPa despite the relatively small difference between the critical pressures of pure CO 2 and pure ethanol, p c = 7.4 and 6.3 MPa, respectively. The large temperature range of that curve (304 K to 515 K) is simply due to the difference between the critical temperatures of the pure substances. Note that the composition varies along the critical line such that the pressure-composition pair that is focussed on does not cross the two-phase region.

Intradiffusion coefficients
The dynamical crossover between the supercritical high density and the low density regimes across the Widom line has been related to the presence of inflection points or extrema of some transport coefficients. 3 Therefore, the analysis of the diffusion behavior of the mixture along the regarded isobar offers insight into the underlying transition dynamics. The temperature dependence of the intradiffusion coefficient of CO 2 is sigmoidal, with an initially rather weak temperature dependence followed by a strong stepwise increase at temperatures between T = 320 and 330 K, cf. Fig. 4. This curve is similar to that observed for the enthalpy and shows an inflection point at the Widom line, i.e. at T B 323 K. The random propagation of molecules is not only dependent on the thermodynamic state point, but is also strongly affected by other factors like molecular size and polarity. Thus, within a mixture, the components generally propagate with different velocity distributions leading to different intradiffusion coefficients. In fact, the intradiffusion coefficient of the bulkier ethanol molecules is on average 40% lower than that of the smaller CO 2 molecules. Although both curves show a sigmoidal behavior, the stepwise increase of the intradiffusion coefficient of ethanol is less pronounced than the one of CO 2 , which can be linked to the presence of local density inhomogeneities caused by self-association. Further, the inflection point   of the curve is observed at the transition between the super-and subcritical regimes T B 328 K.

Microscopic structure
To elucidate the relevant microscopic structure aspects, the center-of-mass radial distribution functions of the ethanolethanol g EtOH-EtOH (r), CO 2 -ethanol g CO2-EtOH (r) and CO 2 -CO 2 g CO2-CO2 (r) pair interactions were analyzed, cf. Fig. 5 and Fig. S5 of the ESI. † The main peak of g EtOH-EtOH (r) with a maximum located at r B 4.3 Å is related to ethanol self-association through hydrogen bonding. This first peak is followed by a well defined shoulder, which is more pronounced at the lowest studied temperature and indicates an overlap with the second solvation shell. The location of the main peak does not change significantly along the studied isobar, suggesting that ethanol remains structured at short intermolecular distances. On the other hand, g CO2-CO2 (r) shows a well defined first peak followed by a considerably smaller second peak, which becomes weaker with increasing temperature and disappears at T B 330 K. This reduction of the long-range structure is expected because of the transition from supercritical gas-like to compressed gas conditions. The observed differences in peak intensity and the relatively small main peaks observed in the radial distribution functions at lower temperatures are mainly a result from statistical standardization, i.e. more ethanol or CO 2 molecules can be found in the far range of the simulation volume. 71 The CO 2 -ethanol pair interaction, shown in Fig. S5 in the ESI, † exhibits a relatively small peak with a shoulder located between r B 3.0 and 6.3 Å, which partially lies within the range of hydrogen bonding interactions and suggests the occurrence of relatively strong CO 2 -ethanol association. 32 However, because of the rather small peak magnitude, these interactions are expected to occur rather sporadically. A quantification of the local inhomogeneities can be achieved by monitoring the average coordination numbers for the first coordination shell defined as N xÀy ¼ 4pr y Ð rc 0 r 2 g xÀy ðrÞdr. Therein, x stands for the central interaction site surrounded by interaction sites of type y, r y is the bulk density of interaction sites of type y, g xÀy (r) represents the radial distribution function of the pair involved in the running average number calculation and r c is the radius of the first coordination shell, i.e. the location of the first minimum of the regarded radial distribution function g xÀy (r). The average coordination number N CO2-CO2 shows a similar behavior as the bulk density. At temperatures between T = 315 and 330 K, the coordination number decreases rapidly to approximately one third of its initial value and remains more stable in the compressed gas region, cf. Fig. 6. The strong decrease of the coordination number is mirrored in the observed stepwise increase of the CO 2 intradiffusion coefficient. Similarly, the inflection point of the coordination number was also found to be located at the Widom line, i.e. at T B 323 K. The average coordination number N EtOH-EtOH , indicating the amount of ethanol self-association, shows a stronger temperature dependence in the supercritical liquid-like region. At T = 305 K, each ethanol molecule is associated on average with 1.4 alcohol molecules, but after an increase of 10 K in temperature this value is reduced to 1.23. This implies a reduction of about 12% in the hydrogen bonded structures and explains the related increase of the ethanol intradiffusion coefficient between  View Article Online T = 305 and 315 K. In the region between 315 and 328 K, the decrease of the average ethanol-ethanol coordination number is much less pronounced than that of CO 2 , suggesting the presence of enhanced ethanol hydrogen bonded structures that are disrupted to a lesser extent by the strong density reduction. These rather stable ethanol hydrogen bonded structures are linked to mixture microheterogeneities observed in the domain influenced by the Widom line and the super-to subcritical transition. The relatively small values of the thermodynamic factor as well as the milder sigmoidal increase of the ethanol intradiffusion coefficient can also be explained with these microscopic structures. At temperatures above T B 328 K, the coordination number decreases almost stepwise, suggesting a steady breakup of the hydrogen bonded structures with temperature in the compressed gas region. The presence of ethanol self-association and microheterogeneities can be visually corroborated when simulation snapshots are analyzed. Although a snapshot represents only one microstate of a molecular system, in case of mixtures with associating components, a single microstate may well represent all possible microstates, since they are permutations of the segregation patterns. 72 Fig. 7 shows snapshots of simulation volumes for four temperatures. In spite of the low alcohol content, ethanol molecules tend to self-associate and form hydrogen bonded networks. At low temperatures up to the Widom line, most of the ethanol molecules are part of clusters that form segregated domains. As the temperature is increased, the ethanol clusters become smaller and are more uniformly distributed in the simulation volume. This observation corresponds to the observed values of the intradiffusion coefficient of ethanol and the relatively low values of the thermodynamic factor.

Mutual diffusion coefficients
The Maxwell-Stefan diffusion coefficient calculated from the Onsager phenomenological coefficients 73 of the studied mixture is shown in Fig. 4. As can be seen for temperatures below T = 328 K, the Maxwell-Stefan diffusivity curve runs almost parallel to and above the one of CO 2 intradiffusion. Higher values of the Maxwell-Stefan diffusion coefficient with respect to the intradiffusion coefficient are related to cluster formation due to alcohol self-association. 74 In the regime following the transition to the compressed gas, the Maxwell-Stefan diffusion coefficient becomes lower than the CO 2 intradiffusion coefficient, which is in line with the strong reduction of the ethanol hydrogen bonded structures indicated by the decrease of the average ethanol-ethanol coordination number. The inflection point of the Maxwell-Stefan diffusivity curve is located in the vicinity of the Widom line, i.e. at T B 324 K.
The Maxwell-Stefan diffusion coefficient obtained from EMD simulation can be straightforwardly related to the Fick diffusion coefficient D through the so-called thermodynamic factor G, i.e. D = G , which is a measure of the mixture's nonideality. Because of the nature of the CO 2 + ethanol mixture and the presence of microheterogeneities, the sampled thermodynamic factor reaches relatively low values G B 0.45. The expected minimum of the thermodynamic factor in the proximity of the Widom line is predicted clearly by the employed equations of state. The Kirkwood-Buff integration results predict a weaker minimum shifted to higher temperatures, cf. Fig. 8. For the sake of consistency and because of the rather large differences in G obtained from the different equations of state, the thermodynamic factor from Kirkwood-Buff integration was employed here to calculate the Fick diffusion coefficient. Since the values of the sampled thermodynamic factor do Fig. 7 Snapshots of the present EMD simulations at selected temperatures. The CO 2 solvent molecules were graphically removed to reveal clustering among ethanol molecules. not show a strong variation in the range of the studied thermodynamic conditions, the temperature dependence of the sampled Fick diffusion coefficient is similar to that observed for the Maxwell-Stefan diffusion coefficient, cf. Fig. 9. However, its inflection point is shifted to a higher temperature, i.e. T B 329 K, which can be explained with the changes observed for the ethanol-ethanol average coordination number near this temperature. In general, values of the Fick diffusion coefficient sampled by EMD are comparable with present measurements within the statistical uncertainties, cf. Table 3. However, the experimental values exhibit a smoother temperature dependence with an inflection point located at T B 327 K. The comparably lower values of the predicted Fick diffusion coefficient in the supercritical liquid-like region, i.e. at T o 320 K, could be due to an overestimation of ethanol selfassociation by molecular simulation. Approaching the Widom line, where density variations are large, simulation and experimental results agree quite well up to T B 332 K, i.e. right after the transition into the subcritical region. At the highest studied temperatures, the Fick diffusion coefficient predicted by EMD is overestimated when compared with the experimental data. This suggests that the actual breakup of ethanol hydrogen bonded structures in the compressed gas region occurs more gradually than predicted by EMD simulation. It should be noted that contrary to EMD simulations, where the mixture concentration is exactly known, the experimental concentration can only be estimated. Therefore, the differences between simulation and experimental values can partially be explained with the difference in ethanol concentration, which has been proven to have a strong influence on the Fick diffusion coefficient even in the dilution limit of scCO 2 mixtures. 66,75 This significant influence of concentration on the Fick diffusion coefficient is mainly related to the strong concentration dependence of the thermodynamic factor given by the proximity of the critical point. 66 A direct comparison of present experimental results with data from previous studies 36,37 is not possible because of they were measured for different state points.
The Stokes-Einstein based equations by Wilke-Chang, 76 Sassiat, 77 Tyn-Calus, 78 Scheibel, 76 Reddy-Doraiswamy 79 and Lai-Tan 80 as well as the free volume based equations by Catchpole and King, 81 He and Yu 82 as well as Funazukuri et al. 83 were tested with respect to their ability to predict the Fick diffusion coefficient at infinite dilution of the CO 2 + ethanol mixture. However, none of these equations was able to capture neither qualitatively nor quantitatively the present experimental Fick diffusion coefficient values, cf. Fig. S7 of the ESI. † The predictive He and Yu 82 equation yields the best results among the tested equations, however, its agreement with the present experimental results is quite poor in the supercritical region, cf. Fig. 9(a). The strong overestimation of the Fick diffusion coefficient in regions with a higher density, where ethanol selfassociation plays a decisive role for molecular mobility, as well as the better agreement with experiments in the compressed gas region, indicate a failure of these models to adequately consider the hydrogen bonded structures in the mixture.

Shear viscosity
The importance of considering information on hydrogen bonding dynamics for the prediction of transport properties can be clearly  observed when the shear viscosity is analyzed. Fig. 9(b) shows the strong difference between the shear viscosity of pure CO 2 and that of the mixture. In the liquid-like region, where intermolecular forces dominate the viscous effects, the shear viscosity of the mixture exhibits an increment of approximately 70% from the value of pure CO 2 . This enhancement can be ascribed to the presence of microheterogeneities caused by hydrogen bonding networks despite the small amount of ethanol present in the mixture. Here, ethanol self-association is expected to play the main role in the formation of microscopic structures because of the low values of the CO 2 -ethanol coordination number, i.e. N CO 2 -EtOH o 0.3. As higher temperatures are reached, the momentum transfer due to molecular thermal motion makes an increasingly important contribution to the shear viscosity and hydrogen bonds break, reducing the viscosity enhancement to less than 15% in the compressed gas region. Further, the shear viscosity curve of the mixture shows an inflection point located at T B 320 K, which cannot be observed for pure CO 2 . Viscosity and diffusion curves exhibit an opposing behavior, as expected from the Stokes-Einstein equation.

Conclusions
A study on the dynamic behavior of ethanol diluted in -CO 2 in the temperature range from T = 304 to 343 K along the isobar p = 10 MPa was conducted employing complementary approaches, i.e. experiment and molecular simulation. Along this path, the studied mixture goes through a transition in the supercritical region by crossing the Widom line and from super-to subcritical regimes. The Fick diffusion coefficient was measured with the Taylor dispersion technique, while several EMD simulation methods were employed. Shear viscosity, intra-and Maxwell-Stefan diffusion coefficients were sampled with equilibrium molecular dynamics, employing rigid, non-polarizable force fields based on Lennard-Jones sites and superimposed point charges or quadrupoles, and the Green-Kubo formalism. The thermodynamic factor was calculated with Kirkwood-Buff integration. In this way, the Fick diffusion coefficient was determined by EMD simulations consistently on the basis of the selected force fields. In order to validate the force fields employed to describe the mixture, the VLE of the mixture was calculated at three temperatures and compared successfully with experimental data from the literature. Further, to gain insight into the microscopic structure of the mixture, center-of-mass radial distribution functions as well as averaged coordination numbers were investigated.
In the region where the crossover between the liquid-like and gas-like regimes occurs, several thermophysical properties become very sensitive to temperature and pressure variations and can therefore change drastically along supercritical paths. Along the studied isobar, the presence of inflection points in the temperature dependence of density, enthalpy, shear viscosity, intra-and transport diffusion coefficients were reported. Mostly, the inflection points were observed at temperatures between T B 321 and B324 K, a range that specifies the region influenced by the Widom line. Moreover, the calculated values of the center-of-mass average coordination numbers up to the first solvation shell also exhibit an inflection point in this region, shedding light on the underlying microscopic structural background of the observed macroscopic property behavior.
The transition between the supercritical gas-like and the subcritical compressed gas regimes at T B 328 K is smooth for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell-Stefan and CO 2 intradiffusion coefficients. However, the intradiffusion coefficient of ethanol and the Fick diffusion coefficient showed an inflection point in the proximity of this transition. It was noticed that the averaged coordination number of ethanol-ethanol pairs, which is an indicator of ethanol self-association, shows an almost stepwise change in this region. This was related to a relatively strong decrease of ethanol self-association when moving from supercritical to subcritical states, which could also be observed by analyzing snapshots of the simulation volumes.
A satisfactory agreement was found between predictive EMD simulation data and experimental results for the Fick diffusion coefficient, especially in the temperature range from the Widom line to the super-to subcritical transition. Both experiment and simulation exhibit a sigmoidal behavior along the studied isobar, however, experiments showed an inflection point at a slightly lower temperature. Possible reasons for the difference between simulation and experimental results were thoroughly discussed.
An analysis of the shear viscosity of the mixture corroborated the strong influence of microheterogeneities given by hydrogen bonded networks on the macroscopic transport properties of the mixture.

Conflicts of interest
There are no conflicts to declare.