Spectral line shape models can successfully reproduce experimental Rayleigh-Brillouin spectra, but they need knowledge about the bulk viscosity ηb. Light scattering involves GHz frequencies, but since ηb is only documented at low frequencies, ηb is usually left as a free parameter, which is determined by a fit of the model to an experimental spectrum. The question is whether models work so well because of this freedom. Moreover, for light scattering in air, spectral models view "air"as an effective molecule. We critically evaluate the use of ηb as a fit parameter by comparing ηb obtained from fits of the Tenti S6 model to the result of Direct Simulation Monte Carlo (DSMC) for a mixture of Nitrogen and Oxygen. These simulations are used to compute light scattering spectra, which are then compared to experiments. The DSMC simulation parameters are cross-checked with a molecular dynamics simulation based on intermolecular potentials. At large values of the uniformity parameter y, y ≈ 4, where the Brillouin contribution to spectra is large, fitted ηb are 20% larger than the ones from DSMC, while the quality of the simulated spectra is comparable to that of the Tenti S6 line shape model. At smaller y, the difference between fitted and simulated ηb can be as large as 100%. We hypothesize the breakdown of the bulk viscosity concept to be the cause of this fallacy.

The number of molecules used in a typical Molecular Dynamics (MD) simulations is orders of magnitude lower than in the thermodynamic limit. It is therefore essential to correct diffusivities computed from Molecular Dynamics simulations for finite-size effects. We present a comprehensive review on finite-size effects of diffusion coefficients by considering self-, Maxwell–Stefan, and Fick diffusion coefficients in pure liquids, as well as binary, ternary, and quaternary mixtures. All finite-size corrections, both analytical and empirical, are discussed in detail. The finite-size effects of rotational and confined diffusion are also briefly discussed.

The translational self-diffusion coefficient and the shear viscosity of water are related by the fractional Stokes–Einstein relation. We report extensive novel molecular dynamics simulations for the self-diffusion coefficient and the shear viscosity of water. The SPC/E and TIP4P/2005 water models are used in the temperature range 220–560 K and at 1 or 1,000 bar. We compute the fractional exponents t, and s that correspond to the two forms of the fractional Stokes–Einstein relation (Formula presented.) and (Formula presented.) respectively. We analyse other available experimental and numerical simulation data. In the current analysis two temperature ranges are considered (above or below 274 K) and in both cases deviations from the Stokes–Einstein relation are observed with different values for the fractional exponents obtained for each temperature range. For temperatures above 274 K, both water models perform comparably, while for temperatures below 274 K TIP4P/2005 outperforms SPC/E. This is a direct result of the ability of TIP4P/2005 to predict water densities more accurately and thus predict more accurately the water self-diffusion coefficient and the shear viscosity.

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.

We present a new plugin for LAMMPS for on-the-fly computation of transport properties (OCTP) in equilibrium molecular dynamics. OCTP computes the self- and Maxwell-Stefan diffusivities, bulk and shear viscosities, and thermal conductivities of pure fluids and mixtures in a single simulation. OCTP is the first implementation in LAMMPS that uses the Einstein relations combined with the order-n algorithm for the efficient sampling of dynamic variables. OCTP has low computational requirements and is easy to use because it follows the native input file format of LAMMPS. A tool for calculating the radial distribution function (RDF) of the fluid beyond the cutoff radius, while taking into account the system size effects, is also part of the new plugin. The RDFs computed from OCTP are needed to obtain the thermodynamic factor, which relates Maxwell-Stefan and Fick diffusivities. To demonstrate the efficiency of the new plugin, the transport properties of an equimolar mixture of water-methanol were computed at 298 K and 1 bar.

The combination of the TraPPE and OPLS/2016 force fields with five water models, TIP3P, SPC/E, OPC, TIP4P/2005 and TIP4P/EW was used to compute mixing enthalpies, excess chemical potentials and activity coefficients of water and methanol. Excess chemical potentials and activity coefficients were computed in an expanded version of the NPT ensemble. We found the best agreement between experimental data for all the computed properties of water–methanol mixtures for the TIP4P/2005-TraPPE force fields. The performance of the spherical cutoff methods in MC and MD simulations was compared to the Ewald summation. The radial distribution functions obtained from the Ewald summation and the Damped-Shifted Force (DSF) method were in excellent agreement. Numerical artifacts appeared at the cutoff radius when the original Wolf method was used to calculate the electrostatic interactions. The calculated excess mixing enthalpies, excess chemical potentials, and activity coefficients of water and methanol obtained from the Wolf method were in good agreement with the DSF method. Our simulation results show that the numerical artifacts of the original Wolf method have little effect for energy calculations in aqueous methanol mixtures.

