Chemistry and Dynamics of Supercritical Carbon Dioxide and Methane in the Slit Pores of Layered Silicates.

ConspectusIn the mid 2010s, high-pressure diffraction and spectroscopic tools opened a window into the molecular-scale behavior of fluids under the conditions of many CO2 sequestration and shale/tight gas reservoirs, conditions where CO2 and CH4 are present as variably wet supercritical fluids. Integrating high-pressure spectroscopy and diffraction with molecular modeling has revealed much about the ways that supercritical CO2 and CH4 behave in reservoir components, particularly in the slit-shaped micro- and mesopores of layered silicates (phyllosilicates) abundant in caprocks and shales. This Account summarizes how supercritical CO2 and CH4 behave in the slit pores of swelling phyllosilicates as functions of the H2O activity, framework structural features, and charge-balancing cation properties at 90 bar and 323 K, conditions similar to a reservoir at ∼1 km depth. Slit pores containing cations with large radii, low hydration energy, and large polarizability readily interact with CO2, allowing CO2 and H2O to adsorb and coexist in these interlayer pores over a wide range of fluid humidities. In contrast, cations with small radii, high hydration energy, and low polarizability weakly interact with CO2, leading to reduced CO2 uptake and a tendency to exclude CO2 from interlayers when H2O is abundant. The reorientation dynamics of confined CO2 depends on the interlayer pore height, which is strongly influenced by the cation properties, framework properties, and fluid humidity. The silicate structural framework also influences CO2 uptake and behavior; for example, smectites with increasing F-for-OH substitution in the framework take up greater quantities of CO2. Reactions that trap CO2 in carbonate phases have been observed in thin H2O films near smectite surfaces, including a dissolution-reprecipitation mechanism when the edge surface area is large and an ion exchange-precipitation mechanism when the interlayer cation can form a highly insoluble carbonate. In contrast, supercritical CH4 does not readily associate with cations, does not react with smectites, and is only incorporated into interlayer slit mesopores when (i) the pore has a z-dimension large enough to accommodate CH4, (ii) the smectite has low charge, and (iii) the H2O activity is low. The adsorption and displacement of CH4 by CO2 and vice versa have been studied on the molecular scale in one shale, but opportunities remain to examine behavioral details in this more complicated, slit-pore inclusive system.

Non-Equilibrium Modeling of Concentration-Driven processes with Constant Chemical Potential Molecular Dynamics Simulations.

ConspectusConcentration-driven processes in solution, i.e., phenomena that are sustained by persistent concentration gradients, such as crystallization and surface adsorption, are fundamental chemical processes. Understanding such phenomena is crucial for countless applications, from pharmaceuticals to biotechnology. Molecular dynamics (MD), both in- and out-of-equilibrium, plays an essential role in the current understanding of concentration-driven processes. Computational costs, however, impose drastic limitations on the accessible scale of simulated systems, hampering the effective study of such phenomena. In particular, due to these size limitations, closed system MD of concentration-driven processes is affected by solution depletion/enrichment that unavoidably impacts the dynamics of the chemical phenomena under study. As a notable example, in simulations of crystallization from solution, the transfer of monomers between the liquid and crystal phases results in a gradual depletion/enrichment of solution concentration, altering the driving force for phase transition. In contrast, this effect is negligible in experiments, given the macroscopic size of the solution volume. Because of these limitations, accurate MD characterization of concentration-driven phenomena has proven to be a long-standing simulation challenge. While disparate equilibrium and nonequilibrium simulation strategies have been proposed to address the study of such processes, the methodologies are in continuous development.In this context, a novel simulation technique named constant chemical potential molecular dynamics (CµMD) was recently proposed. CµMD employs properly designed, concentration-dependent external forces that regulate the flux of solute species between selected subregions of the simulation volume. This enables simulations of systems under a constant chemical drive in an efficient and straightforward way. The CµMD scheme was originally applied to the case of crystal growth from solution and then extended to the simulation of various physicochemical processes, resulting in new variants of the method. This Account illustrates the CµMD method and the key advances enabled by it in the framework of in silico chemistry. We review results obtained in crystallization studies, where CµMD allows growth rate calculations and equilibrium shape predictions, and in adsorption studies, where adsorption thermodynamics on porous or solid surfaces was correctly characterized via CµMD. Furthermore, we will discuss the application of CµMD variants to simulate permeation through porous materials, solution separation, and nucleation upon fixed concentration gradients. While presenting the numerous applications of the method, we provide an original and comprehensive assessment of concentration-driven simulations using CµMD. To this end, we also shed light on the theoretical and technical foundations of CµMD, underlining the novelty and specificity of the method with respect to existing techniques while stressing its current limitations. Overall, the application of CµMD to a diverse range of fields provides new insight into many physicochemical processes, the in silico study of which has been hitherto limited by finite-size effects. In this context, CµMD stands out as a general-purpose method that promises to be an invaluable simulation tool for studying molecular-scale concentration-driven phenomena.

