NLDFT density functional theory

Thus, either type I or type IV isotherms are obtained in sorption experiments on microporous or mesoporous materials. Of course, a material may contain both types of pores. In this case, a convolution of a type I and type IV isotherm is observed. From the amount of gas that is adsorbed in the micropores of a material, the micropore volume is directly accessible (e.g., from t plot of as plot [1]). The low-pressure part of the isotherm also contains information on the pore size distribution of a given material. Several methods have been proposed for this purpose (e.g., Horvath-Kawazoe method) but most of them give only rough estimates of the real pore sizes. Recently, nonlocal density functional theory (NLDFT) was employed to calculate model isotherms for specific materials with defined pore geometries. From such model isotherms, the calculation of more realistic pore size distributions seems to be feasible provided that appropriate model isotherms are available. The mesopore volume of a mesoporous material is also rather easy accessible. Barrett, Joyner, and Halenda (BJH) developed a method based on the Kelvin equation which allows the calculation of the mesopore size distribution and respective pore volume. Unfortunately, the BJH algorithm underestimates pore diameters, especially at... 

Comparison of the pore size distribution determined by the present method with that from the classical methods such as the BJH, the Broekhoff-de Boer and the Saito-Foley methods is shown in Figure 4. Figure 5 shows a close resemblance of the results of our method with those from the recent NLDFT of Niemark et al. [16], and XRD pore diameter for their sample AMI. The results clearly indicate the utility of our method and accuracy comparable to the much more computationally demanding density functional theory. There are several other methods published recently (e. g. [21]), however space limitations do not permit comparison with these results here. It is hoped to discuss these in a future publication. 

Essential progress has been made recently in the area of molecular level modeling of capillary condensation. The methods of grand canonical Monte Carlo (GCMC) simulations [4], molecular dynamics (MD) [5], and density functional theory (DFT) [6] are capable of generating hysteresis loops for sorption of simple fluids in model pores. In our previous publications (see [7] and references therein), we have shown that the non-local density functional theory (NLDFT) with properly chosen parameters of fluid-fluid and fluid-solid intermolecular interactions quantitatively predicts desorption branches of hysteretic isotherms of nitrogen and argon on reference MCM-41 samples with pore channels narrower than 5 nm. 

The non-local density functional theory (NLDFT) with properly chosen parameters of fluid-fluid and fluid-solid intermolecular interactions quantitatively predicts both adsorption and desorption branches of capillary condensation isotherms on MCM-41 materials with the pore sizes from 5 to 10 nm. Both experimental branches can be used for calculating the pore size distributions in this pore size range. However for the samples with smaller pores, the desorption branch has an advantage of being theoretically accurate. Thus, we recommend to use the desorption isotherms for estimating the pore size distributions in mesoporous materials of MCM-41 type, provided that the pore networking effects are absent. 

The nonlocal density functional theory (NLDFT) has been used to characterize mesocaged structures.[164] NLDFT analysis gives accurate information about the cage size, the total meso- and micropore volumes and surface area, and the pore-wall thickness in combination with XRD measurements. Argon- and nitrogen-desorption data on FDU-1 provided evidence that there are two major populations of pore entrances. Argon desorption was superior in providing information about pore connectivity in FDU-1 samples. 

Because of the present availability of commercial software, the nonlocal version of density functional theory (i.e., NLDFT) is now widely used for pore... 

G activated carbon (points) with the fit given by the nonlocal density functional theory (NLDFT) models (line), (b) The pore width distribution for the carbon. 

Figure 7.10 The average pore fluid density in pores of various widths as calculated hy nonlocal density functional theory (NLDFT). Note the periodic nature of the density, with density maxima near the positions of the minima in the distrihutions shown in Figs 7.8(h) and 7.9(h).

The solid line is calculated by the nonlocal density functional theory (NLDFT), which will be described to some detail in Section 11.4. As seen in the figure, the curve obtained with the HK method (dashed fine) correlates with NLDFT much better than that calculated with the Kelvin equation (dash-dotted line). 

Density functional theory is a powerful tool to study many phenomena in physical chemistry and chemical engineering. It was popularized in the early 1960s by a number of authors [72-74]. But it is not until the 1980s that this theory had found widespread appHcations in many interfacial problems. Capillary condensation in pore was systematically studied [75], and the first paper [76] applying this technique to the problem of PSD determination of carbon particle appeared in 1989. This work used a local DFT, and it is now superseded by the NLDFT, which was developed by Tarazona and Evans [77-79]. This is the method that is now widely used in the characterization of pore size distribution. 

This would allow performing accurate PSD calculations using these simple algorithms. Theoretical considerations [13], nonlocal density functional theory (NLDFT) calculations [62, 146], computer simulations [147], and studies of the model adsorbents [63, 88] strongly suggested that the Kelvin equation commonly used to provide a relation between the capillary condensation or evaporation pressure and the pore size underestimates the pore size. 

