Absorption of pressurized methane in normal and supercooled p

Scientific Reports volume 13, Article number: 136 (2023) Cite this article 636 Accesses 1 Altmetric Metrics details Supercooling of liquids leads to peculiarities which are scarcely studied under

Scientific Reports volume 13, Article number: 136 (2023) Cite this article

636 Accesses

1 Altmetric

Metrics details

Supercooling of liquids leads to peculiarities which are scarcely studied under high-pressure conditions. Here, we report the surface tension, solubility, diffusivity, and partial molar volume for normal and supercooled liquid solutions of methane with p-xylene. Liquid bodies of perdeuterated p-xylene (p-C8D10), and, for comparison, o-xylene (o-C8D10), were exposed to pressurized methane (CH4, up to 101 bar) at temperatures ranging 7.0–30.0 °C and observed at high spatial resolution (pixel size 20.3 μm) using a non-tactile neutron imaging method. Supercooling led to the increase of diffusivity and partial molar volume of methane. Solubility and surface tension were insensitive to supercooling, the latter substantially depended on methane pressure. Overall, neutron imaging enabled to reveal and quantify multiple phenomena occurring in supercooled liquid p-xylene solutions of methane under pressures relevant to the freeze-out in the production of liquefied natural gas.

Liquids cooled below their melting points show peculiar physical properties, such as changes in viscosity and diffusivity, which are commonly ascribed to the dynamic heterogeneity on the molecular level1,2,3,4,5,6. Thermal movements of molecules slow down either due to energy barriers or lack of free volume1,7,8,9. Specific interactions, such as hydrogen bonds and interactions among aromatic rings, can be involved. For instance, the formation of ice-like clusters in supercooled water was hypothesized since organic compounds, such as o-, m-, p-xylene, are volatilized (Henry’s constants increase) from their supercooled solutions with water10. Besides water, specific interactions likely influence the properties of other supercooled liquids, such as p-xylene. This compound, which solidifies at 0.84 GPa at room temperature, was reported to form trimers and tetramers at 13.5 GPa11. Thus, the molecular-level heterogeneity of supercooled liquid p-xylene can be expected to occur and to be involved in the spontaneous solidification and hardly predictable properties due to the thermodynamic instability.

The high normal melting point of p-xylene (p-C8H10, 13.25 °C12) causes that this compound can condense, become supercooled and deposit on cold spots in the production of liquefied natural gas (LNG). The other xylene isomers are much less severe, the normal melting points12 are: − 25.17 °C (o-C8H10) and − 47.85 °C (m-C8H10). Besides that, p-xylene is a practical test compound other than water for studying supercooled liquids at rather mild temperatures. Although condensation of the BTEX compounds (benzene, toluene, ethylbenzene, xylenes) in the production of LNG is avoided by their low admitted concentration (< 1 ppm13), this study can bring a better understanding of phenomena occurring on cold spots in the processing of natural gas and broaden the general knowledge.

Recent studies have reported vital experimental data and predictive models describing compositions of coexisting phases13,14 and, importantly, the melting temperature for p-xylene (p-C8H10) in systems comprising of methane (CH4) and p-xylene (p-C8H10)14. Thus, the degree of supercooling can be assessed not only for pure p-xylene but also for its mixtures with methane at elevated pressures. The achievable supercooling generally depends on the chemical composition and surfaces in the experimental device. The supercooling by about 2 °C was reported to induce p-xylene solidification at a cooled copper tip for methane (CH4) and p-xylene (p-C8H10) solution at 4.35 °C and 225 bar14. In contrast to the literature14, in which solidification was induced, we aim at studying the properties of the thermodynamically unstable supercooled liquid.

Experimental data on density, solubility, speed of sound, heat capacity, surface tension and viscosity have been so far reported for several supercooled liquids, chiefly water3,10,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29. With the exception of one investigation on the speed of sound and derived quantities for supercooled water16, the above studies report data observed at pressures near atmospheric or do not report experimental pressure. The common high-throughput methods for studying liquids under high-pressure conditions are the vibrating tube densimetry, pendant drop method, Taylor dispersion method, the method of capillary waves, methods utilizing Raman spectroscopy and Nuclear Magnetic Resonance30,31,32,33,34,35,36,37,38,39,40,41. To our knowledge, no report on their use for supercooled liquids under high pressures is available. We certainly admit that these or other methods can be applied to study the properties of supercooled liquids. For instance, Raman spectroscopy and Nuclear Magnetic Resonance were used for studies on the formation of natural gas hydrate under relevant conditions42,43. As we show in this study, our non-tactile one-pot neutron imaging method44 is applicable for studying systems involving supercooled liquids exposed to pressurized gases.

In this work, methane (CH4) absorption in liquid perdeuterated p-xylene (p-C8D10) and o-xylene (o-C8D10) was studied while this choice of isotopic composition exploits the high neutronic contrast45 between protium (H) and deuterium (D). The influence of the isotopic composition on the physical properties of chemical species has been thoroughly reported for benzene rather than for the xylenes and appears low. For instance, the molar volume (molar mass over density) of perdeuterated benzene (C6D6) differs by less than 0.24% from that of benzene (C6H6) under conditions relevant for this study46; see Supplementary Information (SI) for more discussion. Viscosity, melting point, boiling point and surface tension of perdeuterated benzene differ from those of benzene by 5%47, 1.0 °C48, 0.8 K49 and 2% (− 0.5 mN·m−1)50, respectively. Systematical errors caused by using deuterated xylenes instead of the protium-based (normal) ones are expected to be comparable to those for benzene.

