Molecular Dynamics Simulation of Aqueous Solutions Using Interaction Energy Components: Application to the Solvation Gibbs Energy

Molecular dynamics simulations of aqueous solutions of the solutes acetamide (AcNH₂), acetic acid (AcOH), and acetaldehyde (AcH) were made using Lennard–Jones 12-6-1 potentials to describe the solute–solvent interactions. The Morokuma decomposition scheme and the ESIE solute atomic charges were used to reproduce the exchange, polarization, and electrostatic components of the solute–water interaction energy. A nonlinear perturbation was incorporated into the “slow-growth” technique in order to improve the results for the solvation Gibbs energy that were found to be in agreement with the available experimental and theoretical values.     

Modified PT-LJ-SAFT Density Functional Theory: I. Prediction of surface properties and phase equilibria of non-associating fluids

A new modified version of a Perturbation Density Functional Theory (PT-DFT) based on the Statistical Association Fluid Theory (SAFT) with a Lennard–Jones interaction potential is proposed to model the vapor–liquid phase equilibrium and to predict the interfacial behavior of non-associating hydrocarbon fluids. In the interaction model for the Helmholtz free energy functional the molecules are separated into m spherical segments interacting via a Lennard–Jones potential. The segments form chains of tangent spheres. In the perturbation approximation to Density Functional Theory the interaction potential is split according to WCA and the attractive term to the free energy functional consists of a suitable modification of the perturbation expression. This modification to PT-DFT yields surface tensions for the Lennard–Jones sphere fluid (m = 1.0) which are in perfect agreement with simulation data.The new PT-DFT model combines the high flexibility of the SAFT free energy functional with a modified density functional approach that enables to perform accurate calculations of interfacial properties. To take into account the contributions to surface tension resulting from mesoscale thermal fluctuations a semiempirical model is proposed that allows to correct the microscopic intrinsic surface tension.The model is used to describe the phase equilibrium of lower alkanes and aromatics. The results demonstrate the capability to fit vapor–liquid equilibrium data and to predict very accurately the surface properties of these fluids within the uncertainties of the experimental data.     

The effect of corrugation of pore wall on the thermal diffusion in nanopores by molecular dynamics simulations

In this work we have studied the effect of corrugation on the thermal diffusion (soret effect) in isotopic and non-isotopic fluid mixtures confined in a slit pore. We used a boundary driven non-equilibrium molecular dynamics to simulate thermal diffusion in Lennard–Jones (LJ) binary mixtures confined in structureless Steele 10-4-3 and atomistic Lennard–Jones pore walls. The results showed that for the isotopic mixture thermal diffusion factor for both wall types agrees and the corrugation of the LJ wall has no effect in isotopic mixture. However, for non-isotopic mixture confined in atomistic LJ pore the component with stronger attraction adsorbs more to the wall than the structureless Steele wall. The effect of corrugation of pore wall on the thermal diffusion is noticeable in narrow slit pore and mixture with large difference in molecular attraction parameter of components.     

Sound Velocities of Generalized Lennard-Jones (n − 6) Fluids Near Freezing

In a recent paper [S. Khrapak, Molecules 25, 3498 (2020)], the longitudinal and transverse sound velocities of a conventional Lennard–Jones system at the liquid–solid coexistence were calculated. It was shown that the sound velocities remain almost invariant along the liquid–solid coexistence boundary lines and that their magnitudes are comparable with those of repulsive soft-sphere and hard-sphere models at the fluid–solid phase transition. This implies that attraction does not considerably affect the magnitude of the sound velocities at the fluid–solid phase transition. This paper provides further evidence to this by examining the generalized Lennard–Jones (n − 6) fluids with n ranging from 12 to 7 and demonstrating that the steepness of the repulsive term has only a minor effect on the magnitude of the sound velocities. Nevertheless, these minor trends are identified and discussed.     

On predicting self-diffusion coefficients in fluids

