Modiﬁed Dual-Site Langmuir Adsorption Equilibrium Models from A GCMC Molecular Simulation

: In the modern industrial separation process, the pressure swing adsorption technology is widely used to separate and purify gases due to its low energy consumption, low cost, convenience, reliability, and environmental benignity. The basic elements of the design and application of the pressure swing adsorption process are adsorption isotherms at di ﬀ erent temperatures for adsorbents. The dual-site Langmuir (DSL) adsorption equilibrium model is the mostly used model; however, this model is based on the assumption that the adsorption energy on the surface of an adsorbent is uniform and remains unchanged. Here, a grand canonical Monte Carlo (GCMC) molecular simulation was used to calculate the CO 2 adsorption equilibrium on MIL-101 (Cr) at 298 K. MIL-101 (Cr) was chosen, as it has more a general pore structure with three di ﬀ erent pores. The calculation results showed that the adsorption energies with di ﬀ erent adsorption pressures ﬁtted a normal distribution and the relationship of the average adsorption energies, E with pressures had a linear form described as: E = aP + c. With this relationship, the parameter b = k · exp (E / RT) in the DSL model was modiﬁed to b = k · exp ((aP + c) / RT) , and the modiﬁed DSL model (M-DSL) was used to correlate the adsorption equilibrium data on CO 2 -MIL-101 (Cr), C 2 H 4 -HHPAC, CH 4 -BPL, and CO 2 -H-Mordenite, showing better correlations than those of the DSL model. We also extended the parameter q m in the M-DSL model with the equation q m = k 1 + k 2 T to adsorption equilibrium data for di ﬀ erent temperatures. The obtained model (M-TDSL) was checked with the abovementioned adsorption equilibrium systems. The ﬁtting results also indicated that the M-TDSL model could be used to improve the correlation of adsorption equilibrium data for di ﬀ erent temperatures. The linear relationship between the average adsorption energy and adsorption pressure could be further tested in other adsorption equilibrium models to determine its universality. the equilibrium data on CO 2 -MIL-101 (Cr), C 2 H 4 CH -BPL, and CO 2 to investigate their and generality and provide a basis for improving other adsorption equilibrium models.


