RESEARCH NOTE Molecular simulation of the liquid-liquid equilibria of binary mixtures containing dipolar and non-polar components interacting via the Keesom potential

Quantifying the effect of polar interactions on the behavior of binary mixtures: Phase, interfacial, and excess properties

In this work, polar soft-Statistical Associating Fluid Theory (SAFT) was used in a systematic manner to quantify the influence of polar interactions on the phase equilibria, interfacial, and excess properties of binary mixtures. The theory was first validated with available molecular simulation data and then used to isolate the effect of polar interactions on the thermodynamic behavior of the mixtures by fixing the polar moment of one component while changing the polar moment of the second component from non-polar to either highly dipolar or quadrupolar, examining 15 different binary mixtures. It was determined that the type and magnitude of polar interactions have direct implications on the vapor-liquid equilibria (VLE), resulting in azeotropy for systems of either dipolar or quadrupolar fluids when mixed with non-polar or low polar strength fluids, while increasing the polar strength of one component shifts the VLE to be more ideal. Additionally, excess properties and interfacial properties such as interfacial tension, density profiles, and relative adsorption at the interface were also examined, establishing distinct enrichment in the case of mixtures with highly quadrupolar fluids. Finally, polar soft-SAFT was applied to describe the thermodynamic behavior of binary mixtures of experimental systems exhibiting various intermolecular interactions (non-polar and polar), not only demonstrating high accuracy and robustness through agreement with experimental data but also providing insights into the effect of polarity on the interfacial properties of the studied mixtures. This work proves the value of having an accurate theory for isolating the effect of polarity, especially for the design of ad hoc polar solvents.

Isomolar semigrand ensemble molecular dynamics: Development and application to liquid-liquid equilibria

An extended system molecular dynamics method for the isomolar semigrand ensemble (fixed number of particles, pressure, temperature, and fugacity fraction) is developed and applied to the calculation of liquid-liquid equilibria (LLE) for two Lennard-Jones mixtures. The method utilizes an extended system variable to dynamically control the fugacity fraction xi of the mixture by gradually transforming the identity of particles in the system. Two approaches are used to compute coexistence points. The first approach uses multiple-histogram reweighting techniques to determine the coexistence xi and compositions of each phase at temperatures near the upper critical solution temperature. The second approach, useful for cases in which there is no critical solution temperature, is based on principles of small system thermodynamics. In this case a coexistence point is found by running N-P-T-xi simulations at a common temperature and pressure and varying the fugacity fraction to map out the difference in chemical potential between the two species A and B (mu(A)-mu(B)) as a function of composition. Once this curve is known the equal-distance/equal-area criterion is used to determine the coexistence point. Both approaches give results that are comparable to those of previous Monte Carlo (MC) simulations. By formulating this approach in a molecular dynamics framework, it should be easier to compute the LLE of complex molecules whose intramolecular degrees of freedom are often difficult to properly sample with MC techniques.