Liquids were maintained at a constant temperature and exposed to a methane pressure step at the zero time in cylindrical cells44. Since the cells were axially symmetric (inner diameter 9.0 ± 0.1 mm), the onion-peeling algorithm51 was used to provide the tomographic reconstructions at the central plane of the sample. The overall linear attenuation coefficient of the liquid (Σ) was contributed by the two components, A (CH4) and B (p- or o-C8D10). The Beer-Lambert law thus has the form

Cross-sections (σ) of the pure components were evaluated based on the tomographic reconstruction observed just after the release of pressurized methane to the vessel with sample liquid by assuming negligible evaporation of the liquid and diffusion of methane to the bottom part of the liquid body (Fig. 1). Symbols I, N0 and d are intensity, Avogadro number and length, respectively. The mole concentration (c) and density (\(\rho\)) of pure methane in the gas and supercritical fluid phase were calculated using the Peng-Robinson equation of state52. Evaporation of the xylenes was neglected53,54; see SI for the assessment of the combined systematic uncertainty. For the liquid perdeuterated xylenes, mole concentration and density were calculated from the known equations of state for the pure p- and o-xylene (p- and o-C8H10) liquids46 by assuming that the molar volume of the deuterated and protium-based chemicals are equal (see SI for uncertainty estimation). Cross-sections and densities for the pure components at the studied conditions are listed in Table S1 in SI. In the case of supercooled p-xylene, density was calculated by extrapolating from the region of (normal) liquid. Since the cross-section of p-xylene held constant over all the inspected conditions, such extrapolation provided a meaningful approximation and supercooled p-xylene did not solidify as the solidification of p-xylene is accompanied with the change of density by about 20%12. The level of supercooling was assessed as the difference between the melting temperature and the actual temperature using the available literature data for the melting point of the protium-based p-xylene (p-C8H10) under relevant conditions, see SI and Materials and Methods. It is substantial to note that the liquid swelled during methane absorption, which did not enable for the simplistic use of Beer-Lambert law, Eq. (1).

Tomographic reconstructions of the measuring cell with supercooled perdeuterated p-xylene (p-C8D10) exposed to methane at 1.0 bar (A) and at 100.4 ± 0.2 bar at 7.0 ± 0.2 °C (A)–(D). Time after exposure to pressurized methane, fits of the phase interface shapes (red curves), local level of supercooling (sc), axes and cell radial dimension (blue lines, 2·R = 9.0 ± 0.1 mm), and domains used for the calculation of the linear attenuation coefficient of the pure components (purple boxes) are indicated. Green dashed line is the liquid level at the cell center short after the pressurization. The bottom-left insert shows the body of revolution used for the calculation of the liquid volume.

The shape of the mobile axially symmetric phase interface in gravity (red curves in Fig. 1) was parameterized by solving the Young–Laplace equation in form55,56

The boundary value problem, z'(r = 0) = 0 and z'(r = R) = cot(θ), was solved using the midrich method as in the Maple 2021 software package, r ∈ (0, R). Optimum values of the surface tension (γ) and contact angle at the wall of the cell (θ), and their uncertainty due to random errors (ur, cover factor 2) were calculated using Gauss–Newton and Bonferroni methods57,58, combined systematic uncertainty of surface tension was calculated using the law of uncertainty propagation (see SI). The density difference at the phase interface, ∆ρ, was calculated as follows.

The volume of the liquid (V) was calculated by computing the volume of the solid of revolution (see insert in Fig. 1), which is a generalization of an earlier method utilizing projections (photography in visible light) of phase interfaces in glass tubes without using the reconstruction of the central plane59. The partial molar volume of methane (\(\overline{V}_{{\text{A}}} = \partial V/\partial n_{{\text{A}}}\)) and its uncertainty due to random errors (ur, cover factor 2) were calculated based on the mole amount of absorbed methane (\(n_{{\text{A}}}\)) in the entire liquid body and its volume at fixed T, p, and mole amount of the perdeuterated xylene (\(n_{{\text{B}}}\)), see Fig. 2A and Table S4 in SI. The molar volume (\(V_{{\text{m}}}\)) of the liquid and its density (\(\rho\)) depend on the mole fractions [\(x_{{\text{A}}} = n_{{\text{A}}} /\left( {n_{{\text{A}}} + n_{{\text{B}}} } \right)\), \(x_{{\text{B}}} = 1 - x_{{\text{A}}}\)] and molar masses (M) such that

where the partial molar volume of each perdeuterated liquid xylene was set to its molar volume.

(A) Plot of volume of the liquid in the cell against the mole amount of absorbed methane is shown together with the linear fit and its confidence interval. (B) Time-evolution of the measured radius-averaged relative concentration of methane at the phase interface is shown with an empirical model (the blue curve is the plot of C∞ − C∞⋅exp[− τ / (91.1 s)]). The red curve is the derived model of density, insert shows the relative concentration of methane at the interface with respect to the cell radius for time 1043 s. Both examples are for supercooled perdeuterated p-xylene (p-C8D10) at 7.0 ± 0.2 °C initially exposed to methane at 100.4 ± 0.2 bar.