A Nanotubular Metal-Organic Framework with a Narrow Bandgap from Extended Conjugation*.

A one-dimensional nanotubular metal-organic framework (MOF) [Ni(Cu-H4 TPPA)]⋅2 (CH3 )2 NH2 + (H8 TPPA=5,10,15,20-tetrakis[p-phenylphosphonic acid] porphyrin) constructed by using the arylphosphonic acid H8 TPPA is reported. The structure of this MOF, known as GTUB-4, was solved by using single-crystal X-ray diffraction and its geometric accessible surface area was calculated to be 1102 m2 g-1 , making it the phosphonate MOF with the highest reported surface area. Due to the extended conjugation of its porphyrin core, GTUB-4 possesses narrow indirect and direct bandgaps (1.9 eV and 2.16 eV, respectively) in the semiconductor regime. Thermogravimetric analysis suggests that GTUB-4 is thermally stable up to 400 °C. Owing to its high surface area, low bandgap, and high thermal stability, GTUB-4 could find applications as electrodes in supercapacitors.

Alkali Phosphonate Metal-Organic Frameworks.

A new family of porous metal-organic frameworks (MOFs), namely alkali phosphonate MOFs, is reported. [Na2 Cu(H4 TPPA)]⋅(NH2 (CH3 )2 )2 (GTUB-1) was synthesized using the tetratopic 5,10,15,20-tetrakis[p-phenylphosphonic acid] porphyrin (H8 -TPPA) linker with planar X-shaped geometrical core. GTUB-1 is composed of rectangular void channels with BET surface area of 697 m2 g-1 . GTUB-1 exhibits exceptional thermal stability. The toxicity analysis of the (H8 -TPPA) linker indicates that it is well tolerated by an intestinal cell line, suggesting its suitability for creating phosphonate MOFs for biological applications.

Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling.

The interactions among fluid species such as H2O, CO2, and CH4 confined in nano- and meso-pores in shales and other rocks is of central concern to understanding the chemical behavior and transport properties of these species in the earth's subsurface and is of special concern to geological C-sequestration and enhanced production of oil and natural gas. The behavior of CO2, and CH4 is less well understood than that of H2O. This paper presents the results of a computational modeling study of the partitioning of CO2 and CH4 between bulk fluid and nano- and meso-pores bounded by the common clay mineral montmorillonite. The calculations were done at 323 K and a total fluid pressure of 124 bars using a novel approach (constant reservoir composition molecular dynamics, CRC-MD) that uses bias forces to maintain a constant composition in the fluid external to the pore. This purely MD approach overcomes the difficulties in making stochastic particle insertion-deletion moves in dense fluids encountered in grand canonical Monte Carlo and related hybrid approaches. The results show that both the basal siloxane surfaces and protonated broken edge surfaces of montmorillonite both prefer CO2 relative to CH4 suggesting that methods of enhanced oil and gas production using CO2 will readily displace CH4 from such pores. This preference for CO2 is due to its preferred interaction with the surfaces and extends to approximately 20 Å from them.

Molecular dynamics modeling of carbon dioxide, water and natural organic matter in Na-hectorite.

