Coarse-Grained Model for Perfluorocarbons and Phase Equilibrium Simulation of Perfluorocarbons/CO2 Mixtures

Branching in molecular structure enhancement of solubility in CO2

Most compounds of some 1,000 amu molecular weight (MW) and higher are poorly soluble in carbon dioxide (CO2). Only at very high pressure, there may be mild solubility. This limits the use of CO2 as a solvent and modifications of CO2 properties through additives. We have developed a coarse-grained molecular model to investigate the dependency of the solubility of hydrocarbon oligomers (MW of ∼1,000 amu) in CO2 and on the molecular structure. The coarse-grained model is optimized by the particle swarm optimization algorithm to reproduce density, surface tension, and enthalpy of vaporization of a highly branched hydrocarbon oligomer (poly-1-decene with six repeating units). We demonstrate that branching in molecular structure of oligomers significantly increases solubility in CO2. The branching in molecular structure results in up to 270-time enhancement of solubility in CO2 than an n-alkane with the same MW. The number of structural edges (methyl group) is a key in improved CO2-philicity. The solubility of poly-1-decene with nine repeating units (MW of 1,264.4 amu) is higher in CO2 than poly-1-dodecene with six repeating units (MW of 1,011.93 amu) because it has more structural edges (10 vs. 7). These results shed light on the enhancement of CO2-philicity by altering molecular structure rather than modifying chemical composition in compounds.

Effect of Branching on Mutual Solubility of Alkane-CO2 Systems by Molecular Simulations

Mutual solubilities of hydrocarbon-CO₂ systems are important in a broad range of applications. Experimental data and theoretical understanding of phase behavior of large hydrocarbon molecules and CO₂ are limited. This is especially true in relation to the molecular structure of hydrocarbons when the carbon number exceeds 12. In this work, the continuous fractional component Gibbs ensemble Monte Carlo simulations are used to investigate mutual solubility of different alkane and CO₂ systems and the molecular structure. We investigate the mutual solubility of n-decane, n-hexadecane, n-eicosane, and the corresponding structural isomers in the CO₂-rich and hydrocarbon-rich phase. The focus will be solubility of the heavy normal alkanes and their structural isomers in CO₂. The simulation results are verified by comparing the experimental data when measurements are available. The simulation of phase behavior of the n-decane-CO₂ system agrees with the experiments. We also present simulation results of n-hexadecane-CO₂ and n-eicosane-CO₂ systems away from the critical region partly due to the finite size effect. We establish that solubility of the hydrocarbons in CO₂ is improved by change of the molecular structure in heavier alkanes. The enhanced solubility is limited in decane isomers, but the isomers of hexadecane and eicosane show 2- to 3-time solubility enhancement. The molecular dynamics simulations suggest that the improvement is from a higher coordination number of CO₂ for methyl (CH₃) rather than for methylene (CH₂) groups. This study sets the stage for molecular engineering and synthesis of hydrocarbons that are soluble in CO₂ not only by considering functionality but also by changing the molecular structure. The solubility enhancement is the first step in viscosification of CO₂ which broadens the use of CO₂.

Modeling Solvation in Supercritical CO2

In recent decades, a microscopic understanding of solute-solvent intermolecular interactions has been key to advances in technologies based on supercritical carbon dioxide. In many cases, computational work has provided the impetus for new discoveries, shedding new light on important concepts such as the local structure around the solute in the supercritical medium, the influence of the peculiar properties of the latter on the molecular behavior of dissolved substances and, importantly, CO₂ -philicity. In this Review, the theoretical work that has been relevant to these developments is surveyed and, by presenting some crucial open questions, the possible routes to achieving further progress based on the interplay between theory and experiments is discussed.

Classical force field for hydrofluorocarbon molecular simulations. Application to the study of gas solubility in poly(vinylidene fluoride)

A classical all-atoms force field for molecular simulations of hydrofluorocarbons (HFCs) has been developed. Lennard-Jones force centers plus point charges are used to represent dispersion-repulsion and electrostatic interactions. Parametrization of this force field has been performed iteratively using three target properties of pentafluorobutane: the quantum energy of an isolated molecule, the dielectric constant in the liquid phase, and the compressed liquid density. The accuracy and transferability of this new force field has been demonstrated through the simulation of different thermophysical properties of several fluorinated compounds, showing significant improvements compared to existing models. This new force field has been applied to study solubilities of several gases in poly(vinylidene fluoride) (PVDF) above the melting temperature of this polymer. The solubility of CH4, CO2, H2S, H2, N2, O2, and H2O at infinite dilution has been computed using test particle insertions in the course of a NpT hybrid Monte Carlo simulation. For CH4, CO2, and their mixtures, some calculations beyond the Henry regime have also been performed using hybrid Monte Carlo simulations in the osmotic ensemble, allowing both swelling and solubility determination. An ideal mixing behavior is observed, with identical solubility coefficients in the mixtures and in pure gas systems.

SAFT-γ force field for the simulation of molecular fluids: 2. Coarse-grained models of greenhouse gases, refrigerants, and long alkanes

In the first paper of this series [C. Avendaño, T. Lafitte, A. Galindo, C. S. Adjiman, G. Jackson, and E. A. Müller, J. Phys. Chem. B2011, 115, 11154] we introduced the SAFT-γ force field for molecular simulation of fluids. In our approach, a molecular-based equation of state (EoS) is used to obtain coarse-grained (CG) intermolecular potentials that can then be employed in molecular simulation over a wide range of thermodynamic conditions of the fluid. The macroscopic experimental data for the vapor-liquid equilibria (saturated liquid density and vapor pressure) of a given system are represented with the SAFT-VR Mie EoS and used to estimate effective intermolecular parameters that provide a good description of the thermodynamic properties by exploring a wide parameter space for models based on the Mie (generalized Lennard-Jones) potential. This methodology was first used to develop a simple single-segment CG Mie model of carbon dioxide (CO2) which allows for a reliable representation of the fluid-phase equilibria (for which the model was parametrized), as well as an accurate prediction of other properties such as the enthalpy of vaporization, interfacial tension, supercritical density, and second-derivative thermodynamic properties (thermal expansivity, isothermal compressibility, heat capacity, Joule-Thomson coefficient, and speed of sound). In our current paper, the methodology is further applied and extended to develop effective SAFT-γ CG Mie force fields for some important greenhouse gases including carbon tetrafluoride (CF4) and sulfur hexafluoride (SF6), modeled as simple spherical molecules, and for long linear alkanes including n-decane (n-C10H22) and n-eicosane (n-C20H42), modeled as homonuclear chains of spherical Mie segments. We also apply the SAFT-γ methodology to obtain a CG homonuclear two-segment Mie intermolecular potential for the more challenging polar and asymmetric compound 2,3,3,3-tetrafluoro-1-propene (HFO-1234yf), a novel replacement refrigerant with promising properties. The description of the fluid-phase behavior and the prediction of the other thermophysical properties obtained by molecular simulation using our SAFT-γ CG Mie force fields are found to be of comparable quality (and sometimes superior) to that obtained using the more sophisticated all-atom (AA) and united-atom (UA) models commonly employed in the field. We should emphasize that though the focus of our current work is on simple homonuclear models, the SAFT-γ methodology is based on a group contribution methodology which is naturally suited to the development of more sophisticated heteronuclear models.