The simultaneous measurement of the neutron attenuation and the shape of the liquid (and thus volume) enabled to derive the concentration distributions of both species (Fig. 3A) using Eqs. (1) and (3), constant mole amount of the xylene in the liquid was assumed. This, in turn, enabled to extrapolate the methane concentration in the xylene phase at the phase interface and thus to derive the respective density (Figs. 2B, 3A) needed for the Young–Laplace equation, Eq. (2). Moreover, the z coordinate (depicted in Fig. 1B) was transformed to the B-fixed coordinate \(\xi\) (Fig. 3B) along which mole concentration of the perdeuterated xylene (\({c}_{\mathrm{B}}\)) is fixed.

(A) Experimental relative (C) and mole (c) concentrations of methane, and the derived mole concentration of p-xylene (orange curve) are plotted against the depth coordinate (z) together with a guide for the eye (black curve) used for the extrapolation to the interface. (B) Experimental (blue crosses) and modelled (orange points, Fick’s second law) relative concentrations of methane are plotted against the B-fixed depth coordinate (\(\xi\)) together with the \(\xi /z\) ratio (green curve). Phase interface is indicated by vertical dashed lines, both figures are for supercooled perdeuterated p-xylene (p-C8D10) at 7.0 ± 0.2 °C initially exposed to methane at 100.4 ± 0.2 bar, radius coordinate and time upon pressurization are indicated.

Diffusion of methane in the axially symmetric liquid body (Fig. 1) was modelled using Fick’s second law in cylindrical B-fixed coordinates in form60

where D is diffusivity in the B-fixed reference frame and \(\xi\) ranges from zero to the initial liquid level at a given radius. Concentration at the phase interface (Fig. 2B) and its shape (Fig. 1) were used for the construction of the Dirichlet boundary condition; impermeable walls of the cell were represented by Neumann boundary conditions. We have solved Eq. (4) using an explicit differentiation scheme58, and calculated optimum value and uncertainty due to random errors (ur, cover factor 2) of diffusivity (D) using Gauss–Newton and Bonferroni methods57,58. Fick’s second law, Eq. (4), provided a good approximation of the experimental data (Fig. 3B).

The surface tension of the binary methane solutions with perdeuterated p-xylene (p-C8D10) and o-xylene (o-C8D10) showed a mild dependence on temperature and a strong dependence on (methane) pressure that followed one master trend irrespective of the supercooling and of the actual xylene isomer (Fig. 4, Table S3 in SI). Surface tension measured for the studied perdeuterated xylenes saturated with methane at 1.0 bar followed the correlations from the database for the pure protium-based xylenes12. Interestingly, no influence of the ongoing methane diffusion through the phase interface (Fig. 5A) on the surface tension was detected within the experimental uncertainty. Neither the liquid supercooling nor the actual isomer form of xylene (o- and p-) thus measurably influenced the methane adsorption on the phase interface, which appears insensitive to the concentration gradient in the liquid. Thus, the measurement of interfacial tension among liquid and gas (or supercritical fluid) does not apparently necessitate reaching the phase equilibrium.

(A) Surface tension of p-xylene (p-C8D10) and o-xylene (o-C8D10) plotted against methane pressure and (B) against temperature for a fixed (atmospheric) pressure. Typical uncertainty of surface tension due to random errors (ur, cover factor 2) is 1–2 mN⋅m−1; experimental data are in Table S3 in SI. Curves represent guides for the eye (A) and correlations for protium-based xylenes taken from the database12; the blue part of the curve is the extrapolation for the supercooled liquid at 1.0 bar (B).

(A) Surface tension as a function of time upon exposure to pressurized methane, error bars show uncertainty due to random errors calculated using Bonferroni method57, line shows the average. (B) Henry’s law constant, see Eq. (5), of methane in the studied deuterated xylenes is shown together with the temperature dependence of Henry’s the constant for all data (this work) for perdeuterated p-xylene with the enthalpy of dissolution of − 4.5 kJ⋅mol−1 and datum calculated based on the literature14 using Eq. (5); sc abbreviates supercooled.

Methane solubility in the xylene phase at the interface was expressed using Henry’s law

in which H is the Henry’s law constant; see61 for more discussion. The Henry’s constant for methane and p-xylene showed mild dependences on pressure and temperature, no influence of the liquid supercooling was observed to within the experimental uncertainty (Fig. 5B, Table S4 in SI, rel. uc(\(H\)) ≈ 15%). Our data differed by about 4% from the only available literature datum14 for the Henry’s constant of methane in its solution with p-xylene (p-C8H10) at relevant conditions. By seeing literature data for supercooled water solutions of volatile compounds10, it can be expected that a higher Henry’s constant would be observed for more supercooled solutions of methane and p-xylene (p-C8H10).

Henry’s constant of methane in o-xylene (o-C8D10) was by about 59 bar higher than that in p-xylene (p-C8D10). A similar dependence of the Henry’s constant on the xylene isomerism was reported in the literature62 for several hydrocarbon gases dissolved in o-xylene (o-C8H10) and in p-xylene (p-C8H10).

Diffusivity and partial molar volume of methane both showed positive deviations from the master trends due to the supercooling of p-xylene, while no such deviations were found for control experiments with o-xylene that becomes supercooled at much lower temperatures (Fig. 6, Table S4 in SI). The higher deviations of both quantities were observed at the lowest explored pressure and temperature (average 44.8 bar, 7.0 °C). This suggests that i) supercooling of liquid p-xylene presumably leads to the molecular-level heterogeneity that facilitates diffusion of methane and reduces the free volume accessible to the methane molecules, ii) the dissolved methane disrupts this heterogeneity.

