Seyed Hossein Jamali

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.

Absorption Refrigeration Cycles with Ammonia-Ionic Liquid Working Pairs Studied by Molecular Simulation

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][Tf2N], and [emim][SCN]. As refrigerant NH3 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.

Gibbs ensemble Monte Carlo simulations of multicomponent natural gas mixtures

Abstract 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 (, , , S, and ), and the other containing a mixture of nine components (, , , S, , , , and ) 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.

Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics

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.

Thermodynamic and Transport Properties of Crown-Ethers: Force Field Development and Molecular Simulations

Crown-ethers have recently been used to assemble porous liquids (PLs), which are liquids with permanent porosity formed by mixing bulky solvent molecules (e.g., 15-crown-5 ether) with solvent-inaccessible organic cages. PLs and crown-ethers belong to a novel class of materials, which can potentially be used for gas separation and storage, but their performance for this purpose needs to be assessed thoroughly. Here, we use molecular simulations to study the gas separation performance of crown-ethers as the solvent of porous liquids. The TraPPE force field for linear ether molecules has been adjusted by fitting a new set of torsional potentials to accurately describe cyclic crown-ether molecules. Molecular dynamics (MD) simulations have been used to compute densities, shear viscosities, and self-diffusion coefficients of 12-crown-4, 15-crown-5, and 18-crown-6 ethers. In addition, Monte Carlo (MC) simulations have been used to compute the solubility of the gases CO2, CH4, and N2 in 12-crown-4 and 15-crown-5 ether. The computed properties are compared with available experimental data of crown-ethers and their linear counterparts, i.e., polyethylene glycol dimethyl ethers.