Fick Diffusion Coefficients of Liquid Mixtures Directly Obtained From Equilibrium Molecular Dynamics

Mutual Diffusivities of Mixtures of Carbon Dioxide and Hydrogen and Their Solubilities in Brine: Insight from Molecular Simulations

H2-CO2 mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H2-CO2 mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H2 due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO2, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO2 content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO2-H2-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO2 and H2 in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO2 in NaCl brine decreased in the ternary system compared to the binary CO2-NaCl brine system, the solubility of H2 in NaCl brine increased less in the ternary system compared to the binary H2-NaCl brine system. The cooperative effect of H2-CO2 enhances the H2 solubility while suppressing the CO2 solubility. The water content in the gas phase was found to be intermediate between H2-NaCl brine and CO2-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO2-H2 mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H2-CO2 mixtures.

Mathematical Models for Cholesterol Metabolism and Transport

Cholesterol is an essential component of eukaryotic cellular membranes. It is also an important precursor for making other molecules needed by the body. Cholesterol homeostasis plays an essential role in human health. Having high cholesterol can increase the chances of getting heart disease. As a result of the risks associated with high cholesterol, it is imperative that studies are conducted to determine the best course of action to reduce whole body cholesterol levels. Mathematical models can provide direction on this. By examining existing models, the suitable reactions or processes for drug targeting to lower whole-body cholesterol can be determined. This paper examines existing models in the literature that, in total, cover most of the processes involving cholesterol metabolism and transport, including: the absorption of cholesterol in the intestine; the cholesterol biosynthesis in the liver; the storage and transport of cholesterol between the intestine, the liver, blood vessels, and peripheral cells. The findings presented in these models will be discussed for potential combination to form a comprehensive model of cholesterol within the entire body, which is then taken as an in-silico patient for identifying drug targets, screening drugs, and designing intervention strategies to regulate cholesterol levels in the human body.

Characterizing the Diffusion and Rheological Properties of Aged Asphalt Binder Rejuvenated with Bio-Oil Based on Molecular Dynamic Simulations and Laboratory Experimentations

Soybean-derived bio-oil is one of the vegetable-based oils that is gaining the most interest for potential use in the rejuvenation of aged asphalt binders. This laboratory study was conducted to characterize and quantify the diffusion and rheological properties of bio-oil-rejuvenated aged asphalt binder (BRAA) using soybean oil. In the study, the chemical structure of the soybean oil was comparatively characterized using an element analyzer (EA), gel permeation chromatography (GPC), and a Fourier infrared (FTIR) spectrometer, respectively. Based on the chemical structure of the bio-oil, BRAA molecular models were built for computing the diffusion parameters using molecular dynamic simulations. Likewise, a dynamic shear rheometer (DSR) test device was used for measuring and quantifying the rheological properties of the aged asphalt binder rejuvenated with 0%, 1%, 2%, 3%, 4%, and 5% soybean oil, respectively. The laboratory test results indicate that bio-oil could potentially improve the diffusion coefficients and phase angle of the aged asphalt binder. Similarly, the corresponding decrease in the complex shear modulus has a positive effect on the low-temperature properties of BRAA. For a bio-oil dosage 4.0%, the diffusion coefficients of the BRAA components are 1.52 × 10⁻⁸, 1.33 × 10⁻⁸, 3.47 × 10⁻⁸, 4.82 × 10⁻⁸ and 3.92 × 10⁻⁸, respectively. Similarly, the corresponding reduction in the complex shear modulus from 1.27 × 10⁷ Pa to 4.0 × 10⁵ Pa suggests an improvement in the low-temperature properties of BRAA. Overall, the study contributes to the literature on the potential use of soybean-derived bio-oil as a rejuvenator of aged asphalt binders.

How sensitive are physical properties of choline chloride-urea mixtures to composition changes: Molecular dynamics simulations and Kirkwood-Buff theory

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.

Generalized Form for Finite-Size Corrections in Mutual Diffusion Coefficients of Multicomponent Mixtures Obtained from Equilibrium Molecular Dynamics Simulation