Molecular dynamics (MD) modeling of systems containing a Na-exchanged smectite clay (hectorite) and model natural organic matter (NOM) molecules along with pure H2O, pure CO2, or a mixture of H2O and CO2 provides significant new insight into the molecular scale interactions among silicate surfaces, dissolved cations and organic molecules, H2O and CO2 relevant to geological C-sequestration strategies. The simulations for systems containing H2O show the following results; (1) Na(+) does not bridge between NOM molecules and the clay surface at protonation states comparable to near neutral pH conditions. (2) In systems without CO2 the NOM molecules retain charge balancing cations and drift away from the silicate surface. (3) In systems containing both H2O and CO2, the NOM molecules adopt equilibrium positions at the H2O-CO2 interface with the more hydrophilic structural elements facing the H2O and the more hydrophobic ones facing the CO2. In systems with only CO2, NOM and Na(+) ions are pinned to the clay surface with the hydrophilic structural elements of the NOM pointed toward the clay surface. Dynamically, in systems with only CO2, Na(+) diffusion is nearly eliminated, and in systems with a thin water film on the clay surface diffusion perpendicular the surface is greatly reduced relative to the system with bulk water. Energetically, the results for the systems with only H2O show that hydration of the net charge neutral Na-NOM molecule outweighs the sum of its Coulombic and dispersive interactions with the net charge-neutral Na-clay particle and the interactions of the water molecules with the hydrophobic structural elements of the NOM. The aggregation of NOM molecules in solution appears to be driven not by Na(+) bridging between the molecules but by hydrophobic interactions between them. In contrast, for the systems with only CO2 the interaction between the Na-NOM molecules and the CO2 is outweighed by the interaction of NOM with the clay particle. With both H2O and CO2 present, the energetic interactions leading to the hydration of the Na-clay surface and the hydrophilic structural elements of the Na-NOM molecule and the hydrophobic interactions between the CO2 and the hydrophobic aromatic and aliphatic structural elements of the NOM can both be satisfied, leading to the Na-NOM molecules migrating away from the surface and residing at the H2O-CO2 interface. The MD results suggest some alternative explanations for the previously observed (23)Na NMR behavior of Na-hectorite at elevated temperatures and CO2 pressures.

Synergistic interaction of paclitaxel and curcumin with cyclodextrin polymer complexation in human cancer cells.

The use of cytotoxic chemotherapic agents is the most common method for the treatment of metastatic cancers. Poor water solubility and low efficiency of chemotherapic agents are among the major hurdles of effective chemotherapy treatments. Curcumin and paclitaxel are well-known chemotherapic agents with poor water solubility and undesired side effects. In this study, a novel drug nanocarrier system was formulated by encapsulating curcumin and paclitaxel in poly(ß-cyclodextrin triazine) (PCDT) for the therapy of four cancer models; ovarian, lung, prostate, and breast cancer. Cell viability and colony formation assays revealed enhanced curcumin cytotoxicity upon complexation. Annexin V apoptotic studies showed that the PCDT complexation improved curcumin induced apoptosis in human ovarian cancer cell lines A2780 and SKOV-3, human nonsmall cell lung carcinoma cell line H1299, and human prostate cancer line DU-145, while no significant effect was observed with paclitaxel/PCDT complexation. The bioactivity of combining curcumin and paclitaxel was also investigated. A synergism was found between curcumin and paclitaxel, particularly when complexed with PCDT on A2780, SKOV-3, and H1299 cancer cell lines.

The effect of SO2 on CO2 capture in zeolitic imidazolate frameworks.

The CO2 capture performance of four zeolitic imidazolate frameworks (ZIFs), ZIF-10, ZIF-71, ZIF-68 and ZIF-69, was investigated in the presence of SO2 using molecular simulations. Our results show that SO2 is preferentially adsorbed over CO2. However, it also shows a cooperative effect on CO2 adsorption. In ZIF-10 and ZIF-71, the SO2 presence was found to have almost no effect on CO2 adsorption. In ZIF-68 and ZIF-69, although the two ZIFs adsorbed large amounts of SO2, the decrease in CO2 uptake was relatively small due to the cooperative interaction between the SO2 and CO2 molecules.

Metal-organic framework materials with ultrahigh surface areas: is the sky the limit?

We have synthesized, characterized, and computationally simulated/validated the behavior of two new metal-organic framework (MOF) materials displaying the highest experimental Brunauer-Emmett-Teller (BET) surface areas of any porous materials reported to date (~7000 m(2)/g). Key to evacuating the initially solvent-filled materials without pore collapse, and thereby accessing the ultrahigh areas, is the use of a supercritical CO(2) activation technique. Additionally, we demonstrate computationally that by shifting from phenyl groups to "space efficient" acetylene moieties as linker expansion units, the hypothetical maximum surface area for a MOF material is substantially greater than previously envisioned (~14600 m(2)/g (or greater) versus ~10500 m(2)/g).