Introduction
To date, with the development of adsorption theory and the improvement of adsorption separation technology, the pressure swing adsorption (PSA) technology has been widely used in petrochemical industry, biomedicine, energy, and environmental protection. In the modern industrial separation process, the PSA technology is used to separate and purify gases due to its low energy consumption, low cost, convenience, reliability, and environmental benignity [1][2][3][4]. For example, the PSA process improved by Augelletti et al. was a feasible process to obtain an almost complete separation of biogas components, and a recovery of more than 99% of methane was achieved [5]. The PSA technology also plays an important role in the field of liquid separation. The PSA process could be used to isolate individual bioactive compounds from complex plant extracts and recover solvents, which open new ways for the fractionation of various valuable natural compounds [6,7]. In addition, PSA mathematical models were developed to assist with the process optimization, reducing the need for excessive experimentation [8]. The basic elements of design and application of the PSA process are adsorption isotherms at different temperatures for adsorbents. To date, adsorption equilibrium models that have been used to describe adsorption isotherms have mainly been empirical or semiempirical formulas, which have some disadvantages. Among these models, the dual-site Langmuir (DSL) adsorption equilibrium model is the mostly used [9][10][11][12].
The DSL model is based on the Langmuir model, which assumes that there are two kinds of adsorption sites with different adsorption energies, and the adsorption behavior at each site is consistent with three assumptions of the Langmuir model: (1) the adsorption of gas molecules on a solid surface forms a monolayer; (2) the adsorption energy on the surface of an adsorbent is uniform and remains unchanged; and (3) there is no interaction between adsorbate molecules [13][14][15]. These assumptions make the DSL model suitable for most systems and applicable under different temperatures and pressures. For example, the adsorption module in Aspenplus, the chemical engineering simulation software, uses the DSL model to reflect temperature changes to correlate various adsorption isotherms [16][17][18]. However, the applicability of the DSL model remains limited. The results it provides on the regression correlation effect of the adsorption equilibrium data on materials with a complex pore structure are insufficiently accurate, and it cannot well describe the adsorption equilibrium characteristics of low-temperature and high-temperature ranges and the high-pressure ranges of each temperature range. The main reason for this is the underlying assumption that the adsorption energy on the surface of an adsorbent is uniform and constant [19,20]. In fact, in the actual adsorption process, the characteristic change in adsorption energy occurs via a complex process that is related to the adsorption pressure. Under most conditions, it cannot be simply assumed that it is constant. To date, few studies have focused on the characteristics of adsorption energy with changes in pressure, and, at present, the research on the characteristics of adsorption energy at the micro level is insufficiently deep.
Adsorption energy is a microscopic physical quantity, which needs to be obtained by molecular simulation. Molecular simulation is widely used to calculate the adsorption properties of various new materials, such as metal-organic frameworks (MOFs). MOFs with extremely large pore volumes and high surface areas are highly potential candidates as adsorbents on gas adsorption. For example, Park et al. examined upper bounds on CO 2 swing capacity existing via molecular simulation of a large collection of MOFs, which enabled the formulation of material selection criteria that allow for the efficient identification of materials with high CO 2 swing capacity [21]. In this study, a grand canonical Monte Carlo (GCMC) molecular simulation was used to calculate the CO 2 adsorption equilibrium on a MOF material MIL-101 (Cr) [22][23][24] to explore characteristics of changes in adsorption energy. There are two reasons for the selection of MIL-101 (Cr) as an adsorbent in this study: (1) MIL-101 (Cr) has three different types of pores, which has microporous and mesoporous characteristics, and the structure is relatively complex, so it is more meaningful to calculate its adsorption energy changes, and the results have universal applicability for the adsorption energy changes of other simple systems such as activated carbon and molecular sieves; (2) MIL-101 (Cr) has a high adsorption capacity for CO 2 , and the adsorption capacity varies greatly under different pressures, so it is helpful to research the characteristic changes of adsorption energy with pressure during the adsorption process. The results that were obtained at the molecular simulation level were used to modify the DSL model and the TDSL model. The modified models were used to correlate the adsorption equilibrium data on CO 2 -MIL-101 (Cr), C 2 H 4 -HHPAC, CH 4 -BPL, and CO 2 -H-Mordenite to investigate their accuracy and generality and provide a basis for improving other adsorption equilibrium models.

Experimental
In the past, the studies of adsorption equilibrium models were often limited by data on the adsorption equilibrium in different temperature and pressure ranges. Therefore, the experimental data on ethylene's adsorption equilibrium in powder-activated carbon (HHPAC) were obtained in a temperature range from 223.15 to 313.15 K and a pressure range from 0 to 1.5 MPa.

Determination of Adsorption Isotherms
HHPAC was pretreated in a vacuum at 250 • C for 6 h. The N 2 adsorption desorption isotherm was measured at 77 K. As shown in Table 1, the specific surface area and pore size were calculated by the Brunauer-Emmett-Teller (BET) equation. HHPAC was loaded into the sample cell of the full-automatic volumetric high-pressure adsorption instrument. Then, the sample cell was put into the low-temperature constant-temperature tank. The target temperature of the low-temperature constant-temperature tank and the corresponding parameters of adsorbates were set. The adsorption balance of C 2 H 4 of HHPAC was measured. The adsorption equilibriums at 223. 15

Model of MIL-101 (Cr)
The unit cell framework of MIL-101 (Cr) (as shown in Figure 1) was obtained from the Cambridge Crystallographic Data Centre (CCDC). Material Studio was used to optimize the structure of a crystal cell and convert the original crystal cell to a rhombohedral crystal cell. The volume of the transformed cell was only one-sixth of that of the original cell, but it was still able to accurately express the structural information on MIL-101 (Cr) and reduce the number of required calculations in the molecular simulation. The specific surface area, pore volume, and porosity of MIL-101 (Cr) were calculated by the simulation program RASPA. The structural properties of MIL-101 (Cr) are summarized in Table 2.

Experimental
In the past, the studies of adsorption equilibrium models were often limited by data on the adsorption equilibrium in different temperature and pressure ranges. Therefore, the experimental data on ethylene's adsorption equilibrium in powder-activated carbon (HHPAC) were obtained in a temperature range from 223.15 to 313.15 K and a pressure range from 0 to 1.5 MPa.

