Sorption Hysteresis: A Statistical Thermodynamic Fluctuation Theory.

Hysteresis is observed commonly in sorption isotherms of porous materials. Still, there has so far been no unified approach that can both model hysteresis and assess its underlying energetics. Standard approaches, such as capillary condensation and isotherms based on interfacial equations of state, have not proved to be up to the task. Here, we show that a statistical thermodynamic approach can achieve the following needs simultaneously: (i) showing why adsorption and desorption transitions may be sharp yet continuous; (ii) providing a simple (analytic) isotherm equation for hysteresis branches; (iii) clarifying the energetics underlying sorption hysteresis; and (iv) providing macroscopic and nanoscopic perspectives to understanding hysteresis. This approach identifies the two pairs of parameters (determinable by fitting experimental data) that are required to describe the hysteresis: the free energy per molecule within the pore clusters and the cluster size in the pores. The present paper focuses on providing mechanistic insights to IUPAC hysteresis types H1, H2(a), and H2(b) and can also be applied to the isotherm types IV and V.

Actual Amount Adsorbed as Estimated from the Surface Excess Isotherm.

The amount of adsorption at equilibrium is commonly used for reporting solid/solution isotherms, despite the admonishment by the International Union of Pure and Applied Chemistry (IUPAC) against equating the surface excess (i.e., the measurable quantity for sorption, signifying the competitive sorption of adsorbate and solvent) with the actual amount adsorbed. The consensus, more generally stated, is that the surface excess cannot be divided into individual isotherms for sorbate and solvent unless simplifying model assumptions are introduced. Here we show, contrary to the IUPAC report, that there exists a simple method for assigning the total isotherm to the sorbate's actual amount adsorbed and to the individual solute isotherm. This requires a combination of isotherm and volumetric measurements. For dilute sorbates, we establish criteria to show if the total isotherm is dominated by the amount of sorption at the interface, in agreement with the common assumption in the practical literature. In the absence of the volume data, we propose an approximate yet more versatile method based on the specific surface area to carry out order-of-magnitude analysis to examine whether the actual amount adsorbed dominates surface excess. Application of our methods to the adsorption of sodium decyl sulfate on polystyrene latex, malachite green on activated carbons, and thiophenes on a metal-organic framework all demonstrated the dominance of the actual amount adsorbed, significantly simplifying isotherm analysis in terms of the underlying interactions (i.e., surface-sorbate and net self-interactions at the interface), eliminating the need for excess surface quantities. Analysis of fully miscible solvent-sorbate isotherms (e.g., the mixtures of organic solvents adsorbed on mesoporous silica and carbonaceous adsorbents) indicates the contributions from both sorbate and solvent isotherms.

Sorption from Solution: A Statistical Thermodynamic Fluctuation Theory.

Given an experimental solid/solution sorption isotherm, how can we gain insight into the underlying sorption mechanism on a molecular basis? Classifying sorption isotherms, for both completely and partially miscible solvent/sorbate systems, has been useful, yet the molecular foundation of these classifications remains speculative. Isotherm models, developed predominantly for solid/gas sorption, have been adapted to solid/solution isotherms, yet how their parameters should be interpreted physically has long remained ambiguous. To overcome the inconclusiveness, we establish in this paper a universal theory that can be used for interpreting and modeling solid/solution sorption. This novel theory shares the same theoretical foundation (i.e., the statistical thermodynamic fluctuation theory) not only with solid/gas sorption but also with solvation in liquid solutions and solution nonidealities. The key is the Kirkwood-Buff χ parameter, which quantifies the net self-interaction (i.e., solvent-solvent and sorbate-sorbate interactions minus solvent-sorbate interaction) via the Kirkwood-Buff integral in the same manner as the solvation theory and, unlike the Flory χ, is not limited to the lattice model. We will demonstrate that the Kirkwood-Buff χ is the key not only to isotherm classification but also to generalizing our recent statistical thermodynamic gas (vapor) isotherm, which is capable of fitting most of the solid/solution isotherm types.

