Prediction of Transport Properties by Molecular Simulation: Methanol and Ethanol and Their Mixture

Experimental and computational study of binary mixture ethanolamine and 2-amino-2-methyl-1-propanol

The present work involves experimental and computational investigations into the density of pure and mixed states of ethanolamine (ET) and 2-amino-2-methyl-1-propanol (AMP) under a pressure of 1 atm and temperatures ranging from 293.15 K to 333.15 K The density data were used to derive the excess molar volume, thermal expansion coefficient, and isothermal coefficient of pressure excess molar enthalpy. The Redlich-Kister equation was employed to calculate the excess molar and its accompanying coefficients. In the gas phase, density functional theory (DFT) was utilized to explore the most stable structures of ET … ET, AMP … AMP, and the ET … AMP mixture. Molecular dynamics simulation (MD) was used to calculate the structural properties of these mixtures in the liquid phase. Radial distribution function (RDFs) combined distribution function (CDF) and spatial distribution function (SDF) in different mole fractions calculated in the liquid phase. The intramolecular and intermolecular interactions of ethanolamine and AMP were obtained using the radial distribution function in different molar fractions. It was found that the ethanolamine molecule has a greater tendency to form intramolecular hydrogen bonds, while the AMP molecule has a greater tendency to form intermolecular hydrogen bonds.

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.

Influence of Ethanol Parametrization on Diffusion Coefficients Using OPLS-AA Force Field

Molecular dynamics simulations employing the all-atom optimized potential for liquid simulations (OPLS-AA) force field were performed for determining self-diffusion coefficients (D11) of ethanol and tracer diffusion coefficients (D12) of solutes in ethanol at several temperature and pressure conditions. For simulations employing the original OPLS-AA diameter of ethanol's oxygen atom (σOH), calculated and experimental diffusivities of protic solutes differed by more than 25%. To correct this behavior, the σOH was reoptimized using the experimental D12 of quercetin and of gallic acid in liquid ethanol as benchmarks. A substantial improvement of the calculated diffusivities was found by changing σOH from its original value (0.312 nm) to 0.306 nm, with average absolute relative deviations (AARD) of 3.71% and 4.59% for quercetin and gallic acid, respectively. The new σOH value was further tested by computing D12 of ibuprofen and butan-1-ol in liquid ethanol with AARDs of 1.55% and 4.81%, respectively. A significant improvement was also obtained for the D11 of ethanol with AARD = 3.51%. It was also demonstrated that in the case of diffusion coefficients of non-polar solutes in ethanol, the original σOH=0.312 nm should be used for better agreement with experiment. If equilibrium properties such as enthalpy of vaporization and density are estimated, the original diameter should be once again adopted.

A New Force Field for OH- for Computing Thermodynamic and Transport Properties of H2 and O2 in Aqueous NaOH and KOH Solutions

The thermophysical properties of aqueous electrolyte solutions are of interest for applications such as water electrolyzers and fuel cells. Molecular dynamics (MD) and continuous fractional component Monte Carlo (CFCMC) simulations are used to calculate densities, transport properties (i.e., self-diffusivities and dynamic viscosities), and solubilities of H2 and O2 in aqueous sodium and potassium hydroxide (NaOH and KOH) solutions for a wide electrolyte concentration range (0-8 mol/kg). Simulations are carried out for a temperature and pressure range of 298-353 K and 1-100 bar, respectively. The TIP4P/2005 water model is used in combination with a newly parametrized OH- force field for NaOH and KOH. The computed dynamic viscosities at 298 K for NaOH and KOH solutions are within 5% from the reported experimental data up to an electrolyte concentration of 6 mol/kg. For most of the thermodynamic conditions (especially at high concentrations, pressures, and temperatures) experimental data are largely lacking. We present an extensive collection of new data and engineering equations for H2 and O2 self-diffusivities and solubilities in NaOH and KOH solutions, which can be used for process design and optimization of efficient alkaline electrolyzers and fuel cells.

Non-monotonic composition dependence of the breakdown of Stokes-Einstein relation for water in aqueous solutions of ethanol and 1-propanol: explanation using translational jump-diffusion approach