Determination of Adsorption Isotherms
HHPAC was pretreated in a vacuum at 250 °C for 6 h. The N2 adsorption desorption isotherm was measured at 77 K. As shown in Table 1, the specific surface area and pore size were calculated by the Brunauer-Emmett-Teller (BET) equation. HHPAC was loaded into the sample cell of the full-automatic volumetric high-pressure adsorption instrument. Then, the sample cell was put into the low-temperature constant-temperature tank. The target temperature of the low-temperature constant-temperature tank and the corresponding parameters of adsorbates were set. The adsorption balance of C2H4 of HHPAC was measured. The adsorption equilibriums at 223. 15 Sample Specific surface Area SBET (cm 2 · g −1 ) Pore Volume (cm 3 · g −1 ) HHPAC 750 0.55

Model of MIL-101 (Cr)
The unit cell framework of MIL-101 (Cr) (as shown in Figure 1) was obtained from the Cambridge Crystallographic Data Centre (CCDC). Material Studio was used to optimize the structure of a crystal cell and convert the original crystal cell to a rhombohedral crystal cell. The volume of the transformed cell was only one-sixth of that of the original cell, but it was still able to accurately express the structural information on MIL-101 (Cr) and reduce the number of required calculations in the molecular simulation. The specific surface area, pore volume, and porosity of MIL-101 (Cr) were calculated by the simulation program RASPA. The structural properties of MIL-101 (Cr) are summarized in Table 2.

Force Field
Since the main purpose of this study is to determine the characteristics of changes in the adsorption energy in the adsorption process, the interactions and energy distributions between adsorbates and adsorbents were explored. When selecting a force field, interactions between adsorbates can be ignored to simplify calculations. On this basis, the Lennard-Jones (LJ) potential energy model was selected to describe the van der Waals force in the adsorption process, and the Ewald summation method was used to describe the electrostatic force. The parameters of cross-interactions between different atoms were calculated using the Lorentz−Berthelot mixing rules. The potential parameter values for CO 2 molecules were taken from TraPPE, and the potential parameter values for the MIL-101 (Cr) framework were taken from the universal force field (UFF). According to the literature [25,26], the values were adjusted in order to accurately describe the CO 2 adsorption process on MIL-101 (Cr). The force field parameters are shown in Table 3

Electric Charge
When calculating the CO 2 adsorption, the charge of MIL-101 (Cr) cannot be ignored because of its strong quadrupole moment. The DMol3 module in Materials Studio was used to calculate the atomic charge. Because MIL-101 (Cr) has an ordered symmetry structure, the intercepted cluster was used as a model for calculating the charge (as shown in Figure 2), and the truncated part was calculated by the methyl saturation method [27][28][29]. The charges of the adsorbents and the adsorbates are shown in Table 4.

Force Field
Since the main purpose of this study is to determine the characteristics of changes in the adsorption energy in the adsorption process, the interactions and energy distributions between adsorbates and adsorbents were explored. When selecting a force field, interactions between adsorbates can be ignored to simplify calculations. On this basis, the Lennard-Jones (LJ) potential energy model was selected to describe the van der Waals force in the adsorption process, and the Ewald summation method was used to describe the electrostatic force. The parameters of cross-interactions between different atoms were calculated using the Lorentz−Berthelot mixing rules. The potential parameter values for CO2 molecules were taken from TraPPE, and the potential parameter values for the MIL-101 (Cr) framework were taken from the universal force field (UFF). According to the literature [25,26], the values were adjusted in order to accurately describe the CO2 adsorption process on MIL-101 (Cr). The force field parameters are shown in Table 3   Table 3. Force field parameters.

Electric Charge
When calculating the CO2 adsorption, the charge of MIL-101 (Cr) cannot be ignored because of its strong quadrupole moment. The DMol3 module in Materials Studio was used to calculate the atomic charge. Because MIL-101 (Cr) has an ordered symmetry structure, the intercepted cluster was used as a model for calculating the charge (as shown in Figure 2), and the truncated part was calculated by the methyl saturation method [27][28][29]. The charges of the adsorbents and the adsorbates are shown in Table 4.