Understanding Sorption Mechanisms Directly from Isotherms.

Currently, more than 100 isotherm models coexist for the six IUPAC isotherm types. However, no mechanistic insights can be reached when several models, each claiming a different mechanism, fit an experimental isotherm equally well. More frequently, popular isotherm models [such as the site-specific models like Langmuir, Brunauer-Emmett-Teller (BET), and Guggenheim-Anderson-de Boer (GAB)] have been applied to real and complex systems that break their basic assumptions. To overcome such conundrums, we establish a universal approach to model all isotherm types, attributing the difference to the sorbate-sorbate and sorbate-surface interactions in a systematic manner. We have generalized the language of the traditional sorption models (such as the monolayer capacity and the BET constant) to the model-free concepts of partitioning and association coefficients that can be applied across the isotherm types. Through such a generalization, the apparent contradictions, caused by applying the site-specific models alongside with cross-sectional area of sorbates for the purpose of surface area determination, can be eliminated straightforwardly.

Cooperativity in Sorption Isotherms.

We present a general theory of cooperativity in sorption isotherms that can be applied to sorbent/gas and sorbent/solution isotherms and is valid even when sorbates dissolve into or penetrate the sorbent. Our universal foundation, based on the principles of statistical thermodynamics, is the excess number of sorbates (around a probe sorbate), which can capture the cooperativities of sigmoidal and divergent isotherms alike via the ln-ln gradient of an isotherm (the excess number relationship). The excess number relationship plays a central role in deriving isotherm equations. Its combination with the characteristic relationship (i.e., a succinct summary of the sorption mechanism via the dependence of excess number on interfacial coverage or sorbate activity) yields a differential equation whose solution is an isotherm equation. The cooperative isotherm equations for convergent and divergent cooperativities derived from this novel method can be applied to fit experimental data traditionally fitted via various isotherm models, with a clear statistical thermodynamic interpretation of their parameters..

Long term phase separation dynamics in liquid crystal-enriched microdroplets obtained from binary fluid mixtures.

The dynamics of long term phase separation in binary liquid mixtures remains a subject of fundamental interest. Here, we study a binary liquid mixture, where the minority phase is confined to a liquid crystal (LC)-rich droplet, by investigating the evolution of size, defect and mesogen alignment over time. We track the binary liquid mixture evolving towards equilibrium by visualising the configuration of the liquid crystal droplet through polarisation microscopy. We compare our experimental findings with computational simulations and elucidate differences between bulk phases and confined droplets based on the respective thermodynamics of phase separation. Our work provides insights on how phase transitions on the microscale can deviate from bulk phase diagrams with relevance to other material systems, such as the liquid-liquid phase separation of polymer and protein solutions.

Surface Area Estimation: Replacing the Brunauer-Emmett-Teller Model with the Statistical Thermodynamic Fluctuation Theory.

Surface area estimation using the Brunauer-Emmett-Teller (BET) analysis has been beset by difficulties. The BET model has been applied routinely to systems that break its basic assumptions. Even though unphysical results arising from force-fitting can be eliminated by the consistency criteria, such a practice, in turn, complicates the simplicity of the linearized BET plot. We have derived a general isotherm from the statistical thermodynamic fluctuation theory, leading to facile isotherm fitting because our isotherm is free of the BET assumptions. The reinterpretation of the monolayer capacity and the BET constant has led to a statistical thermodynamic generalization of the BET analysis. The key is Point M, which is defined as the activity at which the sorbate-sorbate excess number at the interface is at its minimum (i.e., the point of strongest sorbate-sorbate exclusion). The straightforwardness of identifying Point M and the ease of fitting by the statistical thermodynamic isotherm have been demonstrated using zeolite 13X and a Portland cement paste. The adsorption at Point M is an alternative for the BET monolayer capacity, making the BET model and its consistency criteria unnecessary. The excess number (i) replaces the BET constant as the measure of knee sharpness and monolayer coverage, (ii) links macroscopic (isotherms) to microscopic (simulation), and (iii) serves as a measure of sorbate-sorbate interaction as a signature of sorption cooperativity in porous materials. Thus, interpretive clarity and ease of analysis have been achieved by a statistical thermodynamic generalization of the BET analysis.

