Hubbry Logo
AdsorptionAdsorptionMain
Open search
Adsorption
Community hub
Adsorption
logo
8 pages, 0 posts
0 subscribers
Be the first to start a discussion here.
Be the first to start a discussion here.
Adsorption
Adsorption
from Wikipedia
Brunauer, Emmett and Teller's model of multilayer adsorption is a random distribution of molecules on the material surface.

Adsorption is the adhesion[1] of atoms, ions or molecules from a gas, liquid or dissolved solid to a surface.[2] This process creates a film of the adsorbate on the surface of the adsorbent. This process differs from absorption, in which a fluid (the absorbate) is dissolved by or permeates a liquid or solid (the absorbent).[3] While adsorption does often precede absorption, which involves the transfer of the absorbate into the volume of the absorbent material, alternatively, adsorption is distinctly a surface phenomenon, wherein the adsorbate does not penetrate through the material surface and into the bulk of the adsorbent.[4] The term sorption encompasses both adsorption and absorption, and desorption is the reverse of sorption.

IUPAC definition

adsorption: An increase in the concentration of a dissolved substance at the interface of a condensed and a liquid phase due to the operation of surface forces. Adsorption can also occur at the interface of a condensed and a gaseous phase.[5]

Like surface tension, adsorption is a consequence of surface energy. In a bulk material, all the bonding requirements (be they ionic, covalent or metallic) of the constituent atoms of the material are fulfilled by other atoms in the material. However, atoms on the surface of the adsorbent are not wholly surrounded by other adsorbent atoms and therefore can attract adsorbates. The exact nature of the bonding depends on the details of the species involved, but the adsorption process is generally classified as physisorption (characteristic of weak van der Waals forces) or chemisorption (characteristic of covalent bonding). It may also occur due to electrostatic attraction.[6][7] The nature of the adsorption can affect the structure of the adsorbed species. For example, polymer physisorption from solution can result in squashed structures on a surface.[8]

Adsorption is present in many natural, physical, biological and chemical systems and is widely used in industrial applications such as heterogeneous catalysts,[9][10] activated charcoal, capturing and using waste heat to provide cold water for air conditioning and other process requirements (adsorption chillers), synthetic resins, increasing storage capacity of carbide-derived carbons and water purification.[11] Adsorption, ion exchange and chromatography are sorption processes in which certain adsorbates are selectively transferred from the fluid phase to the surface of insoluble, rigid particles suspended in a vessel or packed in a column. Pharmaceutical industry applications, which use adsorption as a means to prolong neurological exposure to specific drugs or parts thereof,[citation needed] are lesser known.

The word "adsorption" was coined in 1881 by German physicist Heinrich Kayser (1853–1940).[12]

Isotherms

[edit]

The adsorption of gases and solutes is usually described through isotherms, that is, the amount of adsorbate on the adsorbent as a function of its pressure (if gas) or concentration (for liquid phase solutes) at constant temperature. The quantity adsorbed is nearly always normalized by the mass of the adsorbent to allow comparison of different materials. To date, 15 different isotherm models have been developed.[13]

Freundlich

[edit]

The first mathematical fit to an isotherm was published by Freundlich and Kuster (1906) and is a purely empirical formula for gaseous adsorbates:

where is the mass of adsorbate adsorbed, is the mass of the adsorbent, is the pressure of adsorbate (this can be changed to concentration if investigating solution rather than gas), and and are empirical constants for each adsorbent–adsorbate pair at a given temperature. The function is not adequate at very high pressure because in reality has an asymptotic maximum as pressure increases without bound. As the temperature increases, the constants and change to reflect the empirical observation that the quantity adsorbed rises more slowly and higher pressures are required to saturate the surface.

Langmuir

[edit]

Irving Langmuir was the first to derive a scientifically based adsorption isotherm in 1918.[14] The model applies to gases adsorbed on solid surfaces. It is a semi-empirical isotherm with a kinetic basis and was derived based on statistical thermodynamics. It is the most common isotherm equation to use due to its simplicity and its ability to fit a variety of adsorption data. It is based on four assumptions:

  1. All of the adsorption sites are equivalent, and each site can only accommodate one molecule.
  2. The surface is energetically homogeneous, and adsorbed molecules do not interact.
  3. There are no phase transitions.
  4. At the maximum adsorption, only a monolayer is formed. Adsorption only occurs on localized sites on the surface, not with other adsorbates.

These four assumptions are seldom all true: there are always imperfections on the surface, adsorbed molecules are not necessarily inert, and the mechanism is clearly not the same for the first molecules to adsorb to a surface as for the last. The fourth condition is the most troublesome, as frequently more molecules will adsorb to the monolayer; this problem is addressed by the BET isotherm for relatively flat (non-microporous) surfaces. The Langmuir isotherm is nonetheless the first choice for most models of adsorption and has many applications in surface kinetics (usually called Langmuir–Hinshelwood kinetics) and thermodynamics.

Langmuir suggested that adsorption takes place through this mechanism: , where A is a gas molecule, and S is an adsorption site. The direct and inverse rate constants are k and k−1. If we define surface coverage, , as the fraction of the adsorption sites occupied, in the equilibrium we have:

or

where is the partial pressure of the gas or the molar concentration of the solution. For very low pressures , and for high pressures .

The value of is difficult to measure experimentally; usually, the adsorbate is a gas and the quantity adsorbed is given in moles, grams, or gas volumes at standard temperature and pressure (STP) per gram of adsorbent. If we call vmon the STP volume of adsorbate required to form a monolayer on the adsorbent (per gram of adsorbent), then , and we obtain an expression for a straight line:

Through its slope and y intercept we can obtain vmon and K, which are constants for each adsorbent–adsorbate pair at a given temperature. vmon is related to the number of adsorption sites through the ideal gas law. If we assume that the number of sites is just the whole area of the solid divided into the cross section of the adsorbate molecules, we can easily calculate the surface area of the adsorbent. The surface area of an adsorbent depends on its structure: the more pores it has, the greater the area, which has a big influence on reactions on surfaces.

If more than one gas adsorbs on the surface, we define as the fraction of empty sites, and we have:

Also, we can define as the fraction of the sites occupied by the j-th gas:

where i is each one of the gases that adsorb.

Note:

1) To choose between the Langmuir and Freundlich equations, the enthalpies of adsorption must be investigated.[15] While the Langmuir model assumes that the energy of adsorption remains constant with surface occupancy, the Freundlich equation is derived with the assumption that the heat of adsorption continually decrease as the binding sites are occupied.[16] The choice of the model based on best fitting of the data is a common misconception.[15]

2) The use of the linearized form of the Langmuir model is no longer common practice. Advances in computational power allowed for nonlinear regression to be performed quickly and with higher confidence since no data transformation is required.

BET

[edit]

Often molecules do form multilayers, that is, some are adsorbed on already adsorbed molecules, and the Langmuir isotherm is not valid. In 1938 Stephen Brunauer, Paul Emmett, and Edward Teller developed a model isotherm that takes that possibility into account. Their theory is called BET theory, after the initials in their last names. They modified Langmuir's mechanism as follows:

A(g) + S ⇌ AS,
A(g) + AS ⇌ A2S,
A(g) + A2S ⇌ A3S and so on.
Langmuir (blue) and BET (red) isotherms

The derivation of the formula is more complicated than Langmuir's (see links for complete derivation). We obtain:

where x is the pressure divided by the vapor pressure for the adsorbate at that temperature (usually denoted ), v is the STP volume of adsorbed adsorbate, vmon is the STP volume of the amount of adsorbate required to form a monolayer, and c is the equilibrium constant K we used in Langmuir isotherm multiplied by the vapor pressure of the adsorbate. The key assumption used in deriving the BET equation that the successive heats of adsorption for all layers except the first are equal to the heat of condensation of the adsorbate.

The Langmuir isotherm is usually better for chemisorption, and the BET isotherm works better for physisorption for non-microporous surfaces.

Kisliuk

[edit]
Two adsorbate nitrogen molecules adsorbing onto a tungsten adsorbent from the precursor state around an island of previously adsorbed adsorbate (left) and via random adsorption (right)

In other instances, molecular interactions between gas molecules previously adsorbed on a solid surface form significant interactions with gas molecules in the gaseous phases. Hence, adsorption of gas molecules to the surface is more likely to occur around gas molecules that are already present on the solid surface, rendering the Langmuir adsorption isotherm ineffective for the purposes of modelling. This effect was studied in a system where nitrogen was the adsorbate and tungsten was the adsorbent by Paul Kisliuk (1922–2008) in 1957.[17] To compensate for the increased probability of adsorption occurring around molecules present on the substrate surface, Kisliuk developed the precursor state theory, whereby molecules would enter a precursor state at the interface between the solid adsorbent and adsorbate in the gaseous phase. From here, adsorbate molecules would either adsorb to the adsorbent or desorb into the gaseous phase. The probability of adsorption occurring from the precursor state is dependent on the adsorbate's proximity to other adsorbate molecules that have already been adsorbed. If the adsorbate molecule in the precursor state is in close proximity to an adsorbate molecule that has already formed on the surface, it has a sticking probability reflected by the size of the SE constant and will either be adsorbed from the precursor state at a rate of kEC or will desorb into the gaseous phase at a rate of kES. If an adsorbate molecule enters the precursor state at a location that is remote from any other previously adsorbed adsorbate molecules, the sticking probability is reflected by the size of the SD constant.

These factors were included as part of a single constant termed a "sticking coefficient", kE, described below:

As SD is dictated by factors that are taken into account by the Langmuir model, SD can be assumed to be the adsorption rate constant. However, the rate constant for the Kisliuk model (R') is different from that of the Langmuir model, as R' is used to represent the impact of diffusion on monolayer formation and is proportional to the square root of the system's diffusion coefficient. The Kisliuk adsorption isotherm is written as follows, where θ(t) is fractional coverage of the adsorbent with adsorbate, and t is immersion time:

Solving for θ(t) yields:

Adsorption enthalpy

[edit]

Adsorption constants are equilibrium constants, therefore they obey the Van 't Hoff equation:

As can be seen in the formula, the variation of K must be isosteric, that is, at constant coverage. If we start from the BET isotherm and assume that the entropy change is the same for liquefaction and adsorption, we obtain

that is to say, adsorption is more exothermic than liquefaction.

Single-molecule explanation

[edit]

The adsorption of ensemble molecules on a surface or interface can be divided into two processes: adsorption and desorption. If the adsorption rate wins the desorption rate, the molecules will accumulate over time giving the adsorption curve over time. If the desorption rate is larger, the number of molecules on the surface will decrease over time. The adsorption rate is dependent on the temperature, the diffusion rate of the solute (related to mean free path for pure gas), and the energy barrier between the molecule and the surface. The diffusion and key elements of the adsorption rate can be calculated using Fick's laws of diffusion and Einstein relation (kinetic theory). Under ideal conditions, when there is no energy barrier and all molecules that diffuse and collide with the surface get adsorbed, the number of molecules adsorbed at a surface of area on an infinite area surface can be directly integrated from Fick's second law differential equation to be:[18]

where is the surface area (unit m2), is the number concentration of the molecule in the bulk solution (unit #/m3), is the diffusion constant (unit m2/s), and is time (unit s). Further simulations and analysis of this equation[19] show that the square root dependence on the time is originated from the decrease of the concentrations near the surface under ideal adsorption conditions. Also, this equation only works for the beginning of the adsorption when a well-behaved concentration gradient forms near the surface. Correction on the reduction of the adsorption area and slowing down of the concentration gradient evolution have to be considered over a longer time.[20] Under real experimental conditions, the flow and the small adsorption area always make the adsorption rate faster than what this equation predicted, and the energy barrier will either accelerate this rate by surface attraction or slow it down by surface repulsion. Thus, the prediction from this equation is often a few to several orders of magnitude away from the experimental results. Under special cases, such as a very small adsorption area on a large surface, and under chemical equilibrium when there is no concentration gradience near the surface, this equation becomes useful to predict the adsorption rate with debatable special care to determine a specific value of in a particular measurement.[19]

The desorption of a molecule from the surface depends on the binding energy of the molecule to the surface and the temperature. The typical overall adsorption rate is thus often a combined result of the adsorption and desorption.

Quantum mechanical – thermodynamic modelling for surface area and porosity

[edit]

Since 1980 two theories were worked on to explain adsorption and obtain equations that work. These two are referred to as the chi hypothesis, the quantum mechanical derivation, and excess surface work (ESW).[21] Both these theories yield the same equation for flat surfaces:

where U is the unit step function. The definitions of the other symbols is as follows:

where "ads" stands for "adsorbed", "m" stands for "monolayer equivalence" and "vap" is reference to the vapor pressure of the liquid adsorptive at the same temperature as the solid sample. The unit function creates the definition of the molar energy of adsorption for the first adsorbed molecule by:

The plot of adsorbed versus is referred to as the chi plot. For flat surfaces, the slope of the chi plot yields the surface area. Empirically, this plot was noticed as being a very good fit to the isotherm by Michael Polanyi[22][23][24] and also by Jan Hendrik de Boer and Cornelis Zwikker[25] but not pursued. This was due to criticism in the former case by Albert Einstein and in the latter case by Brunauer. This flat surface equation may be used as a "standard curve" in the normal tradition of comparison curves, with the exception that the porous sample's early portion of the plot of versus acts as a self-standard. Ultramicroporous, microporous and mesoporous conditions may be analyzed using this technique. Typical standard deviations for full isotherm fits including porous samples are less than 2%.

Notice that in this description of physical adsorption, the entropy of adsorption is consistent with the Dubinin thermodynamic criterion, that is the entropy of adsorption from the liquid state to the adsorbed state is approximately zero.

Adsorbents

[edit]

Characteristics and general requirements

[edit]
Activated carbon is used as an adsorbent

Adsorbents are used usually in the form of spherical pellets, rods, moldings, or monoliths with a hydrodynamic radius between 0.25 and 5 mm. They must have high abrasion resistance, high thermal stability and small pore diameters, which results in higher exposed surface area and hence high capacity for adsorption. The adsorbents must also have a distinct pore structure that enables fast transport of the gaseous vapors.[26] Most industrial adsorbents fall into one of three classes:

  • Oxygen-containing compounds – Are typically hydrophilic and polar, including materials such as silica gel, limestone (calcium carbonate)[27] and zeolites.
  • Carbon-based compounds – Are typically hydrophobic and non-polar, including materials such as activated carbon and graphite.
  • Polymer-based compounds – Are polar or non-polar, depending on the functional groups in the polymer matrix.

Silica gel

[edit]
Silica gel adsorber for NO2, Fixed Nitrogen Research Laboratory, ca.1930s

Silica gel is a chemically inert, non-toxic, polar and dimensionally stable (< 400 °C or 750 °F) amorphous form of SiO2. It is prepared by the reaction between sodium silicate and acetic acid, which is followed by a series of after-treatment processes such as aging, pickling, etc. These after-treatment methods results in various pore size distributions.

Silica is used for drying of process air (e.g. oxygen, natural gas) and adsorption of heavy (polar) hydrocarbons from natural gas.

Zeolites

[edit]

Zeolites are natural or synthetic crystalline aluminosilicates, which have a repeating pore network and release water at high temperature. Zeolites are polar in nature.

They are manufactured by hydrothermal synthesis of sodium aluminosilicate or another silica source in an autoclave followed by ion exchange with certain cations (Na+, Li+, Ca2+, K+, NH4+). The channel diameter of zeolite cages usually ranges from 2 to 9 Å. The ion exchange process is followed by drying of the crystals, which can be pelletized with a binder to form macroporous pellets.

Zeolites are applied in drying of process air, CO2 removal from natural gas, CO removal from reforming gas, air separation, catalytic cracking, and catalytic synthesis and reforming.

Non-polar (siliceous) zeolites are synthesized from aluminum-free silica sources or by dealumination of aluminum-containing zeolites. The dealumination process is done by treating the zeolite with steam at elevated temperatures, typically greater than 500 °C (930 °F). This high temperature heat treatment breaks the aluminum-oxygen bonds and the aluminum atom is expelled from the zeolite framework.

Activated carbon

[edit]

The term "adsorption" itself was coined by Heinrich Kayser in 1881 in the context of uptake of gases by carbons.[28]

Activated carbon is a highly porous, amorphous solid consisting of microcrystallites with a graphite lattice, usually prepared in small pellets or a powder. It is non-polar and cheap. One of its main drawbacks is that it reacts with oxygen at moderate temperatures (over 300 °C).

Activated carbon nitrogen isotherm showing a marked microporous type I behavior

Activated carbon can be manufactured from carbonaceous material, including coal (bituminous, subbituminous, and lignite), peat, wood, or nutshells (e.g., coconut). The manufacturing process consists of two phases, carbonization and activation.[29][30] The carbonization process includes drying and then heating to separate by-products, including tars and other hydrocarbons from the raw material, as well as to drive off any gases generated. The process is completed by heating the material over 400 °C (750 °F) in an oxygen-free atmosphere that cannot support combustion. The carbonized particles are then "activated" by exposing them to an oxidizing agent, usually steam or carbon dioxide at high temperature. This agent burns off the pore blocking structures created during the carbonization phase and so, they develop a porous, three-dimensional graphite lattice structure. The size of the pores developed during activation is a function of the time that they spend in this stage. Longer exposure times result in larger pore sizes. The most popular aqueous phase carbons are bituminous based because of their hardness, abrasion resistance, pore size distribution, and low cost, but their effectiveness needs to be tested in each application to determine the optimal product.

Activated carbon is used for adsorption of organic substances[31] and non-polar adsorbates and it is also usually used for waste gas (and waste water) treatment. It is the most widely used adsorbent since most of its chemical (e.g. surface groups) and physical properties (e.g. pore size distribution and surface area) can be tuned according to what is needed.[32] Its usefulness also derives from its large micropore (and sometimes mesopore) volume and the resulting high surface area. Recent research works reported activated carbon as an effective agent to adsorb cationic species of toxic metals from multi-pollutant systems and also proposed possible adsorption mechanisms with supporting evidences.[33]

Water adsorption

[edit]

The adsorption of water at surfaces is of broad importance in chemical engineering, materials science and catalysis. Also termed surface hydration, the presence of physically or chemically adsorbed water at the surfaces of solids plays an important role in governing interface properties, chemical reaction pathways and catalytic performance in a wide range of systems. In the case of physically adsorbed water, surface hydration can be eliminated simply through drying at conditions of temperature and pressure allowing full vaporization of water. For chemically adsorbed water, hydration may be in the form of either dissociative adsorption, where H2O molecules are dissociated into surface adsorbed -H and -OH, or molecular adsorption (associative adsorption) where individual water molecules remain intact [34]

Adsorption solar heating and storage

[edit]

The low cost ($200/ton) and high cycle rate (2,000 ×) of synthetic zeolites such as Linde 13X with water adsorbate has garnered much academic and commercial interest recently for use for thermal energy storage (TES), specifically of low-grade solar and waste heat. Several pilot projects have been funded in the EU from 2000 to the present (2020).[citation needed] The basic concept is to store solar thermal energy as chemical latent energy in the zeolite. Typically, hot dry air from flat plate solar collectors is made to flow through a bed of zeolite such that any water adsorbate present is driven off. Storage can be diurnal, weekly, monthly, or even seasonal depending on the volume of the zeolite and the area of the solar thermal panels. When heat is called for during the night, or sunless hours, or winter, humidified air flows through the zeolite. As the humidity is adsorbed by the zeolite, heat is released to the air and subsequently to the building space. This form of TES, with specific use of zeolites, was first taught by John Guerra in 1978.[35]

Carbon capture and storage

[edit]

Typical adsorbents proposed for carbon capture and storage are zeolites and MOFs.[36] The customization of adsorbents makes them a potentially attractive alternative to absorption. Because adsorbents can be regenerated by temperature or pressure swing, this step can be less energy intensive than absorption regeneration methods.[37] Major problems that are present with adsorption cost in carbon capture are: regenerating the adsorbent, mass ratio, solvent/MOF, cost of adsorbent, production of the adsorbent, lifetime of adsorbent.[38]

In sorption enhanced water gas shift (SEWGS) technology a pre-combustion carbon capture process, based on solid adsorption, is combined with the water gas shift reaction (WGS) in order to produce a high pressure hydrogen stream.[39] The CO2 stream produced can be stored or used for other industrial processes.[40]

Protein and surfactant adsorption

[edit]

Protein adsorption is a process that has a fundamental role in the field of biomaterials. Indeed, biomaterial surfaces in contact with biological media, such as blood or serum, are immediately coated by proteins. Therefore, living cells do not interact directly with the biomaterial surface, but with the adsorbed proteins layer. This protein layer mediates the interaction between biomaterials and cells, translating biomaterial physical and chemical properties into a "biological language".[41] In fact, cell membrane receptors bind to protein layer bioactive sites and these receptor-protein binding events are transduced, through the cell membrane, in a manner that stimulates specific intracellular processes that then determine cell adhesion, shape, growth and differentiation. Protein adsorption is influenced by many surface properties such as surface wettability, surface chemical composition [42] and surface nanometre-scale morphology.[43] Surfactant adsorption is a similar phenomenon, but utilising surfactant molecules in the place of proteins.[44]

Adsorption chillers

[edit]
A schematic diagram of an adsorption chiller: (1) heat is lost through evaporation of refrigerant, (2) refrigerant vapour is adsorbed onto the solid medium, (3) refrigerant is desorbed from the solid medium section not in use, (4) refrigerant is condensed and returned to the start, (5) & (6) solid medium is cycled between adsorption and desorption to regenerate it.

Combining an adsorbent with a refrigerant, adsorption chillers use heat to provide a cooling effect. This heat, in the form of hot water, may come from any number of industrial sources including waste heat from industrial processes, prime heat from solar thermal installations or from the exhaust or water jacket heat of a piston engine or turbine.

Although there are similarities between adsorption chillers and absorption refrigeration, the former is based on the interaction between gases and solids. The adsorption chamber of the chiller is filled with a solid material (for example zeolite, silica gel, alumina, active carbon or certain types of metal salts), which in its neutral state has adsorbed the refrigerant. When heated, the solid desorbs (releases) refrigerant vapour, which subsequently is cooled and liquefied. This liquid refrigerant then provides a cooling effect at the evaporator from its enthalpy of vaporization. In the final stage the refrigerant vapour is (re)adsorbed into the solid.[45] As an adsorption chiller requires no compressor, it is relatively quiet.

Portal site mediated adsorption

[edit]

Portal site mediated adsorption is a model for site-selective activated gas adsorption in metallic catalytic systems that contain a variety of different adsorption sites. In such systems, low-coordination "edge and corner" defect-like sites can exhibit significantly lower adsorption enthalpies than high-coordination (basal plane) sites. As a result, these sites can serve as "portals" for very rapid adsorption to the rest of the surface. The phenomenon relies on the common "spillover" effect (described below), where certain adsorbed species exhibit high mobility on some surfaces. The model explains seemingly inconsistent observations of gas adsorption thermodynamics and kinetics in catalytic systems where surfaces can exist in a range of coordination structures, and it has been successfully applied to bimetallic catalytic systems where synergistic activity is observed.

In contrast to pure spillover, portal site adsorption refers to surface diffusion to adjacent adsorption sites, not to non-adsorptive support surfaces.

The model appears to have been first proposed for carbon monoxide on silica-supported platinum by Brandt et al. (1993).[46] A similar, but independent model was developed by King and co-workers[47][48][49] to describe hydrogen adsorption on silica-supported alkali promoted ruthenium, silver-ruthenium and copper-ruthenium bimetallic catalysts. The same group applied the model to CO hydrogenation (Fischer–Tropsch synthesis).[50] Zupanc et al. (2002) subsequently confirmed the same model for hydrogen adsorption on magnesia-supported caesium-ruthenium bimetallic catalysts.[51] Trens et al. (2009) have similarly described CO surface diffusion on carbon-supported Pt particles of varying morphology.[52]

Adsorption spillover

[edit]

In the case catalytic or adsorbent systems where a metal species is dispersed upon a support (or carrier) material (often quasi-inert oxides, such as alumina or silica), it is possible for an adsorptive species to indirectly adsorb to the support surface under conditions where such adsorption is thermodynamically unfavorable. The presence of the metal serves as a lower-energy pathway for gaseous species to first adsorb to the metal and then diffuse on the support surface. This is possible because the adsorbed species attains a lower energy state once it has adsorbed to the metal, thus lowering the activation barrier between the gas phase species and the support-adsorbed species.

Hydrogen spillover is the most common example of an adsorptive spillover. In the case of hydrogen, adsorption is most often accompanied with dissociation of molecular hydrogen (H2) to atomic hydrogen (H), followed by spillover of the hydrogen atoms present.

The spillover effect has been used to explain many observations in heterogeneous catalysis and adsorption.[53]

Polymer adsorption

[edit]

Adsorption of molecules onto polymer surfaces is central to a number of applications, including development of non-stick coatings and in various biomedical devices. Polymers may also be adsorbed to surfaces through polyelectrolyte adsorption.

In viruses

[edit]

Adsorption is the first step in the viral life cycle. The next steps are penetration, uncoating, synthesis (transcription if needed, and translation), and release. The virus replication cycle, in this respect, is similar for all types of viruses. Factors such as transcription may or may not be needed if the virus is able to integrate its genomic information in the cell's nucleus, or if the virus can replicate itself directly within the cell's cytoplasm.

[edit]

The game of Tetris is a puzzle game in which blocks of 4 are adsorbed onto a surface during game play. Scientists have used Tetris blocks "as a proxy for molecules with a complex shape" and their "adsorption on a flat surface" for studying the thermodynamics of nanoparticles.[54][55]

See also

[edit]

References

[edit]

Further reading

[edit]
[edit]
Revisions and contributorsEdit on WikipediaRead on Wikipedia
from Grokipedia
Adsorption is a surface characterized by the accumulation of atoms, ions, or molecules (the adsorbate) from a gas, , or dissolved solid phase onto the surface of a solid or (the adsorbent), resulting in an increased concentration at the interface due to surface forces. This process differs fundamentally from absorption, where the adsorbate penetrates into the bulk of the adsorbent, as adsorption is confined to the surface layers. Adsorption can be classified into two primary types based on the nature of the interaction forces: and . , or physical adsorption, involves weak intermolecular forces such as van der Waals interactions, leading to reversible binding that is typically multilayered and occurs at lower temperatures. In contrast, , or chemical adsorption, entails stronger valence forces akin to those in chemical bonding, resulting in irreversible or semi-reversible coverage that often requires higher energies and is more selective. The extent and mechanism of adsorption are often described by isotherms, such as the Langmuir model, which assumes monolayer adsorption on a homogeneous surface with no lateral interactions between adsorbates, providing a foundational equation for predicting equilibrium uptake: θ=Kp1+Kp\theta = \frac{K p}{1 + K p}, where θ\theta is the fractional surface coverage, KK is the equilibrium constant, and pp is the partial pressure of the adsorbate. Other models, including Freundlich and BET (Brunauer-Emmett-Teller), extend this to multilayer adsorption and porous materials, with the IUPAC classifying isotherms into six types (I-VI) based on adsorbent porosity and adsorbate-adsorbent affinity. Adsorption plays a critical role in numerous industrial and environmental applications, including gas purification (e.g., removal of impurities in ), water treatment (e.g., heavy metal ion removal using ), (e.g., surface reactions in heterogeneous catalysts), and for separation processes. In , is essential for activating reactants on solid surfaces, while aids in characterizing porous materials' surface area via adsorption at 77 . These applications leverage adsorption's efficiency in selective binding, contributing to processes like solvent recovery in chemical industries and control.

Fundamentals

Definition and Distinction from Absorption

Adsorption is defined as the adhesion of atoms, ions, or molecules from a gas, liquid, or dissolved solid to a surface, resulting in the formation of an adsorbate layer due to interactions with surface energy. This surface phenomenon occurs at the interface between the adsorbate (the substance being adsorbed) and the adsorbent (the material providing the surface), leading to an increase in concentration at that boundary. The term "adsorption" was coined in 1881 by German physicist Heinrich Kayser to describe the uptake of gases by solids like carbon, building on earlier observations of gas-solid interactions. Early studies, including Kayser's work, highlighted how gases accumulate on surfaces without penetrating the bulk, laying the foundation for understanding interfacial accumulation. A key distinction from absorption lies in the location and nature of the process: adsorption is strictly a surface-based where the adsorbate forms a or multilayer film on the external or internal surfaces of the adsorbent, whereas absorption involves the uniform penetration and distribution of a substance throughout the volume of another material. For example, in , gas molecules adsorb onto the surface of a solid catalyst to facilitate reactions, remaining localized at the interface, while in contrast, a solute like salt is absorbed into the bulk volume of during dissolution, becoming fully integrated into the . This surface specificity in adsorption allows for reversible processes in many cases, unlike the often more permanent incorporation in absorption. The driving force behind adsorption is primarily the high surface free energy of the adsorbent, which creates an energetic imbalance that promotes the attachment of adsorbates to minimize the overall free energy of the . Solids and liquids inherently possess excess due to unbalanced intermolecular forces at their boundaries, encouraging the accumulation of foreign to stabilize these interfaces. This role of surface energy underscores adsorption's importance in applications like purification and , where controlled surface interactions are essential.

Types of Adsorption

Adsorption processes are primarily classified into two categories—physisorption and chemisorption—based on the strength and nature of the interactions between the adsorbate molecules and the surface atoms of the adsorbent. This classification provides the foundation for understanding adsorption mechanisms, with involving weaker, non-covalent forces and featuring stronger, covalent or . Physisorption occurs through weak van der Waals forces, resulting in binding energies typically ranging from 1 to 40 kJ/mol. This process is reversible, allowing adsorbates to desorb easily upon changes in conditions, and it can form multilayers as molecules accumulate beyond the initial monolayer. Physisorption is favored at low temperatures, where thermal energy is insufficient to disrupt the weak interactions but sufficient for molecules to approach the surface without significant barriers. In contrast, involves the formation of chemical bonds between the adsorbate and surface atoms, with binding energies generally spanning 50 to 400 kJ/mol. It is often irreversible under ambient conditions or requires activation to reverse, limiting it to a due to the strong, site-specific bonding that saturates available surface sites. Chemisorption typically requires higher temperatures to overcome any barrier, making it more prevalent in catalytic processes or at elevated operational conditions. The key differences between and are summarized in the following table:
Aspect
Binding Energy (kJ/mol)1–4050–400
Activation EnergyLow or negligibleOften high (requires thermal activation)
SpecificityNon-specific (occurs on any surface)Site-specific (depends on surface chemistry)
ReversibilityHighly reversibleIrreversible or activated desorption
These distinctions are drawn from established analyses. Certain processes exhibit hybrid characteristics, such as dissociative chemisorption, where diatomic molecules like dissociate upon adsorption to form atomic species that bind chemically to the surface. A prominent example is the dissociative chemisorption of H₂ on surfaces, such as or , which is crucial in . The prevalence of versus is influenced by external factors including , , and surface chemistry. Lower temperatures and higher pressures promote physisorption by enhancing molecular proximity without needing to surmount energy barriers, while higher temperatures and compatible surface sites favor chemisorption. Surface chemistry plays a pivotal role, as chemisorption demands specific reactive sites, whereas physisorption occurs more universally across diverse surfaces.

Driving Forces and Surface Interactions

Adsorption occurs spontaneously when the change, ΔG\Delta G, for the process is negative, governed by the relation ΔG=ΔHTΔS\Delta G = \Delta H - T\Delta S, where ΔH\Delta H represents the enthalpy change from adsorbate-surface , TT is the absolute , and Δ[S](/page/Glossaryofcurling)\Delta [S](/page/Glossary_of_curling) is the change often negative due to loss of translational freedom of the adsorbate. This thermodynamic favorability arises primarily from exothermic interactions that outweigh entropic penalties, enabling the adsorbate to accumulate at the interface despite the general tendency toward uniform distribution in the bulk phase. At equilibrium, the process balances these energetic gains with the disorder reduction, with lower temperatures typically enhancing adsorption by minimizing the TΔST\Delta S term. The key intermolecular forces driving adsorption include van der Waals interactions—comprising dispersion forces between nonpolar entities, Keesom dipole-dipole alignments, and induced effects—along with electrostatic attractions such as ion- and ion-induced couplings. bonding provides directional specificity, particularly for molecules with polar groups like -OH or -NH, while covalent bonding dominates in stronger interactions where electron sharing occurs between adsorbate and surface atoms. These forces collectively lower the system's , with van der Waals typically contributing 5–40 kJ/mol and hydrogen bonds 10–40 kJ/mol in scenarios. Surface heterogeneity significantly influences interaction strength, as real surfaces feature active sites—such as coordinatively unsaturated atoms on metal oxides or catalysts—defects like steps and vacancies, and porous structures that increase accessible area and trap adsorbates. These non-uniform features create localized potential wells deeper than on ideal flat surfaces, enhancing binding at specific locales while weaker interactions occur elsewhere, leading to site-specific adsorption energies varying by tens of kJ/mol. , in materials like zeolites or activated carbons, further amplifies this by confining adsorbates in nanopores, promoting multilayer effects through repeated interactions. From a kinetic perspective, adsorption begins with the of adsorbate molecules or ions to the surface, governed by Fickian laws where flux is proportional to concentration gradients, followed by the adsorption step with a rate constant reflecting barrier crossing for attachment. Desorption, the reverse , involves overcoming energies typically higher for chemisorbed , resulting in rates that follow Arrhenius dependence on ; overall, equilibrium is reached when adsorption and desorption rates equalize. of adsorbed , often via hopping between sites, also plays a role in redistributing coverage, with barriers influenced by the same heterogeneity factors. Representative examples illustrate these principles: polar molecules such as adsorb preferentially on polar surfaces like silica through strong hydrogen bonding and interactions, achieving coverages up to monolayers at ambient conditions, whereas nonpolar species like interact weakly via dispersion forces on hydrophobic , requiring cryogenic temperatures for significant uptake. Physisorption relies predominantly on these weaker van der Waals and electrostatic forces, contrasting with chemisorption's covalent character.

Adsorption Isotherms and Models

Langmuir Isotherm

The Langmuir isotherm, developed by in 1918, describes the adsorption of molecules onto a solid surface under conditions of coverage and equilibrium. This model assumes that adsorption occurs on a fixed number of identical and uniformly distributed sites on the surface, with each site accommodating only one adsorbate molecule, leading to a saturated at high pressures. Additionally, the model posits no lateral interactions between adsorbed molecules and dynamic equilibrium between adsorption and desorption rates. The derivation begins with the law of mass action applied to the adsorption equilibrium. Consider the reversible process where gas molecules (A) adsorb onto vacant surface sites (S) to form an adsorbed layer (AS):
\ceA(g)+S<=>AS\ce{A(g) + S <=> AS}
The forward rate is proportional to the gas pressure PP and the fraction of vacant sites (1θ)(1 - \theta), while the reverse rate depends on θ\theta. At equilibrium, these rates balance, yielding the fractional surface coverage θ\theta as:
θ=KP1+KP\theta = \frac{K P}{1 + K P}
where KK is the equilibrium constant related to the adsorption-desorption rate constants.
For practical applications, the isotherm is expressed in terms of the amount adsorbed qq (e.g., in moles per unit mass of adsorbent) relative to the maximum capacity qmq_m:
q=qmKP1+KPq = \frac{q_m K P}{1 + K P}
This form highlights the approach to saturation as PP increases. To fit experimental data, linear transformations are used, such as the reciprocal plot:
1q=1qm+1qmKP\frac{1}{q} = \frac{1}{q_m} + \frac{1}{q_m K P}
which yields a straight line when plotting 1/q1/q versus 1/P1/P, allowing estimation of qmq_m and KK.
In gas-solid adsorption systems, the Langmuir isotherm finds key applications, such as modeling the surface coverage in catalytic processes like ammonia synthesis, where the Langmuir-Hinshelwood mechanism describes the dissociative adsorption of and on iron catalysts. Despite its foundational role, the Langmuir model has limitations, as it assumes ideal adsorption on homogeneous surfaces and fails to account for multilayer formation or site heterogeneity, which are common in real systems.

Freundlich Isotherm

The Freundlich isotherm is an empirical model for adsorption, originally developed by Herbert Freundlich in to describe the uptake of organic compounds, such as dyes, from aqueous solutions onto surfaces. This model captures the non-ideal behavior observed in heterogeneous systems, where adsorption does not follow the uniform site assumptions of more theoretical approaches. The isotherm is mathematically expressed as
q=KfP1/nq = K_f P^{1/n}
where qq is the amount of adsorbate per unit mass of adsorbent at equilibrium, PP is the equilibrium (or concentration for solutions), KfK_f is the Freundlich capacity factor indicating adsorption strength, and 1/n1/n is the heterogeneity factor with n>1n > 1 for favorable adsorption processes. A linearized form,
logq=logKf+1nlogP\log q = \log K_f + \frac{1}{n} \log P
facilitates parameter estimation by plotting logq\log q against logP\log P, yielding a straight line whose slope and intercept provide 1/n1/n and KfK_f, respectively.
The parameter 1/n1/n reflects surface heterogeneity, with values between 0 and 1 indicating normal on energetically diverse sites and values greater than 1 suggesting effects or chemisorption-like . This empirical flexibility allows the model to approximate a distribution of adsorption site energies, tying into broader surface interaction concepts. Compared to ideal models assuming homogeneous surfaces, the Freundlich isotherm offers advantages in fitting experimental data for real adsorbents at low coverages, where heterogeneity dominates. However, its limitation lies in the absence of a maximum adsorption capacity, causing predicted [q](/page/Q)[q](/page/Q) to rise indefinitely at high pressures, which does not align with physical saturation. In practice, it excels in applications like dye removal from ; for instance, adsorption on follows the Freundlich model with 1/n0.421/n \approx 0.42 and Kf1.23K_f \approx 1.23 mg/g (L/mg)^{1/n} at 20°C, demonstrating effective heterogeneous binding.

BET Isotherm and Multilayer Adsorption

The Brunauer-Emmett-Teller () theory represents a significant extension of the Langmuir isotherm to account for multilayer adsorption on solid surfaces. Unlike the single-layer limitation of the Langmuir model, BET assumes that adsorption can form multiple layers, with the first layer bound to the surface by stronger interactions (often chemisorption-like), while subsequent layers are held by weaker van der Waals forces akin to in the liquid phase, allowing for an infinite number of layers in theory. This model was developed by Stephen Brunauer, Paul Hugh Emmett, and in 1938, published in the Journal of the American Chemical Society, where it was applied to explain gas adsorption behaviors observed experimentally on materials like iron catalysts. The theory has since become foundational for characterizing porous materials, particularly in measuring surface area and pore structure through experiments using gases such as at 77 K. The derivation of the BET isotherm builds directly on the Langmuir approach by treating adsorption as a series of equilibrium processes across layers. It begins by considering the surface sites: empty sites adsorb gas to form the first layer, while molecules in the first layer serve as "sites" for the second layer, and so on, with the uppermost layer adsorbing and desorbing like a bulk liquid. The rate of adsorption onto a given layer equals the rate of desorption from that layer at equilibrium. For the first layer, the adsorption energy E1E_1 differs from the energy ELE_L of subsequent layers. Summing the coverages over all layers—using Langmuir-like expressions for each—yields the total adsorbed amount after algebraic manipulation and approximation for infinite layers. The resulting BET isotherm equation for the amount adsorbed qq (typically in moles per unit mass) at relative pressure P/P0P/P_0 is: q=qmCP(P0P)[1+(C1)PP0]q = \frac{q_m C P}{(P_0 - P) \left[1 + (C - 1) \frac{P}{P_0}\right]} where qmq_m is the capacity, P0P_0 is the of the adsorbate, and CC is a constant related to the difference in adsorption energies: C=exp(E1ELRT)C = \exp\left(\frac{E_1 - E_L}{RT}\right), with RR the and TT the . Values of C>1C > 1 indicate stronger binding in the first layer compared to , leading to a sigmoid-shaped isotherm. A linear form of the equation is often used for practical analysis: Pq(P0P)=1qmC+C1qmCP0P\frac{P}{q (P_0 - P)} = \frac{1}{q_m C} + \frac{C - 1}{q_m C P_0} P Plotting the left-hand side versus P/P0P/P_0 yields a straight line in the appropriate range, from which qmq_m and CC are obtained via and intercept. The BET model enables surface area determination by leveraging qmq_m. The total surface area SS (in m²/g) is calculated as S=qm×NA×σS = q_m \times N_A \times \sigma, where NAN_A is Avogadro's number (6.022 × 10²³ mol⁻¹) and σ\sigma is the cross-sectional area of the adsorbed (e.g., 0.162 nm² for N₂). This method assumes uniform coverage in the and is routinely applied to estimate the accessible surface in powders, catalysts, and . In porosity measurement, analysis of isotherms reveals multilayer formation indicative of meso- and macropores, complementing techniques like mercury porosimetry. Despite its utility, the isotherm has limitations, particularly becoming invalid near saturation where P/P0P/P_0 approaches 1, as the equation predicts infinite adsorption, diverging from real behavior. It is most applicable to IUPAC Type II isotherms, characteristic of nonporous or macroporous adsorbents showing unrestricted multilayer growth after completion, and Type III isotherms, where weak interactions lead to initial multilayer adsorption without a distinct at coverage. For Type I (microporous) or Type IV/V (mesoporous with ), modifications or alternative models are often needed. The parameter CC provides qualitative insight into differences between layers, with higher CC reflecting greater E1ELE_1 - E_L.

Thermodynamic and Molecular Aspects

Adsorption Enthalpy and Energetics

Adsorption processes are inherently exothermic, releasing upon the binding of adsorbate molecules to a surface, which reflects the strength of the adsorbate-surface interactions. The of adsorption, denoted as ΔH_ads, quantifies this heat exchange and is crucial for understanding the of the process. It can be distinguished as integral enthalpy, representing the average heat released for adsorbing a certain amount of molecules from an empty surface, or differential enthalpy, which measures the incremental heat for adding one more at a given coverage. The differential enthalpy is particularly important as it varies with surface coverage, often decreasing as sites become occupied due to lateral interactions among adsorbates. The differential enthalpy of adsorption can be derived from the Clausius-Clapeyron equation applied to adsorption isotherms at constant coverage θ: ΔHdiff=RT2(lnPT)θ\Delta H_{\text{diff}} = -RT^2 \left( \frac{\partial \ln P}{\partial T} \right)_{\theta} where R is the gas constant, T is temperature, and P is the equilibrium pressure. This equation links the temperature dependence of pressure at fixed coverage to the heat of adsorption. Common measurement techniques include direct calorimetry, which captures the heat released during adsorption under controlled conditions to yield isosteric heats, and indirect methods such as van't Hoff plots constructed from adsorption isotherms at varying temperatures. Calorimetry provides precise, real-time data but requires specialized equipment, while van't Hoff analysis offers accessibility from equilibrium data alone. Typical values of adsorption enthalpy differ markedly between physisorption and chemisorption. Physisorption, involving weak van der Waals forces, exhibits enthalpies in the range of 20-50 kJ/mol, whereas chemisorption, characterized by formation, ranges from 80-250 kJ/mol, reflecting stronger interactions. These values often show coverage dependence, with initial adsorption on high-energy sites yielding higher |ΔH_ads| that diminishes at higher coverages due to site heterogeneity and adsorbate-adsorbate repulsions. In activated adsorption processes, an enthalpy-entropy compensation effect frequently occurs, where higher adsorption enthalpies correlate with more negative entropy changes due to restricted molecular freedom on the surface, yet the net free energy remains favorable. This trade-off influences the temperature dependence of adsorption rates and equilibrium constants across different systems. The exothermic nature of adsorption drives its spontaneity (negative ΔG) under ambient conditions but poses challenges for reversibility, as regeneration of the adsorbent requires supplying sufficient energy—often via heating or pressure reduction—to overcome the binding enthalpy and desorb the molecules. In multilayer adsorption, the first layer typically exhibits higher enthalpies than subsequent layers, which approach the heat of liquefaction of the adsorbate. A representative example is the of N₂ on iron catalysts used in synthesis, where the differential is approximately -146 kJ/mol, derived from temperature-programmed desorption studies, highlighting the strong dissociative binding that underpins catalytic efficiency but necessitates high-temperature regeneration.

Single-Molecule Adsorption Mechanisms

Single-molecule adsorption mechanisms describe the probabilistic and dynamic processes by which individual adsorbate molecules interact with a surface, providing a microscopic foundation for macroscopic adsorption behavior. These mechanisms emphasize the nature of adsorption events, where molecules may approach, bind, , or desorb based on kinetic rates influenced by surface sites and barriers. Unlike ensemble-averaged models, single-molecule perspectives reveal transient states and diffusion pathways that bridge individual events to collective surface coverage. A foundational kinetic framework for single-molecule adsorption is provided by Langmuir kinetics, which models the rate of adsorption as proportional to the gas pressure PP and the fraction of unoccupied sites (1θ)(1 - \theta), while desorption depends on the occupied fraction θ\theta. The net rate equation is given by dθdt=kads(1θ)Pkdesθ,\frac{d\theta}{dt} = k_{\text{ads}} (1 - \theta) P - k_{\text{des}} \theta, where kadsk_{\text{ads}} and kdesk_{\text{des}} are the adsorption and desorption rate constants, respectively. At equilibrium, dθ/dt=0d\theta/dt = 0, yielding the Langmuir isotherm θ=KP1+KP\theta = \frac{K P}{1 + K P} with K=kads/kdesK = k_{\text{ads}}/k_{\text{des}}. This model assumes non-activated adsorption onto discrete sites, capturing the competition between impinging molecules and site availability for isolated adsorbates. Transition state theory (TST) further elucidates these kinetics by incorporating activation barriers for adsorption and desorption. In TST, the rate constants are expressed as k=νexp(Ea/RT)k = \nu \exp(-E_a / RT), where ν\nu is the pre-exponential factor related to the attempt frequency, EaE_a is the activation energy barrier, RR is the gas constant, and TT is temperature. For adsorption, EaE_a may be near zero for non-activated processes, but desorption typically involves overcoming the adsorption energy plus any additional barriers, leading to thermally activated escape from the bound state. This framework highlights how enthalpy barriers influence single-molecule residence times on the surface. Precursor states play a critical role in many adsorption pathways, representing transient, weakly bound physisorbed intermediates that molecules occupy before transitioning to strongly bound chemisorbed states. These states, often stabilized by van der Waals forces at a distance of several angstroms from the surface, allow mobile exploration of sites and enhance sticking probabilities by delaying direct . For instance, in CO adsorption on metal surfaces, the precursor facilitates chemisorption by lowering the effective barrier through energy dissipation during surface accommodation. Direct observation of single-molecule dynamics has been enabled by techniques such as scanning tunneling microscopy (STM) and , which visualize hopping and diffusion events on surfaces. STM studies reveal that adsorbed molecules undergo random walks via thermally activated jumps between adjacent sites, with diffusion coefficients scaling as Dexp(Ed/RT)D \propto \exp(-E_d / RT), where EdE_d is the diffusion barrier. AFM complements this by probing mechanical interactions during motion, confirming sub-angstrom displacements in real time. These methods have quantified hopping rates on the order of seconds to minutes at for light adsorbates. Coverage effects modulate single-molecule behavior, transitioning from isolated diffusion at low coverages to constrained motion at higher ones. At low θ\theta (<0.1 monolayer), attractive adsorbate-adsorbate interactions promote nucleation and island formation, where molecules cluster into compact domains to minimize edge energy. As coverage increases, repulsive interactions dominate, favoring uniform distributions or ordered overlayers that suppress further island growth. This evolution influences overall kinetics, as island edges exhibit distinct binding and diffusion properties compared to terrace sites. A representative example is CO adsorption on the Pt(111) surface, where single-molecule dynamics illustrate these mechanisms. At low coverages, isolated CO molecules bind atop Pt atoms and diffuse via hopping with barriers around 0.6 eV, forming transient islands due to lateral attractions. As coverage approaches 0.5 monolayer, repulsive dipole interactions stabilize the ordered c(4×2) phase, reducing mobility and altering desorption rates. STM observations confirm precursor-mediated adsorption, with molecules entering physisorbed states before chemisorbing, bridging individual events to the equilibrium overlayer structure.

Quantum Mechanical Modeling of Surfaces

Quantum mechanical modeling of surfaces plays a crucial role in understanding adsorption processes at the atomic level, particularly for porous materials where classical approaches fall short in capturing electronic and quantum effects. Density Functional Theory (DFT) is widely employed to compute binding energies and determine preferred adsorption sites on surfaces, providing insights into adsorbate-surface interactions. For instance, DFT calculations reveal that adsorption energies on transition metal surfaces can vary significantly depending on the functional used, with typical values ranging from -1 to -2 eV for common adsorbates like CO or O on metals. The B3LYP hybrid functional, often augmented with dispersion corrections, is particularly effective for modeling these interactions in catalytic systems, as it balances accuracy in describing both covalent and van der Waals forces. Such computations enable the prediction of site selectivity, where adsorbates preferentially bind to high-coordination sites like atop or bridge positions on metal surfaces. Complementing DFT, Grand Canonical Monte Carlo (GCMC) simulations in the grand canonical ensemble are used to predict adsorption isotherms within confined pores, accounting for ensemble averages of particle number under varying chemical potential. These simulations model the filling of slit-shaped or cylindrical pores, showing how adsorbate density increases nonlinearly with pressure due to pore geometry effects. For carbonaceous materials, GCMC integrates with DFT-derived potentials to simulate gas uptake, reproducing experimental isotherms for gases like Kr or CH4 in micropores with deviations under 10%. Non-Local Density Functional Theory (NLDFT) extends classical models like the BET isotherm by incorporating nonlocal correlations in fluid density profiles, offering more accurate pore size distributions (PSD) for microporous materials with pores below 2 nm. NLDFT kernels, calibrated against simulated isotherms, deconvolute experimental N2 or Ar adsorption data to yield PSDs that reveal heterogeneity in amorphous carbons or zeolites, improving surface area estimates by up to 20% over BET in the micropore regime. Applications of these methods are prominent in modeling CO2 adsorption in metal-organic frameworks (MOFs), where DFT identifies open metal sites as high-affinity locations with binding energies around -0.3 to -0.5 eV for CO2. Quantum-derived potentials from DFT feed into GCMC simulations to predict selective CO2 uptake in flexible MOFs, aiding the design of materials for carbon capture with capacities exceeding 5 mmol/g at ambient conditions. However, challenges arise from quantum effects in confinement, such as zero-point energy contributions that alter adsorption energetics in zeolites; for hydrogen isotopes, these effects lead to inverse kinetic isotopic sieving, where heavier D2 diffuses faster than H2 at low temperatures due to reduced zero-point vibrations in tight pores. In zeolites, quantum delocalization can shift binding energies by 0.1-0.2 eV, complicating predictions for light gases. Recent advances post-2020 incorporate machine learning (ML) potentials trained on DFT data to enable large-scale simulations of adsorption in complex porous systems, bypassing the computational cost of ab initio methods. These ML models, such as neural network potentials, accurately reproduce adsorption isotherms and diffusion in MOFs for gases like H2 or CO2, with errors below 5% relative to DFT benchmarks, facilitating high-throughput screening of thousands of structures. As of 2025, further progress includes multi-feature deep learning frameworks for predicting CO adsorption energies and hybrid computational fluid dynamics-ML models for simulating adsorption processes in porous media. For instance, ML-enhanced GCMC has optimized H2 storage in MOFs by predicting binding modes and temperature-dependent distributions in pores.

Adsorbents and Materials

General Characteristics of Adsorbents

Adsorbents are materials designed to capture molecules or ions on their surfaces through physical or chemical interactions, with key properties including high specific surface area, appropriate porosity, selectivity, and stability under operational conditions. The specific surface area, often measured via the Brunauer-Emmett-Teller (BET) method, typically exceeds 500 m²/g for effective adsorbents, enabling greater interaction sites for adsorbates. Porosity is classified by the International Union of Pure and Applied Chemistry (IUPAC) into micropores (<2 nm), mesopores (2-50 nm), and macropores (>50 nm), which determine the accessibility and of adsorbates within the material. Selectivity refers to the preferential adsorption of target species over others in a mixture, influenced by pore size matching and surface chemistry. Stability encompasses mechanical integrity, thermal resistance, and chemical inertness, ensuring the material withstands repeated cycles without degradation. Regeneration of adsorbents is essential for reusability and involves methods such as thermal swing adsorption (TSA), where heat desorbs the adsorbate; , which uses differentials for desorption; and chemical regeneration, employing solvents or reactive agents to displace bound . These techniques balance energy efficiency and material longevity, with TSA suitable for heat-stable adsorbents and PSA for rapid, low-energy processes. Adsorption capacity is quantified through static and dynamic metrics, where static capacity measures equilibrium uptake in a , while dynamic capacity assesses performance in flow-through setups, often 50-70% of static values due to kinetic limitations. The working capacity, a critical indicator for cyclic operations, is defined as the difference between adsorption loading (qadsq_{ads}) and desorption loading (qdesq_{des}), expressed as Δq=qadsqdes\Delta q = q_{ads} - q_{des}, highlighting the reversible fraction available for reuse. The behavior of adsorbents is characterized by IUPAC-classified isotherms: Type I for microporous materials with saturation; Type II for nonporous or macroporous surfaces with multilayer formation; Type III for weak interactions on nonporous substrates; and Type IV for mesoporous structures exhibiting due to . These classifications guide based on expected adsorption mechanisms. Economic considerations in adsorbent selection include production and costs, operational lifetime (often hundreds of cycles), and environmental impacts such as and waste generation during synthesis or disposal. Lifecycle assessments emphasize sustainable sourcing to minimize carbon footprints. Standardized testing ensures reliable capacity evaluation, with ASTM methods like D3860 for aqueous adsorptive capacity and D5160 for gas-phase dynamic adsorption providing protocols for breakthrough volume and saturation limits under controlled conditions.

Activated Carbon

Activated carbon is produced from carbonaceous precursors such as wood, coal, or agricultural residues through a two-step process involving followed by . occurs under inert atmospheres at temperatures around 400-600°C to form a char, while develops either physically or chemically. Physical employs or as oxidizing agents at 700-1000°C, etching the carbon surface to create pores. Chemical , often using (KOH) or , proceeds at 800-1000°C and allows greater control over pore size distribution, yielding higher at lower temperatures compared to physical methods. The structure of consists of featuring disordered, graphene-like sheets arranged in a turbostratic configuration, with cross-linking and defects contributing to its . Micropores, defined as those smaller than 2 nm, dominate the pore volume in many s, accounting for up to 80-90% of the total and enabling high adsorption capacities for small molecules. This hierarchical pore network, including mesopores (2-50 nm) and some macropores (>50 nm), arises from the selective removal of carbon atoms during , resulting in a highly interconnected void space. Key properties of include exceptionally high specific surface areas, often exceeding 2000 m²/g and reaching up to 3000 m²/g in optimized forms, which underpin its adsorption efficacy. Its non-polar, hydrophobic surface exhibits strong affinity for organic compounds through π-π electron donor-acceptor interactions between the aromatic sheets and adsorbate molecules, making it particularly effective for non-polar and weakly polar organics. However, unmodified shows limited adsorption for highly polar gases or ions due to weak electrostatic interactions, necessitating surface modifications for enhanced selectivity. To tailor for specific adsorbates, impregnation with metal oxides or chemicals is commonly applied; for instance, impregnation with caustic soda or enhances (H₂S) removal by promoting , achieving capacities up to 140 kg H₂S per m³ of carbon compared to 10-20 kg for untreated material. In applications, is widely used for air purification, where it captures volatile organic compounds and odors in HVAC systems and industrial emissions. It also plays a critical role in gold recovery from leaching solutions via the carbon-in-pulp process, adsorbing -cyanide complexes with efficiencies over 99% under optimized conditions. Limitations include its reduced performance for polar gases without functionalization and challenges in regeneration, which typically requires thermal treatment at 800-950°C, incurring significant energy costs of 1-2 GJ per ton and contributing to operational expenses.

Silica Gel and Zeolites

Silica gel consists of dioxide (SiO₂) produced through the sol-gel process, involving and of alkoxide precursors like tetraethoxysilane to form a , followed by to yield a porous structure. The step, often conducted at elevated temperatures around 110–180°C, removes solvent and stabilizes the network while preserving high . This material typically achieves a surface area of 300–800 m²/g, enabling extensive adsorption sites. Its hydrophilic nature arises from surface (Si-OH) groups, which form hydrogen bonds with polar adsorbates like . Zeolites are crystalline aluminosilicates featuring a three-dimensional framework of tetrahedra that creates cage-like pores with uniform diameters of 0.3–1 nm. They are synthesized via hydrothermal methods, where sources of silica and alumina, such as and , are mixed with structure-directing agents like tetrapropylammonium hydroxide and heated under autogenous pressure at 100–200°C for hours to days, yielding topologies like . For instance, forms through and in alkaline media, resulting in microporous channels that dictate selective access. Key properties of zeolites include shape-selective adsorption, exemplified by type 5A zeolite, which features 5 Å pores that permit entry of linear n-paraffins while excluding branched isomers due to steric constraints. Additionally, zeolites support , where extra-framework cations like Na⁺ or Ca²⁺ in the lattice can be replaced with other ions, altering acidity and adsorption affinity. , in contrast, enables reversible water adsorption up to 40 wt% at ambient conditions, driven by multilayer on its silanol-rich surface. Both materials serve as drying agents and catalysts, with excelling in broad-polarity adsorption for processes and zeolites in precise molecular sieving for catalytic cracking. However, zeolites face limitations from hydrothermal , where exposure to at high temperatures (>500°C) can hydrolyze Al-O-Si bonds, leading to framework collapse and reduced performance. In comparison, offers greater versatility for polar solutes across a wide range, whereas zeolites provide superior selectivity for size- and shape-based separations in non-aqueous environments.

Applications in Environmental and Energy Systems

Water Adsorption and Purification

Adsorption plays a crucial role in by selectively removing contaminants and enabling through the capture of molecules while excluding salts and impurities. In processes, adsorbents like metal-organic frameworks (MOFs) and zeolites exploit their porous structures to adsorb from saline solutions under low pressure, producing pure upon desorption. This mechanism relies on the high affinity of these materials for , driven by bonding and within their micropores, allowing selective uptake that leaves salts behind. Key processes include fixed-bed adsorption systems, where water flows through columns packed with adsorbents to target specific contaminants such as and . For instance, iron oxide-impregnated in fixed-bed columns effectively removes (III) from , with breakthrough curves indicating saturation after processing several bed volumes depending on initial concentrations. Similarly, Al³⁺-pretreated low-silica synthetic zeolites in fixed-bed setups achieve high fluoride uptake through and surface complexation, demonstrating equilibrium capacities exceeding 10 mg/g under neutral pH conditions. Prominent adsorbents for these applications include modified activated carbons tailored for organic pollutants and silica gels for humidity regulation. Surface-modified activated carbons, such as those treated with or iron impregnation, enhance adsorption of dissolved natural and pesticides via increased surface polarity and pore accessibility, removing up to 90% of targeted organics from . Silica gels, known for their strong water affinity due to polar groups, are widely used in dehumidification to control moisture in purification systems, adsorbing up to 40% of their weight in at relative above 60%. In adsorption desalination (AD) cycles, systems alternate between adsorption and desorption phases to produce , often achieving capacities of 100 m³/day in pilot-scale installations using or beds. Advanced MOFs like Al-fumarate exhibit superior performance, with water production rates of 23.5 m³/tonne of adsorbent per day and efficient regeneration over multiple cycles via low-energy heating. Breakthrough curves in these cycles typically show stable operation until 80-90% saturation, followed by regeneration to restore 95% of initial capacity after 50-100 cycles. Despite these advances, challenges persist in scaling up adsorption systems, including mineral scaling from salt precipitation and that clogs pores and reduces efficiency by up to 50% over time. Recent post-2020 developments in membranes, incorporating MOFs or zeolites into thin-film matrices, address these issues by enhancing antifouling properties through hydrophilic surfaces and dispersion, achieving improvements of 20-30% while minimizing leaching.

Carbon Capture and Storage

Adsorption-based methods play a crucial role in (CCS) by selectively separating CO2 from gas streams using solid sorbents, offering advantages over liquid absorption in terms of lower energy requirements for regeneration and reduced . These techniques are particularly suited for point-source emissions and (DAC), where CO2 concentrations vary from dilute flue gases (3-15% in post-combustion scenarios) to ultra-low levels (0.04% in ambient air). Key processes include post-combustion capture, which targets CO2 from power plant exhaust using amine-impregnated solids or zeolites, and pre-combustion capture, which involves shifting to separate CO2 from hydrogen-rich streams via (PSA) with zeolites. In post-combustion capture, PSA and vacuum swing adsorption (VSA) cycles utilize zeolites like 13X, which exhibit high CO2 adsorption capacities of up to 3-5 mmol/g at 298 and 0.15 bar due to strong quadrupolar interactions with the framework's cationic sites, enabling >90% CO2 recovery in multi-bed systems. Pre-combustion processes employ similar PSA cycles with zeolites such as 5A or 13X to handle higher CO2 partial pressures (15-40%), achieving purities exceeding 95% and recoveries of 85-90% through pressure equalization and steps that minimize use. Working capacity, defined as the difference between adsorption and desorption loadings, is optimized in these cycles to 1-2 mmol/g, balancing selectivity and regeneration efficiency. Advanced adsorbents like metal-organic frameworks (MOFs), exemplified by , enhance selectivity through open metal sites that coordinate with CO2's moment, yielding capacities exceeding 4 mmol/g at 273 K and 1 bar, with isosteric heats of adsorption around 30-40 kJ/mol for chemisorption-enhanced uptake. Amine-impregnated activated carbons and zeolites further improve performance in post-combustion settings by forming carbamates with CO2, achieving working capacities of 2-3 mmol/g even under humid conditions. Langmuir isotherm models briefly illustrate this selectivity, where CO2/N2 separation factors reach 50-100 for optimized MOFs and zeolites at low pressures. Integration with geological storage is demonstrated in DAC pilots, such as ' facilities operational since 2017, which employ temperature-vacuum swing adsorption (TVSA) with amine-functionalized sorbents to capture CO2 from ambient air at scales up to 36,000 tons per year—as in the Mammoth plant operational since 2024—followed by mineralization in basaltic formations. Recent developments include amine-functionalized carbons, like polyethyleneimine (PEI)-loaded variants, which maintain >2 mmol/g capacity in humid flue gases (10-15% H2O) by leveraging water to facilitate formation, reducing selectivity loss to <20% compared to dry conditions. Economically, adsorption-based CCS faces challenges from regeneration energy penalties of 1-2 GJ/ton CO2, but targets aim for cost reductions to $100/ton by 2030 through scalable PSA/TSA cycles and low-cost sorbents like derivatives, potentially capturing 100-200 Mt CO2 annually from industrial sources. These advancements prioritize modular designs for , though humidity tolerance and sorbent stability remain key hurdles for widespread deployment.

Adsorption-Based Heating, Cooling, and Solar Systems

Adsorption-based systems utilize the reversible adsorption and desorption of , typically , onto solid sorbents such as or zeolites to achieve heating, cooling, or , offering an environmentally benign alternative to vapor-compression cycles by leveraging low-grade sources. These systems operate on thermodynamic cycles where input drives desorption, releasing the refrigerant vapor for and subsequent to produce cooling or heating effects. Primarily driven by or , they avoid the use of ozone-depleting substances like chlorofluorocarbons (CFCs), aligning with goals in building HVAC and applications. In adsorption chillers, silica gel-water pairs are commonly employed due to the high affinity of for and its thermal stability, enabling operation with driving heat temperatures between 80°C and 150°C from or solar sources. The (COP) for these systems typically ranges from 0.5 to 0.7, with optimized designs achieving up to 0.697 under controlled conditions. For instance, commercial units like the SorTech ACS08, a 7.5 kW silica gel-water adsorption , demonstrate practical deployment in buildings for , providing chilled water at 7-12°C. Solar heating applications leverage zeolite-water systems for seasonal thermal energy storage, where zeolites' microporous structure allows high water uptake at low relative , storing through adsorption during off-peak periods like nighttime or winter, and releasing it via desorption with solar during the day. This process exploits the exothermic adsorption to deliver space heating, with systems achieving energy densities up to 200-300 kWh/m³ in modified zeolites impregnated with hygroscopic salts for enhanced capacity. 13X, for example, has shown average energy storage densities of 129 kWh/m³ during desorption in open sorption setups. The basic adsorption cycle consists of four main components: an adsorber bed containing the , an where () evaporates at low to absorb and produce cooling, a condenser to liquefy the vapor under higher , and an expansion valve to return the to the . During the adsorption phase, the cooled bed draws vapor from the , lowering the pressure and enabling ; input then desorbs the , completing the cycle in an intermittent manner. Advanced multi-stage cycles, such as two- or four-bed configurations with recovery, enable continuous operation by staggering adsorption and desorption phases across beds, improving efficiency by 20-30% over single-bed systems. Performance metrics for adsorption chillers include specific cooling power (SCP), which measures cooling output per unit mass of adsorbent and can reach 50-200 W/kg in silica gel-water systems, with enhanced designs using heat and mass recovery achieving up to 250 W/kg. These values establish the scalability for medium-sized applications, such as , where SCP influences system footprint and cost-effectiveness. Key advantages of these systems include their reliance on natural refrigerants like , eliminating CFCs and minimizing , alongside seamless integration with intermittent sources like solar or for reduced consumption. However, challenges persist, such as bulky adsorber beds due to slow and kinetics, leading to larger footprints compared to vapor-compression units, and sensitivity to intermittent solar input requiring auxiliary heating for consistent performance. Recent advancements post-2022 focus on hybrid photovoltaic-thermal (PVT) integration with adsorption cooling, where PVT collectors supply both electrical power and low-grade to drive the cycle, enhancing overall to 40-50% in combined output for off-grid applications. Studies in regions like the MENA area have demonstrated fuel savings of up to 60% in institutional buildings through such hybrids, addressing via thermal storage.

Applications in Biological and Chemical Systems

Protein and Surfactant Adsorption

Protein adsorption at interfaces, particularly on hydrophobic surfaces, often results in the formation of irreversible multilayers, where proteins undergo conformational changes leading to denaturation and partial unfolding. This process is driven primarily by hydrophobic interactions and electrostatic forces, with the protein's native structure destabilizing upon contact to maximize surface interactions. A key dynamic in protein adsorption from complex mixtures, such as , is the Vroman effect, where initially adsorbed low-molecular-weight proteins are sequentially displaced by higher-affinity, larger proteins over time due to competitive binding. This effect, first observed in studies of blood-biomaterial interactions, highlights the time-dependent evolution of the adsorbed layer, influencing . In contrast, surfactant adsorption typically forms well-ordered at interfaces like air-water or solid-liquid, reducing and enabling applications in stabilization of colloids. , being amphiphilic molecules, adsorb with hydrophobic tails oriented away from the aqueous phase, lowering the (CMC) and promoting . The Gibbs adsorption isotherm quantifies this relationship, linking surface excess concentration to changes in : Γ=1RTdγdlnC\Gamma = -\frac{1}{RT} \frac{d\gamma}{d \ln C}, where Γ\Gamma is the surface excess, γ\gamma is , RR is the , TT is , and CC is bulk concentration. This thermodynamic framework has been validated for both ionic and non-ionic , providing insights into adsorption isotherms like Langmuir models for coverage. Mechanisms of adsorption for both proteins and surfactants involve a balance of electrostatic attractions, van der Waals forces, and hydrophobic effects, with pH, , and surface charge modulating outcomes. For instance, (HSA), a common model protein, adsorbs readily onto biomaterials like titanium oxides via hydrophobic domains, forming a conditioning film that mediates subsequent cellular responses. In mixed systems, can compete with or enhance protein adsorption, altering layer stability through co-adsorption or displacement. Applications of controlled protein and adsorption span biomedical fields, including systems where protein coatings on nanoparticles improve targeting and reduce . In biosensors, selective protein immobilization enables sensitive detection of analytes, while surfactant monolayers facilitate functionalization for electrochemical assays. However, uncontrolled adsorption leads to challenges like in membranes, where protein layers increase resistance and reduce efficiency in processes. Kinetically, adsorption often follows Langmuir-like models, with initial to the interface followed by attachment barriers, achieving equilibrium in seconds to minutes at low concentrations. Protein adsorption, however, is predominantly -limited, with arrival rates governed by bulk concentration and molecular weight, leading to slower multilayer buildup over minutes to hours; activation energies for unfolding add complexity. Predicting adsorption in complex fluids remains challenging due to multifactorial influences like shear flow and competitive species, complicating model scalability. Recent advances in with dissipation (QCM-D) measurements provide real-time insights into layer thickness, , and reversibility, revealing trapped water in protein films and aiding in antifouling design. For example, QCM-D has quantified adsorption on polymers, showing dissipation increases indicative of soft, hydrated layers that resist simple desorption.

Polymer Adsorption and Spillover Effects

Polymer adsorption onto solid surfaces typically results in complex chain conformations characterized by the "trains, loops, and tails" model, where trains consist of consecutive segments in direct contact with the surface, loops are unbound segments extending into the solution between trains, and tails are unbound end segments protruding from the surface. This model, developed through mean-field lattice theories, predicts that the fraction of trains increases with stronger adsorption , while loops and tails dominate in weaker binding regimes, leading to a diffuse adsorbed layer thickness on the order of the polymer's . The thermodynamics of polymer adsorption are often described using extensions of Flory-Huggins theory, which quantifies the free energy change through contributions from segmental interactions (via the Flory-Huggins parameter χ), entropy loss upon binding, and surface energy terms. In good solvents (low χ), chains adopt more extended conformations with longer loops and tails to maximize solvation, whereas in poor solvents (high χ), flatter profiles with more trains prevail; molecular weight influences this by scaling the layer thickness roughly as N^{1/2} in θ-solvents, where N is the degree of polymerization. Desorption often exhibits hysteresis due to kinetic barriers in chain reconfiguration and bridging between surfaces, requiring higher solvent quality or temperature for complete reversal compared to adsorption. For instance, polyethylene oxide (PEO) adsorption on silica surfaces demonstrates these effects, with surface diffusion coefficients for chain segments on the order of 10^{-11} cm²/s, facilitating loop formation and exchange dynamics. Applications of polymer adsorption leverage these conformational features; in flocculation processes, extended tails and loops enable bridging between colloidal particles to promote aggregation, as seen in with . In lubrication, adsorbed polymer layers, such as PEO on metal oxides, provide steric stabilization and reduce coefficients by up to 50% through sheared loop deformation. Recent advances in single-chain using (AFM) since 2015 have directly visualized these conformations, revealing transitions from to states for isolated chains on or silica, with loop lengths matching scaling predictions. Spillover effects, distinct from bulk polymer adsorption, involve the migration of atomic species across catalytic interfaces, particularly atomic from metal sites to supports. In prototypical systems like on alumina, H₂ dissociates on Pt nanoparticles to form mobile H atoms that spillover onto the Al₂O₃ support via , enhancing overall capacity by factors of up to 2-10 compared to metal-only sites. However, the phenomenon of hydrogen spillover remains controversial, with ongoing debates regarding its mechanisms, extent, and contributions in certain systems. Mechanisms are governed by surface diffusion coefficients, typically 10^{-8} to 10^{-10} cm²/s at 300-500 K, influenced by support and defect sites that facilitate reverse spillover during reaction. This phenomenon, first evidenced in isotopic exchange experiments, plays a key role in by extending reach beyond metal particles. In applications, spillover on Pt-supported catalysts improves performance by dissociating H₂ at Pt sites and migrating atomic H to carbon or oxide supports, reducing overpotentials by approximately 20-70 mV in polymer electrolyte membrane s (PEMFCs) and aiding in mitigating Pt poisoning. For example, on Pt/C electrodes, spillover creates a of reactive H, boosting oxidation reaction rates at low loadings. These effects underscore spillover's potential utility in energy systems, where controlled migration enhances efficiency without increasing content.

Adsorption in Viruses and Portal-Mediated Processes

In viral infections, adsorption represents the initial and critical step where viral capsid proteins bind to specific host cell receptors, facilitating entry and subsequent replication. For instance, the spike (S) protein of mediates adsorption to the (ACE2) receptor on human cells, triggering conformational changes that enable membrane fusion and viral genome release. This (RBD) interaction is characterized by high affinity, with mutations like N501Y enhancing binding strength and contributing to increased transmissibility. Similarly, in HIV-1, the envelope glycoprotein gp120 initiates adsorption by engaging the receptor, inducing a conformational shift that exposes coreceptor-binding sites such as or CXCR4. Mechanisms of viral adsorption often involve multivalent interactions, where multiple binding sites on the viral surface engage host receptors simultaneously to achieve stable attachment, and pH-dependent modulation that fine-tunes affinity during transit through endosomal compartments. In HIV-1, gp120 binding to is strengthened at lower levels, as alters electrostatic interactions and stabilizes the complex, a process modeled computationally to reveal strain-specific variations in founder versus chronic infections. These dynamics underscore adsorption's role in overcoming cellular barriers, with multivalency amplifying beyond single-site affinities. Beyond host interaction, adsorption features prominently in viral packaging, particularly in bacteriophages where portal proteins form a dodecameric channel at the procapsid vertex to facilitate DNA translocation into the . In phages like T7 and ϕ29, the portal acts as a DNA sensor, coordinating motor proteins to "adsorb" and package the under high , ensuring efficient assembly. Single-virus tracking via () microscopy has quantified these kinetics, revealing adsorption rates and dwell times for individual particles, as demonstrated in studies of and entry where binding precedes fusion. Applications of adsorption principles extend to antiviral therapies and biotechnology. Adsorption inhibitors, such as peptide-based entry blockers targeting gp120-CD4 interfaces, disrupt viral attachment and have informed designs for broad-spectrum antivirals against enveloped viruses. In gene therapy, engineered viral vectors like adeno-associated viruses exploit controlled adsorption to deliver therapeutic genes, with biomaterials enhancing targeting by modulating surface interactions. Recent studies from 2020-2025 on COVID-19 have explored nanomaterials for inactivation; for example, graphene oxide nanosheets adsorb SARS-CoV-2 via spike protein interactions, disrupting the viral envelope and achieving near-complete inactivation in vitro, while titanium dioxide surfaces bind and degrade the virus through photocatalytic effects.

References

Add your contribution
Related Hubs
User Avatar
No comments yet.