A series of experimental and simulation studies examined the validity of the Stokes-Einstein relationship (SER) of water in binary water/alcohol mixtures of different mixture compositions. These studies revealed a strong non-monotonic composition dependence of the SER with maxima at the specific alcohol mole fraction where the non-idealities of the thermodynamic and transport properties are observed. The translational jump-diffusion (TJD) approach elucidated the breakdown of the SER in pure supercooled water as caused by the jump translation of molecules. The breakdown of SER in the supercooled water/methanol binary mixture was successfully explained using the same TJD approach. To further generalize the picture, here we focus on the non-monotonic composition dependence of SER breakdown of water in two water/alcohol mixtures (water/ethanol and water/propanol) for a broad temperature range. In agreement with previous studies, maximum breakdown of SER is observed for the mixture with alcohol mole fraction x = 0.2. Diffusion of the water molecules at the maximum SER breakdown point is largely contributed by jump-diffusion. The residual-diffusion, obtained by subtracting the jump-diffusion from the total diffusion, approximately follows the SER for different compositions and temperatures. We also performed hydrogen (H-)bond dynamics and observed that the contribution of jump-diffusion is proportional to the total free energy of activation of breaking all H-bonds that exist around a molecule. This study, therefore, suggests that the more a molecule is trapped by H-bonding, the more likely it is to diffuse through the jump-diffusion mechanism, eventually leading to an increasing degree of SER breakdown.

Thermal conductivity temperature dependence of water confined in nanoporous silicon

Recently, it has been shown that high density nanoconfined water was the reason of the important enhancement of the effective thermal conductivity up to a factor of 50% of a nanoporous silicon filled with water. In this work, using molecular dynamics simulations, we further investigate the role of the temperatureT(from 285 to 360 K) on the thermal conductivity enhancement of nanohybrid porous silicon and water system. Furthermore, by studying and analysing several structural and dynamical parameters of the nanoconfined water, we give physical insights of the observed phenomena. Upon increasing the temperature of the system, the thermal conductivity of the hybrid system increases reaching a maximum forT= 300 K. With this article, we prove the existence of new heat flux channels between a solid matrix and a nanoconfined liquid, with clear signatures both in the radial distribution function, mean square displacements, water molecules orientation, hydrogen bond networks and phonon density of states.

Molecular dynamic simulation studies of adsorption and diffusion behaviors of methanol and ethanol through ZSM-5 zeolite

Due to the importance of synthesis gas's entire conversion to methanol, the separation of methanol from unconverted synthesis gas is an industrial challenge. In this work, the influence of temperature, guest molecules concentrations (methanol and ethanol), and acid site density (Si/Al) of zeolites on the diffusion of methanol and ethanol, pure and binary mixture (80% methanol and 20% ethanol) in silicalite-1 and HZSM-5 (Si/Al = 47 and 23) were studied by using of the COMPASS force-field molecular dynamics method. Also, the adsorption of pure methanol and ethanol and binary mixture through these zeolites has been studied by using the Grand Canonical Monte Carlo (GCMC) method. The calculated adsorption rate and isosteric heat of adsorption for ethanol are lower and higher than methanol, respectively. The results of the binary mixture show that HZSM-5 (Si/Al = 23) has the lowest adsorption selectivity and most diffusion selectivity. The calculated diffusion coefficients of methanol and ethanol guest molecules decreased with rising guest molecule concentration and Si/Al-ratios. The effect of both agents was investigated by analysis of mean square displacement (MSD) and RDF diagram.

Diffusion in multicomponent aqueous alcoholic mixtures