A modification of an existing correlative equation for self-diffusion coefficients is presented. A free-volume theory for liquids has been corrected to reproduce diffusivities of hard-sphere and Lennard–Jones fluids from low-density limit to melting points, and has been applied to correlate a wide database of non-polar, polar, quantic and hydrogen-bonding substances, although it is unable to fit helium, water and hydrogen fluoride. The adjustable parameters are generalized with available fluid properties, and the three resulting predictive formulas present lower deviations than other correlative equations used in a predictive way for non-quantic and non-associated fluids.     

Shear viscosity of Stockmayer fluid: Application of integral equations method to Vesovic–Wakeham scheme

In our previous works, we applied the integral equations method to calculate transport properties of nonpolar fluids such as Lennard–Jones (12-6) fluid [R. Khordad, Physica A 387 (2008) 4519, M.M. Papari, R. Khordad, Z. Akbari, Physica A 388 (2009) 585]. The present work is a continuation of our studies on transport properties of polar fluids. We use the Stockmayer potential and examine theoretically the viscosity and pressure of several refrigerant mixtures such as R125 + R143a, HFC-125 + HFC-134a, HFC-125 + HFC-32, and HFC-134a + HFC-32. We solve numerically the Ornstein–Zernike (OZ) equation using the hypernetted-chain approximation (HNC) for binary fluid mixtures and obtain the pair correlation functions. Finally, the density and temperature dependence of shear viscosity and pressure are studied using Vesovic–Wakeham method and compared with experimental results. According to the results obtained from the present work reveals that the integral equations method is suitable for predicting the pressure and shear viscosity of this class of fluids.     

The van der Waals one-fluid model for viscosity in Lennard–Jones fluids: Influence of size and energy parameters

This work aims to estimate the limitations of the van der Waals one-fluid (vdW1) approximation in the prediction of the viscosity of Lennard–Jones (LJ) mixtures. To do so, non-equilibrium molecular dynamics simulations have been performed. Results on mixtures have been compared to those of their equivalent pure fluids (in the sense of the vdW1 model). Several systems (146 configurations) are studied, which are composed of binary and ternary mixtures in various thermodynamic states and for different combining rules. In a first step, deviations induced separately by LJ molecular parameters (size or energy) have been analyzed. It is shown that the vdW1 model is well designed for the energy parameter in every configuration. On the contrary, for the size parameter, deviations induced by this one-fluid approach are shown to be large in dense state and low temperature systems. In a second step, the coupling effects of the LJ size and energy parameters with the mass are studied. It appears that an accurate one-fluid approximation for viscosity should involve a coupling between the mass and the size parameters in its formulation (which is not the case with the vdW1 model) but not between the mass and the energy.     

Thermodynamics of 1-alkanol+linear alkanoate mixtures

1-Alkanol?+?linear alkanoate mixtures have been investigated in the framework of the DISQUAC model. The interaction parameters for the OH/COO contacts are reported. The quasichemical parameters are independent of the mixture compounds. The dispersive parameters change with the molecular structure of the components. The same behaviour is observed for the OH/CO (carbonyl) and OH/OCOO (carbonate) contacts. DISQUAC represents well the molar excess Gibbs energies, coordinates of azeotropes and molar excess enthalpies. Using binary parameters only, DISQUAC improves meaningfully predictions on this property from the UNIFAC model for 1-alkanol?+?linear alkanoate?+?hydrocarbon systems. In contrast, the Nitta–Chao and the DISQUAC models yield similar results for the thermodynamic properties of the binary and ternary mixtures considered. 1-Alkanol?+?linear alkanoate mixtures are characterized by strong dipolar interactions between like molecules. In 1-alkanol?+?CH₃COO(CH₂) _{u ?1}CH₃ systems, dipole–dipole interactions between ester molecules are more important for u?≤?7. For u?≥?8, the more important contribution to the excess molar enthalpy comes from the disruption of the alkanol–alkanol interactions. For systems containing a polar compound such as alkanone, alkanoate or linear organic carbonate, dipolar interactions increase in the order: alkanone?<?alkanoate?<?carbonate.     

Comparison of perturbation theory and mean spherical approximation for polar fluids and ion–dipole mixtures based on molecular simulation data