Cooperative Sorption on Heterogeneous Surfaces.

Heterogeneous adsorbents, those composed of multiple surface and pore types, can result in stepwise isotherms that have been difficult to model. The complexity of these systems has often led to appealing to empirical equations without physical insights, unrealistic assumptions with many parameters, or applicability limited to a particular class of isotherms. Here, we present a statistical thermodynamic approach to model stepwise isotherms, those consisting of either an initial rise followed by a sigmoid or multiple sigmoidal steps, founded on the rigorous statistical thermodynamic theory of sorption. Our only postulates are (i) the finite ranged nature of the interface and (ii) the existence of several different types of microscopic interfacial subsystems that act independently in sorption. These two postulates have led to the superposition scheme of simple surface (i.e., Langmuir type) and cooperative isotherms. Our approach has successfully modeled the adsorption on micro-mesoporous carbons, gate-opening adsorbents, and hydrogen-bonded organic frameworks. In contrast to the previous models that start with a priori assumptions on sorption mechanisms, the advantages of our approach are that it can be applied universally under the above two postulates and that all of the fitting parameters can be interpreted with statistical thermodynamics, leading to clear insights on sorption mechanisms.

The impact of size and shape in the performance of hydrotropes: a case-study of alkanediols.

Inspired by the recently proposed cooperative mechanism of hydrotropy, where water molecules mediate the aggregation of hydrotrope around the solute, this work studies the impact of apolar volume and polar group position on the performance of hydrotropes. To do so, the ability of two different families of alkanediols (1,2-alkanediols and 1,n-alkanediols) to increase the aqueous solubility of syringic acid is initially investigated. Interestingly, it is observed that in the dilute region (low hydrotrope concentration), the relative position of the hydroxyl groups of the alkanediols does not impact their performance. Instead, their ability to increase the solubility of syringic acid correlates remarkably well with the size of their alkyl chains. However, this is not the case for larger hydrotrope concentrations, where 1,2-alkanediols are found to perform, in general, better than 1,n-alkanediols. These seemingly contradictory findings are reconciled using theoretical and experimental techniques, namely the cooperative model of hydrotropy and chemical environment probes (Kamlet-Taft and pyrene polarity scales). It is found that the number of hydrotropes aggregated around a solute molecule does not increase linearly with the apolar volume of the former, reaching a maximum instead. This maximum is discussed in terms of competing solute-hydrotrope and hydrotrope-hydrotrope interactions. The results suggest that hydrotrope self-aggregation is more prevalent in 1,n-alkanediols, which negatively impacts their performance as hydrotropes. The results reported in this work support the cooperative model of hydrotropy and, from an application perspective, show that hydrotropes should be designed taking into consideration not only their apolar volume but also their ability to stabilize their self-aggregation in water, which negatively impacts their performance as solubility enhancers.

Assessing the hydrotropic effect in the presence of electrolytes: competition between solute salting-out and salt-induced hydrotrope aggregation.

Water solubility enhancement is a long-standing challenge in a multitude of chemistry-related fields. Hydrotropy is a simple and efficient method to improve the solubility of hydrophobic molecules in aqueous media. However, the mechanism behind this phenomenon remains controversial. Herein the impact of salt doping on the hydrotropy phenomenon is determined experimentally using the ionic liquid (IL) 1-butyl-3-methylimidazolium chloride ([C4mim]Cl) as a hydrotope and vanillin as a solute. Hydrophobic interactions were found to be central to the aggregation of the hydrotrope around the solute, and seem to drive hydrotropy. Furthermore, 1H-NMR analysis indicates that hydrotrope-solute interactions present a degree of site-specificity. The addition of chloride salts in the presence of higher IL concentrations promotes a greater relative decrease of the vanillin solubility than in the corresponding system without the IL. This was assigned to the negative impact of increased hydrotrope pre-aggregation in the presence of inorganic salts. The results were rationalised using statistical thermodynamics through which hydrotrope aggregation prior to solute addition is shown to be detrimental to the hydrotropic effect, seemingly confirming solute-induced clustering of the hydrotrope to be the predominant mechanism of hydrotropy.