Simulation Methods
The adsorption characteristics of CO 2 were determined using Monte Carlo (MC) simulations in a grand canonical ensemble (µVT) with the RASPA simulation program for adsorption properties of porous materials [30][31][32]. MC is a method, which combines the probability distribution principle of statistical mechanics and constructs different molecular configurations through the random movement of particles in a system to obtain the statistical thermodynamic information of the system, and it can be applied to a variety of systems. In the grand canonical ensemble, the chemical potential, temperature, and volume were kept constant, while the numbers of particles in the system were changed. The ensemble consisted of a large number of open systems, each of which was in an internal equilibrium state and could exchange matter (particles) and energy with other boxes. Therefore, the characteristics of the grand canonical ensemble are especially suitable for gas adsorption simulation. Chemical potentials were converted to pressures with the Peng−Robinson equation of state (PR-EOS). The form of the PR-EOS was shown as follows: where P is the adsorption pressure, R is the gas constant, T is adsorption temperature, V is the volume of the gas, a is a parameter related to the temperature obtained by RASPA, and b is a parameter related to the chemical potential obtained by RASPA.
The adsorption isotherms that are made by different adsorption capacities can be obtained at different pressures. In general, the experimental measurement represents the excess adsorbed amount, and the GCMC simulation result represents the absolute adsorbed amount. In this paper, the absolute adsorbed amounts that were obtained by GCMC simulations were converted to excess adsorbed amounts in order to perform a comparison with the experimental data, and the amount of excess adsorbed molecules N ex was calculated using Equation (2): where N ab is the amount of absolute adsorbed molecules, V p is the pore volume of the adsorbent, and ρ is the density of the adsorbate calculated using the PR-EOS at a given temperature and pressure.
In the simulation, a rigid structure model was adopted for the adsorbate molecules and the MIL-101 (Cr) framework. That is to say, the bond length, bond angle, and atomic model of the adsorbate molecules did not change during the adsorption process, and the atom framework of the adsorbent MIL-101 (Cr) did not shift during the adsorption process. The numbers of unit cells of the frameworks that were adopted in the simulation were 1 × 1 × 1, and periodic boundary conditions were applied in all three dimensions. In this work, a total of 2 × 10 6 steps were used. Half of these steps were used for equilibration, and the remaining steps were used for subsequent calculations and result statistics. For the CO 2 molecules, four types of moves (random insertion, deletion, translation, and rotation) were used, and every possible move was given equal probability. In this work, the selected adsorption pressure ranged from 1 to 300 kPa. The adsorption capacities at 298 K under different pressures were calculated and compared with the experimental data to verify the reliability of the simulation results. Then, the pressures of 10, 20, 50, 100, 200, and 300 kPa were used to calculate the characteristics and distributions of adsorption energy. In addition, Since the MC method is a random sampling method to obtain statistical averages, there may be calculation errors, so the calculation of adsorption balance data of each pressure point was repeated for three times, and the standard deviations of adsorption capacity were calculated to judge the accuracy of the method and data.

