OPLS-AA Force Field Research Articles

Critical assessment of diffusion coefficients of benzene and its derivatives in supercritical carbon dioxide: Experimental data, molecular dynamics simulations, and modeling

The tracer diffusion coefficient (D12) is a critical transport property essential for research, design and optimization of industrial processes. However, experimental determination of D12 is often challenging, time-consuming, and frequently exhibits significant variability across different systems and conditions. In this study, a comprehensive assessment of D12 for benzene, phenol, chlorobenzene, anisole, and the three isomers of dichlorobenzene in supercritical carbon dioxide (▪) is presented. The available experimental data was meticulously compiled and analyzed, comparing different data points from various authors to assess consistency and accuracy. The experimental values were then compared with predictions from classical molecular dynamics (MD) simulations carried out in this work, in which three distinct ▪ force fields (EPM2, TraPPE, and Zhu et al.) and the OPLS-AA force field for the solute molecules were tested in combination with NVT and NPT ensembles. Among the evaluated combinations of force field and ensemble, the Zhu et al./NPT pair demonstrated excellent agreement with experimental data, achieving an overall deviation of only 5.63 %.Additionally, the performance of alternative approaches for estimating D12, including phenomenological and semi-empirical models—Wilke-Chang, Lai-Tan, Tracer Liu-Silva-Macedo and He-Yu equations—and machine learning (ML) algorithms was evaluated, highlighting their comparative accuracy while accounting for their inherent limitations. Overall, the ML model proved highly effective, with a deviation of only 2.42 % when applied to compounds within the training dataset. However, for compounds outside this dataset, the reliability, flexibility and accuracy of MD simulations, as demonstrated in this study, strongly suggest their application for D12 estimation in ▪/benzene-derivative systems, with possible additional valuable insights from the analysis of MD trajectories.

Study of molecular arrangement and density estimation of soybean oil biodiesel-diesel blends employing molecular dynamic simulation

Blends of diesel–biodiesel are a crucial alternative for petrodiesel substitution. However, biodiesel faces challenges such as high density and kinematic viscosity, which can affect the properties and performance of the blends. Based on this, the investigation of the impact of adding biodiesel into diesel fuel at percentages of 12%, 15%, 18%, and 21% using molecular dynamic simulations is presented. The density of the blends was experimentally measured and calculated with good agreement, which indicates a favorable description of the interactions by the OPLS-AA force field. Obtained results suggest that alkanes and isoalkanes showed higher fluctuations and were most sensitive to the variation of the percentage of biodiesel in diesel. The other classes of molecules did not show significant variations in interactions, suggesting that π-π interactions are unaffected by the biodiesel addition. The spatial distribution functions revealed that diesel molecules are organized homogeneously around methyl oleate and methyl linoleate, whereas aliphatic hydrocarbon and aromatic molecules are in distinct regions. The diffusion coefficients (Di) showed that the diesel molecules possess higher values, and the biodiesel addition increases the Di of the blend due to the methyl esters carbon chains, which present lower diffusion than diesel molecules. Finally, physical parameters and simulations revealed that the biodiesel addition into diesel fuel in the proportion up to 21% did not negatively impact the density and viscosity limits established by standards ASTM, EN, and ANP Resolution.

Atomistic molecular dynamics simulation and COSMO-SAC approach for enhanced 1,3-propanediol extraction with imidazolium-based ionic liquids.