A comprehensive study on various internal energies, pressures and chemical potentials for the pure dipolar hard sphere fluids, pure Stockmayer fluids, the Lennard–Jones and Stockmayer mixtures, the Stockmayer and Stockmayer mixtures and the ion–dipole mixtures is reported based on the perturbation theory (PT) and mean spherical approximation (MSA). Compared with the results of molecular simulations, it is shown that the PT is superior to MSA in most cases.     

Influence of polyvinylalcohol on the solvation volumes of FeCl<Subscript>3</Subscript> and CoCl<Subscript>2</Subscript> in aqueous ethanol solutions

The effect of polyvinylalcohol (PVA) on molar, Van-der Waals and the electrostriction volumes of FeCl₃ and CoCl₂ in pure H₂O and 50% ethanol (EtOH)–water mixtures was studied. Excess volumes were also calculated and their values were discussed in terms of Pierotti theory (Scaled Particle Theory). Different theoretical energy values were calculated to explain the following parameters: volumes, dispersion, induction energies and 6-12 Lennard–Jones parameters. The increase and the decrease in excess volumes of FeCl₃, CoCl₂ in the absence and in the presence of polyvinyl alcohol in both aqueous and 50% EtOH–H₂O mixtures were also discussed in terms of the increase or decrease of different interaction energies.     

Dielectric Relaxation of Some Rigid Molecules in Viscous Media

Abstract

The dielectric relaxation times(τ) have been evaluated from the measured values of dielectric constant(?′) and loss (?′′) of four rigid molecules in paraffin oil + benzene mixtures of varying viscosity extending from moderate to high values. Two of the rigid systems are nearly spherical in nature. An attempt has been made to find out the suitable viscosity term describing the dipolar rotation by applying the measured parameters to the earlier models as proposed by Debye, Hill and Kalman. Hill's model which describes the motion by a mutual viscosity (η_AB) term has been found to yield negative values showing its inadequacy for the systems rotating in high viscosity media. A comparison of Kalman's equation with Debye shows that the former gives better results for the region investigated both in high and low viscosity media. Even the systems having spherical symmetry give more satisfactory results when applied to Kalman's equation. In the process, the thermodynamical parameters and dipole moments of the molecules have also been evaluated and discussed.     

Towards understanding the effect of electrostatic interactions on the density of ionic liquids

In order to have a better understanding on the electrostatic contribution to the thermodynamic property of ionic liquids (ILs), a two-parameter equation of state (EOS) is developed on the basis of hard sphere perturbation theory by accounting for the dispersion interaction with Cotterman et al.’s EOS for L-J fluid and electrostatic interaction with mean spherical approximation (MSA) approach. The EOS is applicable for the density correlation of molecular liquids, and the resulting parameters, viz. Lennard–Jones dispersive parameter ?/k and soft-core diameter σ, can be used to predict the density of molecular mixtures and the corresponding ILs. The results indicate that the density of IL is always about 10% higher than the corresponding stoichiometric molecular mixture with which the IL is produced as an ionic adduct, for example, IL 1-methyl-3-methylimidazolium dimethylphosphate ([MMIM][DMP]) versus equimolar mixture of 1-methylimidazole (MIM) and trimethylphosphate (TMP). Furthermore, the density enhancement of ILs with respect to their corresponding stoichiometric molecular mixtures can be well represented by the electrostatic contribution among ionic species involved.     

Influence of a Surface Microheterogeneity on the Disjoining Pressure of Adsorbed Lennard–Jones Fluid. Computer Simulation

Computer simulation of the adsorption of Lennard–Jones fluid in slitlike micropores with structured walls was performed by the Monte Carlo and molecular dynamic methods. The influence of pore width, temperature, pressure, and surface structure on the disjoining pressure was considered. Structural and energy characteristics of the surface affect significantly the value of the disjoining pressure. Among the surface structural characteristics, the important role is played by the number and localization of vacancy-type defects.     

Molecular dynamics simulation of acetamide solvation using interaction energy components: Application to structural and energy properties