Reaction Optimization for Greener Chemistry with a Comprehensive Spreadsheet Tool.

Green chemistry places an emphasis on safer chemicals, waste reduction, and efficiency. Processes should be optimized with green chemistry at the forefront of decision making, embedded into research at the earliest stage. To assist in this endeavor, we present a spreadsheet that can be used to interpret reaction kinetics via Variable Time Normalization Analysis (VTNA), understand solvent effects with linear solvation energy relationships (LSER), and calculate solvent greenness. With this information, new reaction conditions can be explored in silico, calculating product conversions and green chemistry metrics prior to experiments. The application of this tool was validated with literature case studies. Reaction performance was predicted and then confirmed experimentally for examples of aza-Michael addition, Michael addition, and an amidation. The combined analytical package presented herein permits a thorough examination of chemical reactions, so that the variables that control reaction chemistry can be understood, optimized, and made greener for research and education purposes.

Sorption: A Statistical Thermodynamic Fluctuation Theory.

Can the sorption mechanism be proven by fitting an isotherm model to an experiment? Such a question arises because (i) multiple isotherm models, with different assumptions on sorption mechanisms, often fit an experimental isotherm equally well, (ii) some isotherm models [such as Brunauer-Emmett-Teller (BET) and Guggenheim-Anderson-de Boer (GAB)] fit experimental isotherms that do not satisfy the underlying assumptions of the model, and (iii) some isotherms (such as Oswin and Peleg) are empirical equations that do not have a well-defined basis on sorption mechanisms. To overcome these difficulties, we propose a universal route of elucidating the sorption mechanism directly from an experimental isotherm, without an isotherm model, based on the statistical thermodynamic fluctuation theory. We have shown that how sorbate-sorbate interaction depends on activity is the key to understanding the sorption mechanism. Without assuming adsorption sites and planar layers, an isotherm can be derived, which contains the Langmuir, BET, and GAB models as its special cases. We have constructed a universal approach applicable to adsorption and absorption, solid and liquid sorbents, and vapor and liquid sorbates and demonstrated its efficacy using the humidity sorption isotherm of sucrose from both the solid and liquid sides.

Cooperative Sorption on Porous Materials.

The functional shape of a sorption isotherm is determined by underlying molecular interactions. However, doubts have been raised on whether the sorption mechanism can be understood in principle from analyzing sorption curves via a range of competing models. We have shown recently that it is possible to translate a sorption isotherm to the underlying molecular interactions via rigorous statistical thermodynamics. The aim of this paper is to fill the gap between the statistical thermodynamic theory and analyzing experimental sorption isotherms, especially of microporous and mesoporous materials. Based on a statistical thermodynamic approach to interfaces, we have derived a cooperative isotherm, as a generalization of the Hill isotherm and our cooperative solubilization model, without the need for assumptions on adsorption sites, layers, and pore geometry. Instead, the statistical characterization of sorbates, such as the sorbate-interface distribution function and the sorbate number distribution, as well as the existence of statistically independent units of the interface, underlies the cooperative sorption isotherm. Our isotherm can be applied directly to literature data to reveal a few key system attributes that control the isotherm: the cooperative number of sorbates and the free energy of transferring sorbates from the saturated vapor to the interface. The sorbate-sorbate interaction is quantified also via the Kirkwood-Buff integral and the excess numbers.

Temperature Dependence of Sorption.