1,3-Propanediol (1,3-PDO) is a key chemical in various industries, including pharmaceuticals and material sciences, and is projected to see significant market growth. However, the current challenges in its downstream processing, particularly in terms of cost and efficiency, highlight the need for innovative solutions. Our study delves into using ionic liquids (ILs) as a potential alternative, aiming to address these critical separation challenges more sustainably and efficiently. In this study, we utilized molecular dynamics (MD) simulations and the COSMO-SAC to examine 1,3-propanediol (1,3-PDO) extraction using four imidazolium-based ionic liquids with 1-butyl-3-methylimidazolium [Bmim] cation and with different anions bis(pentafluoroethanesulfonyl)imide [NPF2]-, bis(trifluoromethylsulfonyl)imide [NTF2]-, thiocyanate [SCN]-, and trifluoromethanesulfonate [TFO]-. Molecular dynamics simulations, incorporating analysis of radial distribution functions (RDF) and spatial distribution functions (SDF), revealed that [Bmim][SCN] and [Bmim][TFO] exhibit enhanced interactions with 1,3-PDO. Notably, [Bmim][SCN] formed the most hydrogen bonds, averaging 1.639 per molecule, due to its coordinating [SCN]- anion. This was in contrast to the fewer hydrogen bonds formed by non-coordinating anions in [Bmim][NPF2] and [Bmim][NTF2]. In ternary systems, [Bmim][SCN] and [Bmim][TFO] demonstrated superior selectivity for 1,3-PDO extraction compared to the other ionic liquids, with selectivity values around 29. These findings, supported by COSMO-SAC predictive modeling, highlight the potential of [Bmim][SCN] as a promising candidate for 1,3-PDO extraction, emphasizing the importance of anion selection in optimizing ionic liquid properties for this application. In our study, we employed MD simulations, incorporating the OPLS-AA force field, and COSMO-SAC to investigate the extraction of 1,3-PDO using imidazolium-based ionic liquids: [Bmim][NTF2], [Bmim][NPF2], [Bmim][SCN], and [Bmim][TFO]. The MD simulations were conducted using LAMMPS software, focusing on elucidating the RDF, SDF, and hydrogen bonding. Analysis of the distribution coefficient (β) and selectivity (S) for the ternary mixture was also conducted. These aspects of the simulation were analyzed using TRAVIS and VMD software. Additionally, the COSMO-SAC model was employed to determine the activity coefficients of 1,3-PDO in the ionic liquids, with molecular optimization conducted using Gaussian16 and sigma profile calculations performed using COSMO-SAC.

Extraction of Mechanical Parameters via Molecular Dynamics Simulation: Application to Polyimides.

Polyimides feature a vast number of industrial applications due to their high thermal stability and insulation properties. These polymers exhibit an exceptional combination of thermal stability and mechanical toughness, which allows the semiconductor industry to use them as a mechanical stress buffer. Here, we perform all-atom molecular dynamics (MD) simulations for such materials to assess their predictive capability with respect to their mechanical properties. Specifically, we demonstrate that the OPLS-AA force field can be used to successfully describe an often-used polyimide (i.e., Kapton®) with respect to its Young's modulus and Poisson's ratio. Two different modes to extract these mechanical properties from MD simulations are presented. In particular, our continuous deformation mode simulations almost perfectly replicate the results from real-world experimental data and are in line with predictions using other MD force fields. Our thorough investigation of Kapton® also includes an analysis of the anisotropy of normal stresses, as well as the effect of simulation properties on the predicted Young's moduli. Furthermore, the polyimide pyromellitic dianhydride/2-(4-aminophenyl)-1H-benzimidazole-5-amine (PMDA-BIA) was investigated to draw a more thorough picture of the usability of the OPLS-AA force field for polyimides.

Counter-ion adsorption and electrostatic potential in sodium and choline dodecyl sulfate micelles - a molecular dynamics simulation study.

Choline-based surfactants are interesting both from the practical point of view to obtaining environmental-friendly surfactants as well as from the theoretical side since the interactions between the choline and surfactants can help to understand self-assembly phenomena in deep eutectic solvents. Although no significant change was noticed in the micelle size and shape due to the exchange of the sodium counter-ion by choline in our simulations, the adsorption of the choline cation over the micelle surface is stronger than the adsorption of the sodium, which leads to a reduction of the exposed surface area of the micelle and remarkable effects over the electrostatic potential. The choline neutralizes the surface charge of the surfactant better than sodium; however, this is partially compensated by a stronger water orientation around the SDS micelle. The balance between the contributions from the surfactant, the counter-ion, and water to the electrostatic potential leads to a complex pattern with alternate regions of positive and negative potential at the micelle/water interface which can be important to the incorporation of other charged species at the micelle surface as well as for the interaction between micelles in solution. To evaluate the effects of the counter-ion substitution, micelles of sodium dodecyl sulfate (SDS) and choline dodecyl sulfate (ChDS) were studied and compared by means of molecular dynamics simulations in aqueous solution. In both cases, the simulations started from pre-assembled micelles with 60 dodecyl sulfate ions and 240-ns simulations were performed at NPT ensemble at T = 323.15K and P = 1bar using the Gromacs software with the OPLS-AA force field to describe dodecyl sulfate and choline, Åqvist parameters for sodium, and SPC model for water molecules.