Adsorption Capacity and Adsorption Energy
To verify the reliability of simulation results, the simulated adsorption isotherms of CO 2 on MIL-101 (Cr) at 298 K were compared to experimental data [33,34]. As shown in Figure 3, the simulated adsorption isotherms were slightly higher than the experimental values, and the deviation from the experimental value increased with the increase of pressure. These phenomena provide reasonable explanations: (1) the model of MIL-101 (Cr) that was used in the calculations contained no crystal defects and no impurities. However, actual MIL-101 (Cr) may have defects, and its pores may contain residual terephthalic acid, which can result in a relatively low CO 2 adsorption capacity; (2) when the pressure was low, the CO 2 adsorption capacity on MIL-101 (Cr) was also low, CO 2 molecules only occupied a small amount of space, and the residual impurities in the pores had little effect on the adsorption capacity. With the increase of pressure, the adsorption capacity also increased, and CO 2 molecules occupied more and more space in the pores. The existence of impurities hindered the adsorption of CO 2 , resulting in the increase of deviation between the experimental and simulated values. In addition, Table 5 shows the standard deviations of adsorption capacity. The value of the standard deviations of adsorption capacity are all less than 0.2, which indicated that the results could be repeated and the MC method could be applied to this system. Therefore, this simulation method, which is based on an MC simulation in a grand canonical ensemble, was utilized for further study, and the simulation results were considered to be accurate and reliable.    When calculating the adsorption isotherm, the adsorption energy reflecting the strength of the interaction between an adsorbent and an adsorbate can also be obtained. In this process, the adsorption energy was not constant. As summarized in Table 5, the absolute value of the adsorption energy decreased, as the pressure increased, indicating that the interaction between the adsorbent and the adsorbate decreased as the pressure and adsorption capacity increased. Figure 4 shows the characteristics and distributions of the adsorption energy. The data points in Figure 4 represent the distribution points of the adsorption energy characteristics determined by GCMC molecular simulations, and the curve represents the curve fitted by a normal distribution model in the mathematical software MATLAB to correlate all data points. It can be seen in Figure 4 that the normal distribution model had a good regression effect for the data points. Under different pressures, the values of the distributions of adsorption energy were different, but they all obeyed a normal distribution. The purpose of this study is to explore the relationship between adsorption energy and pressure, so the average value of the adsorption energy was used in subsequent work and is shown as the peak value of the normal distribution ( Figure 4). As the pressure increased, the peak value of the normal distribution shifted to the right, and the distribution area of adsorption energy narrowed, that is, the absolute average value of adsorption energy decreased. This result indicated that the adsorption was gradually weakened, which was also consistent with the change in adsorption energy shown in Table 5. From the above results, it can be seen that the average adsorption energy was correlated with the adsorption pressure, that is, there was a linear relationship between them that can be expressed as follows: where E is the average adsorption energy, P is the adsorption pressure, and a and c are the coefficients of the linear relationship between the adsorption energy and the adsorption pressure.
Appl. Sci. 2020, 10, x 8 of 17 adsorption energy was correlated with the adsorption pressure, that is, there was a linear relationship between them that can be expressed as follows: E aP c  ,  Figure 5 shows the fitting correlation curve of the average adsorption energy data points under different pressures using Equation (3). It can be seen that Equation (3) can well describe the relationship between the adsorption energy and pressure. It also can be seen from Table 6 that the determination coefficient (R 2 ) for linear regression was greater than 0.98 and the residual sum of squared estimate of errors (SSE) was only 0.12. In conclusion, there was a linear relationship between the adsorption energy and the adsorption pressure, which can be written as Equation (3).