A method is proposed for calculating the shear viscosity of a liquid from finite-size effects of self-diffusion coefficients in Molecular Dynamics simulations. This method uses the difference in the self-diffusivities, computed from at least two system sizes, and an analytic equation to calculate the shear viscosity. To enable the efficient use of this method, a set of guidelines is developed. The most efficient number of system sizes is two and the large system is at least four times the small system. The number of independent simulations for each system size should be assigned in such a way that 50%-70% of the total available computational resources are allocated to the large system. We verified the method for 250 binary and 26 ternary Lennard-Jones systems, pure water, and an ionic liquid ([Bmim][Tf2N]). The computed shear viscosities are in good agreement with viscosities obtained from equilibrium Molecular Dynamics simulations for all liquid systems far from the critical point. Our results indicate that the proposed method is suitable for multicomponent mixtures and highly viscous liquids. This may enable the systematic screening of the viscosities of ionic liquids and deep eutectic solvents.

Mutual diffusion coefficients can be successfully predicted with models based on the Darken equation. However, Darken-based models require composition-dependent self-diffusion coefficients which are rarely available. In this work, we present a predictive model for composition-dependent self-diffusion coefficients (also called tracer diffusion coefficients or intradiffusion coefficients) in nonideal binary liquid mixtures. The model is derived from molecular dynamics simulation data of Lennard-Jones systems. A strong correlation between nonideal diffusion effects and the thermodynamic factor is observed. We extend the model by McCarty and Mason (Phys. Fluids 1960, 3, 908-922) for ideal binary gas mixtures to predict the composition-dependent self-diffusion coefficients in nonideal binary liquid mixtures. Our new model is a function of the thermodynamic factor, the self-diffusion coefficients at infinite dilution, and the self-diffusion coefficients of the pure substances, which are readily available. We validate our model with experimental data of 9 systems. For both Lennard-Jones systems and experimental data, the accuracy of the predicted self-diffusion coefficients is improved by a factor of 2 compared to the correlation of McCarty and Mason. Thus, our new model significantly expands the practical applicability of Darken-based models for the prediction of mutual diffusion coefficients.

For absorption refrigeration, it has been shown that ionic liquids have the potential to replace conventional working pairs. Due to the huge number of possibilities, conducting lab experiments to find the optimal ionic liquid is infeasible. Here, we provide a proof-of-principle study of an alternative computational approach. The required thermodynamic properties, i.e., solubility, heat capacity, and heat of absorption, are determined via molecular simulations. These properties are used in a model of the absorption refrigeration cycle to estimate the circulation ratio and the coefficient of performance. We selected two ionic liquids as absorbents: [emim][Tf₂N], and [emim][SCN]. As refrigerant NH₃ was chosen due to its favorable operating range. The results are compared to the traditional approach in which parameters of a thermodynamic model are fitted to reproduce experimental data. The work shows that simulations can be used to predict the required thermodynamic properties to estimate the performance of absorption refrigeration cycles. However, high-quality force fields are required to accurately predict the cycle performance.

Vapour–liquid equilibrium (VLE) and volumetric data of multicomponent mixtures are extremely important for natural gas production and processing, but it is time consuming and challenging to experimentally obtain these properties. An alternative tool is provided by means of molecular simulation. Here, Monte Carlo (MC) simulations in the Gibbs ensemble are used to compute the VLE of multicomponent natural gas mixtures. Two multicomponent systems, one containing a mixture of six components ((Formula presented.), (Formula presented.), (Formula presented.), (Formula presented.)S, (Formula presented.)(Formula presented.) and (Formula presented.)(Formula presented.)), and the other containing a mixture of nine components ((Formula presented.), (Formula presented.), (Formula presented.), (Formula presented.)S, (Formula presented.)(Formula presented.), (Formula presented.)(Formula presented.), (Formula presented.)(Formula presented.), (Formula presented.)(Formula presented.) and (Formula presented.)(Formula presented.)) are simulated. The computed VLE from the MC simulations is in good agreement with available experimental data and the GERG-2008 equation of state modelling. The results show that molecular simulation can be used to predict properties of multicomponent systems relevant for the natural gas industry. Guidelines are provided to setup Gibbs ensemble simulations for multicomponent systems, which is a challenging task due to the increased number of degrees of freedom.

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.

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.