The system-size dependence of computed mutual diffusion coefficients of multicomponent mixtures is investigated, and a generalized correction term is derived. The generalized finite-size correction term was validated for the ternary molecular mixture chloroform/acetone/methanol as well as 28 ternary LJ systems. It is shown that only the diagonal elements of the Fick matrix show system-size dependency. The finite-size effects of these elements can be corrected by adding the term derived by Yeh and Hummer (J. Phys. Chem. B2004, 108, 15873–15879). By performing an eigenvalue analysis of the finite-size effects of the matrix of Fick diffusivities we show that the eigenvector matrix of Fick diffusivities does not depend on the size of the simulation box. Only eigenvalues, which describe the speed of diffusion, depend on the size of the system. An analytic relation for finite-size effects of the matrix of Maxwell–Stefan diffusivities was developed. All Maxwell–Stefan diffusivities depend on the system size, and the required correction depends on the matrix of thermodynamic factors.

Mutual and Thermal Diffusivities as well as Fluid-Phase Equilibria of Mixtures of 1-Hexanol and Carbon Dioxide

This work contributes to an improved understanding of the fluid-phase behavior and diffusion processes in mixtures of 1-hexanol and carbon dioxide (CO₂) at temperatures around the upper critical end point (UCEP) of the system. Raman spectroscopy and dynamic light scattering were used to determine the composition at saturation conditions as well as Fick and thermal diffusivities. An acceleration of the Fick diffusive process up to CO₂ mole fractions of about 0.2 was found, followed by a strong slowing-down approaching vapor-liquid-liquid equilibrium or critical conditions. The acceleration of the Fick diffusive process vanished at temperatures much higher than the UCEP. Experimental Fick diffusivity data were compared with predictions from equilibrium molecular dynamics simulations and excess Gibbs energy calculations using interaction parameters from the literature. Both theoretical methods were not able to predict that the thermodynamic factor is equal to zero at the spinodal composition, stressing the need for new methodologies under such conditions. Thus, new sets of temperature-dependent interaction parameters were developed for the nonrandom two-liquid model, which improve the prediction of the Fick diffusion coefficient considerably. The link between the Fick diffusion coefficient and the nonrandomness of the liquid phases is also discussed.

Diffusion in Binary Aqueous Solutions of Alcohols by Molecular Simulation

Based on the molecular dynamics method, the calculations for diffusion coefficients were carried out in binary aqueous solutions of three alcohols: ethanol, isopropanol, and tert-butanol. The intermolecular potential TIP4P/2005 was used for water; and five force fields were analyzed for the alcohols. The force fields providing the best accuracy of calculation were identified based on a comparison of the calculated self-diffusion coefficients of pure alcohols with the experimental data for internal (Einstein) diffusion coefficients of alcohols in solutions. The temperature and concentration dependences of the interdiffusion coefficients were determined using Darken’s Equation. Transport (Fickian) diffusion coefficients were calculated using a thermodynamic factor determined by the non-random two-liquid (NRTL) and Willson models. It was demonstrated that for adequate reproduction of the experimental data when calculating the transport diffusion coefficients, the thermodynamic factor has to be 0.64. Simple approximations were obtained, providing satisfactory accuracy in calculating the concentration and temperature dependences of the transport diffusion coefficients in the studied mixtures.

Multicomponent mutual diffusion in the warm, dense matter regime

We present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.

Optimizing Nonbonded Interactions of the OPLS Force Field for Aqueous Solutions of Carbohydrates: How to Capture Both Thermodynamics and Dynamics

Knowledge on thermodynamic and transport properties of aqueous solutions of carbohydrates is of great interest for process and product design in the food, pharmaceutical, and biotechnological industries. Molecular simulation is a powerful tool to calculate these properties, but current classical force fields cannot provide accurate estimates for all properties of interest. The poor performance of the force fields is mainly observed for concentrated solutions, where solute–solute interactions are overestimated. In this study, we propose a method to refine force fields, such that solute–solute interactions are more accurately described. The OPLS force field combined with the SPC/Fw water model is used as a basis. We scale the nonbonded interaction parameters of sucrose, a disaccharide. The scaling factors are chosen in such a way that experimental thermodynamic and transport properties of aqueous solutions of sucrose are accurately reproduced. Using a scaling factor of 0.8 for Lennard-Jones energy parameters (ϵ) and a scaling factor of 0.95 for partial atomic charges (q), we find excellent agreement between experiments and computed liquid densities, thermodynamic factors, shear viscosities, self-diffusion coefficients, and Fick (mutual) diffusion coefficients. The transferability of these optimum scaling factors to other carbohydrates is verified by computing thermodynamic and transport properties of aqueous solutions of d-glucose, a monosaccharide. The good agreement between computed properties and experiments suggests that the scaled interaction parameters are transferable to other carbohydrates, especially for concentrated solutions.

Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics

A method is proposed for calculating the shear viscosity of a liquid from finite-size effects of self-diffusion coefficients in Molecular Dynamics simulations. This method uses the difference in the self-diffusivities, computed from at least two system sizes, and an analytic equation to calculate the shear viscosity. To enable the efficient use of this method, a set of guidelines is developed. The most efficient number of system sizes is two and the large system is at least four times the small system. The number of independent simulations for each system size should be assigned in such a way that 50%–70% of the total available computational resources are allocated to the large system. We verified the method for 250 binary and 26 ternary Lennard-Jones systems, pure water, and an ionic liquid ([Bmim][Tf₂N]). The computed shear viscosities are in good agreement with viscosities obtained from equilibrium Molecular Dynamics simulations for all liquid systems far from the critical point. Our results indicate that the proposed method is suitable for multicomponent mixtures and highly viscous liquids. This may enable the systematic screening of the viscosities of ionic liquids and deep eutectic solvents.

Application of Reverse Nonequilibrium Molecular Dynamics to the Calculation of the Mutual Diffusion Coefficient of Alkane Mixtures

In a recent publication, a reverse nonequilibrium molecular dynamics (RNEMD) method was presented for computing the mutual diffusion coefficient of liquid mixtures. A concentration gradient and a subsequent mass flux are induced in the system by suitably exchanging molecules in different regions. The algorithm has been successfully tested on Lennard-Jones mixtures and molecular fluid mixtures with molecules having the same number of particles. In this work, a modification is made to the RNEMD method to determine the mutual diffusion coefficient of binary liquid mixtures with molecules having different sizes and masses. To migrate molecules of a different type, the splitting method has been used in this work. Investigation of the resulting steady-state mass fraction profile allows the evaluation of the mutual diffusion coefficient. For validation, the mutual diffusion coefficients of ethane-propane and ethane-pentane liquid mixtures at different compositions and temperatures have been obtained using this method. The mutual diffusion coefficients obtained from the RNEMD simulations are within the error bars of values obtained by equilibrium molecular dynamics for the identical model and conditions. The excess energy released due to the exchange of molecules is efficiently removed by strongly coupling a local thermostat in the region around the insertion point. There is no heating of the analysis region.

Interplay of structure and diffusion in ternary liquid mixtures of benzene + acetone + varying alcohols

The Fick diffusion coefficient matrix of ternary mixtures containing benzene + acetone + three different alcohols, i.e., methanol, ethanol, and 2-propanol, is studied by molecular dynamics simulation and Taylor dispersion experiments. Aiming to identify common features of these mixtures, it is found that one of the main diffusion coefficients and the smaller eigenvalue do not depend on the type of alcohol along the studied composition path. Two mechanisms that are responsible for this invariant behavior are discussed in detail, i.e., the interplay between kinetic and thermodynamic contributions to Fick diffusion coefficients and the presence of microscopic heterogeneities caused by hydrogen bonding. Experimental work alone cannot explain these mechanisms, while present simulations on the molecular level indicate structural changes and uniform intermolecular interactions between benzene and acetone molecules in the three ternary mixtures. The main diffusion coefficients of these ternary mixtures exhibit similarities with their binary subsystems. Analyses of radial distribution functions and hydrogen bonding statistics quantitatively evidence alcohol self-association and cluster formation, as well as component segregation. Furthermore, the excess volume of the mixtures is analyzed in the light of intermolecular interactions, further demonstrating the benefits of the simultaneous use of experiment and simulation. The proposed framework for studying diffusion coefficients of a set of ternary mixtures, where only one component varies, opens the way for further investigations and a better understanding of multicomponent diffusion. The presented numerical results may also give an impulse to the development of predictive approaches for multicomponent diffusion.

Influence of Liquid Structure on Fickian Diffusion in Binary Mixtures of n-Hexane and Carbon Dioxide Probed by Dynamic Light Scattering, Raman Spectroscopy, and Molecular Dynamics Simulations