(A) Diffusivity and (B) partial molar volume of methane for its binary solutions with perdeuterated p-xylene (p-C8D10) and o-xylene (o-C8D10). Typical uncertainties of diffusivity and partial molar volume due to random error (ur, cover factor 2) are 0.3·10–9 m2 s−1 and 7 cm3 mol−1, respectively. Curves in (A) represent the Wilke–Chang model61,63 with parameters from the database12 (viscosity was temperature-extrapolated in the case of p-xylene supercooled at 1.0 bar) and with the association factors adjusted to 1.5. Green line in (B) shows mean over all non-supercooled systems (this work), literature data for the partial molar volume of methane in n-hexane at infinite dilution are shown64. Experimental data are listed in Table S4 in SI, sc abbreviates supercooled.

The supercooling of liquid mixtures of p-xylene (p-C8D10) and methane (CH4) led to the peculiar increase of the methane partial molar volume and diffusivity under pressures relevant to the p-xylene freeze-out in the production of liquefied natural gas (7.0–30.0 °C, 1.0–101.1 bar). Thus, methane diffused more readily in the supercooled solutions, which also showed higher swelling than the normal ones. Systems involving o-xylene (o-C8D10) and not showing supercooling were studied for reference. Surface tension was influenced by temperature and pressure and was sensitive neither to the liquid supercooling nor to the isomerism of xylene (ortho, para), no impact of supercooling on the methane solubility was discerned. Our inherently non-tactile neutron imaging method enabled to observe the supercooled liquid bodies and to derive information on their composition and shape at a high spatial resolution (pixel size 20.3 μm). The provision of experimental data at conditions inherently relevant to applications can be fruitful for the community studying supercooled liquids on the molecular level.

The neutron imaging experiments were conducted at the NEUTRA beamline65 at Paul Scherrer Institute at the measuring position No. 2 (L/D = 365) by MIDI-box detector system using a 30 µm-thick Gd2O2S:Tb scintillator screen and a sCMOS camera (Andor Neo) fitted with a 100-mm objective (Zeiss Makro-Planar). This experimental arrangement yielded the image of 2560 (W) × 2160 (H) pixels in size with an isotropic pixel pitch of 20.3 µm. The mean distance of the sample to the detector has been equal to 23 mm. The spatial resolution of the resulting images is therefore estimated to be better than 80 µm (based on the intrinsic detector resolution and the beam geometrical unsharpness).The acquisition scheme of the neutron radiographies consisted of several (usually seven) series of 50 images each of the 10 s acquisition time for each investigated system. The raw data were subject to the open beam and the black body corrections66,67. For the evaluation of the data from the first two respective series, ten data points were provided as an average of 10 images having the respective time stamp of the average time of the respective 10 images; for the latter series, the entire 50 images were averaged into a single data point having the time stamp of the average of the 50 images.

The experimental setup was described in our previous study44 and consisted of a pair of equivalent axially symmetric titanium measuring cells from which one was used for this study were placed in a duralumin block maintained at a constant temperature to within ± 0.2 °C using a Julabo F12-MA water circulator, the temperature was measured using a platinum resistance thermometer (Greissinger GMH 3710), the setup was placed in a duralumin box purged with nitrogen to avoid moisture condensation on the outer parts of the measuring setup. The cells were rinsed with acetone, vacuumed (< 0.01 Pa, Leybold D4B), and twice washed with the fresh sample liquid prior to the filling. One cell contained p-xylene (p-C8D10), the other o-xylene (o-C8D10) or a heterogeneous system consisting of water (the bottom phase) and p-xylene (p-C8D10); data for the latter system will be published elsewhere. Chemicals and gases were used as obtained and are listed in Table 1. The interior of the apparatus was then purged, and the liquids were separately bubbled first with nitrogen and second with methane at atmospheric pressure. Experiments included the steep change of methane pressure from atmospheric to a fixed value, which was maintained constant using a pressure reducer RSD 1 (Siad) to within 0.2 bar and sensed using a PXM409-175BAV (Omega) with a DP41-B control unit (Omega). Experiments were terminated by bubbling the liquid with methane at atmospheric pressure, which was done at 15.0 °C in the case of experiments below the normal melting point of p-xylene.

The steep change of methane pressure in the cell inherently causes a shift of the melting temperature of the initially pure liquid p-xylene (p-C8D10), which then changes as methane diffuses into the liquid. The degree of supercooling thus changes as methane diffuses through the liquid (Fig. 1 and Table S3 in SI). The melting temperature of perdeuterated p-xylene (p-C8D10) and p-xylene (p-C8H10) were measured by immersing two ampules with the solidified compounds at atmospheric pressure (rest was air) into the water bath of the Julabo F12-MA water circulator, the temperature was measured using a platinum resistance thermometer (Greissinger GMH 3710). The normal (at 1.0 bar) melting temperature of perdeuterated p-xylene (p-C8D10) was 14.1 ± 0.2 °C and that of p-xylene (p-C8H10) was 13.3 ± 0.2 °C; the latter agrees with the value from the literature for p-xylene (p-C8H10), 13.25 °C12. Hence, the melting temperatures of perdeuterated p-xylene (p-C8D10) were considered to be by higher 0.8 °C than those for p-xylene (p-C8H10)14 also in the case of the solutions with methane. The pressure dependence of the melting point of pure p-xylene (p-C8H10) was taken from the literature68, see Table S3 in SI for details.

