Diffusion in Multicomponent Liquids: From Microscopic to Macroscopic Scales

Measurement of Diffusion of Atmospheric Gases in a Liquid Perfluoro Compound by Means of an Optical Technique

We report on accurate measurements of Fickian diffusion coefficients in binary mixtures consisting of hydrofluoroether (a perfluoro compound of methoxy-nonafluorobutane or HFE-7100) with dissolved atmospheric gases CO₂, N₂, and O₂ in the limit of an infinite dilution of the gas. We show that the use of optical digital interferometry (ODI) allows the determination of diffusion coefficients of dissolved gases with relatively small standard uncertainties for this class of experiments. In addition, we illustrate the ability of an optical approach to determine the gas concentration. We compare the capacity of four mathematical models, singly used in the literature, to obtain diffusion coefficients by applying them to the processing of a large amount of experimental data. We quantify their systematic errors and standard uncertainties. The temperature behavior of the diffusion coefficients measured in the range of 10 to 40 °C is consistent with the behavior of the same gases in other solvents available in the literature.

Diffusion in multicomponent aqueous alcoholic mixtures

The Fick diffusion coefficient matrix of the highly associating quaternary mixture water + methanol + ethanol + 2-propanol as well as its ternary and binary subsystems is analyzed with molecular dynamics simulation techniques. Three of the ternary subsystems are studied in this sense for the first time. The predictive capability of the employed force fields, which were sampled with the Green–Kubo formalism and Kirkwood–Buff integration, is confirmed by comparison with experimental literature data on vapor-liquid equilibrium, shear viscosity and Fick diffusion coefficient, wherever possible. A thorough analysis of the finite size effects on the simulative calculation of diffusion coefficients of multicomponent systems is carried out. Moreover, the dependence of the Fick diffusion coefficient matrix on the velocity reference frame and component order is analyzed. Their influence is found to be less significant for the main matrix elements, reaching a maximum variation of 19%. The large differences found for the cross elements upon variation of the reference frame hinder a straightforward interpretation of the Fick diffusion coefficient matrix with respect to the presence of diffusive coupling effects.

Generalized Form for Finite-Size Corrections in Mutual Diffusion Coefficients of Multicomponent Mixtures Obtained from Equilibrium Molecular Dynamics Simulation

The system-size dependence of computed mutual diffusion coefficients of multicomponent mixtures is investigated, and a generalized correction term is derived. The generalized finite-size correction term was validated for the ternary molecular mixture chloroform/acetone/methanol as well as 28 ternary LJ systems. It is shown that only the diagonal elements of the Fick matrix show system-size dependency. The finite-size effects of these elements can be corrected by adding the term derived by Yeh and Hummer (J. Phys. Chem. B2004, 108, 15873–15879). By performing an eigenvalue analysis of the finite-size effects of the matrix of Fick diffusivities we show that the eigenvector matrix of Fick diffusivities does not depend on the size of the simulation box. Only eigenvalues, which describe the speed of diffusion, depend on the size of the system. An analytic relation for finite-size effects of the matrix of Maxwell–Stefan diffusivities was developed. All Maxwell–Stefan diffusivities depend on the system size, and the required correction depends on the matrix of thermodynamic factors.

Diffusion of methane in supercritical carbon dioxide across the Widom line

Diffusion of methane diluted in supercritical carbon dioxide is studied by experiment and molecular simulation in the temperature range from 292.55 to 332.85 K along the isobars 9.0, 12.5 and 14.7 MPa. Measurements of the Fick diffusion coefficient are carried out with the Taylor dispersion technique. Molecular dynamics simulation and the Green-Kubo formalism are employed to obtain Fick, Maxwell-Stefan and intradiffusion coefficients as well as shear viscosity. The obtained diffusion coefficients are on the order of 10-8 m2/s. The composition, temperature and density dependence of diffusion is analyzed. The Fick diffusion coefficient of methane in carbon dioxide shows an anomaly in the near-critical region. This behavior can be attributed to the crossing of the so-called Widom line, where the supercritical fluid goes through a transition between liquid-like and gas-like states. Further, several classical equations are tested on their ability to predict this behavior and it is found that equations that explicitly include the density are better suited to predict the sharp variation of the diffusion coefficient near the critical region predicted by molecular simulation.

Do ternary liquid mixtures exhibit negative main Fick diffusion coefficients?

Measured main Fick diffusion coefficients are throughout positive. However, they may appear to be negative after transformation to the molar reference frame, if the excess volume is significant and the experimental uncertainties of the cross diffusion coefficients are large.

Optimizing Nonbonded Interactions of the OPLS Force Field for Aqueous Solutions of Carbohydrates: How to Capture Both Thermodynamics and Dynamics