Biohybrid Membrane Formation by Directed Insertion of Aquaporin into a Solid-State Nanopore.

Biohybrid nanopores combine the durability of solid-state nanopores with the precise structure and function of biological nanopores. Particular care must be taken to control how biological nanopores adapt to their surroundings once they come into contact with the solid-state nanopores. Two major challenges are to precisely control this adaptability under dynamic conditions and provide predesigned functionalities that can be manipulated for engineering applications. In this work, we report on the computational design of a distinctive class of biohybrid active membrane layers, built from the directed-insertion of an aquaporin-incorporated lipid nanodisc into a model alkyl-functionalized silica pore. We show that in an aqueous environment when a pressure difference exists between the two sides of the solid-state nanopore, the preferential interactions between the hydrocarbon tail of the lipid molecules that surround the aquaporin protein and the alkyl group functionalizing the interior surface of the silica nanopore enable the insertion of the aquaporin-incorporated lipid shell into the nanopore by forcing out the water molecules. The same preferential interactions are responsible for the structural stability of the inserted aquaporin-incorporated lipid shell as well as the water sealing properties of the lipid-alkyl interface. We further show that the aquaporin protein stabilized in the alkyl-functionalized silica nanopore preserves its biological structure and function in both pure and saline water, and, remarkably, its water permeability is equal to the one measured in the biological environment. The designed biohybrid membrane could pave the way for the development of durable transformative devices for water filtration.

Tuning the Transport Properties of Gases in Porous Graphene Membranes with Controlled Pore Size and Thickness.

Porous graphene membranes have emerged as promising alternatives for gas-separation applications due to their atomic thickness enabling ultrahigh permeance, but they suffer from low gas selectivity. Whereas decreasing the pore size below 3 nm is expected to increase the gas selectivity due to molecular sieving, it is rather challenging to generate a large number of uniform small pores on the graphene surface. Here, a pore-narrowing approach via gold deposition onto porous graphene surface is introduced to tune the pore size and thickness of the membrane to achieve a large number of small pores. Through the systematic approach, the ideal combination is determined as pore size below 3 nm, obtained at the thickness of 100 nm, to attain high selectivity and high permeance. The resulting membrane shows a H2 /CO2 separation factor of 31.3 at H2 permeance of 2.23 × 105 GPU (1 GPU = 3.35 × 10-10  mol s-1 m-2 Pa-1 ), which is the highest value reported to date in the 105 GPU permeance range. This result is explained by comparing the predicted binding energies of gas molecules with the Au surface, -5.3 versus -21 kJ mol-1 for H2 and CO2 , respectively, increased surface-gas interactions and molecular-sieving effect with decreasing pore size.

Evaluation of the BET method for determining surface areas of MOFs and zeolites that contain ultra-micropores.

The BET analysis is commonly used for determining surface areas of metal-organic frameworks (MOFs) and zeolites that contain "ultra-micropores" (<7 A) even though it is often stated that the BET surface areas obtained for such small pores are not really meaningful in an absolute sense. In this study, nitrogen and argon isotherms in MOFs and zeolites (most of them having ultra-micropores) were predicted by grand canonical Monte Carlo (GCMC) simulations and used as pseudoexperimental data to evaluate the BET method for these structures. The BET surface areas calculated from the simulated nitrogen and argon isotherms agree well with the accessible surface areas obtained directly from crystal structures in a geometric fashion. However, this was only true when the BET analysis was performed using the appropriate pressure range based on published "consistency criteria"; the BET analysis underestimates the surface areas of all the selected MOFs and zeolites if it is done in the "standard" BET pressure range. Moreover, the BET theory was shown to work well for a zeolite with a highly heterogeneous surface. These results validate the application of the BET method for determining surface areas of MOFs and zeolites even for materials having ultra-micropores or heterogeneous pores.

Does Confinement Enable Methane Hydrate Growth at Low Pressures? Insights from Molecular Dynamics Simulations.

