<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1372060</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2024.1372060</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The adsorption and desorption processes of organic working fluids R1234yf, R134a, R32 in MOF-5 and Co-MOF-74: a molecular simulation study</article-title>
<alt-title alt-title-type="left-running-head">Chen</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2024.1372060">10.3389/fenrg.2024.1372060</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Chen</surname>
<given-names>Xinping</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2398610/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
</contrib-group>
<aff>
<institution>College of Artificial Intelligence and Big Data</institution>, <institution>Chongqing College of Electronic Engineering</institution>, <addr-line>Chongqing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2259061/overview">Chao Su</ext-link>, China Jiliang University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2612088/overview">Tao Wang</ext-link>, Southwest University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2637917/overview">Bo Song</ext-link>, Xi&#x2019;an Jiaotong University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2637932/overview">Jiang Pan</ext-link>, China Jiliang University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Xinping Chen, <email>202321001@cqcet.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>27</day>
<month>02</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1372060</elocation-id>
<history>
<date date-type="received">
<day>17</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>02</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Chen.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Chen</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The combination of nanoporous materials with organic working fluids holds the promise of further enhancing the performance of refrigerants based thermodynamics cycles. In this study, the adsorption and desorption properties of several organic refrigerants, e.g., R1234yf, R134a, R32, and their mixtures in metal-organic framework materials MOF-5 and Co-MOF-74 are investigated via molecular dynamics methods. The results indicate that the adsorption capacity is inversely proportional to the temperature during the adsorption process, and the adsorption capacity of the R1234yf/R32 mixture (molar ratio 2:1) is higher than that of the corresponding pure working fluids. The desorption amount, desorption regeneration rate, and desorption heat are directly proportional to the temperature. The interactions between different molecular atoms in the mixed working fluids promote the desorption process.</p>
</abstract>
<kwd-group>
<kwd>metal-organic frameworks</kwd>
<kwd>refrigerants</kwd>
<kwd>nanofluid</kwd>
<kwd>adsorption</kwd>
<kwd>molecular dynamics</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Sustainable Energy Systems</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Energy serves as the foundation for the development of modern society. And its utilization reflects the level of scientific and economic strength of a country (<xref ref-type="bibr" rid="B4">Carbon, 2021</xref>). Therefore, the energy crisis and environmental pollution have drawn a lot of attentions. Thus, enhancing the efficiency of fossil fuel utilization and enlarging the scale of renewable energy are proposed to maintain the sustainable development of the world (<xref ref-type="bibr" rid="B6">Chen et al., 2019</xref>; <xref ref-type="bibr" rid="B35">Xiao et al., 2022</xref>).</p>
<p>So far, the thermodynamic cycles are the main approach for energy conversion. Organic Rankine Cycle (ORC) uses organic refrigerants as working fluids to recover the low-grade energies, e.g., industrial waste heat, geothermal energy and solar energy, which have a very good application prospects (<xref ref-type="bibr" rid="B29">Vivian et al., 2015</xref>; <xref ref-type="bibr" rid="B37">Zhang et al., 2019</xref>; <xref ref-type="bibr" rid="B23">Miao et al., 2020</xref>; <xref ref-type="bibr" rid="B30">Wang et al., 2024</xref>). However, the efficiency of ORC is relatively low because of the low temperature of heating source. Since the working fluid plays as the energy carrier in the thermodynamic cycles, the properties optimization of organic refrigerants could result in the efficiency enhancement of ORC (<xref ref-type="bibr" rid="B28">Su et al., 2017</xref>; <xref ref-type="bibr" rid="B32">Wang et al., 2020</xref>; <xref ref-type="bibr" rid="B18">Li et al., 2022</xref>). One practical method is using the zeotropic mixture as working fluids, which have extensively reported (<xref ref-type="bibr" rid="B33">Wang et al., 2019</xref>; <xref ref-type="bibr" rid="B22">McLinden and Huber, 2020</xref>; <xref ref-type="bibr" rid="B41">Zhu et al., 2021</xref>). Besides, <xref ref-type="bibr" rid="B21">McGrail et al. (2013)</xref> proposed to use metal-organic frameworks (MOFs) to enhance the properties of refrigerants. This is the so called metal organic heat carriers (MOHCs). The adsorption energy is stored in MOFs at the condensation stage and then into thermal energy by desorbing the fluid molecules from MOFs during heating process.</p>
<p>MOFs are porous materials with periodic network structures formed by the self-assembly of metal ions or clusters and organic ligands through coordination (<xref ref-type="bibr" rid="B15">Kong et al., 2013</xref>). MOFs exhibit high specific surface area, large pore volume, and tunable structures, making them highly promising for applications in gas storage, heterogeneous catalysis, adsorption-based energy storage and other fields (<xref ref-type="bibr" rid="B9">Furukawa et al., 2013</xref>).</p>
<p>
<xref ref-type="bibr" rid="B26">Qasem et al. (2018)</xref> investigated the adsorption of CO<sub>2</sub> by MOF-5 and MOF-177 under pressures ranging from 5 to 50&#xa0;bar. The results showed that MOFs could store more CO<sub>2</sub> while consuming less energy compared to traditional storage systems. <xref ref-type="bibr" rid="B31">Wang et al. (2018)</xref> studied the adsorption and energy storage of R161 and similar refrigerants in MOF-5. It was found that fluorine atoms in the organic working pairs enhanced their adsorption in MOFs. <xref ref-type="bibr" rid="B13">Hu et al. (2018a)</xref>, <xref ref-type="bibr" rid="B14">Hu et al. (2018b)</xref>, <xref ref-type="bibr" rid="B12">Hu et al. (2019)</xref> investigated thermal energy storage properties of several refrigerants in MOF-74 and MOF-5. The results showed that the energy storage enhancement ratios of R1234yf, R1234ze and R134a with Mg-MOF-74 nanoparticles are higher than those of other M-MOF-74 (M &#x3d; Co, Ni, Zn) materials.</p>
<p>In recent years, with the rapid development of computer technology, molecular simulation techniques have demonstrated their unique advantages in screening high-quality adsorbents for energy storage working pairs (<xref ref-type="bibr" rid="B38">Zhang et al., 2022</xref>; <xref ref-type="bibr" rid="B16">Lei et al., 2023</xref>). Currently, research on organic working pairs in nano-porous materials mostly focuses on the discussion of adsorption mechanisms and analysis of storage properties (<xref ref-type="bibr" rid="B40">Zhou et al., 2019</xref>; <xref ref-type="bibr" rid="B17">Li et al., 2020</xref>; <xref ref-type="bibr" rid="B39">Zheng et al., 2020</xref>; <xref ref-type="bibr" rid="B10">Garc&#xed;a et al., 2021</xref>). The reports on desorption processes and regeneration properties of MOHCs are few. Besides, for different heat source environments, there are multiple choices of working fluids for ORC. With the growing awareness of environmental protection, the fourth-generation refrigerants with low global warming potential (GWP) and zero ozone depletion potential (ODP) will gradually replace the existing third-generation refrigerants. Among them, R1234yf, R134a, R32 are currently the hot topics of research (<xref ref-type="bibr" rid="B8">Fouad and Vega, 2018</xref>; <xref ref-type="bibr" rid="B20">Liu et al., 2019</xref>). MOF-5 (<xref ref-type="bibr" rid="B2">Annapureddy et al., 2014</xref>), as a fundamental and classic MOF structure, exhibits good thermal stability and chemical stability at around 300&#xb0;C, making it suitable for the application environment of low-grade energy in ORC systems. MOF-74 (<xref ref-type="bibr" rid="B5">Chaemchuen et al., 2018</xref>), being one of the MOFs with the highest unsaturated metal site density, can enhance its adsorption performance by adjusting the types of metal ions. Among the MOF-74 series, Co-MOF-74 stands out as the material with the maximum adsorption heat, gaining widespread attention.</p>
<p>Consequently, in this work, the metal-organic frameworks materials of MOF-5 and Co-MOF-74 are employed along with organic working fluids R32, R1234yf, R134a and their mixtures to form corresponding adsorption working pairs. Molecular dynamics simulations are conducted to investigate the desorption and regeneration behaviors of the organic working fluids within the MOFs. This work also aims to provide useful insight in adsorption refrigeration.</p>
</sec>
<sec sec-type="methods" id="s2">
<title>2 Methodology</title>
<sec id="s2-1">
<title>2.1 Simulation models</title>
<p>This study simulated the adsorption and desorption regeneration processes of pure refrigerants (R1234yf, R134a, R32) and their mixtures (<xref ref-type="table" rid="T1">Table 1</xref>) within MOF-5 and Co-MOF-74.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>The molar ratio of mixed refrigerants.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Mixed refrigerants</th>
<th align="center">molar ratio</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">R1234yf &#x2b; R32</td>
<td align="center">1:3</td>
</tr>
<tr>
<td align="center">R1234yf &#x2b; R32</td>
<td align="center">1:1</td>
</tr>
<tr>
<td align="center">R1234yf &#x2b; R32</td>
<td align="center">3:1</td>
</tr>
<tr>
<td align="center">R134a &#x2b; R32</td>
<td align="center">1:3</td>
</tr>
<tr>
<td align="center">R134a &#x2b; R32</td>
<td align="center">1:1</td>
</tr>
<tr>
<td align="center">R134a &#x2b; R32</td>
<td align="center">3:1</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The spatial structures of MOF-5 (Fig. a) and Co-MOF-74 (Fig. b) were shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. The MOF-5 adopts a 2 &#xd7; 2 &#xd7; 2 supercell structure (X: 5.1664&#xa0;nm, Y: 5.1664&#xa0;nm, Z: 5.1664&#xa0;nm). It consists of 1536&#xa0;C atoms, 768&#xa0;H atoms, 832 O atoms, and 256 Zn atoms. Besides, the Co-MOF-74 utilizes a 1 &#xd7; 1 &#xd7; 4 supercell structure (X: 2.6136&#xa0;nm, Y: 2.6136&#xa0;nm, Z: 2.7768&#xa0;nm). It contains 288&#xa0;C atoms, 72&#xa0;H atoms, 216 O atoms, and 72 Co atoms. The atomic partial charges for the MOF materials were calculated by the Extended Charge Equilibration (EQeq) method (<xref ref-type="bibr" rid="B34">Wilmer et al., 2012</xref>). The interatomic interactions were described by the 12&#x2013;6 Lennard-Jones (LJ) potential (<xref ref-type="bibr" rid="B36">Xu et al., 2019</xref>). Specific force field parameters were listed in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Structures of MOFs: <bold>(A)</bold> MOF-5 2 &#xd7; 2&#xd7;2 supercell; <bold>(B)</bold> Co-MOF-74 1 &#xd7; 1&#xd7;4 supercell.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g001.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>The force field parameters used for MOF-5 and Co-MOF-74.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Parameter</th>
<th align="center">MOF-Zn</th>
<th align="center">MOF-C</th>
<th align="center">MOF-O</th>
<th align="center">MOF-H</th>
<th align="center">MOF-Co</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<italic>&#x3c3;</italic> (&#xc5;)</td>
<td align="center">2.462</td>
<td align="center">3.431</td>
<td align="center">3.118</td>
<td align="center">2.571</td>
<td align="center">2.560</td>
</tr>
<tr>
<td align="center">
<italic>&#x3b5;</italic>/<italic>k</italic>
<sub>B</sub> (K)</td>
<td align="center">62.40</td>
<td align="center">52.84</td>
<td align="center">30.19</td>
<td align="center">22.14</td>
<td align="center">7.045</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The molecular structures of the three selected organic working fluids in this study are illustrated in <xref ref-type="fig" rid="F2">Figure 2</xref>. A full atomistic force field model, proposed by Raabe and Maginn (<xref ref-type="bibr" rid="B27">Raabe and Maginn, 2010</xref>), was employed to investigate these organic working fluids. The potential energy of interparticle interactions can be described by Eq. <xref ref-type="disp-formula" rid="e1">1</xref>:<disp-formula id="equ1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>U</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>&#x3c6;</mml:mi>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c6;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>&#x3c7;</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>&#x3c7;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="e1">
<mml:math id="m2">
<mml:mrow>
<mml:mfenced open="" close="}" separators="|">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mi>i</mml:mi>
</mml:munder>
</mml:mstyle>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>12</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>6</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Spatial structures of R1234yf, R134a and R32.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g002.tif"/>
</fig>
<p>Furthermore, the cutoff distance for intermolecular interactions was set to 12&#xa0;&#xc5;. The non-bonded interaction parameters between different types of atoms were calculated by the Lorentz-Berthelot combination rules (<xref ref-type="bibr" rid="B7">Delhommelle and Milli&#xc9;, 2001</xref>). The particle-particle/particle-mesh (PPPM) method (<xref ref-type="bibr" rid="B11">Hockney et al., 1974</xref>) was employed to calculate the long-range Coulomb interactions.</p>
<p>The initial model for adsorption simulation is shown in <xref ref-type="fig" rid="F3">Figure 3A</xref> (R1234yf/MOF-5). The model consisted of a 2 &#xd7; 2 &#xd7; 2 supercell of MOF-5 and two liquid membrane of R1234yf. The simulation box size was X: 53.355&#xa0;&#xc5;, Y: 53.355&#xa0;&#xc5;, Z: 181.134&#xa0;&#xc5;. Each liquid membrane consisted of 648 R1234yf molecules distributed on both sides of the simulation box with dimensions of X: 53.355&#xa0;&#xc5;, Y: 53.355&#xa0;&#xc5; and Z: 34.80&#xa0;&#xc5;. The MOF-5 supercell was positioned in the center of the system along the Z direction and accounts for 25% mass fraction in the entire system. Periodic boundary conditions were applied in all three directions for the simulation box. The initial model for desorption regeneration was obtained by removing the unabsorbed working fluids outside the MOF from the model at the end of adsorption at 300&#xa0;K, as shown in <xref ref-type="fig" rid="F3">Figure 3B</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>The initial model of adsorption and desorption of the R1234yf/MOF-5 working pair system. <bold>(A)</bold> adsorption; <bold>(B)</bold> desorption.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g003.tif"/>
</fig>
</sec>
<sec id="s2-2">
<title>2.2 Simulation details</title>
<p>In this study, the molecular dynamics simulation processes were conducted by the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) (<xref ref-type="bibr" rid="B25">Plimpton, 1995</xref>). The visualization of the models and research processes were performed with the OVITO software (<xref ref-type="bibr" rid="B1">Alexander, 2010</xref>). The thermophysical data for the organic working fluids were obtained from the National Institute of Standards and Technology (NIST) (<xref ref-type="bibr" rid="B42">webbook</xref>). The entire simulation processes were carried out under the NVT canonical ensemble. And the adsorption temperatures were set at 300&#xa0;K, 330&#xa0;K, 360&#xa0;K, and 390&#xa0;K, and the desorption temperatures were set at 330&#xa0;K, 360&#xa0;K, and 390&#xa0;K. The Berendsen thermostat method (<xref ref-type="bibr" rid="B3">Berendsen et al., 1984</xref>) was utilized to control the system temperature. The time step was set at 1 fs with adsorption and desorption simulations performed at each temperature for 5&#xa0;ns. Atom information was output every 10,000 steps.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Analysis of the adsorption process</title>
<p>To validate the accuracy of the computational model utilized in this study, the adsorption capacity of R134a in MOF-5 was simulated. The adsorption capacity of R134a in MOF-5 at 330&#xa0;K of present model is about 1.442&#xa0;g/g while that of previous work (<xref ref-type="bibr" rid="B19">Liu et al., 2022</xref>) is about 1.405&#xa0;g/g, yielding a relative error of 2.5%. Thus, the utilization of the computational model in investigating the adsorption and desorption processes of organic working fluids in metal-organic framework materials is deemed feasible in this study.</p>
<sec id="s3-1">
<title>3.1 Adsorption capacity</title>
<p>The molecules of working fluid located in the region of MOFs are regarded as the adsorbed particles, which is used for calculating the adsorption capacities, m<sub>adsorption</sub>, and its specific calculation formula (<xref ref-type="disp-formula" rid="e2">2</xref>) is as follows:<disp-formula id="e2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2219;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mi mathvariant="bold-italic">O</mml:mi>
<mml:mi mathvariant="bold-italic">F</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>The adsorption capacities of pure working fluids, including R1234yf, R134a, R32, as well as their mixed counterparts in MOF-5, are statistically analyzed and are presented in <xref ref-type="fig" rid="F4">Figure 4</xref>. The results show that a decrease in adsorption capacity with increasing temperature for both pure and mixed working fluids, which agree with previous related works (<xref ref-type="bibr" rid="B13">Hu et al., 2018a</xref>; <xref ref-type="bibr" rid="B14">Hu et al., 2018b</xref>; <xref ref-type="bibr" rid="B12">Hu et al., 2019</xref>; <xref ref-type="bibr" rid="B19">Liu et al., 2022</xref>). This can be attributed to the increase in thermodynamic energy of the working fluid and the enhancement of molecular thermal motion caused by temperature increase, leading to facilitating desorption of organic working fluids from the surface of porous materials. In terms of adsorption capacity for pure working fluids, the ranking order is R134a &#x3e; R1234yf &#x3e; R32.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Adsorption of pure working fluids and their mixtures in MOF-5. <bold>(A)</bold> pure working fluids, <bold>(B)</bold> R1234yf/R32, <bold>(C)</bold> R134a/R32.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g004.tif"/>
</fig>
<p>And the analysis of adsorption data for mixed refrigerants reveals that in the adsorption of mixed refrigerants, the adsorption capacities of each component are lower than the pure refrigerants. This is attributed to the shared occupancy of the MOF-5 internal space by the mixed refrigerants. Among the three organic refrigerants, R1234yf and R134a have similar sizes, both larger than R32. The similar phenomenon is also can be found in Li et al. work (<xref ref-type="bibr" rid="B17">Li et al., 2020</xref>). Due to the smaller pore size of MOF-5, R1234yf and R134a occupy the majority of the adsorption sites, which results in a significant decrease in the adsorption capacity of R32. As the quantity of R1234yf and R134a increases, the decrease in the adsorption capacity of R32 becomes more pronounced. Additionally, the adsorption capacity of R1234yf in R1234yf/R32 is higher than that of R134a in R134a/R32, which is contrary to the results obtained for pure refrigerants. This is due to the stronger adsorption selectivity of R1234yf towards R32 compared to R134a during the adsorption of mixed refrigerants. The observation is further supported by the higher adsorption capacity of R32 in R134a/R32 compared to R1234yf/R32. When the molar ratio of R1234yf to R32 is 2:1 in R1234yf/R32, the adsorption capacity exceeds that of the individual pure refrigerant components. This is attributed to the utilization of the MOF-5 pore space through the combination of the large molecular structure of R1234yf and the small molecular structure of R32 at the current molar ratio, which is resulting in superior performance compared to pure refrigerants.</p>
<p>The adsorption amounts of organic refrigerants in Co-MOF-74 are shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. Similar to MOF-5, the adsorption amounts of organic refrigerants decrease with the increase of temperature. The relationship among the adsorption amounts of pure refrigerants is as follows: R134a &#x3e; R1234yf &#x3e; R32.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Adsorption of pure working fluids and their mixtures in Co-MOF-74. <bold>(A)</bold> pure working fluids, <bold>(B)</bold> R1234yf/R32, <bold>(C)</bold> R134a/R32.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g005.tif"/>
</fig>
<p>The analysis of the adsorption data for mixed refrigerants reveals that in the adsorption of mixed refrigerants, the adsorption capacities of each component are lower than those of pure refrigerants. This can be attributed to the shared occupancy of the Co-MOF-74 internal space by the mixed refrigerants. In the case of R1234yf/R32 and R134a/R32 mixed refrigerants with a molar ratio of 2:1, both R1234yf and R134a exhibit an increase in adsorption amount with increasing temperature. This is due to the enhanced molecular thermal motion at higher temperatures, which makes it easier for the smaller-sized R32 molecules to desorb from the porous material, which creates vacant adsorption sites for a small amount of R1234yf and R134a molecules to enter.</p>
<p>In the R134a/R32 mixed refrigerant system, the adsorption amount of R134a is higher than that of R1234yf in the R1234yf/R32 mixed refrigerant. This indicates that during the adsorption of mixed refrigerants in Co-MOF-74, R134a exhibits stronger adsorption selectivity towards R32 compared to R1234yf. When the molar ratio of R1234yf to R32 is 2:1 in the R1234yf/R32 mixed refrigerant, the adsorption amount exceeds that of the individual pure refrigerant components. This is attributed to the effective utilization of the Co-MOF-74 pore space through the combination of the large molecular structure of R1234yf and the small molecular structure of R32 at the current molar ratio, which results in superior performance compared to pure refrigerants.</p>
<p>Next, the adsorption capacity of organic working fluids in two different metal-organic framework (MOF) materials, MOF-5 and Co-MOF-74, was analyzed. It was observed that, for different organic working fluids, the adsorption capacity followed the trend of MOF-5&#x3e;Co-MOF-74. This can be attributed to the distinct structural characteristics of the two MOF materials. MOF-5, with its three-dimensional network structure, exhibits a larger specific surface area and higher porosity compared to Co-MOF-74, which possesses a three-dimensional structure with one-dimensional pore channels. A larger specific surface area and higher porosity result in more available adsorption sites per unit mass of the MOFs, leading to a higher adsorption capacity of organic working fluids.</p>
<p>Furthermore, in the adsorption of mixed working fluids in MOF-5, it was observed that the adsorption capacity of R1234yf in the R1234yf/R32 mixture was higher than that of R134a in the R134a/R32 mixture. However, in Co-MOF-74, the opposite trend was observed. This discrepancy can be attributed to the smaller pore size of MOF-5, where R1234yf, being the largest molecule among the three working fluids, occupies a majority of the available space, resulting in a higher adsorption capacity. On the other hand, Co-MOF-74 has a larger pore size, and the size of the organic working fluid molecules has a lesser impact on the adsorption capacity. As a result, the pure working fluid with the highest adsorption capacity, R134a, exhibits a higher adsorption capacity in Co-MOF-74. Additionally, the smaller pore size of MOF-5 also prevents the phenomenon of increasing adsorption capacity of mixed working fluid components with rising temperature, which is observed in Co-MOF-74.</p>
</sec>
<sec id="s3-2">
<title>3.2 Adsorption heat</title>
<p>Adsorption heat is the thermal effect generated during the adsorption process of a substance and serves as an important indicator of adsorption characteristics. During the adsorption process, gas or liquid molecules move towards the surface of the porous medium, leading to a significant decrease in their thermodynamic energy. This energy is converted into surface energy, releasing heat, which is known as adsorption heat. The adsorption heat can be calculated using the following equation (<xref ref-type="disp-formula" rid="e3">3</xref>) (<xref ref-type="bibr" rid="B24">Niu et al., 2019</xref>):<disp-formula id="e3">
<mml:math id="m4">
<mml:mrow>
<mml:mi mathvariant="bold">&#x394;</mml:mi>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the energy of the stable configuration formed between the organic adsorbate and the metal-organic framework material. <inline-formula id="inf2">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the energy of the metal-organic framework material before adsorption, and <inline-formula id="inf3">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the energy of the organic adsorbate before adsorption. <inline-formula id="inf4">
<mml:math id="m8">
<mml:mrow>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the adsorption heat. Since the adsorption process typically releases heat, a larger <inline-formula id="inf5">
<mml:math id="m9">
<mml:mrow>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> indicates a greater amount of heat released during the adsorption process, resulting in a more stable configuration after adsorption.</p>
<p>The adsorption heat of organic working fluids in MOF-5 is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. At about 330&#xa0;K, the adsorption heat suddenly increases, which is attributed to a phase transition of the organic working fluid at this temperature, leading to an increase in <inline-formula id="inf6">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Apart from the phase transition point, as the temperature rises, the adsorption of organic working fluids per unit mass of the MOF decreases, and the adsorption heat also decreases accordingly. The magnitude of adsorption heat for pure working fluids in MOF-5 follows the order: R134a &#x3e; R32 &#x3e; R1234yf. The higher the proportion of fluorine atoms, the greater the adsorption heat, possibly due to the stronger interaction forces between the fluorine atoms and other atoms during adsorption.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Adsorption heat of pure and mixed working fluids in MOF-5. <bold>(A)</bold> pure working fluids, <bold>(B)</bold> mixed working fluids.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g006.tif"/>
</fig>
<p>The adsorption heat curve of mixed working fluids is more complex compared to pure working fluids, which is due to the presence of multiple components in the mixture. During adsorption, although the adsorption amount decreases with the increase of temperature, the changes in composition vary due to different adsorption selectivity. Overall, the adsorption heat of the R134a/R32 mixed working fluid is higher than that of the R1234yf/R32 mixed working fluid. This is partly because the adsorption heat of pure R134a is greater than that of R1234yf. Additionally, R1234yf exhibits stronger adsorption selectivity in MOF-5, resulting in a lower adsorption amount of R32 in the R1234yf/R32 mixed working fluid compared to the adsorption amount of R32 in the R134a/R32 mixed working fluid.</p>
<p>The adsorption heats of organic adsorbates in Co-MOF-74 are shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. Similar to MOF-5, there is a sharp transition point due to phase change in the adsorption heat at 330&#xa0;K, and the overall adsorption heat decreases with increasing temperature. The order of adsorption heats for pure adsorbates in Co-MOF-74 is R32 &#x3e; R134a &#x3e; R1234yf, which is different from that in MOF-5 where the order is R134a &#x3e; R32 &#x3e; R1234yf. This is partly because the mass fractions of fluorine in R32 (73.07%) and R134a (74.50%) are similar, and partly because the difference in adsorption between R32 and R134a in Co-MOF-74 is smaller than that in MOF-5.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Adsorption heat of pure and mixed working fluids in Co-MOF-74. <bold>(A)</bold> pure working fluids. <bold>(B)</bold> mixed working fluids.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g007.tif"/>
</fig>
<p>The adsorption heat of R1234yf/R32 mixed refrigerant is higher than that of R134a/R32 mixed refrigerant, because the adsorption of R32 in R1234yf/R32 mixed refrigerant is greater than that in R134a/R32 mixed refrigerant. In addition, in both types of mixed refrigerants, the order of adsorption heats is 1:2&#xa0;M ratio &#x3e;1:1&#xa0;M ratio &#x3e;2:1&#xa0;M ratio. That is, the adsorption heat of mixed refrigerants is proportional to the mass fraction of R32.</p>
<p>Then, the comparative analysis of the adsorption heats of organic adsorbates in two types of MOF materials reveals that for different organic adsorbates, the order of adsorption heats is consistently MOF-5 &#x3e; Co-MOF-74. Therefore, it can be concluded that for the same adsorbate at the same temperature, its adsorption heat in a unit mass of MOF is directly proportional to the specific surface area and porosity of the MOF material. The adsorption heat increases when the organic adsorbate undergoes phase transition, and generally decreases with increasing temperature. On the other hand, for mixed refrigerants, although the adsorption capacity decreases with increasing temperature, the variation in composition differs due to different adsorption selectivity. This complexity results in a more intricate adsorption heat curve. Overall, for mixed refrigerants with the same composition, the adsorption heat increases with an increase in the mass fraction of the component with higher adsorption heat.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Analysis of the desorption process</title>
<sec id="s4-1">
<title>4.1 Desorption regeneration rate</title>
<p>After completing the adsorption process in MOFs, the accumulated organic refrigerants within them need to undergo a desorption regeneration process to achieve the recycling of the adsorbent. Therefore, the regenerability of the refrigerant in MOFs is a key indicator during the desorption process. Due to variations in the initial adsorption amounts among different adsorption models, it is not appropriate to simply use the desorption amount as a measure of the regenerability of organic refrigerants in MOFs. In this study, the desorption regeneration rate is employed to evaluate the regenerability of organic refrigerants. The desorption regeneration rate is defined as the ratio of the desorption amount to the original adsorption amount, which can be calculated using Eq. (<xref ref-type="disp-formula" rid="e4">4</xref>):<disp-formula id="e4">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">R</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn mathvariant="bold">1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">T</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mn mathvariant="bold">100</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <inline-formula id="inf7">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf8">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the moles of working fluids in the MOF at adsorption and desorption temperature, respectively.</p>
<p>The desorption regeneration rates of organic refrigerants and their mixtures in MOF-5 are shown in <xref ref-type="fig" rid="F8">Figure 8</xref>. As the temperature increases, the desorption regeneration rate also increases. The reason for this phenomenon is that with the temperature rising, the thermodynamic energy of the organic refrigerant molecules increases, which makes them more likely to desorb from the MOF. The desorption regeneration rates of pure refrigerants are in the following order: R32 &#x3e; R1234yf &#x3e; R134a. This is due to the smaller size of R32 molecules which makes them easier to desorb from the MOF. R1234yf and R134a have similar molecular sizes, but R134a exhibits stronger interactions with MOF-5, which results in a slightly lower desorption regeneration rate compared to R1234yf.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Desorption regeneration rates of pure and mixed working fluids in MOF-5. <bold>(A)</bold> desorption and regeneration rates of each organic working fluid, <bold>(B)</bold> R1234yf/R32, <bold>(C)</bold> R134a/R32.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g008.tif"/>
</fig>
<p>Among mixtures of the same composition but different molar ratios, the desorption regeneration rate is directly proportional to the mass fraction of R32. The order of desorption regeneration rates is as follows: molar ratio (1:2) &#x3e; molar ratio (1:1) &#x3e; molar ratio (2:1). For mixtures with the same molar ratio but different components, the order of desorption regeneration rates is R134a/R32 &#x3e; R1234yf/R32. This is because the adsorption amount of R32 varies in different mixtures. During adsorption, R1234yf exhibits stronger selectivity for R32 compared to R134a. The adsorption amount of R32 in R134a/R32 is greater than that in R1234yf/R32, which results in a higher desorption regeneration rate for R134a/R32 compared to R1234yf/R32.</p>
<p>Analyzing the desorption regeneration rates of each component in the mixture, it is observed that the desorption regeneration rate of R32 is not only greater than the other components but also higher than the desorption regeneration rate of the pure refrigerant. This finding indicates that the intermolecular interactions between different molecular species in the mixture play a promoting role in the desorption process.</p>
<p>The desorption regeneration rates of organic refrigerants in Co-MOF-74 are shown in <xref ref-type="fig" rid="F9">Figure 9</xref>. As the temperature increases, the desorption regeneration rate also increases. The order of desorption regeneration rates of pure refrigerants is the same as that in MOF-5, namely, R32 &#x3e; R1234yf &#x3e; R134a. For mixtures with the same composition but different molar ratios, the desorption regeneration rate is directly proportional to the mass fraction of R32. The order of desorption regeneration rates is molar ratio (1:2) &#x3e; molar ratio (1:1) &#x3e; molar ratio (2:1). For mixtures with the same molar ratio but different components, the order of desorption regeneration rates is R134a/R32 &#x3e; R1234yf/R32.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Desorption regeneration rates of pure and mixed working fluids in Co-MOF-74. <bold>(A)</bold> desorption and regeneration rates of each organic working fluid, <bold>(B)</bold> R1234yf/R32, <bold>(C)</bold> R134a/R32.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g009.tif"/>
</fig>
<p>Besides, it is found that R32 which has a smaller molecular structure exhibits a higher desorption regeneration rate than the other components in the mixture. Due to the intermolecular interactions between different molecular species in the mixture, the desorption regeneration rates of each component in the mixture are higher than those of their corresponding pure refrigerants. And for the same organic working fluid, the desorption regeneration rate follows the order Co-MOF-74 &#x3e; MOF-5. This is attributed to the different structures of the MOFs. Co-MOF-74, as a one-dimensional porous metal-organic framework material with a three-dimensional structure, has larger pore size and simpler internal structure compared to MOF-5, making it easier for organic working fluids to desorb.</p>
</sec>
<sec id="s4-2">
<title>4.2 Desorption heat</title>
<p>Desorption heat refers to the thermal effect generated during the desorption process of a working fluid. Here, desorption heat is defined as the energy released per unit mass of MOF during the desorption process of the organic working fluid, starting from the system state when the adsorption of the working fluid ends at 300&#xa0;K. This can be obtained by calculating the energy difference before and after the desorption process using a simulation model, the specific calculation method is shown in Eq. (<xref ref-type="disp-formula" rid="e5">5</xref>):<disp-formula id="e5">
<mml:math id="m14">
<mml:mrow>
<mml:mo>&#x2206;</mml:mo>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">d</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">o</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">b</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where <inline-formula id="inf9">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the energy of the stable configuration formed by the adsorption working pairs. <inline-formula id="inf10">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the energy after desorption of MOF, <inline-formula id="inf11">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the energy of working fluids after desorption.</p>
<p>The desorption heats of pure adsorbates and their mixed refrigerants in MOF-5 are shown in <xref ref-type="fig" rid="F10">Figure 10</xref>. It can be observed that as the temperature increases, both the desorption heats of pure adsorbates and mixed refrigerants increase. This is because at higher temperatures, more organic fluid adsorbates are released from the MOF, resulting in the release of more heat. For pure adsorbates, the order of desorption heats is R32 &#x3e; R134a &#x3e; R1234yf, which is inversely proportional to the molecular volume of the organic adsorbates.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Desorption heat of pure and mixed working fluids in MOF-5. <bold>(A)</bold> pure working fluids, <bold>(B)</bold> mixed working fluids.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g010.tif"/>
</fig>
<p>For mixed refrigerants with the same composition but different molar ratios, the relationship between desorption heats follows the order: molar ratio (1:2) &#x3e; molar ratio (1:1) &#x3e; molar ratio (2:1). In other words, for mixed refrigerants with the same composition, the desorption heat increases as the mass fraction of R32 in the mixed refrigerant increases. As for mixed refrigerants with the same molar ratio but different components, the relationship between desorption heats follows the order: R134a &#x3e; R1234yf, which similarly aligns with the conclusion that desorption heat is inversely proportional to the molecular volume of the organic adsorbates.</p>
<p>The desorption heats of pure refrigerants and their mixed refrigerants in Co-MOF-74 are shown in <xref ref-type="fig" rid="F11">Figure 11</xref>. As the temperature increases, the desorption heat also increases. In the desorption process of pure refrigerants, R32 exhibits the highest desorption heat, reaching 262.27&#xa0;kJ/kg at 390&#xa0;K, while R1234yf exhibits the lowest desorption heat, with a value of 157.31&#xa0;kJ/kg at 390&#xa0;K.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Desorption heat of pure and mixed working fluids in Co-MOF-74. <bold>(A)</bold> pure working fluids, <bold>(B)</bold> mixed working fluids.</p>
</caption>
<graphic xlink:href="fenrg-12-1372060-g011.tif"/>
</fig>
<p>For mixed refrigerants with the same composition but different molar ratios, the desorption heat increases as the mass fraction of R32 in the mixed refrigerant increases, following the order: molar ratio (1:2) &#x3e; molar ratio (1:1) &#x3e; molar ratio (2:1). As for mixed refrigerants with the same molar ratio but different components, the relationship between desorption heats follows the order: R1234yf &#x3e; R134a.</p>
<p>Besides, for the same working fluid, the desorption heat in MOF-5 is greater than that in Co-MOF-74. Hence, desorption heat exhibits a similar relationship to adsorption amount and heat, increasing with higher specific surface area and porosity. The order of desorption heats for pure working fluids in both MOF materials is as follows: R32 &#x3e; R134a &#x3e; R1234yf. For mixed working fluids with the same composition, the desorption heat increases with the mass fraction of R32 in the mixture. In Co-MOF-74, the desorption heat varies significantly with changes in molar ratio for mixed working fluids.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>This work investigates the adsorption and desorption regeneration processes of pure refrigerants R1234yf, R134a, R32, and their mixtures in MOF-5 and Co-MOF-74 through simulation. The following conclusions can be drawn,<list list-type="simple">
<list-item>
<p>(1) The adsorption capacity of organic refrigerants is inversely proportional to temperature. The mixed refrigerant R1234yf/R32 (with a molar ratio of 2:1) fully utilizes the space within the MOF, resulting in a higher adsorption capacity compared to the individual pure components.</p>
</list-item>
<list-item>
<p>(2) When organic refrigerants undergo phase transitions, their adsorption heat suddenly increases due to the increased thermodynamic potential of the refrigerant itself. In the absence of phase transitions, the adsorption heat decreases with increasing temperature.</p>
</list-item>
<list-item>
<p>(3) With the increase in temperature, the thermal motion capability of molecules enhances, resulting in higher desorption capacity, desorption regeneration rate, and desorption heat of organic refrigerants in the MOF material.</p>
</list-item>
<list-item>
<p>(4) Co-MOF-74 with larger pore size exhibits a higher desorption regeneration rate than MOF-5, while MOF-5 has a higher desorption heat compared to Co-MOF-74.</p>
</list-item>
<list-item>
<p>(5) R32 is more easily desorbed from the MOF material with smaller molecular structure. The interactions between atoms of different molecules in the mixed refrigerant play a promoting role in the desorption process.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>XC: Writing&#x2013;original draft.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alexander</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Visualization and analysis of atomistic simulation data with OVITO&#x2013;the Open Visualization Tool</article-title>. <source>Model. Simul. Mater. Sci. Eng.</source> <volume>18</volume> (<issue>1</issue>), <fpage>015012</fpage>. <pub-id pub-id-type="doi">10.1088/0965-0393/18/1/015012</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Annapureddy</surname>
<given-names>H. V. R.</given-names>
</name>
<name>
<surname>Motkuri</surname>
<given-names>R. K.</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>P. T.</given-names>
</name>
<name>
<surname>Truong</surname>
<given-names>T. B.</given-names>
</name>
<name>
<surname>Thallapally</surname>
<given-names>P. K.</given-names>
</name>
<name>
<surname>McGrail</surname>
<given-names>B. P.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>Computational studies of adsorption in metal organic frameworks and interaction of nanoparticles in condensed phases</article-title>. <source>Mol. Simul.</source> <volume>40</volume> (<issue>7-9</issue>), <fpage>571</fpage>&#x2013;<lpage>584</lpage>. <pub-id pub-id-type="doi">10.1080/08927022.2013.829224</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Berendsen</surname>
<given-names>H. J. C.</given-names>
</name>
<name>
<surname>Postma</surname>
<given-names>J. P. M.</given-names>
</name>
<name>
<surname>van Gunsteren</surname>
<given-names>W. F.</given-names>
</name>
<name>
<surname>DiNola</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Haak</surname>
<given-names>J. R.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Molecular dynamics with coupling to an external bath</article-title>. <source>J. Chem. Phys.</source> <volume>81</volume> (<issue>8</issue>), <fpage>3684</fpage>&#x2013;<lpage>3690</lpage>. <pub-id pub-id-type="doi">10.1063/1.448118</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<collab>Carbon</collab> (<year>2021</year>). <article-title>Project team on the strategy and pathway for peaked carbon emissions and carbon neutrality. Analysis of a peaked carbon emission pathway in China toward carbon neutrality</article-title>. <source>Engineering</source> <volume>7</volume> (<issue>12</issue>), <fpage>1673</fpage>&#x2013;<lpage>1677</lpage>. <pub-id pub-id-type="doi">10.1016/j.eng.2021.10.003</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chaemchuen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Klomkliang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Yusubov</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Verpoort</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Tunable metal&#x2013;organic frameworks for heat transformation applications</article-title>. <source>Nanomaterials</source> <volume>8</volume> (<issue>9</issue>), <fpage>661</fpage>. <pub-id pub-id-type="doi">10.3390/nano8090661</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Dynamic analysis and control strategies of Organic Rankine Cycle system for waste heat recovery using zeotropic mixture as working fluid</article-title>. <source>Energy Convers. Manag.</source> <volume>192</volume>, <fpage>321</fpage>&#x2013;<lpage>334</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2019.04.049</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Delhommelle</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Milli&#xc9;</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Inadequacy of the Lorentz-Berthelot combining rules for accurate predictions of equilibrium properties by molecular simulation</article-title>. <source>Mol. Phys.</source> <volume>99</volume> (<issue>8</issue>), <fpage>619</fpage>&#x2013;<lpage>625</lpage>. <pub-id pub-id-type="doi">10.1080/00268970010020041</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fouad</surname>
<given-names>W. A.</given-names>
</name>
<name>
<surname>Vega</surname>
<given-names>L. F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Next generation of low global warming potential refrigerants: thermodynamic properties molecular modeling</article-title>. <source>AIChE J.</source> <volume>64</volume> (<issue>1</issue>), <fpage>250</fpage>&#x2013;<lpage>262</lpage>. <pub-id pub-id-type="doi">10.1002/aic.15859</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Furukawa</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Cordova</surname>
<given-names>K. E.</given-names>
</name>
<name>
<surname>O&#x2019;Keeffe</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yaghi</surname>
<given-names>O. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The chemistry and applications of metal-organic frameworks</article-title>. <source>Science</source> <volume>341</volume> (<issue>6149</issue>), <fpage>1230444</fpage>. <pub-id pub-id-type="doi">10.1126/science.1230444</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garc&#xed;a</surname>
<given-names>E. J.</given-names>
</name>
<name>
<surname>Bahamon</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Vega</surname>
<given-names>L. F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Systematic search of suitable metal&#x2013;organic frameworks for thermal energy-storage applications with low global warming potential refrigerants</article-title>. <source>ACS Sustain. Chem. Eng.</source> <volume>9</volume> (<issue>8</issue>), <fpage>3157</fpage>&#x2013;<lpage>3171</lpage>. <pub-id pub-id-type="doi">10.1021/acssuschemeng.0c07797</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hockney</surname>
<given-names>R. W.</given-names>
</name>
<name>
<surname>Goel</surname>
<given-names>S. P.</given-names>
</name>
<name>
<surname>Eastwood</surname>
<given-names>J. W.</given-names>
</name>
</person-group> (<year>1974</year>). <article-title>QUIET HIGH-RESOLUTION COMPUTER MODELS OF A PLASMA</article-title>. <source>J. Of Comput. Phys.</source> <volume>14</volume> (<issue>2</issue>), <fpage>148</fpage>&#x2013;<lpage>158</lpage>. <pub-id pub-id-type="doi">10.1016/0021-9991(74)90010-2</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Thermal energy storage of R1234yf/MOF-5 and R1234ze(Z)/MOF-5 nanofluids: a molecular simulation study</article-title>. <source>Energy Procedia</source> <volume>158</volume>, <fpage>4604</fpage>&#x2013;<lpage>4610</lpage>. <pub-id pub-id-type="doi">10.1016/j.egypro.2019.01.870</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2018a</year>). <article-title>Molecular simulation of thermal energy storage of mixed CO2/IRMOF-1 nanoparticle nanofluid</article-title>. <source>Int. J. Heat Mass Transf.</source> <volume>125</volume>, <fpage>1345</fpage>&#x2013;<lpage>1348</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijheatmasstransfer.2018.04.162</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2018b</year>). <article-title>Thermal energy storage of R1234yf, R1234ze, R134a and R32/MOF-74 nanofluids: a molecular simulation study</article-title>. <source>Materials</source> <volume>11</volume> (<issue>7</issue>), <fpage>1164</fpage>. <pub-id pub-id-type="doi">10.3390/ma11071164</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Swisher</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Smit</surname>
<given-names>B.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>Mapping of functional groups in metal-organic frameworks</article-title>. <source>Science</source> <volume>341</volume> (<issue>6148</issue>), <fpage>882</fpage>&#x2013;<lpage>885</lpage>. <pub-id pub-id-type="doi">10.1126/science.1238339</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lei</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Xi</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Enhancing selective adsorption of CO2 through encapsulating FeTPPs into Cu-BTC</article-title>. <source>Chem. Eng. J.</source> <volume>461</volume>, <fpage>141977</fpage>. <pub-id pub-id-type="doi">10.1016/j.cej.2023.141977</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Molecular simulation of energy storage of R1234yf, R1234ze(z), R32, and their mixtures in Co-MOF-74 materials</article-title>. <source>Chin. Sci. Bulletin-Chinese</source> <volume>65</volume> (<issue>7</issue>), <fpage>633</fpage>&#x2013;<lpage>640</lpage>. <pub-id pub-id-type="doi">10.1360/tb-2019-0435</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Prediction of critical properties and boiling point of fluorine/chlorine-containing refrigerants</article-title>. <source>Int. J. Refrig.</source> <volume>143</volume>, <fpage>28</fpage>&#x2013;<lpage>36</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijrefrig.2022.06.024</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Molecular simulation of energy storage properties of R32, R134A and R1234YF in MOF-5 AND MOF-177</article-title>. <source>Int. J. Mod. Phys. B</source> <volume>36</volume> (<issue>01</issue>), <fpage>2250006</fpage>. <pub-id pub-id-type="doi">10.1142/s0217979222500060</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Investigation on the condensation process of HFO refrigerants by molecular dynamics simulation</article-title>. <source>J. Mol. Liq.</source> <volume>288</volume>, <fpage>111034</fpage>. <pub-id pub-id-type="doi">10.1016/j.molliq.2019.111034</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McGrail</surname>
<given-names>B. P.</given-names>
</name>
<name>
<surname>Thallapally</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Blanchard</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Nune</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jenks</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Dang</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Metal-organic heat carrier nanofluids</article-title>. <source>Nano Energy</source> <volume>2</volume> (<issue>5</issue>), <fpage>845</fpage>&#x2013;<lpage>855</lpage>. <pub-id pub-id-type="doi">10.1016/j.nanoen.2013.02.007</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McLinden</surname>
<given-names>M. O.</given-names>
</name>
<name>
<surname>Huber</surname>
<given-names>M. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>(R)Evolution of refrigerants</article-title>. <source>J. Chem. Eng. Data</source> <volume>65</volume> (<issue>9</issue>), <fpage>4176</fpage>&#x2013;<lpage>4193</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jced.0c00338</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Selection criteria of zeotropic mixtures for subcritical organic Rankine cycle based on thermodynamic and thermo-economic analysis</article-title>. <source>Appl. Therm. Eng.</source> <volume>180</volume>, <fpage>115837</fpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2020.115837</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Niu</surname>
<given-names>F. F.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Pang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A first-principles study: adsorption of small gas molecules on GeP3 monolayer</article-title>. <source>Surf. Sci.</source> <volume>684</volume>, <fpage>37</fpage>&#x2013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1016/j.susc.2019.02.008</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Plimpton</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Fast parallel algorithms for short-range molecular dynamics</article-title>. <source>J. Comput. Phys.</source> <volume>117</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1006/jcph.1995.1039</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qasem</surname>
<given-names>N. A. A.</given-names>
</name>
<name>
<surname>Ben-Mansour</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Habib</surname>
<given-names>M. A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>An efficient CO2 adsorptive storage using MOF-5 and MOF-177</article-title>. <source>Appl. Energy</source> <volume>210</volume>, <fpage>317</fpage>&#x2013;<lpage>326</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2017.11.011</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Raabe</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Maginn</surname>
<given-names>E. J.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Molecular modeling of the Vapor&#x2212;Liquid equilibrium properties of the alternative refrigerant 2,3,3,3-Tetrafluoro-1-propene (HFO-1234yf)</article-title>. <source>J. Phys. Chem. Lett.</source> <volume>1</volume> (<issue>1</issue>), <fpage>93</fpage>&#x2013;<lpage>96</lpage>. <pub-id pub-id-type="doi">10.1021/jz900070h</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Su</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Recent advances in modeling the vapor-liquid equilibrium of mixed working fluids</article-title>. <source>Fluid Phase Equilibria</source> <volume>432</volume>, <fpage>28</fpage>&#x2013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1016/j.fluid.2016.10.016</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vivian</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Manente</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lazzaretto</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A general framework to select working fluid and configuration of ORCs for low-to-medium temperature heat sources</article-title>. <source>Appl. Energy</source> <volume>156</volume>, <fpage>727</fpage>&#x2013;<lpage>746</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2015.07.005</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Off-design performance evaluation of thermally integrated pumped thermal electricity storage systems with solar energy</article-title>. <source>Energy Convers. Manag.</source> <volume>301</volume>, <fpage>118001</fpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2023.118001</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Energy storage analysis of a mixed r161/MOF-5 nanoparticle nanofluid based on molecular simulations</article-title>. <source>Materials</source> <volume>11</volume> (<issue>5</issue>), <fpage>848</fpage>. <pub-id pub-id-type="doi">10.3390/ma11050848</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Huo</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Selection principle of working fluid for organic Rankine cycle based on environmental benefits and economic performance</article-title>. <source>Appl. Therm. Eng.</source> <volume>178</volume>, <fpage>115598</fpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2020.115598</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Huo</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Carbon footprint analysis of organic rankine cycle system using zeotropic mixtures considering leak of fluid</article-title>. <source>J. Clean. Prod.</source> <volume>239</volume>, <fpage>118095</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2019.118095</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="web">
<collab>webbook</collab> <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://webbook.nist.gov/chemistry/">https://webbook.nist.gov/chemistry/</ext-link>
</comment>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wilmer</surname>
<given-names>C. E.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>K. C.</given-names>
</name>
<name>
<surname>Snurr</surname>
<given-names>R. Q.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>An extended charge equilibration method</article-title>. <source>J. Of Phys. Chem. Lett.</source> <volume>3</volume> (<issue>19</issue>), <fpage>2506</fpage>&#x2013;<lpage>2511</lpage>. <pub-id pub-id-type="doi">10.1021/jz3008485</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiao</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Life cycle assessment of the solar thermal power plant integrated with air-cooled supercritical CO2 Brayton cycle</article-title>. <source>Renew. Energy</source> <volume>182</volume>, <fpage>119</fpage>&#x2013;<lpage>133</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2021.10.001</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>G. J.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Molecular simulations on CO2 adsorption and adsorptive separation in fullerene impregnated MOF-177, MOF-180 and MOF-200</article-title>. <source>Comput. Mater. Sci.</source> <volume>168</volume>, <fpage>58</fpage>&#x2013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1016/j.commatsci.2019.05.039</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Multi-factor evaluation method for the assessment of trans-critical organic Rankine cycle with low GWP fluids</article-title>. <source>Energy Procedia</source> <volume>158</volume>, <fpage>1626</fpage>&#x2013;<lpage>1631</lpage>. <pub-id pub-id-type="doi">10.1016/j.egypro.2019.01.380</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Huo</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Shale gas transport through the inorganic cylindrical and conical nanopores: a density gradient driven molecular dynamics</article-title>. <source>Int. J. Heat Mass Transf.</source> <volume>183</volume>, <fpage>122126</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijheatmasstransfer.2021.122126</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Barpaga</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Trump</surname>
<given-names>B. A.</given-names>
</name>
<name>
<surname>Shetty</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Bhattacharya</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Molecular insight into fluorocarbon adsorption in pore expanded metal&#x2013;organic framework analogs</article-title>. <source>J. Am. Chem. Soc.</source> <volume>142</volume> (<issue>6</issue>), <fpage>3002</fpage>&#x2013;<lpage>3012</lpage>. <pub-id pub-id-type="doi">10.1021/jacs.9b11963</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Energy storage analysis of UIO-66 and water mixed nanofluids: an experimental and theoretical study</article-title>. <source>Energies</source> <volume>12</volume> (<issue>13</issue>), <fpage>2521</fpage>. <pub-id pub-id-type="doi">10.3390/en12132521</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Thermodynamic analysis and parametric optimization of ejector heat pump integrated with organic Rankine cycle combined cooling, heating and power system using zeotropic mixtures</article-title>. <source>Appl. Therm. Eng.</source> <volume>194</volume>, <fpage>117097</fpage>. <pub-id pub-id-type="doi">10.1016/j.applthermaleng.2021.117097</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>