Experimental data are listed in Supplementary Information.

Ediger, M. D. & Harrowell, P. Perspective: Supercooled liquids and glasses. J. Chem. Phys. 137, 080901. https://doi.org/10.1063/1.4747326 (2012).

Article ADS Google Scholar

Kob, W., Donati, C., Plimpton, S. J., Poole, P. H. & Glotzer, S. C. Dynamical Heterogeneities in a supercooled Lennard-Jones liquid. Phys. Rev. Lett. 79, 2827–2830. https://doi.org/10.1103/PhysRevLett.79.2827 (1997).

Article ADS Google Scholar

Debenedetti, P. G. Supercooled and glassy water. J. Phys. Condens. Matter 15, R1669–R1726. https://doi.org/10.1088/0953-8984/15/45/r01 (2003).

Article ADS Google Scholar

Scala, A., Starr, F. W., La Nave, E., Sciortino, F. & Stanley, H. E. Configurational entropy and diffusivity of supercooled water. Nature 406, 166–169. https://doi.org/10.1038/35018034 (2000).

Article ADS Google Scholar

Sosso, G. C. et al. Crystal nucleation in liquids: Open questions and future challenges in molecular dynamics simulations. Chem. Rev. 116, 7078–7116. https://doi.org/10.1021/acs.chemrev.5b00744 (2016).

Article Google Scholar

Gallo, P. et al. Advances in the study of supercooled water. Eur. Phys. J. E 44, 143. https://doi.org/10.1140/epje/s10189-021-00139-1 (2021).

Article Google Scholar

Del Gado, E., Ilg, P., Kröger, M. & Öttinger, H. C. Nonaffine deformation of inherent structure as a static signature of cooperativity in supercooled liquids. Phys. Rev. Lett. 101, 095501. https://doi.org/10.1103/PhysRevLett.101.095501 (2008).

Article ADS Google Scholar

Mosayebi, M., Del Gado, E., Ilg, P. & Öttinger, H. C. Deformation of inherent structures to detect long-range correlations in supercooled liquids. J. Chem. Phys. 137, 024504. https://doi.org/10.1063/1.4732859 (2012).

Article ADS Google Scholar

Morineau, D., Dosseh, G., Pellenq, R. J. M., Bellissent-Funel, M. C. & Alba-Simionesco, C. Thermodynamic and structural properties of fragile glass-forming toluene and meta-xylene: Experiments and Monte-Carlo simulations. Mol. Simul. 20, 95–113. https://doi.org/10.1080/08927029708024170 (1997).

Article Google Scholar

Sieg, K., Starokozhev, E., Schmidt, M. U. & Püttmann, W. Inverse temperature dependence of Henry’s law coefficients for volatile organic compounds in supercooled water. Chemosphere 77, 8–14. https://doi.org/10.1016/j.chemosphere.2009.06.028 (2009).

Article ADS Google Scholar

Bai, Y. et al. Pressure-induced crystallization and phase transformation of para-xylene. Sci. Rep. 7, 5321. https://doi.org/10.1038/s41598-017-05639-9 (2017).

Article ADS Google Scholar

Design Institute for Physical Properties, Sponsored by AIChE, DIPPR Project 801 - Full Version, Retrieved from https://app.knovel.com/hotlink/toc/id:kpDIPPRPF7/dippr-project-801-full/dippr-project-801-full (January 30, 2020).

Stringari, P. et al. Toward an optimized design of the LNG production process: Measurement and modeling of the solubility limits of p-xylene in methane and methane + ethane mixtures at low temperature. Fluid Phase Equilib. 556, 113406. https://doi.org/10.1016/j.fluid.2022.113406 (2022).

Article Google Scholar

Siahvashi, A. et al. Solubility of p-xylene in methane and ethane and implications for freeze-out at LNG conditions. Exp. Therm. Fluid Sci. 105, 47–57. https://doi.org/10.1016/j.expthermflusci.2019.03.010 (2019).

Article Google Scholar

Trinh, E. H. & Ohsaka, K. Measurement of density, sound velocity, surface tension, and viscosity of freely suspended supercooled liquids. Int. J. Thermophys. 16, 545–555. https://doi.org/10.1007/BF01441920 (1995).

Article ADS Google Scholar

Lin, C. W. & Trusler, J. P. M. The speed of sound and derived thermodynamic properties of pure water at temperatures between (253 and 473) K and at pressures up to 400 MPa. J. Chem. Phys. 136, 094511. https://doi.org/10.1063/1.3688054 (2012).

Article ADS Google Scholar

Floriano, M. A. & Angell, C. A. Surface tension and molar surface free energy and entropy of water to − 272 °C. J. Phys. Chem. 94, 4199–4202. https://doi.org/10.1021/j100373a059 (1990).

Article Google Scholar

Vinš, V. et al. Possible anomaly in the surface tension of supercooled water: New experiments at extreme supercooling down to − 314 °C. J. Phys. Chem. Lett. 11, 4443–4447. https://doi.org/10.1021/acs.jpclett.0c01163 (2020).

Article Google Scholar

Sorensen, C. M. Densities and partial molar volumes of supercooled aqueous solutions. J. Chem. Phys. 79, 1455–1461. https://doi.org/10.1063/1.445908 (1983).