Theoretical study on formation mechanism of acetic acid associating configurations and their distributions under saturated conditions.

The vapor-liquid equilibrium (VLE) properties of acetic acid systems generally behave strong non-ideality due to the associating interaction among acetic acid molecules. Theoretical study of the associating mechanism will provide guidance for the VLE property prediction, which is crucial for the designing on the separation process of the acetic acid systems. In this work, the association conformers and their distribution on acetic acid molecules in saturated gas and liquid phase were firstly studied. The proportions of the acetic acid monomer and multimers were obtained, which will contribute to the foundation for the vapor-liquid equilibrium simulations. The association mechanism on acetic acid molecules was then investigated by comparing among the structures and non-bonded interaction energies of different dimers. The structure of the cyclic dimer containing two OC-HO hydrogen bonds, may be found probably when acetic acid molecules approached. Electronic properties of different acetic acid dimers showed that the electrons around carbonyl oxygen atoms were deflected by the attraction of hydrogen atoms in the other molecule, which polarized the acetic acid molecules when the hydrogen bonds between acetic acid molecules were formed, providing theoretical basis for the polarized acetic acid molecular model. In this work, the molecular dynamics (MD) simulations and DFT calculations were conducted through the software GROMACS and Gaussian 09, respectively. For the MD simulations, the OPLS-AA force field was used as the atomic force field, with the cubic simulation cells constructed by Packmol program. For the DFT calculations, the M06-2X functional was employed for the optimization of the associating structures with the 6-311G** basis sets. Hydrogen bonding energies of dimers were corrected for the basis set superposition error (BSSE) and the deformation energies of monomers. Furthermore, the energy decomposition analysis was conducted at DFT/M06-2X/def-tzp level by the ADF software, and the wave function analysis was conducted by the Multiwfn software including the atom in molecule (AIM) topology analysis, the electronic potential analysis, and the electron density difference analysis.

Investigation of structural, dynamics, and dielectric properties of an aqueous potassium fluoride system at various concentrations by molecular dynamics simulations

Potassium-ion-based batteries have emerged as promising alternatives to traditional lithium-ion batteries for energy storage systems due to their affordability, wide accessibility and comparable chemical characteristics to lithium. This study employs molecular dynamics simulations to explore the physical phenomena of potassium fluoride in aqueous solutions. The interatomic interactions were defined using the OPLS-AA force field, while the SPC/E water model and ions were represented as charged Lennard?Jones particles. The simulations were conducted across concentrations ranging from 0.1 to 1.0 mol kg-1. The insights derived from this investigation provide valuable understanding into the behaviour of KF electrolytes and their potential utility in energy storage systems. A comprehensive comprehension of the impact of KF electrolyte concentration on structural, dynamic and dielectric properties is pivotal for the design and optimization of potassium-ion batteries, as well as other electrochemical devices leveraging KF-based electrolytes. This research significantly contributes to the ongoing endeavours aimed at developing efficient and economically viable energy storage solutions that transcend the confines of traditional lithium-ion batteries.

Enthalpies and entropies of hydration from Monte Carlo simulations.

The changes in free energy, enthalpy, and entropy for transfer of a solute from the gas phase into solution are the fundamental thermodynamic quantities that characterize the solvation process. Owing to the development of methods based on free-energy perturbation theory, computation of free energies of solvation has become routine in conjunction with Monte Carlo (MC) statistical mechanics and molecular dynamics (MD) simulations. Computation of the enthalpy change and by inference the entropy change is more challenging. Two methods are considered in this work corresponding to direct averaging for the solvent and solution and to computing the temperature derivative of the free energy in the van't Hoff approach. The application is for neutral organic solutes in TIP4P water using long MC simulations to improve precision. Definitive results are also provided for pure TIP4P water. While the uncertainty in computed free energies of hydration is ca. 0.05 kcal mol-1, it is ca. 0.4 kcal mol-1 for the enthalpy changes from either van't Hoff plots or the direct method with sampling for 5 billion MC configurations. Partial molar volumes of hydration are also computed by the direct method; they agree well with experimental data with an average deviation of 3 cm3 mol-1. In addition, the results permit breakdown of the errors in the free energy changes from the OPLS-AA force field into their enthalpic and entropic components. The excess hydrophobicity of organic solutes is enthalpic in origin.

