1
|
Shagurin A, Miannay FA, Kiselev MG, Jedlovszky P, Affouard F, Idrissi A. Widom Line in Supercritical Water in Terms of Changes in Local Structure: Theoretical Perspective. J Phys Chem Lett 2024; 15:5831-5837. [PMID: 38787641 DOI: 10.1021/acs.jpclett.4c01142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Performing molecular dynamics simulations with the TIP4P/2005 water model along 9 isobars (from 175 to 375 bar) in the temperature range between 300 and 1100 K, we have found that the loci of the extrema in the rate of change of specific structural properties can be used to define purely structure-based Widom lines. We have examined several parameters that describe the local structure of water, such as the tetrahedral arrangement, nearest neighbor distance, local density around the molecules, and the size of the largest dense domain. The last two parameters were determined using the Voronoi polyhedral and density-based spatial clustering of applications with noise methods, respectively. By analyzing the moments of the associated distributions, we show that along a given isobar, the temperature at which we observe a maximum in the fluctuation, the rate of change of the average values, or in the skewness values unambiguously determines the Widom line that is in agreement with the experimentally detected, thermodynamic response function-based ones.
Collapse
Affiliation(s)
- Artem Shagurin
- University of Lille, CNRS UMR 8516 -LASIRe - Laboratoire Avancé de Spectroscopie pour les Interactions la Réactivité et l'environnement, 59000 Lille, France
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 - UMET - Unité Matériaux et Transformations, F-59000 Lille, France
- Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Ivanovo, 153045 Russia
| | - Francois A Miannay
- University of Lille, CNRS UMR 8516 -LASIRe - Laboratoire Avancé de Spectroscopie pour les Interactions la Réactivité et l'environnement, 59000 Lille, France
| | - Michael G Kiselev
- Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Ivanovo, 153045 Russia
| | - Pal Jedlovszky
- Department of Chemistry, Eszterházy Károly Catholic University, Leányka u. 6, 3300 Eger, Hungary
| | - Frederic Affouard
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 - UMET - Unité Matériaux et Transformations, F-59000 Lille, France
| | - Abdenacer Idrissi
- University of Lille, CNRS UMR 8516 -LASIRe - Laboratoire Avancé de Spectroscopie pour les Interactions la Réactivité et l'environnement, 59000 Lille, France
| |
Collapse
|
2
|
Saric D, Guevara-Carrion G, Gaponenko Y, Shevtsova V, Vrabec J. Diffusion of hydrocarbons diluted in supercritical carbon dioxide. Sci Rep 2023; 13:16107. [PMID: 37752219 PMCID: PMC10522683 DOI: 10.1038/s41598-023-42892-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Accepted: 09/15/2023] [Indexed: 09/28/2023] Open
Abstract
Mutual diffusion of six hydrocarbons (methane, ethane, isobutane, benzene, toluene or naphthalene) diluted in supercritical carbon dioxide ([Formula: see text]) is studied by molecular dynamics simulation near the Widom line, i.e., in the temperature range from 290 to 345 K along the isobar 9 MPa. The [Formula: see text] + aromatics mixtures are additionally sampled at 10 and 12 MPa and an experimental database with Fick diffusion coefficient data for those systems is provided. Taylor dispersion experiments of [Formula: see text] with benzene, toluene, n-dodecane and 1,2,3,4-tetrahydronaphthalene are conducted along the [Formula: see text] 10 MPa isobar. Maxwell-Stefan and Fick diffusion coefficients are analyzed, together with the thermodynamic factor that relates them. It is found that the peculiar behavior of the Fick diffusion coefficient of some [Formula: see text] mixtures in the extended critical region is a consequence of the thermodynamic factor minimum due to pronounced clustering on the molecular scale. Further, the strong dependence of the Fick diffusion coefficient on the molecular mass of the solute as well as the breakdown of the Stokes-Einstein relation near the Widom line are confirmed. Eleven correlations for the prediction of the Fick diffusion coefficient of [Formula: see text] mixtures are assessed. An alternative two-step approach for the prediction of the infinite dilution Fick diffusion coefficient of supercritical [Formula: see text] mixtures is proposed. It requires only the state point in terms of temperature and pressure (or density) as well as the molecular solute mass as input parameters. First, entropy scaling is applied to estimate the self-diffusion coefficient of [Formula: see text]. Subsequently, this coefficient is used to determine the infinite dilution Fick diffusion coefficient of the mixture, based on the finding that these two diffusion coefficients exhibit a linear relationship, where the slope depends only on the molecular solute mass.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany
| | | | - Yury Gaponenko
- MRC, CP-165/62, Université libre de Bruxelles (ULB), Ave. F.D. Roosevelt 50, B-1050, Brussels, Belgium
| | - Valentina Shevtsova
- Fluid Mechanics Group, Faculty of Engineering, Mondragon University, 20500, Mondragon, Spain
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, 48009, Bilbao, Spain
| | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany.
| |
Collapse
|
3
|
Fomin YD. Boiling line and near-critical maxima of propane-nitrogen mixtures. Phys Rev E 2022; 106:064102. [PMID: 36671079 DOI: 10.1103/physreve.106.064102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 11/09/2022] [Indexed: 06/17/2023]
Abstract
It is well known that some thermodynamic quantities demonstrate maxima in the vicinity of a critical point. The lines of these maxima in the density-temperature or pressure-temperature planes are called "Widom lines." The behavior of Widom lines of one-component fluids has already been well studied in a number of papers by different authors. However, up to now the understanding of Widom lines in binary mixtures is still lacking. In this paper we study the boiling curve and the near-critical maxima of mixtures of nitrogen and propane by means of molecular dynamics simulation. We calculate the boiling curves and estimate the critical temperatures in a set of concentrations from pure nitrogen to pure propane. The influence of the composition of the mixture on the Widom lines of the system is evaluated. We find that the mixture of propane and nitrogen behaves as a type I mixture in the van Konynenburg-Scott classification, i.e., when the concentration is changed, the critical point and the corresponding Widom lines continuously shift in the density-temperature plane.
Collapse
Affiliation(s)
- Yu D Fomin
- Vereshchagin Institute of High Pressure Physics, Russian Academy of Sciences, Kaluzhskoe shosse 14, Troitsk, Moscow 108840, Russia
| |
Collapse
|
4
|
Kalabiński J, Drozd-Rzoska A, Rzoska SJ. New insight into 3-picoline-deuterium oxide (D 2O) mixtures of limited miscibility with the lower critical consolute temperature. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2022; 45:73. [PMID: 36074188 PMCID: PMC9458598 DOI: 10.1140/epje/s10189-022-00224-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/20/2022] [Accepted: 08/05/2022] [Indexed: 06/15/2023]
Abstract
Coexistence curves in mixtures of limited miscibility with the lower critical consolute temperature (LCT), of 3-picoline with deuterium oxide (D2O), and D2O/H2O have been determined. They were tested with respect to the order parameter, the diameter of the binodal, and coordinates of the critical consolute point (critical temperature [Formula: see text] and critical concentration [Formula: see text]). Studies were carried out using the innovative method based on the analysis of relative volumes occupied by coexisting phases, yielding high-resolution data. The clear violation of the Cailletet-Mathew law of rectilinear diameter for the LCT mixtures of limited miscibility is evidenced. For the order parameter, the new distortion-sensitive analysis method yielded the evidence for the model-value of the order parameter critical exponent [Formula: see text], up to ca. 1 K from [Formula: see text]. Finally, the simple & easy method for determining the critical concentration by testing relative volumes of coexisting phases (or alternatively fractional meniscus heights, h) is presented. The significance of the invariant value [Formula: see text] is highlighted. The appearance of the milky & bluish critical opalescence is also shown.
Collapse
Affiliation(s)
- Jakub Kalabiński
- Institute of High Pressure Physics Polish Academy of Sciences, ul. Sokołowska 29/37, 01-142, Warsaw, Poland.
| | - Aleksandra Drozd-Rzoska
- Institute of High Pressure Physics Polish Academy of Sciences, ul. Sokołowska 29/37, 01-142, Warsaw, Poland
| | - Sylwester J Rzoska
- Institute of High Pressure Physics Polish Academy of Sciences, ul. Sokołowska 29/37, 01-142, Warsaw, Poland
| |
Collapse
|
5
|
Widom line of supercritical CO2 calculated by equations of state and molecular dynamics simulation. J CO2 UTIL 2022. [DOI: 10.1016/j.jcou.2022.102075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
6
|
widom and extrema lines as CRITERIA for OPTIMIZING operating conditions IN supercritical processes. J Supercrit Fluids 2022. [DOI: 10.1016/j.supflu.2022.105587] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
7
|
Servis MJ, Nayak S, Seifert S. The pervasive impact of critical fluctuations in liquid-liquid extraction organic phases. J Chem Phys 2021; 155:244506. [PMID: 34972370 DOI: 10.1063/5.0074995] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Overall, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.
Collapse
Affiliation(s)
- Michael J Servis
- Argonne National Laboratory, Chemical Sciences and Engineering Division, Lemont, Illinois 60439, USA
| | - Srikanth Nayak
- Argonne National Laboratory, Chemical Sciences and Engineering Division, Lemont, Illinois 60439, USA
| | - Soenke Seifert
- Argonne National Laboratory, X-ray Science Division, Lemont, Illinois 60439, USA
| |
Collapse
|
8
|
The hybrid ergodic lattice gas model for critical fluids and the molecular nature of the critical point. J Supercrit Fluids 2021. [DOI: 10.1016/j.supflu.2021.105505] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
9
|
The Anomalous Behavior of Thermodynamic Parameters in the Three Widom Deltas of Carbon Dioxide-Ethanol Mixture. Int J Mol Sci 2021; 22:ijms22189813. [PMID: 34575970 PMCID: PMC8472178 DOI: 10.3390/ijms22189813] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2021] [Revised: 09/07/2021] [Accepted: 09/08/2021] [Indexed: 11/23/2022] Open
Abstract
Using molecular dynamics, we demonstrated that in the mixture of carbon dioxide and ethanol (25% molar fraction) there are three pronounced regions on the p-T diagram characterized by not only high-density fluctuations but also anomalous behavior of thermodynamic parameters. The regions are interpreted as Widom deltas. The regions were identified as a result of analyzing the dependences of density, density fluctuations, isobaric thermal conductivity, and clustering of a mixture of carbon dioxide and ethanol in a wide range of pressures and temperatures. Two of the regions correspond to the Widom delta for pure supercritical carbon dioxide and ethanol, while the third region is in the immediate vicinity of the critical point of the binary mixture. The origin of these Widom deltas is a result of the large mixed linear clusters formation.
Collapse
|
10
|
Xu J, Wang Y, Ma X. Phase distribution including a bubblelike region in supercritical fluid. Phys Rev E 2021; 104:014142. [PMID: 34412334 DOI: 10.1103/physreve.104.014142] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Accepted: 07/02/2021] [Indexed: 11/07/2022]
Abstract
Pseudoboiling in supercritical fluid (SF) has been paid great attention in recent years. Available works mainly focus on thermodynamics analysis. Fewer studies were reported on the spatial time phase distribution. Here, SF is investigated in a multiphase fluid framework using molecular dynamics (MD) simulations. A simulation box contains 10 976 argon atoms, with periodic boundary conditions applied on all the box surfaces. Pressure and temperature are well controlled. Based on MD simulation results, an onset pseudoboiling temperature T^{-} and a termination pseudoboiling temperature T^{+} are defined using the neighboring molecules method, the radial distribution function method, and the two-body excess entropy method. The two transition temperatures divide the whole phase diagram into three regimes of liquidlike, two-phase-like (TPL), and gaslike, and the MD determined T^{-} and T^{+} well matched the thermodynamics-determined values. In the TPL regime, nanovoids are observed to have two distinct characteristics: (1) Particles are sparsely distributed to have gas density inside the void, but are densely populated to have liquid density outside the void. (2) Voids have a curved interface. These characteristics are very similar to bubble characteristics in subcritical pressure. Hence, voids in the supercritical state are called "bubblelike" in this paper. Nonlinear dynamics demonstrates chaotic behavior in the TPL regime, similar to the two-phase regime in the subcritical domain. The above findings give strong evidence that SF in the TPL regime consists of a mixture of bubblelike voids and surrounding liquids. Our work highlights the multiphase feature of a SF, hence, the well-established multiphase theory in subcritical pressures can be introduced to handle the complex SF.
Collapse
Affiliation(s)
- Jinliang Xu
- Beijing Key Laboratory of Multiphase Flow and Heat Transfer for Low Grade Energy Utilization, North China Electric Power University, Beijing, 102206, China.,Key Laboratory of Power Station Energy Transfer Conversion and System, North China Electric Power University, Ministry of Education, Beijing, 102206, China
| | - Yan Wang
- Beijing Key Laboratory of Multiphase Flow and Heat Transfer for Low Grade Energy Utilization, North China Electric Power University, Beijing, 102206, China
| | - Xiaojing Ma
- Beijing Key Laboratory of Multiphase Flow and Heat Transfer for Low Grade Energy Utilization, North China Electric Power University, Beijing, 102206, China
| |
Collapse
|
11
|
Servis MJ, Stephenson GB. Mesostructuring in Liquid-Liquid Extraction Organic Phases Originating from Critical Points. J Phys Chem Lett 2021; 12:5807-5812. [PMID: 34137623 DOI: 10.1021/acs.jpclett.1c01429] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Organic phase structure plays an important role in solute extraction energetics and phase behavior of liquid-liquid extraction (LLE) systems. For a binary extractant (amphiphile)/solvent mixture of relevance to LLE, we find that the organic phase mesostructuring is consistent with extractant concentration fluctuations as the compositional isotherm traverses the Widom line above its liquid-liquid critical point. This reveals a different mechanism for the well-documented heterogeneities in LLE organic phases that are typically attributed to micellization.
Collapse
Affiliation(s)
- Michael J Servis
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - G B Stephenson
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
12
|
Servis MJ, Clark AE. Cluster Identification Using Modularity Optimization to Uncover Chemical Heterogeneity in Complex Solutions. J Phys Chem A 2021; 125:3986-3993. [PMID: 33929191 DOI: 10.1021/acs.jpca.0c11320] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Structural heterogeneity is commonly manifested in solutions and liquids that feature competition of different interparticle forces. Identifying and characterizing heterogeneity across different length scales requires multimodal experimental measurement and/or the application of new techniques for the interrogation of atomistic simulation data. Within the latter, the parsing of networks of interparticle interactions (chemical networks) has been demonstrated to be a valuable tool for identifying subensembles of chemical environments. However, chemical networks can adopt a wide variety of topologies that challenge generalizable methods for identifying heterogeneous behavior, and few network analysis algorithms have been proposed for multiscale resolution. In this study, we apply a method of partitioning using the graph theoretic concept of clusters and communities. Using a modularity optimization algorithm, the cluster partition creates subgraphs based on their relative internal and external connectivities. The methodology is tested on two soft matter systems that have significantly different network topologies so as to probe its ability to identify multiple scale features and its generalizability. A binary Lennard-Jones fluid is first examined, where one component causes subgraphs that have high internal network connectivity yet are still connected to the rest of the interparticle network of interactions. The impact of connectivity and edge weighting on the cluster partition is investigated. In the second system, hierarchically organized molecular structures comprised of hydrogen bonded water molecules are identified at a liquid/liquid interface. These structures have a much more sparse network with significantly varied internal connectivity that is a challenge to differentiate from the background hydrogen bonding network of water molecules at the instantaneous interface. The organized macrostructures are effectively isolated from the background network using the cluster partition, and a time-dependent implementation allows us to reveal their reactivity. These studies indicate that cluster partitioning based upon intermolecular network connectivity patterns is broadly generalizable, depending only on user-defined intermolecular connectivity, is operable across different length scales, and is extensible to the study of dynamic phenomena.
Collapse
Affiliation(s)
- Michael J Servis
- Department of Chemistry, Washington State University, Pullman, Washington 99164, United States
| | - Aurora E Clark
- Department of Chemistry, Washington State University, Pullman, Washington 99164, United States.,Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States.,Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
13
|
Vankar H, Rana V, Dey S, Patel H, Jain V. Molecular interaction in binary mixtures of 3-Bromoanisole and methanol: A microwave dielectric relaxation spectroscopy and molecular dynamic simulation study. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.115186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
14
|
Chatwell RS, Guevara-Carrion G, Gaponenko Y, Shevtsova V, Vrabec J. Diffusion of the carbon dioxide-ethanol mixture in the extended critical region. Phys Chem Chem Phys 2021; 23:3106-3115. [PMID: 33491706 DOI: 10.1039/d0cp04985a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effect of traces of ethanol in supercritical carbon dioxide on the mixture's thermodynamic properties is studied by molecular simulations and Taylor dispersion measurements. This mixture is investigated along the isobar p = 10 MPa in the temperature range between T = 304 and 343 K. Along this path, the mixture undergoes two transitions: First, the Widom line is crossed, marking the transition from liquid-like to gas-like conditions. A second transition occurs from the supercritical gas-like domain to a subcritical gas. The Widom line crossover entails inflection points for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell-Stefan and intradiffusion coefficients. On the other hand, the transition between the super- and subcritical regions is found to be generally smooth, an observation that is qualitatively confirmed by experimental Taylor dispersion measurements. Dedicated atomistic simulations show the presence of microheterogeneities due to ethanol self-association along the investigated path, which lead to the mixture's anomalous behavior in its extended critical region.
Collapse
Affiliation(s)
- René Spencer Chatwell
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany.
| | | | - Yuri Gaponenko
- Microgravity Research Center, Université Libre de Bruxelles, 1050 Bruxelles, Belgium
| | - Valentina Shevtsova
- Microgravity Research Center, Université Libre de Bruxelles, 1050 Bruxelles, Belgium
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany.
| |
Collapse
|
15
|
Abdulagatov IM, Skripov PV. Thermodynamic and Transport Properties of Supercritical Fluids: Review of Thermodynamic Properties (Part 1). RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2021. [DOI: 10.1134/s1990793120070192] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
16
|
|
17
|
Mareev E, Semenov T, Lazarev A, Minaev N, Sviridov A, Potemkin F, Gordienko V. Optical Diagnostics of Supercritical CO 2 and CO 2-Ethanol Mixture in the Widom Delta. Molecules 2020; 25:molecules25225424. [PMID: 33228172 PMCID: PMC7699601 DOI: 10.3390/molecules25225424] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2020] [Revised: 11/16/2020] [Accepted: 11/17/2020] [Indexed: 11/16/2022] Open
Abstract
The supercritical CO2 (scCO2) is widely used as solvent and transport media in different technologies. The technological aspects of scCO2 fluid applications strongly depend on spatial–temporal fluctuations of its thermodynamic parameters. The region of these parameters’ maximal fluctuations on the p-T (pressure-temperature) diagram is called Widom delta. It has significant practical and fundamental interest. We offer an approach that combines optical measurements and molecular dynamics simulation in a wide range of pressures and temperatures. We studied the microstructure of supercritical CO2 fluid and its binary mixture with ethanol in a wide range of temperatures and pressures using molecular dynamics (MD) simulation. MD is used to retrieve a set of optical characteristics such as Raman spectra, refractive indexes and molecular refraction and was verified by appropriate experimental measurements. We demonstrated that in the Widom delta the monotonic dependence of the optical properties on the CO2 density is violated. It is caused by the rapid increase of density fluctuations and medium-sized (20–30 molecules) cluster formation. We identified the correlation between cluster parameters and optical properties of the media; in particular, it is established that the clusters in the Widom delta acts as a seed for clustering in molecular jets. MD demonstrates that the cluster formation is stronger in the supercritical CO2-ethanol mixture, where the extended binary clusters are formed; that is, the nonlinear refractive index significantly increased. The influence of the supercritical state in the cell on the formation of supersonic cluster jets is studied using the Mie scattering technique.
Collapse
Affiliation(s)
- Evgenii Mareev
- Institute of Photon Technologies of Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Pionerskaya 2, Troitsk, 108840 Moscow, Russia; (N.M.); (A.S.); (V.G.)
- Faculty of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory bld.1/2, 119991 Moscow, Russia; (T.S.); (F.P.)
- Correspondence:
| | - Timur Semenov
- Faculty of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory bld.1/2, 119991 Moscow, Russia; (T.S.); (F.P.)
- Institute on Laser and Information Technologies—Branch of the Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Svyatoozerskaya 1, Shatura, 140700 Moscow, Russia
| | - Alexander Lazarev
- Faculty of Chemistry, M. V. Lomonosov Moscow State University, Leninskie Gory bld.1/2, 119991 Moscow, Russia;
| | - Nikita Minaev
- Institute of Photon Technologies of Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Pionerskaya 2, Troitsk, 108840 Moscow, Russia; (N.M.); (A.S.); (V.G.)
| | - Alexander Sviridov
- Institute of Photon Technologies of Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Pionerskaya 2, Troitsk, 108840 Moscow, Russia; (N.M.); (A.S.); (V.G.)
| | - Fedor Potemkin
- Faculty of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory bld.1/2, 119991 Moscow, Russia; (T.S.); (F.P.)
| | - Vyacheslav Gordienko
- Institute of Photon Technologies of Federal Scientific Research Centre “Crystallography and Photonics” of Russian Academy of Sciences, Pionerskaya 2, Troitsk, 108840 Moscow, Russia; (N.M.); (A.S.); (V.G.)
- Faculty of Physics, M. V. Lomonosov Moscow State University, Leninskie Gory bld.1/2, 119991 Moscow, Russia; (T.S.); (F.P.)
| |
Collapse
|
18
|
Optimization of process parameters for preparation of polystyrene PM2.5 particles by supercritical antisolvent method using BBD-RSM. Sci Rep 2020; 10:11187. [PMID: 32636462 PMCID: PMC7341738 DOI: 10.1038/s41598-020-67994-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2019] [Accepted: 06/17/2020] [Indexed: 11/18/2022] Open
Abstract
The objective of this study is to optimize the process parameters for preparing polystyrene (PS) PM2.5 particles by supercritical antisolvent (SAS) method. Toluene was selected as the solvent and supercritical carbon dioxide (SC-CO2) was used as the antisolvent. The Box–Behnken design-response surface method was applied to investigate the effect of crystallizer pressure, PS massic concentration, flow ratio of CO2/solution and crystallizer temperature on the size and the distribution of PS particles, systematically. It is found that crystallizer temperature is the most significant variable on the size and the distribution of PS particles, followed by flow ratio of CO2/solution and PS massic concentration, and crystallizer pressure is the slightest significant factor. The particle size increases with the increase of crystallizer temperature. The optimum conditions are obtained as crystallizer pressure 9.8 MPa, PS massic concentration 1.6 wt%, flow ratio of CO2/solution 140 g/g and crystallizer temperature 309 K. Under these conditions, the PS particle with the size of 2.78 μm and a narrow size distribution has been prepared, meeting PM2.5 standard aerosols. The results suggest that it is feasible to produce PM2.5 standard aerosols by SAS.
Collapse
|
19
|
Utilization of CO2 as Cushion Gas for Depleted Gas Reservoir Transformed Gas Storage Reservoir. ENERGIES 2020. [DOI: 10.3390/en13030576] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Underground gas storage reservoirs (UGSRs) are used to keep the natural gas supply smooth. Native natural gas is commonly used as cushion gas to maintain the reservoir pressure and cannot be extracted in the depleted gas reservoir transformed UGSR, which leads to wasting huge amounts of this natural energy resource. CO2 is an alternative gas to avoid this particular issue. However, the mixing of CO2 and CH4 in the UGSR challenges the application of CO2 as cushion gas. In this work, the Donghae gas reservoir is used to investigate the suitability of using CO2 as cushion gas in depleted gas reservoir transformed UGSR. The impact of the geological and engineering parameters, including the CO2 fraction for cushion gas, reservoir temperature, reservoir permeability, residual water and production rate, on the reservoir pressure, gas mixing behavior, and CO2 production are analyzed detailly based on the 15 years cyclic gas injection and production. The results showed that the maximum accepted CO2 concentration for cushion gas is 9% under the condition of production and injection for 120 d and 180 d in a production cycle at a rate of 4.05 kg/s and 2.7 kg/s, respectively. The typical curve of the mixing zone thickness can be divided into four stages, which include the increasing stage, the smooth stage, the suddenly increasing stage, and the periodic change stage. In the periodic change stage, the mixed zone increases with the increasing of CO2 fraction, temperature, production rate, and the decreasing of permeability and water saturation. The CO2 fraction in cushion gas, reservoir permeability, and production rate have a significant effect on the breakthrough of CO2 in the production well, while the effect of water saturation and temperature is limited.
Collapse
|
20
|
Zerón I, Torres-Arenas J, de Jesús E, Ramírez B, Benavides A. Discrete potential fluids in the supercritical region. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111518] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
21
|
Karalis K, Ludwig C, Niceno B. Supercritical water anomalies in the vicinity of the Widom line. Sci Rep 2019; 9:15731. [PMID: 31673024 PMCID: PMC6823507 DOI: 10.1038/s41598-019-51843-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 09/05/2019] [Indexed: 12/21/2022] Open
Abstract
Supercritical water is used in a variety of chemical and industrial applications. As a consequence, a detailed knowledge of the structure-properties correlations is of uttermost importance. Although supercritical water was considered as a homogeneous fluid, recent studies revealed an anomalous behaviour due to nanoscale density fluctuations (inhomogeneity). The inhomogeneity is clearly demarked through the Widom line (maxima in response factions) and drastically affect the properties. In the current study the physical properties of supercritical water have been determined by classical molecular dynamics simulations using a variety of polarized and polarizable interatomic potentials. Their validity which was not available at supercritical conditions has been assessed based on the ability to reproduce experimental data. Overall, the polarized TIP4P/2005 model accurately predicted the properties of water in both liquid-like and gas-like regions. All interatomic potentials captured the anomalous behaviour providing a direct evidence of molecular-scale inhomogeneity.
Collapse
Affiliation(s)
- Konstantinos Karalis
- Laboratory for Scientific Computing and Modelling (LSM), NES Division, Paul Scherrer Institute, 5232, Villigen PSI, Switzerland.
| | - Christian Ludwig
- Laboratory for Bioenergy and Catalysis (LBK), ENE Division, Paul Scherrer Institute, 5232, Villigen PSI, Switzerland
- École Polytechnique Fédérale de Lausanne (EPFL), ENAC IIE GR-LUD, 1015, Lausanne, Switzerland
| | - Bojan Niceno
- Laboratory for Scientific Computing and Modelling (LSM), NES Division, Paul Scherrer Institute, 5232, Villigen PSI, Switzerland.
- Eidgenössische Technische Hochschule Zürich (ETHZ), MAVT-LKE, 8092, Zurich, Switzerland.
| |
Collapse
|
22
|
|
23
|
Aasen A, Hammer M, Ervik Å, Müller EA, Wilhelmsen Ø. Equation of state and force fields for Feynman–Hibbs-corrected Mie fluids. I. Application to pure helium, neon, hydrogen, and deuterium. J Chem Phys 2019. [DOI: 10.1063/1.5111364] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Ailo Aasen
- Department of Energy and Process Engineering, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
- SINTEF Energy Research, NO-7465 Trondheim, Norway
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Morten Hammer
- SINTEF Energy Research, NO-7465 Trondheim, Norway
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Åsmund Ervik
- SINTEF Energy Research, NO-7465 Trondheim, Norway
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Øivind Wilhelmsen
- Department of Energy and Process Engineering, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
- SINTEF Energy Research, NO-7465 Trondheim, Norway
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
24
|
Kryuchkov NP, Mistryukova LA, Brazhkin VV, Yurchenko SO. Excitation spectra in fluids: How to analyze them properly. Sci Rep 2019; 9:10483. [PMID: 31324848 PMCID: PMC6642218 DOI: 10.1038/s41598-019-46979-y] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Accepted: 07/09/2019] [Indexed: 11/25/2022] Open
Abstract
Although the understanding of excitation spectra in fluids is of great importance, it is still unclear how different methods of spectral analysis agree with each other and which of them is suitable in a wide range of parameters. Here, we show that the problem can be solved using a two-oscillator model to analyze total velocity current spectra, while other considered methods, including analysis of the spectral maxima and single mode analysis, yield rough results and become unsuitable at high temperatures and wavenumbers. To prove this, we perform molecular dynamics (MD) simulations and calculate excitation spectra in Lennard-Jones and inverse-power-law fluids at different temperatures, both in 3D and 2D cases. Then, we analyze relations between thermodynamic and dynamic features of fluids at (Frenkel) crossover from a liquid- to gas-like state and find that they agree with each other in the 3D case and strongly disagree in 2D systems due to enhanced anharmonicity effects. The results provide a significant advance in methods for detail analysis of collective fluid dynamics spanning fields from soft condensed matter to strongly coupled plasmas.
Collapse
Affiliation(s)
- Nikita P Kryuchkov
- Bauman Moscow State Technical University, 2nd Baumanskaya street 5, Moscow, 105005, Russia
| | - Lukiya A Mistryukova
- Bauman Moscow State Technical University, 2nd Baumanskaya street 5, Moscow, 105005, Russia
| | - Vadim V Brazhkin
- Institute for High Pressure Physics RAS, Kaluzhskoe shosse, 14, Troitsk, Moscow, 108840, Russia
| | - Stanislav O Yurchenko
- Bauman Moscow State Technical University, 2nd Baumanskaya street 5, Moscow, 105005, Russia.
| |
Collapse
|
25
|
Ploetz EA, Smith PE. Gas or Liquid? The Supercritical Behavior of Pure Fluids. J Phys Chem B 2019; 123:6554-6563. [DOI: 10.1021/acs.jpcb.9b04058] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Elizabeth A. Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid Campus Dr. North, Manhattan, Kansas 66506-0401, United States
| | - Paul E. Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid Campus Dr. North, Manhattan, Kansas 66506-0401, United States
| |
Collapse
|
26
|
Guevara-Carrion G, Ancherbak S, Mialdun A, Vrabec J, Shevtsova V. Diffusion of methane in supercritical carbon dioxide across the Widom line. Sci Rep 2019; 9:8466. [PMID: 31186475 PMCID: PMC6560060 DOI: 10.1038/s41598-019-44687-1] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 05/22/2019] [Indexed: 11/16/2022] Open
Abstract
Diffusion of methane diluted in supercritical carbon dioxide is studied by experiment and molecular simulation in the temperature range from 292.55 to 332.85 K along the isobars 9.0, 12.5 and 14.7 MPa. Measurements of the Fick diffusion coefficient are carried out with the Taylor dispersion technique. Molecular dynamics simulation and the Green-Kubo formalism are employed to obtain Fick, Maxwell-Stefan and intradiffusion coefficients as well as shear viscosity. The obtained diffusion coefficients are on the order of 10-8 m2/s. The composition, temperature and density dependence of diffusion is analyzed. The Fick diffusion coefficient of methane in carbon dioxide shows an anomaly in the near-critical region. This behavior can be attributed to the crossing of the so-called Widom line, where the supercritical fluid goes through a transition between liquid-like and gas-like states. Further, several classical equations are tested on their ability to predict this behavior and it is found that equations that explicitly include the density are better suited to predict the sharp variation of the diffusion coefficient near the critical region predicted by molecular simulation.
Collapse
Affiliation(s)
- Gabriela Guevara-Carrion
- Thermodynamics and Process Engineering, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany
| | - Sergiy Ancherbak
- MRC, CP165/62, Université Libre de Bruxelles, Av. F. D. Roosevelt, 50, B-1050, Brussels, Belgium
| | - Aliaksandr Mialdun
- MRC, CP165/62, Université Libre de Bruxelles, Av. F. D. Roosevelt, 50, B-1050, Brussels, Belgium
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany.
| | - Valentina Shevtsova
- MRC, CP165/62, Université Libre de Bruxelles, Av. F. D. Roosevelt, 50, B-1050, Brussels, Belgium
| |
Collapse
|
27
|
Yoon TJ, Ha MY, Lee WB, Lee YW. A corresponding-state framework for the structural transition of supercritical fluids across the Widom delta. J Chem Phys 2019; 150:154503. [DOI: 10.1063/1.5086467] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Tae Jun Yoon
- School of Chemical and Biological Engineering, Institute of Chemical Processes, Seoul National University, Seoul 08826, South Korea
| | - Min Young Ha
- School of Chemical and Biological Engineering, Institute of Chemical Processes, Seoul National University, Seoul 08826, South Korea
| | - Won Bo Lee
- School of Chemical and Biological Engineering, Institute of Chemical Processes, Seoul National University, Seoul 08826, South Korea
| | - Youn-Woo Lee
- School of Chemical and Biological Engineering, Institute of Chemical Processes, Seoul National University, Seoul 08826, South Korea
| |
Collapse
|
28
|
Takemoto A, Kinugawa K. Quantumness and state boundaries hidden in supercritical helium-4: A path integral centroid molecular dynamics study. J Chem Phys 2018; 149:204504. [DOI: 10.1063/1.5053988] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ayumi Takemoto
- Division of Chemistry, Graduate School of Humanities and Sciences, Nara Women’s University, Nara 630-8506, Japan
| | - Kenichi Kinugawa
- Division of Chemistry, Graduate School of Humanities and Sciences, Nara Women’s University, Nara 630-8506, Japan
| |
Collapse
|
29
|
de Haro ML, Rodríguez-Rivas Á, Yuste SB, Santos A. Structural properties of the Jagla fluid. Phys Rev E 2018; 98:012138. [PMID: 30110748 DOI: 10.1103/physreve.98.012138] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Indexed: 06/08/2023]
Abstract
The structural properties of the Jagla fluid are studied by Monte Carlo (MC) simulations, numerical solutions of integral equation theories, and the (semi-analytical) rational-function approximation (RFA) method. In the latter case, the results are obtained from the assumption (supported by our MC simulations) that the Jagla potential and a potential with a hard core plus an appropriate piecewise constant function lead to practically the same cavity function. The predictions obtained for the radial distribution function, g(r), from this approach are compared against MC simulations and integral equations for the Jagla model, and also for the limiting cases of the triangle-well potential and the ramp potential, with a general good agreement. The analytical form of the RFA in Laplace space allows us to describe the asymptotic behavior of g(r) in a clean way and compare it with MC simulations for representative states with oscillatory or monotonic decay. The RFA predictions for the Fisher-Widom and Widom lines of the Jagla fluid are obtained.
Collapse
Affiliation(s)
- Mariano López de Haro
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México (U.N.A.M.), Temixco, Morelos 62580, Mexico
| | - Álvaro Rodríguez-Rivas
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| | - Santos B Yuste
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| | - Andrés Santos
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| |
Collapse
|
30
|
Yoon TJ, Lee YW. Current theoretical opinions and perspectives on the fundamental description of supercritical fluids. J Supercrit Fluids 2018. [DOI: 10.1016/j.supflu.2017.11.022] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
31
|
Bryk T, Gorelli FA, Mryglod I, Ruocco G, Santoro M, Scopigno T. Behavior of Supercritical Fluids across the "Frenkel Line". J Phys Chem Lett 2017; 8:4995-5001. [PMID: 28945381 DOI: 10.1021/acs.jpclett.7b02176] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The "Frenkel line" (FL), the thermodynamic locus where the time for a particle to move by its size equals the shortest transverse oscillation period, has been proposed as a boundary between recently discovered liquid-like and gas-like regions in supercritical fluids. We report a simulation study of isothermal supercritical neon in a range of densities intersecting the FL. Specifically, structural properties and single-particle and collective dynamics are scrutinized to unveil the onset of any anomalous behavior at the FL. We find that (i) the pair distribution function smoothly evolves across the FL displaying medium-range order, (ii) low-frequency transverse excitations are observed below the "Frenkel frequency", and (iii) the high-frequency shear modulus does not vanish even for low-density fluids, indicating that positive sound dispersion characterizing the liquid-like region of the supercritical state is unrelated to transverse dynamics. These facts critically undermine the definition of the FL and its significance for any relevant partition of the supercritical phase.
Collapse
Affiliation(s)
- T Bryk
- Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine , 1 Svientsitskii Street, UA-79011 Lviv, Ukraine
- Institute of Applied Mathematics and Fundamental Sciences, Lviv Polytechnic National University , UA-79013 Lviv, Ukraine
| | - F A Gorelli
- Istituto Nazionale di Ottica INO-CNR , I-50019 Sesto Fiorentino, Italy
- European Laboratory for Non Linear Spectroscopy, LENS , I-50019 Sesto Fiorentino, Italy
| | - I Mryglod
- Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine , 1 Svientsitskii Street, UA-79011 Lviv, Ukraine
| | - G Ruocco
- Dipartimento di Fisica, Universita di Roma La Sapienza , I-00185 Roma, Italy
- Center for Life Nano Science @Sapienza, Istituto Italiano di Tecnologia , 295 Viale Regina Elena, I-00161 Roma, Italy
| | - M Santoro
- Istituto Nazionale di Ottica INO-CNR , I-50019 Sesto Fiorentino, Italy
- European Laboratory for Non Linear Spectroscopy, LENS , I-50019 Sesto Fiorentino, Italy
| | - T Scopigno
- Dipartimento di Fisica, Universita di Roma La Sapienza , I-00185 Roma, Italy
- Center for Life Nano Science @Sapienza, Istituto Italiano di Tecnologia , 295 Viale Regina Elena, I-00161 Roma, Italy
| |
Collapse
|