The Fick diffusion coefficient matrix of the highly associating quaternary mixture water + methanol + ethanol + 2-propanol as well as its ternary and binary subsystems is analyzed with molecular dynamics simulation techniques. Three of the ternary subsystems are studied in this sense for the first time. The predictive capability of the employed force fields, which were sampled with the Green–Kubo formalism and Kirkwood–Buff integration, is confirmed by comparison with experimental literature data on vapor-liquid equilibrium, shear viscosity and Fick diffusion coefficient, wherever possible. A thorough analysis of the finite size effects on the simulative calculation of diffusion coefficients of multicomponent systems is carried out. Moreover, the dependence of the Fick diffusion coefficient matrix on the velocity reference frame and component order is analyzed. Their influence is found to be less significant for the main matrix elements, reaching a maximum variation of 19%. The large differences found for the cross elements upon variation of the reference frame hinder a straightforward interpretation of the Fick diffusion coefficient matrix with respect to the presence of diffusive coupling effects.

Molecular modelling techniques for predicting liquid-liquid interfacial properties of methanol plus alkane (n-hexane, n-heptane, n-octane) mixtures

In this work, the liquid-liquid interfacial properties of methanol plus n-alkane (n-hexane, n-heptane, n-octane) mixtures are investigated at atmospheric pressure by two complementary molecular modelling techniques; namely, molecular dynamic simulations (MD) and density gradient theory (DGT) coupled with the PC-SAFT (perturbed-chain statistical associating fluid theory) equation of state. Furthermore, two molecular models of methanol are used, which are based on a non-polarisable three site approach. On the one hand, is the original (flexible) TraPPE-UA model force field. On the other hand, is the rigid approximation denoted as OPLS/2016. In both cases, n-alkanes are modelled using the TraPPE-UA model. Simulations are performed using the direct coexistence technique in the ensemble. Special attention is paid to the comparison between the estimations obtained from different methanol models, the available experimental data and theoretical calculations. In all cases, the rigid model is capable of predicting the experimental phase equilibrium and interfacial properties accurately. Unsurprisingly, the methanol-rich density and interfacial tension are overestimated using the TraPPE model combined with Lorentz-Berthelot mixing rules for predicting the mixture behaviour. Accurate comparison between MD and DGT plus PC-SAFT requires consideration of the cross-interactions between individual influence parameters and fitting the βij values. This latter aspect is particularly important because it allows the exploitation of the link between the EOS model and the direct molecular simulation of the corresponding fluid. At the same time, it was demonstrated that the key property defining the interfacial tension value is the absolute concentration of methanol in the methanol-rich phase. This behaviour indicates that there are more hydrogens bonded with each other, and they interact favourably with an increasing number of carbon atoms in the alkane.

Frequency-domain hot-wire sensor and 3D model for thermal conductivity measurements of reactive and corrosive materials at high temperatures

High temperature solids and liquids are becoming increasingly important in next-generation energy and manufacturing systems that seek higher efficiencies and lower emissions. Accurate measurements of thermal conductivity at high temperatures are required for the modeling and design of these systems, but commonly employed time-domain measurements can have errors from convection, corrosion, and ambient temperature fluctuations. Here, we describe the development of a frequency-domain hot-wire technique capable of accurately measuring the thermal conductivity of solid and molten compounds from room temperature up to 800 °C. By operating in the frequency-domain, we can lock into the harmonic thermal response of the material and reject the influence of ambient temperature fluctuations, and we can keep the probed volume below 1 µl to minimize convection. The design of the microfabricated hot-wire sensor, electrical systems, and insulating wire coating to protect against corrosion is covered in detail. Furthermore, we discuss the development of a full three-dimensional multilayer thermal model that accounts for both radial conduction into the sample and axial conduction along the wire and the effect of wire coatings. The 3D, multilayer model facilitates the measurement of small sample volumes important for material development. A sensitivity analysis and an error propagation calculation of the frequency-domain thermal model are performed to demonstrate what factors are most important for thermal conductivity measurements. Finally, we show thermal conductivity measurements including model data fitting on gas (argon), solid (sulfur), and molten substances over a range of temperatures.

A Force Field for Poly(oxymethylene) Dimethyl Ethers (OMEn)