Knowledge on thermodynamic and transport properties of aqueous solutions of carbohydrates is of great interest for process and product design in the food, pharmaceutical, and biotechnological industries. Molecular simulation is a powerful tool to calculate these properties, but current classical force fields cannot provide accurate estimates for all properties of interest. The poor performance of the force fields is mainly observed for concentrated solutions, where solute–solute interactions are overestimated. In this study, we propose a method to refine force fields, such that solute–solute interactions are more accurately described. The OPLS force field combined with the SPC/Fw water model is used as a basis. We scale the nonbonded interaction parameters of sucrose, a disaccharide. The scaling factors are chosen in such a way that experimental thermodynamic and transport properties of aqueous solutions of sucrose are accurately reproduced. Using a scaling factor of 0.8 for Lennard-Jones energy parameters (ϵ) and a scaling factor of 0.95 for partial atomic charges (q), we find excellent agreement between experiments and computed liquid densities, thermodynamic factors, shear viscosities, self-diffusion coefficients, and Fick (mutual) diffusion coefficients. The transferability of these optimum scaling factors to other carbohydrates is verified by computing thermodynamic and transport properties of aqueous solutions of d-glucose, a monosaccharide. The good agreement between computed properties and experiments suggests that the scaled interaction parameters are transferable to other carbohydrates, especially for concentrated solutions.

Composition dependent transport diffusion in non-ideal mixtures from spatially resolved nuclear magnetic resonance spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy is a well-established technique for the measurement of intra-diffusion coefficients. Recently, such information has been used as a basis of predictive models to extrapolate to the Fick diffusion coefficient of liquid mixtures. The present work presents a new approach to directly access the Fick diffusion coefficient by spatially resolved NMR experiments. The Fick diffusion coefficient of the binary mixture TEA/H2O was determined at two temperatures, 283.2 K and 275.2 K. The results are consistent with values previously reported either from optical experiments or predictive Darken-type models developed for this system. The proposed methodology adds high-resolution NMR to the toolbox for the study of the transport diffusion of multicomponent mixtures. It is, however, still limited to mixtures with liquid-liquid equilibrium phase separation.

Interplay of structure and diffusion in ternary liquid mixtures of benzene + acetone + varying alcohols

The Fick diffusion coefficient matrix of ternary mixtures containing benzene + acetone + three different alcohols, i.e., methanol, ethanol, and 2-propanol, is studied by molecular dynamics simulation and Taylor dispersion experiments. Aiming to identify common features of these mixtures, it is found that one of the main diffusion coefficients and the smaller eigenvalue do not depend on the type of alcohol along the studied composition path. Two mechanisms that are responsible for this invariant behavior are discussed in detail, i.e., the interplay between kinetic and thermodynamic contributions to Fick diffusion coefficients and the presence of microscopic heterogeneities caused by hydrogen bonding. Experimental work alone cannot explain these mechanisms, while present simulations on the molecular level indicate structural changes and uniform intermolecular interactions between benzene and acetone molecules in the three ternary mixtures. The main diffusion coefficients of these ternary mixtures exhibit similarities with their binary subsystems. Analyses of radial distribution functions and hydrogen bonding statistics quantitatively evidence alcohol self-association and cluster formation, as well as component segregation. Furthermore, the excess volume of the mixtures is analyzed in the light of intermolecular interactions, further demonstrating the benefits of the simultaneous use of experiment and simulation. The proposed framework for studying diffusion coefficients of a set of ternary mixtures, where only one component varies, opens the way for further investigations and a better understanding of multicomponent diffusion. The presented numerical results may also give an impulse to the development of predictive approaches for multicomponent diffusion.

Finite-Size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics

Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.

Mutual diffusion governed by kinetics and thermodynamics in the partially miscible mixture methanol + cyclohexane

To gain an understanding of the transport and thermodynamic behavior of the highly non-ideal mixture methanol + cyclohexane, three complementary approaches, i.e. experiment, molecular simulation and predictive equations, are employed. The temperature and composition dependence of different diffusion coefficients is studied around the miscibility gap at ambient pressure. On the one hand Fick diffusion coefficients are measured experimentally by interferometric probing and on the other hand Maxwell-Stefan diffusion coefficients and intradiffusion coefficients are sampled by equilibrium molecular dynamics simulation at five temperatures below the upper critical temperature of ∼319 K. The spinodal curve is determined from extrapolation of the experimental Fick diffusion coefficient data and compared to predictions from excess Gibbs energy models. It is found that these models are not capable to correctly describe the activity coefficients over the whole composition range of the studied mixture. Thus, different parameter sets for a modified Wilson model are used for calculations of the thermodynamic factor, which is needed to transform Maxwell-Stefan into Fick diffusion coefficients and vice versa. Further, predictive equations for the Maxwell-Stefan diffusion coefficient, which are based on intradiffusion coefficients, are compared to simulation results. Using different approaches provides a clearer understanding of the relations between kinetic and thermodynamic properties contributing to the diffusion behavior of partially miscible mixtures.

The effect of alcohols as the third component on diffusion in mixtures of aromatics and ketones

The Fick diffusion coefficient matrix of three ternary mixtures composed of an aromatic (benzene), a ketone (acetone) and one of three different alcohols (methanol, ethanol or 2-propanol) is investigated with laboratory and numerical work.