Understanding how sorption depends on temperature on a molecular basis has been made difficult by the coexistence of isotherm models, each assuming a different sorption mechanism and the routine application of planar, multilayer sorption models (such as Brunauer-Emmett-Teller (BET) and Guggenheim-Anderson-de Boer (GAB)) beyond their premises. Furthermore, a common observation that adsorption isotherms measured at different temperatures fall onto a single "characteristic curve" when plotted against the adsorption potential has not been given a clear explanation, due to its ambiguous foundation. Extending our recent statistical thermodynamic fluctuation theory of sorption, we have generalized the classical isosteric theory of sorption into a statistical thermodynamic fluctuation theory and clarified how sorption depends on temperature. We have shown that a characteristic curve exists when sorbate number increment contributes purely energetically to the interface, whereas the correlation between sorbate number and entropy drives the temperature dependence of an isotherm. This theory rationalizes the opposite temperature dependence of water vapor sorption on activated carbons with uniform versus broad pore size distributions and can be applied to moisture sorption on starch gels. The adsorption potential is a convenient variable for sorption in its ability to unify sorbate-sorbate fluctuation and the isosteric thermodynamics of sorption.

Cooperativity in micellar solubilization.

Sudden onset of solubilization is observed widely around or below the critical micelle concentration (CMC) of surfactants. It has also been reported that micellization is induced by the solutes even below CMC and the solubilized solute increases the aggregation number of the surfactant. These observations suggest enhanced cooperativity in micellization upon solubilization. Recently, we have developed a rigorous statistical thermodynamic theory of cooperative solubilization. Its application to hydrotropy revealed the mechanism of cooperative hydrotropy: hydrotrope self-association enhanced by solutes. Here we generalize our previous cooperative solubilization theory to surfactants. We have shown that the well-known experimental observations, such as the reduction of CMC in the presence of the solutes and the increase of aggregation number, are the manifestations of cooperative solubilization. Thus, the surfactant self-association enhanced by a solute is the driving force of cooperativity and a part of a universal cooperative solubilization mechanism common to hydrotropes and surfactants at low concentrations.

Fluctuation adsorption theory: quantifying adsorbate-adsorbate interaction and interfacial phase transition from an isotherm.

How adsorbate-adsorbate interaction determines the functional shape of an adsorption isotherm is an important and challenging question. Many models for the adsorption isotherm have been proposed to answer this question. However, a successful fitting of an isotherm on its own is insufficient to prove the correctness of the model assumptions. Instead, starting from the principles of statistical thermodynamics, we propose how adsorbate-adsorbate interactions can be quantified from an isotherm. This was made possible by extending the key tools of solution statistical thermodynamics to adsorbates at the interface, namely, the Kirkwood-Buff and McMillan-Mayer theories, as well as their relationship to the thermodynamic phase stability condition. When capillary condensation and interfacial phase transition are absent, adsorbate-adsorbate interactions can be quantified from an isotherm using the Kirkwood-Buff integrals, and virial coefficients can yield multiple-body interaction between adsorbates. Such quantities can be obtained directly from the fitting parameters for the well-known isotherm models (e.g., Langmuir, BET). The size of the adsorbate cluster involved in capillary condensation and interfacial phase transition can also be evaluated from the isotherm, which was demonstrated for the adsorption isotherm of water on activated carbons of varying pore sizes from the literature. Signatures of isotherm classifications by IUPAC have been characterized in terms of multiple-body interactions between adsorbates.

Salt-induced LCST-type thermal gelation of methylcellulose: quantifying non-specific interactions via fluctuation theory.

What drives the phase separation of water-soluble polymers in the presence of electrolytes was quantified on a molecular scale via statistical thermodynamic fluctuation theory. Quantifying polymer-water and polymer-salt interactions enabled us to identify the dominant interaction for phase separation. As a model system, the lower critical solution temperature (LCST) type thermal gelation of methylcellulose (MC) in aqueous salt solutions was chosen. The Kirkwood-Buff integrals for intermolecular interactions, calculated from the published calorimetric and volumetric data, showed that (1) the accumulation of salts around MC molecules (favourable interaction between salts and MC) inhibits thermal gelation and the depletion of salts from MC (unfavourable interaction between salts and MC) promotes the gelation, and (2) this salt-MC interaction is the dominant factor (50-100 times stronger than the water-MC interaction). This insight from the fluctuation theory is at odds with the age-old consensus regarding the driving force of thermal gelation: water structure change in the presence of salts induces the promotion or inhibition of thermal gelation. However, our conclusion is founded upon the ability of the fluctuation theory to quantify water-MC and salt-MC interaction independently via the Kirkwood-Buff integrals. Flory-Huggins (FH) theory, on the contrary, could not separate these two interactions owing to the lack of a thermodynamic degree of freedom because the lattice solution is assumed to be fully packed. In addition, the dominant contribution from salt depletion poses difficulty for the χ parameter, which is essentially the difference of contact energies. Our approach, requiring calorimetric and volumetric data alone as input, provides a simple and versatile method towards elucidating the effect of cosolvents on biopolymer phase separation of physiological importance.