Article ADS Google Scholar

Rubio, J. E. F., Baonza, V. G., Taravillo, M., Núñez, J. & Cáceres, M. Dynamic light scattering in liquid and supercooled diphenylmethane. J. Chem. Phys. 120, 1426–1435. https://doi.org/10.1063/1.1628226 (2004).

Article ADS Google Scholar

Brodin, A. et al. Brillouin-scattering study of propylene carbonate: An evaluation of phenomenological and mode coupling analyses. Phys. Rev. E 65, 051503. https://doi.org/10.1103/PhysRevE.65.051503 (2002).

Article ADS MathSciNet Google Scholar

Cicerone, M. T. & Ediger, M. D. Enhanced translation of probe molecules in supercooled o-terphenyl: Signature of spatially heterogeneous dynamics?. J. Chem. Phys. 104, 7210–7218. https://doi.org/10.1063/1.471433 (1996).

Article ADS Google Scholar

Trejo González, J. A., Longinotti, M. P. & Corti, H. R. Diffusion-viscosity decoupling in supercooled glycerol aqueous solutions. J. Phys. Chem. B 119, 257–262. https://doi.org/10.1021/jp509055v (2015).

Article Google Scholar

Corti, H. R., Frank, G. A. & Marconi, M. C. Diffusion−viscosity decoupling in supercooled aqueous trehalose solutions. J. Phys. Chem. B 112, 12899–12906. https://doi.org/10.1021/jp802806p (2008).

Article Google Scholar

Kawasaki, T. & Kim, K. Spurious violation of the Stokes–Einstein–Debye relation in supercooled water. Sci. Rep. 9, 8118. https://doi.org/10.1038/s41598-019-44517-4 (2019).

Article ADS Google Scholar

Dueby, S., Dubey, V. & Daschakraborty, S. Decoupling of translational diffusion from the viscosity of supercooled water: Role of translational jump diffusion. J. Phys. Chem. B 123, 7178–7189. https://doi.org/10.1021/acs.jpcb.9b01719 (2019).

Article Google Scholar

Dubey, V. & Daschakraborty, S. Breakdown of the Stokes-Einstein relation in supercooled water/methanol binary mixtures: Explanation using the translational jump-diffusion approach. J. Phys. Chem. B 124, 10398–10408. https://doi.org/10.1021/acs.jpcb.0c07318 (2020).

Article Google Scholar

Mishima, O. Volume of supercooled water under pressure and the liquid-liquid critical point. J. Chem. Phys. 133, 144503. https://doi.org/10.1063/1.3487999 (2010).

Article ADS Google Scholar

Holten, V., Sengers, J. V. & Anisimov, M. A. Equation of state for supercooled water at pressures up to 400 MPa. J. Phys. Chem. Ref. Data 43, 043101. https://doi.org/10.1063/1.4895593 (2014).

Article ADS Google Scholar

Ukai, T., Kodama, D., Miyazaki, J. & Kato, M. Solubility of methane in alcohols and saturated density at 280.15 K. J. Chem. Eng. Data 47, 1320–1323. https://doi.org/10.1021/je020108p (2002).

Article Google Scholar

Dechoz, J. & Rozé, C. Surface tension measurement of fuels and alkanes at high pressure under different atmospheres. Appl. Surf. Sci. 229, 175–182. https://doi.org/10.1016/j.apsusc.2004.01.057 (2004).

Article ADS Google Scholar

Piszko, M. et al. Diffusivities in binary mixtures of n-hexane or 1-hexanol with dissolved CH4, Ne, Kr, R143a, SF6, or R236fa close to infinite dilution. J. Chem. Eng. Data 66, 2218–2232. https://doi.org/10.1021/acs.jced.1c00084 (2021).

Article Google Scholar

Cadogan, S. P., Mistry, B., Wong, Y., Maitland, G. C. & Trusler, J. P. M. Diffusion coefficients of carbon dioxide in eight hydrocarbon liquids at temperatures between (298.15 and 423.15) K at pressures up to 69 MPa. J. Chem. Eng. Data 61, 3922–3932. https://doi.org/10.1021/acs.jced.6b00691 (2016).

Article Google Scholar

Taib, M. B. M. & Trusler, J. P. M. Diffusion coefficients of methane in methylbenzene and heptane at temperatures between 323 and 398 K at pressures up to 65 MPa. Int. J. Thermophys. 41, 119. https://doi.org/10.1007/s10765-020-02700-0 (2020).

Article ADS Google Scholar

Bellaire, D., Großmann, O., Münnemann, K. & Hasse, H. Diffusion coefficients at infinite dilution of carbon dioxide and methane in water, ethanol, cyclohexane, toluene, methanol, and acetone: A PFG-NMR and MD simulation study. J. Chem. Thermody. 166, 106691. https://doi.org/10.1016/j.jct.2021.106691 (2022).

Article Google Scholar

Hayama, H., Fukuzawa, K., Yasuda, K. & Ohmura, R. Interfacial tension between (methane + ethane + propane) gas mixture and water from 283.2 to 298.2 K under up to 10 MPa. J. Chem. Thermodyn. 108, 71–75. https://doi.org/10.1016/j.jct.2017.01.007 (2017).

Article Google Scholar

