Thermodynamic perturbation theory for multipolar and ionic liquids

Polar soft-SAFT: theory and comparison with molecular simulations and experimental data of pure polar fluids

The consideration of polar interactions is of vital importance for the development of predictive and accurate thermodynamic models for polar fluids, as they govern most of their thermodynamic properties, making them highly non-ideal fluids.

Linear and nonlinear magnetic properties of ferrofluids

Within a high-magnetic-field approximation, employing Ruelle's algebraic perturbation theory, a field-dependent free-energy expression is proposed which allows one to determine the magnetic properties of ferrofluids modeled as dipolar hard-sphere systems. We compare the ensuing magnetization curves, following from this free energy, with those obtained by Ivanov and Kuznetsova [Phys. Rev. E 64, 041405 (2001)] as well as with new corresponding Monte Carlo simulation data. Based on the power-series expansion of the magnetization, a closed expression for the magnetization is also proposed, which is a high-density extension of the corresponding equation of Ivanov and Kuznetsova. From both magnetization equations the zero-field susceptibility expression due to Tani et al. [Mol. Phys. 48, 863 (1983)] can be obtained, which is in good agreement with our MC simulation results. From the closed expression for the magnetization the second-order nonlinear magnetic susceptibility is also derived, which shows fair agreement with the corresponding MC simulation data.

A simulation assessment of the thermodynamics of dense ion-dipole mixtures with polarization

Molecular dynamics (MD) simulations are employed to ascertain the relative importance of various electrostatic interaction contributions, including induction interactions, to the thermodynamics of dense, hot ion-dipole mixtures. In the absence of polarization, we find that an MD-constrained free energy term accounting for the ion-dipole interactions, combined with well tested ionic and dipolar contributions, yields a simple, fairly accurate free energy form that may be a better option for describing the thermodynamics of such mixtures than the mean spherical approximation (MSA). Polarization contributions induced by the presence of permanent dipoles and ions are found to be additive to a good approximation, simplifying the thermodynamic modeling. We suggest simple free energy corrections that account for these two effects, based in part on standard perturbative treatments and partly on comparisons with MD simulation. Even though the proposed approximations likely need further study, they provide a first quantitative assessment of polarization contributions at high densities and temperatures and may serve as a guide for future modeling efforts.

Statistical thermodynamics of fluids with both dipole and quadrupole moments

New Gibbs ensemble simulation data for a polar fluid modeled by a square-well potential plus dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole interactions are presented. This simulation data is used in order to assess the applicability of the multipolar square-well perturbation theory [A. L. Benavides, Y. Guevara, and F. del Río, Physica A 202, 420 (1994)] to systems where more than one term in the multipole expansion is relevant. It is found that this theory is able to reproduce qualitatively well the vapor-liquid phase diagram for different multipolar moment strengths, corresponding to typical values of real molecules, except in the critical region. Hence, this theory is used to model the behavior of substances with multiple chemical bonds such as carbon monoxide and nitrous oxide and we found that with a suitable choice of the values of the intermolecular parameters, the vapor-liquid equilibrium of these species is adequately estimated.

Electrostatic fluctuations in cavities within polar liquids and thermodynamics of polar solvation

We present the results of numerical simulations of the fluctuations of the electrostatic potential and electric field inside cavities created in the fluid of dipolar hard spheres. We found that the thermodynamics of polar solvation changes dramatically when the cavity size becomes about 4-5 times larger than the size of the liquid particle. The range of small cavities can be reasonably understood within the framework of current solvation models. On the contrary, the regime of large cavities is characterized by a significant softening of the cavity interface resulting in a decay of the fluctuation variances with the cavity size much faster than anticipated by both the continuum electrostatics and microscopic theories. For instance, the variance of the electrostatic potential at the cavity center decays with the cavity radius R0 approximately as 1R{0};{4-6} instead of the 1R{0} scaling expected from the standard electrostatics. Our results suggest that cores of nonpolar molecular assemblies in polar liquids lose solvation strength much faster than is traditionally anticipated.

Exp6-polar thermodynamics of dense supercritical water

We introduce a simple polar fluid model for the thermodynamics of dense supercritical water based on a Buckingham (exp-6) core and point dipole representation of the water molecule. The proposed exp6-polar thermodynamics, which is based on ideas originally applied to dipolar hard spheres, performs very well when tested against molecular dynamics simulations. Comparisons of the model predictions with experimental data available for supercritical water yield excellent agreement for the shock Hugoniot, isotherms, and sound speeds, and are also quite good for the self-diffusion constant and relative dielectric constant.

Model energy landscapes of low-temperature fluids: Dipolar hard spheres