Water Networks Can Determine the Affinity of Ligand Binding to Proteins.

Solvent organization is a key but underexploited contributor to the thermodynamics of protein-ligand recognition, with implications for ligand discovery, drug resistance, and protein engineering. Here, we explore the contribution of solvent to ligand binding in the Haemophilus influenzae virulence protein SiaP. By introducing a single mutation without direct ligand contacts, we observed a >1000-fold change in sialic acid binding affinity. Crystallographic and calorimetric data of wild-type and mutant SiaP showed that this change results from an enthalpically unfavorable perturbation of the solvent network. This disruption is reflected by changes in the normalized atomic displacement parameters of crystallographic water molecules. In SiaP's enclosed cavity, relative differences in water-network dynamics serve as a simple predictor of changes in the free energy of binding upon changing protein, ligand, or both. This suggests that solvent structure is an evolutionary constraint on protein sequence that contributes to ligand affinity and selectivity.

Statistical Thermodynamics Unveils How Ions Influence an Aqueous Diels-Alder Reaction.

The kinetics of Diels-Alder (DA) reactions in water has been known to be altered by salts for a long time. Yet the question how salts influence the reaction rate, either as rate-enhancing or rate-reducing additives, has so far remained unresolved. Conflicting hypotheses involve (i) indirect salt contributions through the modulation of internal pressure and (ii) making (or breaking) of the so-called "water-structure" by salts that strengthen (or weaken) the hydrophobic effect. In contrast to the qualitative nature of these hypotheses, here we answer this question quantitatively through a combination of transition state theory and fluctuation adsorption-solvation theory (FAST) using the DA reaction between anthracene-9-carbinol and N-ethylmaleimide as an example. We show that rate enhancement is driven by the salting out of the hydrophobic reactant, while rate-enhancing salts exhibit stronger affinity to the transition state.

Quantifying non-specific interactions via liquid chromatography.

Determinations of solute-cosolute interactions from chromatography have often resulted in problems, such as the "antibinding" (or a negative binding constant) between the solute and micelle in micellar liquid chromatography (MLC) or indeterminacy of salt-ligand binding strength in high-performance affinity chromatography (HPAC). This shows that the stoichiometric binding models adopted in many chromatographic analyses cannot capture the non-specific nature of solvation interactions. In contrast, an approach using statistical thermodynamics handles these complexities without such problems and directly links chromatographic data to, for example, solubility data via a universal framework based on Kirkwood-Buff integrals (KBI) of the radial distribution functions. The chromatographic measurements can now be interpreted within this universal theoretical framework that has been used to rationalize small solute solubility, biomolecular stability, binding, aggregation and gelation. In particular, KBI analysis identifies key solute-cosolute interactions, including excluded volume effects. We present (i) how KBI can be obtained directly from the cosolute concentration dependence of the distribution coefficient, (ii) how the classical binding model, when used solely as a fitting model, can yield the KBIs directly from the literature data, and (iii) how chromatography and solubility measurements can be compared in the unified theoretical framework provided via KBIs without any arbitrary assumptions about the stationary phase. To perform our own analyses on multiple datasets we have used an "app". To aid readers' understanding and to allow analyses of their own datasets, the app is provided with many datasets and is freely available on-line as an open-source resource.