Molecular dynamics simulation of an aqueous solution of acetamide was performed using Lennard–Jones 12-6-1 potentials to describe the solute–solvent interactions, and TIP3P to describe the water–water interactions. The Morokuma decomposition scheme and the ESIE solute atomic charges were used to reproduce the molecular parameters of the solute–water interaction potential. The results showed that the functions that use the EX-PL-DIS-ES interaction model lead to good values of the structural and energy properties (in particular, the hydration shell and the solvation energies) when they are compared with those from using AMBER-derived parameters, and with the available theoretical and experimental data.     

Quantum computation of the properties of helium using two-body and three-body intermolecular potentials: a molecular dynamics study

We have performed the molecular dynamics simulation to obtain energy, pressure, and self-diffusion coefficient of helium at different temperatures and densities using Lennard–Jones (LJ), Hartree–Fock dispersion-Individual damping (HFD-ID) potential, and the HFD-like potential which has been obtained with an inversion of viscosity data at zero pressure supplemented by quantum corrections following the Feynman–Hibbs approach. The contribution of three-body interactions using an accurate simple relationship reported by Wang and Sadus between two-body and three-body interactions has been also involved for non-effective potentials (HFD-ID and HFD-like) in simulation. Our results show a good agreement with corresponding experimental data. A comparison of our simulated results with other molecular simulations using different potentials is also included.     

Molecular simulation of oxygen on supported platinum clusters

Molecular dynamics simulations are performed to study oxygen adsorption on platinum clusters supported on a graphite surface. The Sutton–Chen many-body potential is used for the Pt–Pt interaction, whereas a Steele potential was used to represent the carbon surface. The oxygen–oxygen intramolecular force is modeled by a harmonic oscillator model and other interactions are described by the Lennard–Jones potential. The results indicate an optimum loading of platinum for maximum specific adsorption of oxygen. Adsorption isotherms are constructed and the energies and orientation of adsorbed oxygen are reported. The relevance of this study to electrode processes is discussed.     

Viscosity Prediction for Oxygen,Nitrogen and Their Mixtures at Zero and Moderately Dense Regimes via Semi‐Empirically Based Assessment

The viscosity coefficients for the gaseous states of N₂ and O₂ and their mixtures are determined at zero and moderately density regimes. The Lennard‐Jones 12–6 (LJ 12–6) potential energy function is used as the initial model potential required y the technique. The interaction potential energies from the inversion procedure reproduce the viscosity commensurate to the best measurements. The initial density dependence of gaseous viscosity coefficient according to the Rainwater‐Friend theory, which was given by Najafi et al., has been considered for pure N₂ and pure O₂.     

Modeling the phase behavior of H2S+n-alkane binary mixtures using the SAFT-VR+D approach

A statistical associating fluid theory for potential of variable range has been recently developed to model dipolar fluids (SAFT-VR+D) [Zhao and McCabe, J. Chem. Phys. 2006, 125, 104504]. The SAFT-VR+D equation explicitly accounts for dipolar interactions and their effect on the thermodynamics and structure of a fluid by using the generalized mean spherical approximation (GMSA) to describe a reference fluid of dipolar square-well segments. In this work, we apply the SAFT-VR+D approach to real mixtures of dipolar fluids. In particular, we examine the high-pressure phase diagram of hydrogen sulfide+n-alkane binary mixtures. Hydrogen sulfide is modeled as an associating spherical molecule with four off-center sites to mimic hydrogen bonding and an embedded dipole moment (micro) to describe the polarity of H2S. The n-alkane molecules are modeled as spherical segments tangentially bonded together to form chains of length m, as in the original SAFT-VR approach. By using simple Lorentz-Berthelot combining rules, the theoretical predictions from the SAFT-VR+D equation are found to be in excellent overall agreement with experimental data. In particular, the theory is able to accurately describe the different types of phase behavior observed for these mixtures as the molecular weight of the alkane is varied: type III phase behavior, according to the scheme of classification by Scott and Konynenburg, for the H2S+methane system, type IIA (with the presence of azeotropy) for the H2S+ethane and+propane mixtures; and type I phase behavior for mixtures of H2S and longer n-alkanes up to n-decane. The theory is also able to predict in a qualitative manner the solubility of hydrogen sulfide in heavy n-alkanes.