A united atom force field for the homologous series of the poly(oxymethylene) dimethyl ethers (OMEn), H₃C-O-(CH₂O)_n-CH₃, is presented. OMEn are oxygenates and promising new synthetic fuels and solvents. The molecular geometry of the OMEn, the internal degrees of freedom, and their electrostatic properties were obtained from quantum mechanical calculations. To model repulsion and dispersion, Lennard-Jones parameters were fitted to the experimental liquid densities and vapor pressures of pure OMEn (n = 1-4). The critical properties of OMEn (n = 1-4) were determined from the simulation data. Additionally, the shear viscosity of pure liquid OMEn is evaluated and compared with literature data. Finally, the solubility of CO₂ in OME2, OME3, and OME4 is predicted using a literature model for CO₂ and the Lorentz-Berthelot combining rules. The results agree well with experimental data from the literature.

Molecular dynamics simulations of nano-confined methanol and methanol-water mixtures between infinite graphite plates: Structure and dynamics

Molecular dynamics simulations are used to investigate microscopic structures and dynamics of methanol and methanol-water binary mixture films confined between hydrophobic infinite parallel graphite plate slits with widths, H, in the range of 7-20 Å at 300 K. The initial geometric densities of the liquids were chosen to be the same as bulk methanol at the same temperature. For the two narrowest slit widths, two smaller initial densities were also considered. For the nano-confined system with H = 7 Å and high pressure, a solid-like hexagonal arrangement of methanol molecules arranged perpendicular to the plates is observed which reflects the closest packing of the molecules and partially mirrors the structure of the underlying graphite structure. At lower pressures and for larger slit widths, in the contact layer, the methanol molecules prefer having the C-O bond oriented parallel to the walls. Layered structures of methanol parallel to the wall were observed, with contact layers and additional numbers of central layers depending on the particular slit width. For methanol-water mixtures, simulations of solutions with different composition were performed between infinite graphite slits with H = 10 and 20 Å at 300 K. For the nanoslit with H = 10 Å, in the solution mixtures, three layers of molecules form, but for all mole fractions of methanol, methanol molecules are excluded from the central fluid layer. In the nanopore with H = 20 Å, more than three fluid layers are formed and methanol concentrations are enhanced near the confining plates walls compared to the average solution stoichiometry. The self-diffusion coefficients of methanol and water molecules in the solution show strong dependence on the solution concentration. The solution mole fractions with minimal diffusivity are the same in confined and non-confined bulk methanol-water mixtures.

A Simple Two Dimensional Model of Methanol

Methanol is the simplest alcohol and possible energy carrier because it is easier to store than hydrogen and burns cleaner than fossil fuels. It is a colorless liquid, completely miscible with water and organic solvents and is very hygroscopic. Here, simple two-dimensional models of methanol, based on Mercedes-Benz (MB) model of water, are examined by Monte Carlo simulations. Methanol particles are modeled as dimers formed by an apolar Lennard-Jones disk, mimicking the methyl group, and a sphere with two hydrogen bonding arms for the hydroxyl group. The used models are the one proposed by Hribar-Lee and Dill (Acta Chimica Slovenica, 53:257, 2006.) with the overlapping discs and a new model with tangentially fused dimers. The comparison was done between the models, in connection to the MB water, as well as with experimental results and with new simulations done for 3D models of methanol. Both 2D models show similar trends in structuring and thermodynamics. The difference is the most pronounced at lower temperatures, where the smaller model exhibits spontaneous crystallization, while the larger model shows metastable states. The 2D structural organization represents well the clustering tendency observed in 3D models, as well as in experiments. The models qualitatively agree with the bulk methanol thermodynamic properties like density and isothermal compressibility, however, heat capacity at the constant pressure shows trend more similar to the water behavior. This work on the smallest amphiphilic organic solute provides a simple testing ground to study the competition between polar and non-polar effects within the molecule and physical properties.

1-Butanol as a Solvent for Efficient Extraction of Polar Compounds from Aqueous Medium: Theoretical and Practical Aspects