Natural methane hydrates are estimated to be the largest source of unexploited hydrocarbon fuel. The ideal conditions for methane hydrate formation are low temperatures and high pressures. On the other hand, recent experimental studies suggest that porous materials, thanks to their confinement effects, can enable methane hydrate formation at milder conditions, although there has not been a consensus on this. A number of studies have investigated methane hydrate growth in confinement by employing molecular simulations; however, these were carried out at either very high pressures or very low temperatures. Therefore, the effects of confinement on methane hydrate growth at milder conditions have not yet been elaborated by molecular simulations. In order to address this, we carried out a systematic study by performing molecular dynamics (MD) simulations of methane water systems. Using a direct phase coexistence approach, microsecond-scale MD simulations in the isobaric-isothermal (NPT) ensemble were performed in order to study the behavior of methane hydrates in the bulk and in confined nanospaces of hydroxylated silica pores at external pressures ranging from 1 to 100 bar and a simulation temperature corresponding to a 2 °C experimental temperature. We validated the combination of the TIP4P/ice water and TraPPE-UA methane models in order to correctly predict the behavior of methane hydrates in accordance to their phase equilibria. We also demonstrated that the dispersion corrections applied to short-range interactions lead to artificially induced hydrate growth. We observed that in the confinement of a hydroxylated silica pore, a convex-shaped methane nanobuble forms, and methane hydrate growth primarily takes place in the center of the pore rather than the surfaces where a thin water layer exists. Most importantly, our study showed that in the nanopores methane hydrate growth can indeed take place at pressures which would be too low for the growth of methane hydrates in the bulk.

Modeling of Gas Transport through Polymer/MOF Interfaces: A Microsecond-Scale Concentration Gradient-Driven Molecular Dynamics Study.

Membrane-based separation technologies offer a cost-effective alternative to many energy-intensive gas separation processes, such as distillation. Mixed matrix membranes (MMMs) composed of polymers and metal-organic frameworks (MOFs) have attracted a great deal of attention for being promising systems to manufacture durable and highly selective membranes with high gas fluxes and high selectivities. Therefore, understanding gas transport through these MMMs is of significant importance. There has been longstanding speculation that the gas diffusion behavior at the interface formed between the polymer matrix and MOF particles would strongly affect the global performance of the MMMs due to the potential presence of nonselective voids or other defects. To shed more light on this paradigm, we have performed microsecond long concentration gradient-driven molecular dynamics (CGD-MD) simulations that deliver an unprecedented microscopic picture of the transport of H2 and CH4 as single components and as a mixture in all regions of the PIM-1/ZIF-8 membrane, including the polymer/MOF interface. The fluxes of the permeating gases are computed and the impact of the polymer/MOF interface on the H2/CH4 permselectivity of the composite membrane is clearly revealed. Specifically, we show that the poor compatibility between PIM-1 and ZIF-8, which manifests itself by the presence of nonselective void spaces at their interface, results in a decrease of the H2/CH4 permselectivity for the corresponding composite membrane as compared to the performances simulated for PIM-1 and ZIF-8 individually. We demonstrate that CGD-MD simulations based on an accurate atomistic description of the polymer/MOF composite is a powerful tool for characterization and understanding of gas transport and separation mechanisms in MMMs.

Semiconductive microporous hydrogen-bonded organophosphonic acid frameworks.

Herein, we report a semiconductive, proton-conductive, microporous hydrogen-bonded organic framework (HOF) derived from phenylphosphonic acid and 5,10,15,20-tetrakis[p-phenylphosphonic acid] porphyrin (GTUB5). The structure of GTUB5 was characterized using single crystal X-ray diffraction. A narrow band gap of 1.56 eV was extracted from a UV-Vis spectrum of pure GTUB5 crystals, in excellent agreement with the 1.65 eV band gap obtained from DFT calculations. The same band gap was also measured for GTUB5 in DMSO. The proton conductivity of GTUB5 was measured to be 3.00 × 10-6 S cm-1 at 75 °C and 75% relative humidity. The surface area was estimated to be 422 m2 g-1 from grand canonical Monte Carlo simulations. XRD showed that GTUB5 is thermally stable under relative humidities of up to 90% at 90 °C. These findings pave the way for a new family of organic, microporous, and semiconducting materials with high surface areas and high thermal stabilities.

Rapid and Efficient Removal of Perfluorooctanoic Acid from Water with Fluorine-Rich Calixarene-Based Porous Polymers.