An analytical model of non-Gaussian energy landscape of low-temperature fluids is developed based on the thermodynamics of the fluid of dipolar hard spheres. The entire excitation profile of the liquid, from the high-temperature liquid to the point of ideal-glass transition, has been obtained from Monte Carlo simulations. The fluid of dipolar hard spheres loses stability close to the point of ideal-glass transition transforming via a first-order transition into a columnar liquid phase of dipolar chains locally arranged in a body-centered-tetragonal order. Significant non-Gaussianity of the energy landscape is responsible for narrowing of the distribution of potential energies and energies of inherent structures with decreasing temperature. We suggest that the proposed functionality of the enumeration function is widely applicable to both polar and nonpolar low-temperature liquids.

Predicting the Phase Behavior of Nitrogen + n-Alkanes for Enhanced Oil Recovery from the SAFT-VR Approach: Examining the Effect of the Quadrupole Moment

The phase behavior of nitrogen + n-alkanes is studied within the framework of the statistical associating fluid theory for potentials of variable range (SAFT-VR). The effect of the quadrupole moment of nitrogen on the phase behavior is considered through an extension of the SAFT-VR equation that includes an additional contribution to the Helmholtz free energy due to quadrupolar interactions. A significant improvement in the description of the phase diagram of the binary mixtures of nitrogen with different n-alkanes is obtained with the new approach when compared to predictions from the original SAFT-VR EOS (i.e., without the quadrupolar term). The experimental value for the quadrupole moment of nitrogen is used in the new equation; thus, no additional parameters are employed. Given the nonideal nature of the binary mixtures, a binary interaction parameter is needed to describe the full-phase diagram and high-pressure critical lines of these systems; however, this can be fitted to a single system and successfully used to predict the phase behavior of other binary mixtures without additional fitting. Furthermore, only a single, transferable, cross-energy parameter is required when the quadrupolar term is considered, whereas a cross-range parameter is also needed with the original SAFT-VR approach. The inclusion of the quadrupolar term in the equation of state therefore reduces the need to use effective parameters by explicitly including at the molecular level interactions due to the quadrupole moment.

Perturbed Chain-Statistical Associating Fluid Theory Extended to Dipolar and Quadrupolar Molecular Fluids

The perturbed chain statistical associating fluid theory (PC-SAFT) is extended to polar molecular fluids, namely dipolar and quadrupolar fluids. The extension is based on the perturbation theory for polar fluids by Stell and co-workers. Appropriate expressions are proposed for dipole-dipole, quadrupole-quadrupole, and dipole-quadrupole interactions. Furthermore, induced dipole interactions are calculated explicitly in the model. The new polar PC-SAFT model is relatively complex; for this purpose, a truncated polar PC-SAFT model is proposed using only the leading term in the polynomial expansion for polar interactions. The new model is used for the calculation of thermodynamic properties of various quadrupolar pure fluids. In all cases, the agreement between experimental data and model predictions is very good.

Modeling of the Carbon Dioxide Solubility in Imidazolium-Based Ionic Liquids with the tPC-PSAFT Equation of State

In this work, an equation of state (EoS) is developed to predict accurately the phase behavior of ionic liquid + CO2 systems based on the truncated perturbed chain polar statistical associating fluid theory (tPC-PSAFT) EoS. This EoS accounts explicitly for the dipolar interactions between ionic liquid molecules, the quadrupolar interactions between CO2 molecules, and the Lewis acid-base type of association between the ionic liquid and the CO2 molecules. Physically meaningful model pure-component parameters for ionic liquids are estimated based on literature data. All experimental vapor-liquid equilibrium data are correlated with a single linearly temperature-dependent binary interaction parameter. The ability of the model to describe accurately carbon dioxide solubility in various 1-alkyl-3-methylimidazolium-based ionic liquids with different alkyl chain lengths and different anions at pressures from 0 to 100 MPa and carbon dioxide fractions from 0 to 75 mol % is demonstrated. In all cases, good agreement with experimental data is obtained.

Robust quantification of short echo time1H magnetic resonance spectra using the Padé approximant

Accurate quantification of in vivo short echo time spectra is hampered by the presence of overlapping peaks and a significant baseline. In this work the Padé approximant in conjunction with Monte Carlo simulation is used to extract peak areas from short echo time 1H spectra. We exploit the fact that the Padé approximant is known to model broad non-Lorentzian signals as arbitrary sums of Lorentzian components to separate baseline components from sharper metabolite signals by combining the Padé approximant with Monte Carlo simulation. The simulation results demonstrate that the Padé approximant-Monte Carlo hybrid analysis is able to separate the metabolite signals from the baseline, while a least squares fitting of a time domain model may result in significant bias of the peak area estimations. For the in vivo data the estimates of the peak areas using the Padé approximant and AMARES compare well, with the exception of the NAA peak at 2.02 ppm. We suggest that the discrepancy may be due to the baseline contamination as supported by the simulation results; however, without an in vivo gold standard this remains difficult to demonstrate.