Self-diffusion and shear viscosity of pure 1-alkanol unary system: molecular dynamics simulation and review of experimental data.

Self-diffusion coefficients and shear viscosity coefficients of pure 1-alkanol liquids from methanol to 1-hexanol were predicted using molecular dynamics (MD) simulations. These coefficients have been calculated using the Green-Kubo and Einstein methods at a range of temperatures of 200-330 K with increments of 10 K. Two force fields, TraPPE-UA and OPLS-AA were applied. The predicted results were compared to the experimental data, and the activation energies for self-diffusion and shear viscosity were calculated using the Arrhenius equation. The Stokes-Einstein equation was used to examine its capability in predicting the relationship between self-diffusion and shear viscosity, and the effective hydrodynamic radius was determined using both the experimental data and the results from MD simulations. The TraPPE-UA force field showed better results for the transport properties of methanol, while the OPLS-AA force field performed well for predicting shear viscosity but weakly for self-diffusion, particularly at low temperatures and for 1-alkanol with higher methylene numbers. Using the mean squared displacement method for self-diffusion was found to be more accurate than the Green-Kubo method, while the Green-Kubo method was slightly better for calculating shear viscosity. The Stokes-Einstein equation is valid for pure 1-alkanol liquids with temperature-dependent effective hydrodynamic radius.

Comparison of the Interaction and Structure of Lignin in Pure Systems and in Asphalt Media by Molecular Dynamics Simulations.

Lignin is a class of organic aromatic polymers contributing to the rigidity and strength of plants and has been proposed as a modifier to improve asphalt performance on road pavement. However, contradicting experimental results on the lignin miscibility in asphalt were found from different studies, and lignin has been found to self-assemble in different solutions. Thus, investigating the interaction and microstructure of lignin in asphalt media in molecular detail is necessary. Molecular dynamics (MD) simulations using both the LAMMPS program with the OPLS-aa force field and the NAMD program with the CHARMM force field have been conducted on pure lignin (including lignin monomer, dimer, and polymer with 17 and 31 units) and their mixtures with model asphalt molecules at different temperatures. Consistent results were observed from both programs and force fields in terms of density, hydrogen bonds, diffusion coefficient, radius of gyration, and radial distribution function. Glass transition was observed in the pure lignin systems based on density and diffusion coefficient calculations at different temperatures. Lignin can form intramolecular hydrogen bonds and intermolecular hydrogen bonds with other lignin and 1,7-dimethylnapthalene in the asphalt mixture, which has dependence on temperature and lignin chain length. Correlating the lignin size and chain length using the power-law relationship showed that lignin polymers in pure systems are in quasi-relaxed structures at different temperatures; lignin molecules stay in quasi-relaxed structures in asphalt mixtures at high temperatures but in collapsed structures at low temperatures. Implementing lignin monomer, dimer, and polymer into the model asphalt mixture can improve its density. Although lignin in different chain lengths aggregates in asphalt, lignin can modify the packing between different components in asphalt media at different temperatures. The work suggests that temperature can significantly influence the miscibility of lignin polymer in asphalt and that lignin can function as both a modifier and a resin in asphalt.

Molecular dynamics simulations and FTIR spectroscopy investigations on the hydration, transport, and dielectric properties of the NaF(aq) system at various concentrations