The extraction of polar molecules from aqueous solution is a challenging task in organic synthesis. 1-Butanol has been used sporadically as an eluent for polar molecules, but it is unclear which molecular features drive its efficiency. Here, we employ free energy simulations to study the partitioning of 15 solutes between water and 1-butanol. The simulations demonstrate that the high affinity of polar molecules to the wet 1-butanol phase is associated with its nanostructure. Small inverse micelles of water are able to accommodate polar solutes and locally mimic an aqueous environment. We verify the simulations based on partition coefficients between water and 1-octanol, and include a blind prediction of the water/1-butanol partition coefficient of cyclohexane-1,2-diol. The calculations are in excellent agreement with experiment, reaching root-mean-square deviations below 0.7 kcal/mol. Actual extractions of cyclohexane-1,2-diol from buffer solutions that mimic cell lysates and suspensions in biocatalytic reactions further exemplify our findings. The yields highlight that extractions with 1-butanol can be significantly more efficient than the conventional protocol based on ethyl acetate.

Dynamical Ordering of Hydrogen Molecules Induced by Heat Flux

Achieving a direct nonequilibrium simulation for hydrogen systems has been quite challenging because nuclear quantum effects (NQEs) have to be taken into account. We directly simulated nonequilibrium hydrogen molecules under a temperature gradient with the recently developed nonempirical molecular dynamics method, which describes nonspherical hydrogen molecules with the NQEs. We found dynamical ordering purely induced by heat flux, which should be distinguished from static ordering like orientational alignment, as decelerated translational motions and enhanced intensity of H-H vibrational power spectra despite the little structural ordering. This dynamical ordering, which was enhanced with stronger heat flux while independent of system size, can be regarded as self-solidification of hydrogen molecules for their efficient heat conduction.

Carbon membranes for efficient water-ethanol separation

We demonstrate, on the basis of molecular dynamics simulations, the possibility of an efficient water-ethanol separation using nanoporous carbon membranes, namely, carbon nanotube membranes, nanoporous graphene sheets, and multilayer graphene membranes. While these carbon membranes are in general permeable to both pure liquids, they exhibit a counter-intuitive "self-semi-permeability" to water in the presence of water-ethanol mixtures. This originates in a preferred ethanol adsorption in nanoconfinement that prevents water molecules from entering the carbon nanopores. An osmotic pressure is accordingly expressed across the carbon membranes for the water-ethanol mixture, which agrees with the classic van't Hoff type expression. This suggests a robust and versatile membrane-based separation, built on a pressure-driven reverse-osmosis process across these carbon-based membranes. In particular, the recent development of large-scale "graphene-oxide" like membranes then opens an avenue for a versatile and efficient ethanol dehydration using this separation process, with possible application for bio-ethanol fabrication.

Temperature of maximum density and excess thermodynamics of aqueous mixtures of methanol

In this work, we present a study of representative excess thermodynamic properties of aqueous mixtures of methanol over the complete concentration range, based on extensive computer simulation calculations. In addition to test various existing united atom model potentials, we have developed a new force-field which accurately reproduces the excess thermodynamics of this system. Moreover, we have paid particular attention to the behavior of the temperature of maximum density (TMD) in dilute methanol mixtures. The presence of a temperature of maximum density is one of the essential anomalies exhibited by water. This anomalous behavior is modified in a non-monotonous fashion by the presence of fully miscible solutes that partly disrupt the hydrogen bond network of water, such as methanol (and other short chain alcohols). In order to obtain a better insight into the phenomenology of the changes in the TMD of water induced by small amounts of methanol, we have performed a new series of experimental measurements and computer simulations using various force fields. We observe that none of the force-fields tested capture the non-monotonous concentration dependence of the TMD for highly diluted methanol solutions.

A new intermolecular potential for simulations of methanol: The OPLS/2016 model

