We present a new molecular simulation code, Brick-CFCMC, for performing Monte Carlo simulations using state-of-the-art simulation techniques. The Continuous Fractional Component (CFC) method is implemented for simulations in the NVT/NPT ensembles, the Gibbs Ensemble, the Grand-Canonical Ensemble, and the Reaction Ensemble. Molecule transfers are facilitated by the use of fractional molecules which significantly improve the efficiency of the simulations. With the CFC method, one can obtain phase equilibria and properties such as chemical potentials and partial molar enthalpies/volumes directly from a single simulation. It is possible to combine trial moves from different ensembles. This enables simulations of phase equilibria in a system where also a chemical reaction takes place. We demonstrate the applicability of our software by investigating the esterification of methanol with acetic acid in a two-phase system.@en

Water transport is one of the major factors responsible for moisture damage in asphalt pavements. To study the thermodynamics and kinetics of water transport in bitumen and to uncover microscale mechanisms of moisture-induced damage, molecular dynamics simulations were performed for up to 600 ns for water–bitumen systems with realistic water contents that varied from 0 to 1.76 wt%. Hydrogen bonding interactions and clustering of water molecules at various combinations of temperature and water content were investigated, and their effects on the self-diffusion coefficient of water and bitumen properties are computed and discussed. It is shown that the saturated water concentration in bitumen is small, especially at low temperatures, and additional water molecules tend to form large water clusters via hydrogen bonding, indicating micro-phase separation of the water and bitumen phases inside the simulation box. Hydrogen bonding and water clustering play a crucial role on the magnitude of the self-diffusion coefficient of water. Physical properties of bitumen that include viscosity and cohesive energy are affected by water. The presence of large water clusters is indicative of how degradation in cohesion is observed on the microscale.@en

Deep eutectic solvents (DESs) are considered as green alternatives to room temperature ionic liquids (RTILs), due to their lower-cost synthesis and more environmentally friendly nature. In this work, Monte Carlo (MC) simulations have been used to compute the solubilities of CO2, H2S, CH4, CO, H2, and N2 in choline chloride urea (ChClU) and choline chloride ethylene glycol (ChClEg) DESs. Due to the strong intermolecular interactions of DESs, leading to high viscosities, MC simulations present significant challenges with respect to system equilibration and solute molecule insertions. The Continuous Fractional Component Monte Carlo (CFCMC) method has been used with our open-source code, Brick-CFCMC, to improve molecule insertions and equilibration of the system, and directly compute the excess chemical potential and solubility (in terms of the Henry constant) of the gas molecules in the DESs. Pure DES properties, such as density and radial distribution functions (RDFs), were well reproduced by MC simulations. The solubilities of gases were, however, underestimated by the CFCMC simulations compared to available experimental data from literature. The order of solubilities of the different gases in ChClU at 328 K was obtained as H2S > CO2 > CH4 > H2 > CO > N2, which reasonably agrees with experimental data from literature. The OPLS force field resulted in larger average Henry constants in ChClEg, compared to the GAFF force field, implying the better suitability of the GAFF force field for the calculations. Smaller ionic charge scaling factors were shown to increase the solubilities of the gases in the DESs, but result in lower densities. The differences between the computed Henry constants from MC simulations and experimental data from literature may be caused by the unsuitability of the used force field parameters of the DESs in combination with those of the solute gases. Nonetheless, experimental data from literature is scarce (except for CO2) and in some cases contradictory, which makes the comparison with the computational results difficult.@en

The solubility parameter (SP) of a solvent is a key property that measures the polarity and quantifies the ‘like-dissolves-like’ principle, which is an important rule in chemistry for screening solvents for separation processes. It is challenging to experimentally obtain solubility parameters of non-volatile solvents like ionic liquids (ILs), deep eutectic solvents (DESs), and polymers. Here, Molecular Dynamics (MD) simulations have been used to compute the Hildebrand and Hansen solubility parameters of DESs, which are green solvents with potential applications in many different fields. The results from MD simulations are compared with limited available experimental data and commonly used SP correlations for non-volatile solvents. Very limited information is available in literature for the vapor phase composition of DESs. Solubility parameters are computed based on the vaporization of hydrogen bond donor (HBD) and hydrogen bond acceptor (HBA) components of the DESs as well as clusters, consisting of HBD and HBA components. The relatively large SPs computed from MD indicate that the investigated choline chloride-based DESs are polar solvents. The values of SPs are not significantly affected by temperature. A comparison of vaporization enthalpies of HBD, HBA and clusters from the DES mixture suggests that it is more likely for HBD molecules to vaporize from the DES mixture and dominate the vapor phase.@en