This study contributes to a fundamental understanding of how the liquid structure in a model system consisting of weakly associative n-hexane ( n-C₆H₁₄) and carbon dioxide (CO₂) influences the Fickian diffusion process. For this, the benefits of light scattering experiments and molecular dynamics (MD) simulations at macroscopic thermodynamic equilibrium were combined synergistically. Our reference Fickian diffusivities measured by dynamic light scattering (DLS) revealed an unusual trend with increasing CO₂ mole fractions up to about 70 mol %, which agrees with our simulation results. The molecular impacts on the Fickian diffusion were analyzed by MD simulations, where kinetic contributions related to the Maxwell-Stefan (MS) diffusivity and structural contributions quantified by the thermodynamic factor were studied separately. Both the MS diffusivity and the thermodynamic factor indicate the deceleration of Fickian diffusion compared to an ideal mixture behavior. Computed radial distribution functions as well as a significant blue-shift of the CH stretching modes of n-C₆H₁₄ identified by Raman spectroscopy show that the slowing down of the diffusion is caused by a structural organization in the binary mixtures over a broad concentration range in the form of self-associated n-C₆H₁₄ and CO₂ domains. These networks start to form close to the infinite dilution limits and seem to have their largest extent at a solute-solvent transition point at about 70 mol % CO₂. The current results not only improve the general understanding of mass diffusion in liquids but also serve to develop sound prediction models for Fick diffusivities.

Finite-Size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics

Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.

Thermal, Mutual, and Self-Diffusivities of Binary Liquid Mixtures Consisting of Gases Dissolved in n-Alkanes at Infinite Dilution

In the present study, dynamic light scattering (DLS) experiments and molecular dynamics (MD) simulations were used for the investigation of the molecular diffusion in binary mixtures of liquids with dissolved gases at macroscopic thermodynamic equilibrium. Model systems based on the n-alkane n-hexane or n-decane with dissolved hydrogen, helium, nitrogen, or carbon monoxide were studied at temperatures between 303 and 423 K and at gas mole fractions below 0.06. With DLS, the relaxation behavior of microscopic equilibrium fluctuations in concentration and temperature is analyzed to determine simultaneously mutual and thermal diffusivity in an absolute way. The present measurements document that even for mole gas fractions of 0.007 and Lewis numbers close to 1, reliable mutual diffusivities with an average expanded uncertainty ( k = 2) of 13% can be obtained. By use of suitable molecular models for the mixture components, the self-diffusion coefficient of the gases was determined by MD simulations with an averaged expanded uncertainty ( k = 2) of 7%. The DLS experiments showed that the thermal diffusivity of the studied systems is not affected by the dissolved gas and agrees with the reference data for the pure n-alkanes. In agreement with theory, mutual diffusivities and self-diffusivities were found to be equal mostly within combined uncertainties at conditions approaching infinite dilution of the gas. Our DLS and MD results, representing the first available data for the present systems, reveal distinctly larger mass diffusivities for mixtures containing hydrogen or helium compared to mixtures containing nitrogen or carbon monoxide. On the basis of the broad range of mass diffusivities of the studied gas-liquid systems covering about 2 orders of magnitude from about 10^-9 to 10^-7 m²·s^-1, effects of the solvent and solute properties on the temperature-dependent mass diffusivities are discussed. This contributed to the development of a simple semiempirical correlation for the mass diffusivity of the studied gases dissolved in n-alkanes of varying chain length at infinite dilution as a function of temperature. The generalized expression requiring only information on the kinematic viscosity and molar mass of the pure solvent as well as the molar mass and acentric factor of the solute represents the database from this work and further literature with an absolute average deviation of about 11%.

MD simulation study of the diffusion and local structure of n-alkanes in liquid and supercritical methanol at infinite dilution

The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility. Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low.

The multiscale coarse-graining method. XI. Accurate interactions based on the centers of charge of coarse-grained sites

It is essential to be able to systematically construct coarse-grained (CG) models that can efficiently and accurately reproduce key properties of higher-resolution models such as all-atom. To fulfill this goal, a mapping operator is needed to transform the higher-resolution configuration to a CG configuration. Certain mapping operators, however, may lose information related to the underlying electrostatic properties. In this paper, a new mapping operator based on the centers of charge of CG sites is proposed to address this issue. Four example systems are chosen to demonstrate this concept. Within the multiscale coarse-graining framework, CG models that use this mapping operator are found to better reproduce the structural correlations of atomistic models. The present work also demonstrates the flexibility of the mapping operator and the robustness of the force matching method. For instance, important functional groups can be isolated and emphasized in the CG model.

The reaction enthalpy of hydrogen dissociation calculated with the Small System Method from simulation of molecular fluctuations