In this work, a new rigid-nonpolarizable model of methanol is proposed. The model has three sites, located at the same positions as those used in the OPLS model previously proposed by Jorgensen [J. Phys. Chem. 90, 1276 (1986)]. However, partial charges and the values of the Lennard-Jones parameters were modified by fitting to an adequately selected set of target properties including solid-fluid experimental data. The new model was denoted as OPLS/2016. The overall performance of this model was evaluated and compared to that obtained with other popular models of methanol using a similar test to that recently proposed for water models. In the test, a certain numerical score is given to each model. It was found that the OPLS/2016 obtained the highest score (7.4 of a maximum of 10) followed by L1 (6.6), L2 (6.4), OPLS (5.8), and H1 (3.5) models. The improvement of OPLS/2016 with respect to L1 and L2 is mainly due to an improvement in the description of fluid-solid equilibria (the melting point is only 14 K higher than the experimental value). In addition, it was found that no methanol model was able to reproduce the static dielectric constant and the isobaric heat capacity, whereas the better global performance was found for models that reproduce the vaporization enthalpy once the so-called polarization term is included. Similar conclusions were suggested previously in the analysis of water models and are confirmed here for methanol.

NMR relaxation induced by iron oxide particles: testing theoretical models

Superparamagnetic iron oxide particles find their main application as contrast agents for cellular and molecular magnetic resonance imaging. The contrast they bring is due to the shortening of the transverse relaxation time T 2 of water protons. In order to understand their influence on proton relaxation, different theoretical relaxation models have been developed, each of them presenting a certain validity domain, which depends on the particle characteristics and proton dynamics. The validation of these models is crucial since they allow for predicting the ideal particle characteristics for obtaining the best contrast but also because the fitting of T 1 experimental data by the theory constitutes an interesting tool for the characterization of the nanoparticles. In this work, T 2 of suspensions of iron oxide particles in different solvents and at different temperatures, corresponding to different proton diffusion properties, were measured and were compared to the three main theoretical models (the motional averaging regime, the static dephasing regime, and the partial refocusing model) with good qualitative agreement. However, a real quantitative agreement was not observed, probably because of the complexity of these nanoparticulate systems. The Roch theory, developed in the motional averaging regime (MAR), was also successfully used to fit T 1 nuclear magnetic relaxation dispersion (NMRD) profiles, even outside the MAR validity range, and provided a good estimate of the particle size. On the other hand, the simultaneous fitting of T 1 and T 2 NMRD profiles by the theory was impossible, and this occurrence constitutes a clear limitation of the Roch model. Finally, the theory was shown to satisfactorily fit the deuterium T 1 NMRD profile of superparamagnetic particle suspensions in heavy water.

Effects of molecular structure on microscopic heat transport in chain polymer liquids

In this paper, we discuss the molecular mechanism of the heat conduction in a liquid, based on nonequilibrium molecular dynamics simulations of a systematic series of linear- and branched alkane liquids, as a continuation of our previous study on linear alkane [T. Ohara et al., J. Chem. Phys. 135, 034507 (2011)]. The thermal conductivities for these alkanes in a saturated liquid state at the same reduced temperature (0.7Tc) obtained from the simulations are compared in relation to the structural difference of the liquids. In order to connect the thermal energy transport characteristics with molecular structures, we introduce the new concept of the interatomic path of heat transfer (atomistic heat path, AHP), which is defined for each type of inter- and intramolecular interaction. It is found that the efficiency of intermolecular AHP is sensitive to the structure of the first neighbor shell, whereas that of intramolecular AHP is similar for different alkane species. The dependence of thermal conductivity on different lengths of the main and side chain can be understood from the natures of these inter- and intramolecular AHPs.

Properties of a soft-core model of methanol: an integral equation theory and computer simulation study

Thermodynamic and structural properties of a coarse-grained model of methanol are examined by Monte Carlo simulations and reference interaction site model (RISM) integral equation theory. Methanol particles are described as dimers formed from an apolar Lennard-Jones sphere, mimicking the methyl group, and a sphere with a core-softened potential as the hydroxyl group. Different closure approximations of the RISM theory are compared and discussed. The liquid structure of methanol is investigated by calculating site-site radial distribution functions and static structure factors for a wide range of temperatures and densities. Results obtained show a good agreement between RISM and Monte Carlo simulations. The phase behavior of methanol is investigated by employing different thermodynamic routes for the calculation of the RISM free energy, drawing gas-liquid coexistence curves that match the simulation data. Preliminary indications for a putative second critical point between two different liquid phases of methanol are also discussed.