Carbon dioxide (CO₂) can electrochemically be converted to a range of products including formic acid (HCOOH) and acetic acid (CH₃COOH). The yield of the products in an electrolysis cell depends on the solubility of CO₂ in the (aqueous) mixture. In absence of experimental data, Monte Carlo simulations in the Gibbs ensemble are used to compute the VLE of the binary systems, CO₂-H₂O, CO₂-HCOOH and CO₂-CH₃COOH, and the ternary systems, CO₂-HCOOH-H₂O and CO₂-CH₃COOH-H₂O. In addition, the PC-SAFT equation of state (EoS) is used to model the VLE of these strongly associating mixtures. Both methods correctly predicts the liquid-phase compositions, but the gas-phase compositions are less accurately described. The challenges to model these systems are related to the simultaneous formation of dimers, rings, and chains, which requires accurate force fields and advanced biasing schemes in MC simulations, and association theories that can account for this effect.

Natural gas, synthesis gas, and flue gas typically contain a large number of impurities (e.g., acidic gases), which should be removed to avoid environmental and technological problems, and to meet customer specifications. One approach is to use physical solvents to remove the acidic gases. If no experimental data are available, the solubility data required for designing the sweetening process can be obtained from molecular simulations. Here, Monte Carlo (MC) simulations are used to compute the solubility of the gas molecules, i.e., carbonyl sulfide, carbon disulfide, sulfur dioxide, hydrogen sulfide, methyl mercaptan, carbon dioxide, and methane in the commercial solvents tetraethylene-glycol-dimethyl-ether (Selexol), n-methyl-2-pyrrolidone, propylene carbonate, methanol (Rectisol), and the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([bmim][Tf₂N]). Henry coefficients of the gases in the investigated solvents are obtained from the computed solubilities. The ratio of Henry coefficients is used to compute ideal selectivities of the solvents. The solubilites and selectivities computed from MC simulations are compared with available experimental data. Some guidelines are provided to remove acidic gases using the investigated solvents. Rectisol is the best solvent for acid gas removal, but it should be used at low temperatures. Selexol and the ionic liquid have similar selectivity of sulfur compounds with respect to methane and may be used at elevated pressures and temperatures since both have low vapor pressures. The solubility of carbon disulfide, sulfur dioxide, and methyl mercaptan in these solvents is the highest. Hence, these components can be removed easily prior to hydrogen sulfide, carbonyl sulfide, and carbon dioxide in a pre-absorber.

Monte Carlo (MC) simulations in ensembles with a fixed chemical potential or fugacity, for example the grand-canonical or the osmotic ensemble, are often used to compute phase equilibria. Chemical potentials can be computed either with an equation of state (EoS) or from molecular simulations. The accuracy of the computed chemical potentials depends on the quality of the (critical) parameters used in the EoS and the applied force field in the simulations. We investigated the consistency of both approaches for computing fugacities of the industrially relevant gases CO₂, CH₄, CO, H₂, N₂, and H₂S. The critical temperature (T_c), pressure (P_c), and acentric factors (ω) of these gases are computed from MC simulations in the Gibbs ensemble. The effect of cutoff radius and tail corrections on the computed values of T_c, P_c, and ω is investigated. In addition, MC simulations in the Gibbs ensemble are used to compute the VLE of the 15 possible binary systems comprising the gases CO₂, CH₄, CO, H₂, N₂, and H₂S, and the ternary systems CO₂/CH₄/H₂S and CO₂/CO/H₂. Binary interaction parameters (k_ij) of these natural/synthesis gas mixtures are obtained by fitting the Peng-Robinson (PR) EoS to the binary VLE data from the MC simulations. The computed properties from the MC simulations are compared with the PR EoS, the GERG EoS, and experimental results. The MC results show that including tail corrections in the simulations is crucial to obtain accurate critical properties. The force fields used for the gases can reproduce the fugacities of the gases within 5% of the experimental data. The dew-point curves of all the 15 binaries were predicted correctly by the MC simulations, but the bubble-point curves for the systems H₂/CO, CH₄/H₂, H₂S/N₂, and H₂S/CO significantly deviate from the experiments.