Sodium fluoride has many applications in various fields, such as dentistry, geochemistry, and rechargeable battery technology. Indeed, the determination of the hydration structures, dynamics, and dielectric properties of the NaF electrolyte in aqueous solutions is a crucial element in order to optimize its performance and predict its intervention mechanisms in industrial processes. In this study, we used the molecular dynamics approach to study the NaF(aq) binary system at different concentrations ranging from 0.1 to 1.0 mol.kg−1 using the OPLS-AA force field combined with the extended simple point load water model (SPC/E). Thus, the ions were represented as charged Lennard-Jones particles. The radial distribution functions (RDF) have been calculated to analyze the hydration shell structures of different ion pairs. The microstructures, self-diffusion coefficient, and dielectric constant were determined by molecular dynamics at various concentrations. Moreover, our FTIR spectroscopy data confirm the results obtained from molecular simulations.

A Comparative Analysis of Force Fields of Drug-Like Molecules: Omeprazole, Favipiravir and Amoxicillin by using Bonded and Non-Bonded Potential Distributions

This research aims to compare and analyze the force field parameters of drug-like compounds generated by the online servers CGenFF, SwissParam, and LigParGen. CHARMM27, OPLS-AA, and CHARMM36 force fields of amoxicillin, favipiravir, and omeprazole have been examined by bonded and non-bonded energy distributions of the potential function. The local terms of the force field function are bond length, bond angle, and dihedral potentials, while the non-local terms are van der Waals and pairwise electrostatic interactions. The parabolic curves of bond length potential appear to internally bend because of the larger values of force constants. The bond length and angle energies are lowest at reference values for each of the corresponding atom pairs in the molecules. The higher values of the bending force constant positively deflect the sinusoidal curve from the equilibrium position according to the dihedral angle energy distribution. The Lennard-Jones interaction exhibits the same repellent nature at closer ranges and the same attracting nature at farther ranges whereas a mildly differentiated interaction is seen at the intermediate state as a result of changing well-depth parameters. Depending on the size of the optimized partial charges, it is discovered in electrostatic interaction that the distributions of end-to-end lengths of identical atom pairs are more exponentially skewed from one another. The main cause of inconsistency in the results of MD simulations for the same protein-ligand system is due to the shifted topology and parameters.

Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides.

To assess the salt solubilities of six alkali halides in aqueous systems, we proposed a thermodynamic cycle and an efficient molecular modeling methodology. The Gibbs free energy changes for vaporization, dissociation, and dissolution were calculated using the experimental data of ionic thermodynamic properties obtained from the NBS tables. Additionally, the Marcus' and Tissandier's solvation free energy data for Li+, Na+, K+, Cl-, and Br- ions were compared with the conventional solvation free energies by substituting into our self-consistent thermodynamic cycle. Furthermore, Tissandier's absolute solvation free energy data were used as the training set to refit the Lennard-Jones parameters of OPLS-AA force field for ions. To predict salt solubilities, an assumption of a pseudo-solvent was proposed to characterize the coupling work of a solute with its environment from infinitely diluted to saturated solutions, indicating that the Gibbs energy change of solvation process is a function of ionic strength. Following the self-consistency of the cycle, the newly derived formulas were used to determine the salt solubilities by interpolating the intersection of Gibbs free energy of dissolution and the zero free energy line. The refined ion parameters can also predict the structure and thermodynamic properties of aqueous electrolyte solutions, such as densities, pair correlation functions, hydration numbers, mean activity coefficients, vapor pressures, and the radial dependences of the net charge at 298.15 K and 1 bar. Our method can be used to characterize the solid-liquid equilibria of ions or charged particles in aqueous systems. Furthermore, for highly concentrated strong electrolyte systems, it is essential to introduce accurate water models and polarizable force fields.

Evolution of the hydrogen-bonded network in methanol-water mixtures upon cooling

The hydrogen-bonded structure of methanol – water mixtures is investigated over the entire alcohol concentration range (from xMethanol = 0.1 to 1.0) at several temperatures, from 300 K down to the freezing point of the given mixture. Classical molecular dynamics simulations have been carried out, using the all-atom OPLS-AA force field for methanol and the TIP4P/2005 model for water molecules. Simulation trajectories (‘particle configurations’) obtained have been analyzed, in order to characterize the hydrogen-bonded network in the mixtures. The temperature and concentration dependence of the average hydrogen bond (H-bond) numbers between different types of molecules, the donor/acceptor roles of water and methanol molecules, and hydrogen bond number distributions have been revealed. The topology of the total system, as well as that of the water and methanol subsystems, has been investigated by calculating the cluster size distributions, the number of primitive rings, and ring size and ring type distributions. It has been found that upon cooling, the average number of H-bonded water molecules increases at every concentration and temperature investigated. As far as the connectivity of the hydrogen-bonded network is concerned, the percolation threshold has been shown to be above xM = 0.9 already at room temperature.