Hydrophobic deep eutectic solvents (DESs) have recently gained much attention as water-immiscible solvents for a wide range of applications. However, very few studies exist in which the hydrophobicity of these DESs is quantified. In this work, the interfacial properties of hydrophobic DESs with water were computed at various temperatures using molecular dynamics simulations. The considered DESs were tetrabutylammonium chloride-decanoic acid (TBAC-dec) with a molar ratio of 1:2, thymol-decanoic acid (Thy-dec) with a molar ratio of 1:2, and dl-menthol-decanoic acid (Men-dec) with a molar ratio of 2:1. The following properties were investigated in detail: interfacial tensions, water-in-DES solubilities (and salt-in-water solubilities for TBAC-dec/water), density profiles, and the number densities of hydrogen bonds. Different ionic charge scaling factors were used for TBAC-dec. Thy-dec and Men-dec showed a high level of hydrophobicity with negligible computed water-in-DES solubilities. For charge scaling factors of 0.7 and 1 for the thymol and decanoic acid components of Thy-dec, the computed interfacial tensions of the DESs are in the following order: TBAC-dec (ca. 4 mN m-1) < Thy-dec (20 mN m-1) < Men-dec (26 mN m-1). The two sets of charge scaling factors for Thy-dec did not lead to different density profiles but resulted in considerable differences in the DES/water interfacial tensions due to different numbers of decanoic acid-water hydrogen bonds at the interfaces. Large peaks were observed for the density profiles of (the hydroxyl oxygen of) decanoic acid at the interfaces of all DES/water mixtures, indicating a preferential alignment of the oxygen atoms of decanoic acid toward the aqueous phase. @en

With the emergence of hydrophobic deep eutectic solvents (DESs), the scope of applications of DESs has been expanded to include situations in which miscibility with water is undesirable. Whereas most studies have focused on the applications of hydrophobic DESs from a practical standpoint, few theoretical works exist that investigate the structural and thermodynamic properties at the nanoscale. In this study, Molecular Dynamics (MD) simulations have been performed to model DESs composed of tetraalkylammonium chloride hydrogen bond acceptor and decanoic acid hydrogen bond donor (HBD) at a molar ratio of 1:2, with three different cation chain lengths (4, 7, and 8). After fine-tuning force field parameters, densities, viscosities, self-diffusivities, and ionic conductivities of the DESs were computed over a wide temperature range. The liquid structure was examined using radial distribution functions (RDFs) and hydrogen bond analysis. The MD simulations reproduced the experimental density and viscosity data from the literature reasonably well and were used to predict diffusivities and ionic conductivities, for which experimental data are scarce or unavailable. It was found that although an increase in the cation chain length considerably affected the density and transport properties of the DESs (i.e., yielding smaller densities and slower dynamics), no significant influence was observed on the RDFs and the hydrogen bonds. The self-diffusivities showed the following order for the mobility of the various components: HBD > anion > cation. Strong hydrogen bonds between the hydroxyl and carbonyl groups of decanoic acid and between the hydroxyl group of decanoic acid and chloride were observed to dominate the intermolecular interactions.@en

Despite the widespread acknowledgment that deep eutectic solvents (DESs) have negligible vapor pressures, very few studies in which the vapor pressures of these solvents are measured or computed are available. Similarly, the vapor phase composition is known for only a few DESs. In this study, for the first time, the vapor pressures and vapor phase compositions of choline chloride urea (ChClU) and choline chloride ethylene glycol (ChClEg) DESs are computed using Monte Carlo simulations. The partial pressures of the DES components were obtained from liquid and vapor phase excess Gibbs energies, computed using thermodynamic integration. The enthalpies of vaporization were computed from the obtained vapor pressures, and the results were in reasonable agreement with the few available experimental data in the literature. It was found that the vapor phases of both DESs were dominated by the most volatile component (hydrogen bond donor, HBD, i.e., urea or ethylene glycol), i.e., 100% HBD in ChClEg and 88%-93% HBD in ChClU. Higher vapor pressures were observed for ChClEg compared to ChClU due to the higher volatility of ethylene glycol compared to urea. The influence of the liquid composition of the DESs on the computed properties was studied by considering different mole fractions (i.e., 0.6, 0.67, and 0.75) of the HBD. Except for the partial pressure of ethylene glycol in ChClEg, all the computed partial pressures and enthalpies of vaporization showed insensitivity toward the liquid composition. The activity coefficient of ethylene glycol in ChClEg was computed at different liquid phase mole fractions, showing negative deviations from Raoult’s law.@en