Modification of the DSL Model
The form of the DSL adsorption equilibrium model [3,[35][36][37] was shown as follows: where P is the adsorption equilibrium pressure, q is the adsorption capacity at an equilibrium pressure, qm1 and qm2 represent the saturated adsorption capacities at adsorption sites 1 and 2, respectively, and b1 and b2 represent the Langmuir adsorption equilibrium constants of adsorption site 1 and 2, respectively, which can be written as follows: where ki is a pre-exponential factor, R is the gas constant, and T is the adsorption temperature. Equation (6) can be obtained by combining Equations (3) and (5): By substituting Equation (6)

Modification of the DSL Model
The form of the DSL adsorption equilibrium model [3,[35][36][37] was shown as follows: where P is the adsorption equilibrium pressure, q is the adsorption capacity at an equilibrium pressure, q m1 and q m2 represent the saturated adsorption capacities at adsorption sites 1 and 2, respectively, and b 1 and b 2 represent the Langmuir adsorption equilibrium constants of adsorption site 1 and 2, respectively, which can be written as follows: where k i is a pre-exponential factor, R is the gas constant, and T is the adsorption temperature. Equation (6) can be obtained by combining Equations (3) and (5): By substituting Equation (6) into Equation (4), a new modified DSL model (M-DSL) without the adsorption energy E as a parameter can be obtained. The form of the M-DSL model was described as follows: Because there are many parameters in Equation (7), it can be simplified by combining relevant parameters. The simplification was expressed in Equation (8), which was the simplified form of the M-DSL model: q = q m1 k m1 exp(a m1 P)P 1 + k m1 exp(a m1 P)P + q m2 k m2 exp(a m2 P)P 1 + k m2 exp(a m2 P)P , where k m1 and k m2 are two new parameters that are obtained by multiplying k i by e ci/RT , and a m1 and a m2 are two new coefficients of the linear relationship between the adsorption energy and pressure.

Regression of Adsorption Equilibrium Data with the DSL Model and the M-DSL Model
The DSL model and the M-DSL model were used to regress the adsorption equilibrium data by the nonlinear regression method, and the corresponding model parameters were obtained. The nonlinear regression was realized by the lsqnonlin function in MATLAB, and we used the Levenberg-Marquardt algorithm. The results were calculated with Equation (9), which was written as: where i is the number of experimental data points, N is the total number of experimental data points, q i,exp is the experimental value, and q i,pre is the regression value. The lsqnonlin function in MATLAB is based on the nonlinear least squares method. In the process of regression using this method, the selection of an initial value is very important, different initial values may obtain different results.
The specific operation process is as follows: (1) an appropriate initial value is selected according to experience of regression calculation; (2) the results of the first regression are taken as the initial values of regression; and (3) regression is carried out, until the obtained model parameter value no longer changes. After the operation process is completed, the R 2 of the model reaches the maximum, and the SSE has the minimum value, and the results are the final model parameter values. Therefore, using this method to regress the same system could get the same results, which showed that the results had good repeatability, and the model parameters obtained by this method were considered to be accurate and reasonable. The accuracy of the regression was evaluated by the square sum of the error between the original experimental data and the fitting data, that is, the SSE, which was shown in Equation (10): where i is the number of experimental data points, N is the total number of experimental data points, y i,exp is the experimental value, and y i, pre is the regression value. The closer the SSE is to 0, the closer the predicted values of the model are to the experimental values, and the better the model describes the experimental data. The square of the correlation coefficient (R 2 ) can also be used as an evaluation index of data fitting. It is defined as the ratio of the sum of squares of the fitting deviation to the sum of squares of the actual deviation and can be calculated with Equation (11): where i is the number of experimental data points, N is the total number of experimental data points, y i,exp is the experimental value, y i,pre is the regression value, andȳ is the average number of experimental values used to determine the adsorption capacity. R 2 ranges from 0 to 1. The smaller the value is, the worse the fit between the model and the data is; the closer the value is to 1, the better the fit of the model is. In this study, SSE and R 2 were combined in order to analyze the regression results. We used the adsorption equilibrium data on C 2 H 4 -HHPAC, CO 2 -MIL-101 (Cr), CH 4 -BPL active carbon, and a CO 2 -H-Mordenite molecular sieve [31,32], as determined by experiments, in the DSL model and the M-DSL model for the nonlinear regression. The reason why we selected these data is that they have a wide range of adsorption temperatures and pressures and are commonly used to test adsorption equilibrium models. In addition, the adsorbents used in these systems are commonly used in industry, among which BPL is cylindrical, H-Mordenite is spherical and HHPAC is powder. Therefore, the shapes of adsorbents selected in this study almost contained the shapes of adsorbents commonly used in industry. Similarly, the adsorbates in these systems contained common organic gases such as CH 4 and C 2 H 4 and common inorganic gas CO 2 . Therefore, it could be said that the adsorption systems studied in this paper were sufficiently comprehensive, and the conclusions could be applied to most systems in industry.
Supplementary Figure S1 compares the correlation curves of the DSL model and the M-DSL model for different systems. Tables 7 and 8 show the corresponding parameter values that were obtained by the regression of the adsorption equilibrium data on each system by using the DSL and M-DSL models, respectively. As shown in Figure S1, for the C 2 H 4 -HHPAC system, the DSL model was able to relate the entire adsorption isotherm, and the regression effects were quite different under different temperatures and pressures. From Tables 7 and 8, when the temperature was 253.15, 273.15, and 293.15 K, the DSL model was able to correlate the adsorption equilibrium data, which verified the applicability of the DSL model in the normal temperature range. However, the M-DSL model can better describe the adsorption isotherms of the system at these temperatures. The value of R 2 was greater than 0.998, and the SSE value was less than 0.1. By calculating the change in the SSE, we can see that the regression accuracy of the data with the M-DSL model was approximately 30% higher than that with the DSL model. These results indicated that the M-DSL model inherited the advantages of the DSL model in the normal temperature range and had better applicability. When the temperatures were 223.15 and 233.15 K, the adsorption pressure gradually increased close to the saturated vapor pressure of C 2 H 4 , and, during this time, the adsorption capacity increased rapidly. At the temperature of 313.15 K, the adsorption capacity decreased significantly. These results showed that the DSL model was unable to well correlate the adsorption equilibrium data on C 2 H 4 -HHPAC in the low-temperature range and the high-temperature range. Under these conditions, the R 2 of the DSL model was only 0.98, and the SSE was greater than 0.8. However, the M-DSL model can accurately describe adsorption isotherms in the low-temperature range. The regression accuracy of the M-DSL model to the adsorption equilibrium data was improved by more than 50%, the R 2 was greater than 0.995, and the SSE was less than 0.3. In addition, the M-DSL model had obvious advantages in the high-temperature range. When the temperature was 313.15 K, the correlation curve of the M-DSL model was smoother than that of the DSL model. The R 2 of the M-DSL model increased to 0.999, the SSE decreased to 0.03, and the regression accuracy increased by 94%. We drew similar conclusions for the CH 4 -BPL activated carbon and CO 2 -H-Mordenite systems. The regression accuracies of the M-DSL model and the DSL model were 22% and 8%, respectively. In conclusion, in ranges, in which large changes in adsorption capacity occur, such as high-temperature or low-temperature ranges, the adsorption energy cannot be simply assumed to be constant. The M-DSL model, which takes into account the linear relationship between adsorption energy and adsorption pressure, has a wider application range than the DSL model. It is also able to correlate adsorption equilibrium data well in the low-temperature range, the high-temperature range, and a relatively high-pressure range, where the DSL model is inapplicable.

Adsorbent
Adsorbate Regarding the MIL-101 (Cr)-CO 2 system, because the pore structure of the adsorbent was more complex than that of carbon materials and molecular sieves, the regression effect of the DSL model was poor. The value of R 2 was only 0.98, which cannot produce the expected effect. The regression effect of CO2-MIL-101 (Cr) was greatly improved when using the M-DSL model. The R 2 value was 0.9999, which is close to its upper limit (1), and the SSE was also very small (almost 0). The accuracy of the regression was 98% higher than that of the original DSL model. This result showed that the M-DSL model better correlated the adsorption equilibrium data on systems with a complex structure. Considering the linear relationship between the adsorption energy and pressure in an adsorption system with a complex pore structure can greatly improve the accuracy of the DSL model and improve the universality of the DSL model.

Modification of the TDSL Model That Reflects Temperature Changes
In the model parameter regression, we found that the parameters that were obtained by the DSL model and the M-DSL model were quite different at different temperatures and that the adsorption energy, which was found to have a linear relationship with the pressure obtained by a GCMC simulation, can modify the DSL model at different temperatures and had a good effect. In order to further verify this conclusion, we chose the DSL model that reflects temperature changes (TDSL), which is used in the adsorption module of Aspenplus and has been used in previous research [16][17][18], and attempted to improve its accuracy and universality. The TDSL model is a two-parameter adsorption equilibrium model that is based on the DSL model and takes temperature and pressure as variables together.
It is also based on the three assumptions of the Langmuir model. With this model, the adsorption equilibrium data obtained at different temperatures in the same system can be correlated at the same time. The form of the TDSL model is shown in Equation (12): where P is the adsorption equilibrium pressure, q is the adsorption capacity at the equilibrium pressure, R is the gas constant, T is the adsorption temperature, k 1 , k 2 , k 5 , and k 6 are the coefficients of the linear relationship between the saturated adsorption capacity and temperature in the original DSL model, k 3 and k 7 are the adsorption energies corresponding to different adsorption sites, and k 4 and k 8 are pre-exponential factors.
Based on the conclusion that adsorption energy has a linear relationship with pressure, the TDSL model can be modified to obtain the M-TDSL model by substituting Equation (5) into Equation (12). The form of the modified TDSL (M-TDSL) model is shown in Equation (13): where a 1 , c 1 , a 2 , and c 2 are coefficients of the linear relationship between the adsorption energy and the adsorption pressure. The TDSL model and the obtained M-TDSL model were used to correlate the adsorption equilibrium data of the above systems at different temperatures at the same time.
The results are shown in Tables 9 and 10.  Supplementary Figure S2 shows the correlation curves of the TDSL model and the M-TDSL model for different systems. Tables 9 and 10 show the corresponding parameter values that were obtained by the regression of the adsorption equilibrium data on each system by using the TDSL and M-TDSL models, respectively. Overall, the TDSL model was able to correlate the adsorption equilibrium data on CO 2 -HHPAC well. Its R 2 was greater than 0.99, and its SSE was in a reasonable range. From the physical meaning of the parameters, it can be seen that the values of k 3 and k 7 were the adsorption energies at the two adsorption sites, respectively, and the results also showed that the adsorption energies were −19.98 and −28.99 kJ·mol −1 , respectively, which are in the general range of physical adsorption energies. This indicated that the regression results were reliable and accurate and that the physical meanings of the corresponding parameters were clear and definite. The regression effect of the M-TDSL model on the adsorption equilibrium data was slightly improved as compared to that of the original model, and the accuracy of the regression was increased by 14%. For the CH 4 -BPL and CO 2 -H-Mordenite systems, the TDSL model had a poor correlation, which cannot produce the expected effect. In fact, the M-TDSL model was better than the TDSL model with respect to the fitting correlation effects of these systems. Its R 2 was greater than 0.99, its SSE was significantly smaller than that of the TDSL model, and the regression accuracies of the M-TDSL and TDSL were also improved by 23% and 31%, respectively, compared with that of the original TDSL model. In the M-TDSL model, the two adsorption energies were numerically equal to the values of (c 1 + a 1 P) and (c 2 + a 2 P), both of which were in the general range of physical adsorption energies. In conclusion, the average regression accuracy when using the M-TDSL model was 23% higher than that when using the original TDSL model. Therefore, the obtained M-TDSL model had higher precision and a wider application range than the TDSL model, and the model parameters also had clear physical significance.
In addition, the results also showed that our conclusion-adsorption energy has a linear relationship with pressure-has certain universality across different adsorption equilibrium models. This relationship provides a basis for modifying other forms of adsorption equilibrium models and can be further tested when modifying other adsorption equilibrium models.

Conclusions
In this study, a GCMC molecular simulation was used to calculate the CO 2 adsorption equilibrium on MIL-101 (Cr) at 298 K, and the characteristics and distribution of adsorption energy with pressure were obtained. The calculated results showed that the adsorption energies with different adsorption pressures fitted a normal distribution and the relationship of average adsorption energies with pressures had a linear form. The linear relationship between the adsorption energies and pressure was used to modify the parameters of the DSL model, and the M-DSL model was used to correlate the adsorption equilibrium data on CO 2 -MIL-101 (Cr), C 2 H 4 -HHPAC, CH 4 -BPL, and CO 2 -H-Mordenite. The results showed that the accuracy and generality of the M-DSL model were significantly improved after the modification and the average regression accuracy of the correlation of adsorption equilibrium data on each system was increased by 38%. The M-DSL model can better correlate the adsorption equilibrium data in high-or low-temperature ranges, and its correlation effect in the normal temperature range was also better than that of the DSL model. In addition, the M-DSL model showed a better correlation for systems with a more complex pore structure, such as CO 2 -MIL-101 (Cr). Similar to the DSL model, the linear relationship between the adsorption energy and pressure was also used to modify the TDSL model. The obtained M-TDSL model was checked with the abovementioned adsorption equilibrium systems. The results also indicated that the M-TDSL model could be used to improve the correlation of adsorption equilibrium data obtained at different temperatures. The average accuracy of the correlation for each system was increased by 23%, and the physical meanings of the parameters in the M-DSL model were clear and definite. In the process of practical applications, the modified models would increase the number of required calculations due to their more complex form, which will become a limitation. However, their good accuracy and versatility can make up for this disadvantage. Finally, our conclusion that adsorption energy has a linear relationship with pressure provides a basis for modifying other types of adsorption equilibrium models. The linear relationship between the adsorption energy and pressure could be further tested in other adsorption equilibrium models to determine its universality.