Viscosity evolution of water glycol in deep-sea environment at high pressure and low temperature

Using water glycol as a working medium is an important development direction for deep sea hydraulic drives, its viscosity is an important parameter for hydrodynamic lubrication and seriously influenced by the deep-sea environment. In this paper, Non-Equilibrium Molecular Dynamics method was used to predict the viscosity of water glycol with 0%, 25%, 50%, 75% and 100% water content in the deep-sea environment based on the OPLS-AA force field. The simulation results have an error of less than 10% compared to available experimental data. Detailed analysis of hydrogen bonding and water molecular structure is provided to explain why the viscosity of water decreases with increasing sea depth from 2000 m to 11,000 m. A fitted equation is proposed to calculate the viscosity of water glycol by analyzing the effect of temperature, pressure, and water content on the viscosity of water glycol in deep-sea environment.

A Computational Investigation on Rho-related GTP-binding Protein RhoB through Molecular Modeling and Molecular Dynamics Simulation Study

Background: An indispensable member of the Rho family, RhoB is an isoprenylated small GTPases that modulate the cellular cytoskeletal organization. While DNA gets damaged, it takes part in the neoplastic apoptotic mechanism. In this study, we evaluated the structure of Rho-related GTP-binding protein RhoB due to the unavailability of 3D structure in the protein data bank database. Results: The expected pI value of RhoB was 5.10 (acidic). The target–template alignment was computed using the GMQE value meanwhile 6hxu.1.A from Homo sapiens was selected as the template structure. The Swiss model was exploited to complete the model construction task. The structural compatibility and stability were revealed after a 100ns molecular dynamics simulation using GROMACA employing the OPLS-AA force field. Based on their fluctuating activity and their location between 100 and 110 and 140 and 150, PCA analysis discovered relevant residues. Conclusion: By providing an insight into the biophysical phenomenon of Rho-related GTP-binding protein RhoB inhibitors, this study will assist future investigations addressing the relationship between gene mutation and abnormalities produced by protein Rho-related GTP-binding protein RhoB in apoptotic events.

Solving Chemical Absorption Equilibria using Free Energy and Quantum Chemistry Calculations: Methodology, Limitations, and New Open-Source Software.