Beside classical methods of pore size analysis, there are many advanced methods. Seaton et al. [161] proposed a method based on the mean field theory. Initially this method was less accurate in the range of small pore sizes, but even so it g ve a more realistic -way for evaluation of the pore size distribution than the classical methods based on the Kelvin equation [162]. More rigorous methods based on molecular approaches such as grand canonical Monte Carlo (GCMC) simulations [147, 163-165] and nonlocal density functional theory (NLDFT) [86, 146, 147, 161, 163-169] have been developed and their use for pore size analysis of active carbons is continuously growing. 

The non-local density functional theory (NLDFT) is well established and widely presented in the literature. The distribution of density in a confined pore can be obtained for an open system in which a pore is allowed to exchange mass with the surroundings. From the thermodynamic principle, the density distribution is obtained by minimization of the following grand potential written below for the one-dimensional case [170] ... 

Finn and Monson [139] first tested the predictability of IAS theory for binary systems using the isothermal isobaric Monte Carlo simulation on a single surface. However, this system does not represent real adsorption systems. Tan and Gubbins [140,141] conducted detailed studies on the binary equilibria of the methane-ethane system in slit-shaped micropores using the nonlocal density function theory (NLDFT). The selectivity of ethane to methane was studied in terms of pore width, temperature, pressure, and molar fractions. 

HCP-l,4-benzenedimethanol (HCP-BDM) and HCP-ben l alcohol (HCP-BA) networks synthesized by the Friedel-Crafts self-condensation method have shown high selectivity for CO2 over N2, measured by nitrogen adsorption isotherm at two different temperatures, 273 and 298 Nonlocal density functional theory (NLDFT) calculations confirmed the pore sizes to be below 2 nm for both of the networks. Such small pore sizes, as well as the high ojq gen content in both the polymers, is most probably tbe reason for the strong interactions with polar CO2 rather than N2, making these good candidates for the selective separation of CO2 from N2. 

Note that several variants of the above functional (i.e., Eq. 41) exist. Instead of using WDA for treating as shown in Eq. (43), the EDA can also be appHed which provides comparable accuracy as discussed by Fu et al. (2015a). In addition, if we simply ignore this correlation term F ot and apply WCA scheme for the decomposition of repulsion and attraction in the LJ potential, the above combined functional recovers to the so-called nonlocal density functional theory (NLDFT) initiated by Balbuena and Gubbins (1993) and extensively appfred by Neimark et al. (Landers et al., 2013 Olivier et al., 1994). 

We suggest a model of adsorption in pores with amorphous and microporous solid walls, named the quenched solid non-local density functional theory (QSNLDFT) model. We consider a multicomponent non-local density functional theory (NLDFT), in which the solid is treated as a quenched component with a fixed spatially distributed density. Drawing on several prominent examples, we show that QSNLDFT model produces smooth Isotherms of mono- and polymolecular adsorption, which resemble experimental isotherms on amorphous surfaces. The model reproduces typical behaviors of N2 isotherms on micro- mesoporous materials, such as SBA-15. QSNLDFT model offers a systematic approach to the account for the surface roughness/heterogeneity in pore structure characterization methods. 

Recent progress in the theory of adsorption on porous solids, in general, and in the adsorption methods of pore structure characterization, in particular, has been related, to a large extent, to the application of the density functional theory (DFT) of Inhomogeneous fluids [1]. DFT has helped qualitatively describe and classify the specifics of adsorption and capillary condensation in pores of different geometries [2-4]. Moreover, it has been shown that the non-local density functional theory (NLDFT) with suitably chosen parameters of fluid-fluid and fluid-solid interactions quantitatively predicts the positions of capillary condensation and desorption transitions of argon and nitrogen in cylindrical pores of ordered mesoporous molecular sieves of MCM-41 and SBA-15 types [5,6]. NLDFT methods have been already commercialized by the producers of adsorption equipment for the interpretation of experimental data and the calculation of pore size distributions from adsorption isotherms [7-9]. 

In this paper, we suggest a systematic approach that extends the applicability of NLDFT models to heterogeneous surfaces of amorphous and microporous solids. The main idea is to use a multicomponent NLDFT, in which the solid is treated as one of the components with a fixed spatially distributed density. The model, named quenched solid non-local density functional theory (QSNLDFT), is an extension of the quenehed-annealed DFT model of systems with hard-core interactions recently proposed by Schmidt and coworkers [23,24]. Drawing on several prominent examples, we show that the proposed model produces smooth isotherms in the region of multiplayer adsorption. Moreover, the effects of wall microporosity can be naturally incorporated into the model. Although the parameters of the model have not been yet optimized to describe quantitatively a particular experimental system, the model generates adsorption isotherms which are in qualitative agreement with experimental isotherms of N2 or Ar adsorption on amorphous silica materials. 

The model (2)-(4) is referred to as the quenched solid non-local density functional theory (QSNLDFT). There are several advantages in considering the solid as a quenched component of the system rather than a source of the external field. On the one hand, this approach offers flexibility in the description of the fluid-solid boundary by varying the solid density and the thickness of the diffuse solid surface layer. On the other hand, it retains the main advantage of NLDFT computational efficiency because even a one-dimensional solid density distribution ceui include the effects of surface roughness and heterogeneity. For example, the solid density distribution can be taken from simulations of amorphous silica surfaces [29,30]. 

---