On account of its nonbiodegradable nature and persistence in the environment, perfluorooctanoic acid (PFOA) accumulates in water resources and poses serious environmental issues in many parts of the world. Here, we present the development of two fluorine-rich calix[4]arene-based porous polymers, FCX4-P and FCX4-BP, and demonstrate their utility for the efficient removal of PFOA from water. These materials featured Brunauer-Emmett-Teller (BET) surface areas of up to 450 m2 g-1, which is slightly lower than their nonfluorinated counterparts (up to 596 m2 g-1). FCX4-P removes PFOA at environmentally relevant concentrations with a high rate constant of 3.80 g mg-1 h-1 and reached an exceptional maximum PFOA uptake capacity of 188.7 mg g-1. In addition, it could be regenerated by simple methanol wash and reused without a significant decrease in performance.

Screening of metal-organic frameworks for carbon dioxide capture from flue gas using a combined experimental and modeling approach.

A diverse collection of 14 metal-organic frameworks (MOFs) was screened for CO(2) capture from flue gas using a combined experimental and modeling approach. Adsorption measurements are reported for the screened MOFs at room temperature up to 1 bar. These data are used to validate a generalized strategy for molecular modeling of CO(2) and other small molecules in MOFs. MOFs possessing a high density of open metal sites are found to adsorb significant amounts of CO(2) even at low pressure. An excellent correlation is found between the heat of adsorption and the amount of CO(2) adsorbed below 1 bar. Molecular modeling can aid in selection of adsorbents for CO(2) capture from flue gas by screening a large number of MOFs.

Correction: A cobalt arylphosphonate MOF - superior stability, sorption and magnetism.

Correction for 'A cobalt arylphosphonate MOF - superior stability, sorption and magnetism' by Yunus Zorlu et al., Chem. Commun., 2019, DOI: 10.1039/c8cc09655d.

A cobalt arylphosphonate MOF - superior stability, sorption and magnetism.

We report a novel metal-organic framework (MOF) based on a cobalt arylphosphonate, namely, [Co2(H4-MTPPA)]·3NMP·H2O (1·3NMP·H2O), which was prepared solvothermically from the tetrahedral linker tetraphenylmethane tetrakis-4-phosphonic acid (H8-MTPPA) and CoSO4·7H2O in N-methyl-2-pyrrolidone (NMP). Compound 1 has the highest porosity (BET surface area of 1034 m2 g-1) ever reported for a MOF based on an aryl phosphonic acid linker. The indigo blue crystals of 1·3NMP·H2O are composed of edge-shared eight-membered Co2P2O4 rings, and are thermally very stable up to 500 °C.

Tipping Point for Expansion of Layered Aluminosilicates in Weakly Polar Solvents: Supercritical CO2.

Layered aluminosilicates play a dominant role in the mechanical and gas storage properties of the subsurface, are used in diverse industrial applications, and serve as model materials for understanding solvent-ion-support systems. Although expansion in the presence of H2O is well-known to be systematically correlated with the hydration free energy of the interlayer cation, particularly in environments dominated by nonpolar solvents (i.e., CO2), uptake into the interlayer is not well-understood. Using novel high-pressure capabilities, we investigated the interaction of dry supercritical CO2 with Na-, NH4-, and Cs-saturated montmorillonite, comparing results with predictions from molecular dynamics simulations. Despite the known trend in H2O and that cation solvation energies in CO2 suggest a stronger interaction with Na, both the NH4- and Cs-clays readily absorbed CO2 and expanded, while the Na-clay did not. The apparent inertness of the Na-clay was not due to kinetics, as experiments seeking a stable expanded state showed that none exists. Molecular dynamics simulations revealed a large endothermicity to CO2 intercalation in the Na-clay but little or no energy barrier for the NH4- and Cs-clays. Indeed, the combination of experiment and theory clearly demonstrate that CO2 intercalation of Na-montmorillonite clays is prohibited in the absence of H2O. Consequently, we have shown for the first time that in the presence of a low dielectric constant, gas swelling depends more on the strength of the interaction between the interlayer cation and aluminosilicate sheets and less on that with solvent. The finding suggests a distinct regime in layered aluminosilicate swelling behavior triggered by low solvent polarizability, with important implications in geomechanics, storage, and retention of volatile gases, and across industrial uses in gelling, decoloring, heterogeneous catalysis, and semipermeable reactive barriers.