We developed an open-source chemical reaction equilibrium solver in Python (CASpy, https://github.com/omoultosEthTuDelft/CASpy) to compute the concentration of species in any reactive liquid-phase absorption system. We derived an expression for a mole fraction-based equilibrium constant as a function of excess chemical potential, standard ideal gas chemical potential, temperature, and volume. As a case study, we computed the CO2 absorption isotherm and speciation in a 23 wt % N-methyldiethanolamine (MDEA)/water solution at 313.15 K, and compared the results with available data from the literature. The results show that the computed CO2 isotherms and speciations are in excellent agreement with experimental data, demonstrating the accuracy and the precision of our solver. The binary absorptions of CO2 and H2S in 50 wt % MDEA/water solutions at 323.15 K were computed and compared with available data from the literature. The computed CO2 isotherms showed good agreement with other modeling studies from the literature while the computed H2S isotherms did not agree well with experimental data. The experimental equilibrium constants used as an input were not adjusted for H2S/CO2/MDEA/water systems and need to be adjusted for this system. Using free energy calculations with two different force fields (GAFF and OPLS-AA) and quantum chemistry calculations, we computed the equilibrium constant (K) of the protonated MDEA dissociation reaction. Despite the good agreement of the OPLS-AA force field (ln[K] = -24.91) with the experiments (ln[K] = -23.04), the computed CO2 pressures were significantly underestimated. We systematically investigated the limitations of computing CO2 absorption isotherms using free energy and quantum chemistry calculations and showed that the computed values of μiex are very sensitive to the point charges used in the simulations, which limits the predictive power of this method.

Implementation of solvation free energy framework to predict the vapor-liquid equilibrium behaviors for the water–hydrazine and ethanol–water azeotropic systems

An azeotrope is a special vapor–liquid equilibrium in which the composition of two coexisting phases is equal. Due to its strong deviation from the Raoult’s law, it poses a great challenge to the simple distillation process and the understanding of fundamental molecular mechanisms. In this study, we calculated the coupling work of solute species based on the solvation free energy framework to predict the vapor–liquid equilibrium (VLE) behavior and azeotrope point of mixtures. Two prototypical binary mixtures: the negative water-hydrazine and positive ethanol–water azeotrope were investigated using molecular dynamic simulations and the OPLS-AA force field. All the simulations were performed with the TIP3P water model, because it has a good capability to represent enthalpy properties. The molecular models of hydrazine and ethanol have been re-parameterized to reproduce the density and solvation free energy properties of pure substance. Using free energy perturbation simulation, the coupling work of each substance was studied against liquid mixture composition. The coupling work is also translated to x-y or T-x-y phase diagrams for comparison with available experimental data. The tendency of two azeotropes is that the more volatile components experience the less pronounced coupling work changes in their local environments as the azeotropes form. The coupling work is also divided into van der Waals and Coulomb contributions to investigate the interactions between particles. Lithium Chloride was investigated as an entrainer for the separation of ethanol–water system. The coupling works of binary mixtures were calculated at different salt concentrations to study the salt-out effect for the violate species and salt-in effect for the less violate species. By adding LiCl, it was found that the azeotrope point was efficiently eliminated at a salinity 0.05 mol fraction at 1 bar. All the results show that the simulation methodology of coupling work are reliable for calculating the liquid–vapor equilibrium behavior. This method is also valuable for selecting and designing new extractants for distillation process.

Extension of the TraPPE Force Field for Battery Electrolyte Solvents.

Optimizing electrolyte formulations is key to improving performance of Li-/Na-ion batteries, where transport properties (diffusion coefficient, viscosity) and permittivity need to be predicted as functions of temperature, salt concentration and solvent composition. More efficient and reliable simulation models are urgently needed, owing to the high cost of experimental methods and the lack of united-atom molecular dynamics force fields validated for electrolyte solvents. Here the computationally efficient TraPPE united-atom force field is extended to be compatible with carbonate solvents, optimizing the charges and dihedral potential. Computing the properties of electrolyte solvents, ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), diethyl carbonate (DEC), and dimethoxyethane (DME), we observe that the average absolute errors in the density, self-diffusion coefficient, permittivity, viscosity, and surface tension are approximately 15% of the corresponding experimental values. Results compare favorably to all-atom CHARMM and OPLS-AA force fields, offering computational performance improvement of at least 80%. We further use TraPPE to predict the structure and properties of LiPF6 salt in these solvents and their mixtures. EC and PC form complete solvation shells around Li+ ions, while the salt in DMC forms chain-like structures. In the poorest solvent, DME, LiPF6 forms globular clusters despite DME's higher permittivity than DMC.

The effect of various partial atomic charges on the bulk and liquid/vacuum interface properties of iodobenzene derivatives at their melting points

In this work, various atomic charge schemes including natural bond orbital (NBO), electrostatic potential based (CHELPG), and σ-hole model charges were applied in the OPLS-AA force field to investigate their effects on the thermophysical and structural properties of iodobenzene and its derivatives. Molecular dynamics simulations presented in this work show that the studied structural and thermophysical properties are in good agreement with experiments when the CHELPG charge was coupled with the OPLS-AA force field. Also, the arrangement of iodobenzene derivatives in the liquid phase was investigated via combined radial/angular distribution functions (CDFs) analyses and halogen-bonding theory. The most probable orientation of iodobenzene derivatives at the liquid/vacuum interface was assigned by atom density profile and bivariate orientational distribution maps. For all studied iodobenzene derivatives, benzene rings are oriented such that the iodine atoms tend toward the vacuum phase.