Sachs, W. & Meyn, V. Pressure and temperature dependence of the surface tension in the system natural gas/water principles of investigation and the first precise experimental data for pure methane/water at 25°C up to 46.8 MPa. Colloids Surf. A Physicochem. Eng. Asp. 94, 291–301. https://doi.org/10.1016/0927-7757(94)03008-1 (1995).

Article Google Scholar

Satherley, J., Cooper, D. L. & Schiffrin, D. J. Surface tension, density and composition in the methane-pentane system at high pressure. Fluid Phase Equilib. 456, 193–202. https://doi.org/10.1016/j.fluid.2017.10.023 (2018).

Article Google Scholar

Huang, Y., Wu, C., Chen, Y., Chou, I. M. & Jiang, L. Measurement of diffusion coefficients of hydrogen sulfide in water and brine using in-situ Raman spectroscopy. Fluid Phase Equilib. 556, 113381. https://doi.org/10.1016/j.fluid.2022.113381 (2022).

Article Google Scholar

Pan, Z. & Trusler, J. P. M. Interfacial tensions of systems comprising N2, 7 mass% KI (aq), decane and iododecane at elevated pressures and temperatures. Fluid Phase Equilib. 556, 113364. https://doi.org/10.1016/j.fluid.2021.113364 (2022).

Article Google Scholar

Bellaire, D., Münnemann, K. & Hasse, H. Mutual diffusion coefficients from NMR imaging. Chem. Eng. Sci. 255, 117655. https://doi.org/10.1016/j.ces.2022.117655 (2022).

Article Google Scholar

Davies, S. R., Sloan, E. D., Sum, A. K. & Koh, C. A. In situ studies of the mass transfer mechanism across a methane hydrate film using high-resolution confocal raman spectroscopy. J. Phys. Chem. C 114, 1173–1180. https://doi.org/10.1021/jp909416y (2010).

Article Google Scholar

Moudrakovski, I. L., McLaurin, G. E., Ratcliffe, C. I. & Ripmeester, J. A. Methane and carbon dioxide hydrate formation in water droplets: Spatially resolved measurements from magnetic resonance microimaging. J. Phys. Chem. B 108, 17591–17595. https://doi.org/10.1021/jp0473220 (2004).

Article Google Scholar

Vopička, O. et al. One-pot neutron imaging of surface phenomena, swelling and diffusion during methane absorption in ethanol and n-decane under high pressure. PLoS ONE 15, e0238470. https://doi.org/10.1371/journal.pone.0238470 (2020).

Article Google Scholar

Dunning, J. R., Pegram, G. B., Fink, G. A. & Mitchell, D. P. Interaction of neutrons with matter. Phys. Rev. 48, 265–280. https://doi.org/10.1103/PhysRev.48.265 (1935).

Article ADS Google Scholar

Cibulka, I. & Takagi, T. P−ρ−T data of liquids: Summarization and evaluation. 5. Aromatic hydrocarbons. J. Chem. Eng. Data 44, 411–429. https://doi.org/10.1021/je980278v (1999).

Article Google Scholar

Matsuo, S. & Makita, T. Isotope effect on the viscosity of benzene and cyclohexane mixtures under high pressures. Int. J. Thermophys. 14, 67–77. https://doi.org/10.1007/BF00522662 (1993).

Article ADS Google Scholar

Craven, C. J., Hatton, P. D., Howard, C. J. & Pawley, G. S. The structure and dynamics of solid benzene. I. A neutron powder diffraction study of deuterated benzene from 4 K to the melting point. J. Chem. Phys. 98, 8236–8243. https://doi.org/10.1063/1.464528 (1993).

Article ADS Google Scholar

Zhao, H., Unhannanant, P., Hanshaw, W. & Chickos, J. S. Enthalpies of vaporization and vapor pressures of some deuterated hydrocarbons. Liquid−vapor pressure isotope effects. J. Chem. Eng. Data 53, 1545–1556. https://doi.org/10.1021/je800091s (2008).

Article Google Scholar

Gaines, G. L. & Le Grand, D. G. The surface tension of some deuterated hydrocarbons. Colloids Surf. A 82, 299–300. https://doi.org/10.1016/0927-7757(93)02635-R (1994).

Article Google Scholar

Dasch, C. J. One-dimensional tomography: A comparison of Abel, onion-peeling, and filtered backprojection methods. Appl. Opt. 31, 1146–1152. https://doi.org/10.1364/AO.31.001146 (1992).

Article ADS Google Scholar

Peng, D.-Y. & Robinson, D. B. A new two-constant equation of state. Ind. Eng. Chem. Fundam. 15, 59–64. https://doi.org/10.1021/i160057a011 (1976).

Article Google Scholar

Legret, D., Richon, D. & Renon, H. Vapor-liquid equilibrium of methane-benzene, methane-methylbenzene (toluene), methane-1,3-dimethylbenzene (m-xylene), and methane-1,3,5-trimethylbenzene (mesitylene) at 313.2 K up to the critical point. J. Chem. Eng. Data 27, 165–169. https://doi.org/10.1021/je00028a020 (1982).

Article Google Scholar

Ng, H. J., Huang, S. S. S. & Robinson, D. B. Equilibrium phase properties of selected m-xylene binary systems. m-Xylene-methane and m-xylene carbon dioxide. J. Chem. Eng. Data 27, 119–122. https://doi.org/10.1021/je00028a004 (1982).

Article Google Scholar

