Kirkwood–Buff (KB) integrals provide a connection between microscopic properties and thermodynamic properties of multicomponent fluids. The estimation of KB integrals using molecular simulations of finite systems requires accounting for finite size effects. In the small system method, properties of finite subvolumes with different sizes embedded in a larger volume can be used to extrapolate to macroscopic thermodynamic properties. KB integrals computed from small subvolumes scale with the inverse size of the system. This scaling was used to find KB integrals in the thermodynamic limit. To reduce numerical inaccuracies that arise from this extrapolation, alternative approaches were considered in this work. Three methods for computing KB integrals in the thermodynamic limit from information of radial distribution functions (RDFs) of finite systems were compared. These methods allowed for the computation of surface effects. KB integrals and surface terms in the thermodynamic limit were computed for Lennard–Jones (LJ) and Weeks–Chandler–Andersen (WCA) fluids. It was found that all three methods converge to the same value. The main differentiating factor was the speed of convergence with system size L. The method that required the smallest size was the one which exploited the scaling of the finite volume KB integral multiplied by L. The relationship between KB integrals and surface effects was studied for a range of densities.

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.

Monte Carlo simulation (MC) is combined with equations of state (EoS) to develop a methodology for the calculation of the vapor–liquid equilibrium (VLE) of multicomponent hydrocarbon mixtures with high asymmetry. MC simulations are used for the calculation of the VLE of binary methane mixtures with long n-alkanes, for a wide range of temperatures and pressures, to obtain sufficient VLE data for the consistent fitting of binary interaction parameters (BIPs) for the EoS. The Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), and Perturbed Chain - Statistical Associating Fluid Theory (PC-SAFT) EoS are considered. The ability of each EoS to correlate the VLE data is assessed and the selected ones are used to predict the VLE of multicomponent gas condensate mixtures. MC simulations proved to be very accurate in predicting the VLE in all conditions and mixtures considered. The BIPs regressed from the simulation dataset lead to equally accurate modeling results for multicomponent mixtures, compared to those regressed from experimental data.

We present a detailed overview of classical molecular simulation studies examining the self-diffusion coefficient of water. The self-diffusion coefficient is directly associated with the calculations of tracer or mutual diffusion coefficient of mixtures and, therefore, is a fundamental transport property, essential for an accurate description of mass transfer processes in biological, geological (i.e. energy or environmentally related), and chemical systems. In the current review we explore two distinct research areas. Namely, we discuss the self-diffusion of water in the bulk phase and under confinement. Different aspects that affect the diffusion process, including the molecular models, the system-size effects, the temperature and pressure conditions and the type of confinement are discussed. Finally, possible directions for future research are outlined.

Clathrate hydrates have characteristic properties that render them attractive for a number of industrial applications. Of particular interest are the following two cases: (i) the incorporation of large amounts of gas molecules into the solid structure has resulted in considering hydrates as possible material for the storage/transportation of energy or environmental gases, and (ii) the selective incorporation of guest molecules into the solid structure has resulted in considering hydrates for gas-mixture separations. For the proper design of such industrial applications, it is essential to know accurately a number of thermodynamic, structural and transport properties. Such properties can either be measured experimentally or calculated at different scales that span the molecular scale-up to the continuum scale. By using clathrate hydrates as a particular case study, we demonstrate that performing studies at multiple length scales can be utilised in order to obtain properties that are essential to process design.

An alternative method for calculating partial molar excess enthalpies and partial molar volumes of components in Monte Carlo (MC) simulations is developed. This method combines the original idea of Frenkel, Ciccotti, and co-workers with the recent continuous fractional component Monte Carlo (CFCMC) technique. The method is tested for a system of Lennard–Jones particles at different densities. As an example of a realistic system, partial molar properties of a [NH₃, N₂, H₂] mixture at chemical equilibrium are computed at different pressures ranging from P = 10 to 80 MPa. Results obtained from MC simulations are compared to those obtained from the PC-SAFT Equation of State (EoS) and the Peng–Robinson EoS. Excellent agreement is found between the results obtained from MC simulations and PC-SAFT EoS, and significant differences were found for PR EoS modelling. We find that the reaction is much more exothermic at higher pressures.