We show how we can find the enthalpy of a chemical reaction under non-ideal conditions using the Small System Method to sample molecular dynamics simulation data for fluctuating variables. This method, created with Hill's thermodynamic analysis, is used to find properties in the thermodynamic limit, such as thermodynamic correction factors, partial enthalpies, volumes, heat capacities and compressibility. The values in the thermodynamic limit at (T,V, μj) are then easily transformed into other ensembles, (T,V,Nj) and (T,P,Nj), where the last ensemble gives the partial molar properties which are of interest to chemists. The dissociation of hydrogen from molecules to atoms was used as a convenient model system. Molecular dynamics simulations were performed with three densities; ρ = 0.0052 g cm(-3) (gas), ρ = 0.0191 g cm(-3) (compressed gas) and ρ = 0.0695 g cm(-3) (liquid), and temperatures in the range; T = 3640-20,800 K. The enthalpy of reaction was observed to follow a quadratic trend as a function of temperature for all densities. The enthalpy of reaction was observed to only have a small pressure dependence. With a reference point close to an ideal state (T = 3640 K and ρ = 0.0052 g cm(-3)), we were able to calculate the thermodynamic equilibrium constant, and thus the deviation from ideal conditions for the lowest density. We found the thermodynamic equilibrium constant to increase with increasing temperature, and to have a negligible pressure dependence. Taking the enthalpy variation into account in the calculation of the thermodynamic equilibrium constant, we found the ratio of activity coefficients to be in the order of 0.7-1.0 for the lowest density, indicating repulsive forces between H and H2. This study shows that the compressed gas- and liquid density values at higher temperatures are far from those calculated under ideal conditions. It is important to have a method that can give access to partial molar properties, independent of the ideality of the reacting mixture. Our results show how this can be achieved with the use of the Small System Method.

Partial molar enthalpies and reaction enthalpies from equilibrium molecular dynamics simulation

We present a new molecular simulation technique for determining partial molar enthalpies in mixtures of gases and liquids from single simulations, without relying on particle insertions, deletions, or identity changes. The method can also be applied to systems with chemical reactions. We demonstrate our method for binary mixtures of Weeks-Chandler-Anderson particles by comparing with conventional simulation techniques, as well as for a simple model that mimics a chemical reaction. The method considers small subsystems inside a large reservoir (i.e., the simulation box), and uses the construction of Hill to compute properties in the thermodynamic limit from small-scale fluctuations. Results obtained with the new method are in excellent agreement with those from previous methods. Especially for modeling chemical reactions, our method can be a valuable tool for determining reaction enthalpies directly from a single MD simulation.

Calculation of the chemical potential and the activity coefficient of two layers of CO2 adsorbed on a graphite surface

Thermodynamics of two layers of CO₂ on a graphite surface obtained directly from the simulations and the Small System Method.

Transport diffusivity of propane and propylene inside SWNTs from equilibrium molecular dynamics simulations

The gas transport of two model gases (propane and propylene) inside the single-walled nanotubes (SWNTs) of various diameters was systematically investigated using the molecular dynamics (MD) simulations. The thermodynamic factor can be obtained directly from equilibrium MD simulations following the newly-minted method proposed by Schnell et al. (Chem. Phys. Lett., 2011, 504, 199-201). This process eliminates the need to implement the tedious and challenging Monte Carlo simulations for the adsorption isotherm, from which the thermodynamic factor is usually extracted. The satisfactory agreement between simulation and the literature is found for self-diffusivity, corrected diffusivity and transport diffusivity, as well as for the thermodynamic factor. The ideal selectivity for a propane-propylene mixture through SWNT membranes could be optimized through adjusting the concentration gradient. This method can be readily extended to the binary and multiple-component systems.

Multicomponent diffusion in molten LiCl-KCl: dynamical correlations and divergent Maxwell-Stefan diffusivities

Multicomponent diffusional mechanisms in the ternary LiCl-KCl system are elucidated using the Green-Kubo formalism and equilibrium molecular dynamics simulations. The Maxwell-Stefan (MS) diffusion matrix is evaluated from the Onsager dynamical matrix that contains the diffusion flux correlation functions. From the temporal behavior of the correlation functions, we observe that the Li-Li and Li-Cl ion pairs have a pronounced cage dynamics that remains noticeably strong even at high temperatures. Even though the Onsager coefficients, which are the time integrals of the diffusion flux correlation functions, portray a relatively smooth variation across various compositions and temperatures, we observe a sign change and a divergent-like behavior for the MS diffusivity of the K-Li ion pair at a temperature of ~1100 K for the eutectic composition, and at a KCl mole fraction of ~0.49 at 1043 K. Negative MS diffusivities, while unusual, are however shown to satisfy the nonnegative entropic constraints.