Henriksson, U. & Eriksson, J. C. Thermodynamics of capillary rise: Why is the meniscus curved?. J. Chem. Educ. 81, 150. https://doi.org/10.1021/ed081p150 (2004).

Article Google Scholar

Adamson, A. W. & Gast, A. P. Physical Chemistry of Surfaces. (Wiley, 1997).

Seber, G. A. F. & Wild, C. J. Nonlinear Regression. (John Wiley and Sons, Inc., 2003).

Kubíček, M. Numerické Algoritmy řešení Chemicko-inženýrských úloh. (SNTL, 1983).

Řehák, K., Voňka, P. & Dreiseitlová, J. Revision of the volumetric method for measurements of liquid–liquid equilibria in binary systems. Fluid Phase Equilib. 230, 109–120. https://doi.org/10.1016/j.fluid.2004.12.004 (2005).

Article Google Scholar

Watson, E. B., Wanser, K. H. & Farley, K. A. Anisotropic diffusion in a finite cylinder, with geochemical applications. Geochim. Cosmochim. Acta 74, 614–633. https://doi.org/10.1016/j.gca.2009.10.013 (2010).

Article ADS Google Scholar

Poling, B. E., Prausnitz, J. M. & O’Connell., J. P. The Properties of Gases and Liquids. (McGraw-Hill Professional, 2000).

Miyano, Y. et al. Henry’s law constants and infinite dilution activity coefficients of propane, propene, butane, 2-methylpropane, 1-butene, 2-methylpropene, trans-2-butene, cis-2-butene, 1,3-butadiene, dimethyl ether, chloroethane, and 1,1-difluoroethane in benzene, toluene, o-xylene, m-xylene, p-xylene, and styrene. J. Chem. Eng. Data 52, 291–297. https://doi.org/10.1021/je060395n (2007).

Article Google Scholar

Wilke, C. & Chang, P. Correlation of diffusion coefficients in dilute solutions. AIChE J. 1, 264–270. https://doi.org/10.1002/aic.690010222 (1955).

Article Google Scholar

Böttger, A., Pérez-Salado Kamps, Á. & Maurer, G. Solubility of methane in n-hexane and a petroleum benzine at ambient temperatures. J. Chem. Thermodyn. 99, 97–104. https://doi.org/10.1016/j.jct.2016.03.038 (2016).

Article Google Scholar

Lehmann, E. H., Vontobel, P. & Wiezel, L. Properties of the radiography facility NEUTRA at SINQ and its use as European reference facility. Nondestr. Test. Eval. 16, 191–202. https://doi.org/10.1080/10589750108953075 (2001).

Article Google Scholar

Boillat, P. et al. Chasing quantitative biases in neutron imaging with scintillator-camera detectors: A practical method with black body grids. Opt. Express 26, 15769–15784. https://doi.org/10.1364/OE.26.015769 (2018).

Article ADS Google Scholar

Carminati, C. et al. Implementation and assessment of the black body bias correction in quantitative neutron imaging. PLoS ONE 14, e0210300. https://doi.org/10.1371/journal.pone.0210300 (2019).

Article Google Scholar

Nagaoka, K. & Makita, T. Effect of pressure on the solid-liquid phase equilibria of (carbon tetrachloride + p-xylene) and (carbon tetrachloride+benzene) systems. Int. J. Thermophys. 9, 535–545. https://doi.org/10.1007/BF00503152 (1988).

Article ADS Google Scholar

Blau, B. et al. The swiss spallation neutron source SINQ at Paul Scherrer Institut. Neutron News 20, 5–8. https://doi.org/10.1080/10448630903120387 (2009).

Article Google Scholar

Download references

This work is based on experiments performed at the Swiss spallation neutron source SINQ, Paul Scherrer Institute, Villigen, Switzerland69. Experiments were conducted within beamtime proposal 20200129 at NEUTRA thermal neutron imaging beamline at the Paul Scherrer Institute. Financial support obtained from the Ministry of Education, Youth and Sports of the Czech Republic, specific university research grant A1_FCHI_2020_002, is gratefully acknowledged.

Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28, Prague 6, Czech Republic

Ondřej Vopička, Tereza-Markéta Durďáková & Petr Číhal

Laboratory for Neutron Scattering and Imaging, Paul Scherrer Institut, 5232, Villigen PSI, Switzerland

Pierre Boillat & Pavel Trtik

Electrochemistry Laboratory, Paul Scherrer Institut, 5232, Villigen PSI, Switzerland

Pierre Boillat

You can also search for this author in PubMed Google Scholar

O.V. and P.T. conceived and designed this study, realized the experiment together with T-M.D., P.Č., analyzed data together with T-M.D., P.Č., P.B., prepared tools for data acquisition and analysis together with P.B., and wrote the manuscript. All the authors discussed the results and commented on the manuscript.

Correspondence to Ondřej Vopička or Pavel Trtik.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Vopička, O., Durďáková, TM., Číhal, P. et al. Absorption of pressurized methane in normal and supercooled p-xylene revealed via high-resolution neutron imaging. Sci Rep 13, 136 (2023). https://doi.org/10.1038/s41598-022-27142-6

Download citation

Received: 07 September 2022

Accepted: 27 December 2022

Published: 04 January 2023

DOI: https://doi.org/10.1038/s41598-022-27142-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Think You Know Where Proppant Goes During Fracturing? Based on Recent Testing: Don’t Bet on It

Collision Repair Paint Problem Troubleshooting