Tributylamine (TBuA) and triethylamine (TEtA) are the most commonly used ion pair reagents in ion pair chromatography especially for the analysis of oligonucleotides. In order to improve the understanding of the retention and separation mechanism of oligonucleotides in ion pair chromatography, it is important to understand the retention mechanism and the nature of interaction of these ion pair reagents with the stationary phase in the chromatographic column. Adsorption isotherm is helpful in evaluating such interactions, and subsequently predicting the retention mechanism. Alkylamines are very polar molecules which lack suitable chromophore and are commonly present in charged forms. Therefore, their determination and the subsequent acquisition of their adsorption isotherms using traditional liquid chromatography is very difficult. In this study, we first developed an analytical method for the determination of TBuA and TEtA in a typical chromatographic mobile phase (acetonitrile-water) and then used the same method to acquire the adsorption isotherms for tributylammonium acetate (TBuAA) and triethylammonium acetate (TEtAA). This method started with the conversion of the alkylammonium ions to free neutral forms by treating the sample with a strong base, followed by pentane-mediated extraction and finally the analysis of the extracts using gas chromatography-flame ionization detector (GC-FID). This three-step method was validated for parameters like range, linearity, intra-day and inter-day precision and accuracy, limit of detection and limit of quantitation. For the adsorption isotherms, the C18 column was first equilibrated with the solutions having different concentrations of alkylammonium ions and then stripped with eluent devoid of alkylammonium ions. Several stripping eluents were investigated and it was discovered that the eluent requirement could be decreased by the addition of sodium chloride. The effluents from the stripping phase were collected and analyzed using the developed analytical method to acquire the adsorption data. Under the investigated conditions, adsorption of TBuAA and TEtAA showed type III and type I isotherm behavior respectively. 

This study investigated the influence of pH on the retention of solutes using a mixed-mode column with carboxyl (-COOH) groups acting as weak cation exchanger bonded to the terminal of C18 ligands (C18-WCX column) and a traditional reversed-phase C18 column. First, a model based on electrostatic theory was derived and successfully used to predict the retention of charged solutes (charged, and ionizable) as a function of mobile phase pH on a C18-WCX column. While the Horváth model predicts the pH-dependent retention of ionizable solutes in reversed-phase liquid chromatography (RPLC) solely based on solute ionization, the developed model incorporates the concept of surface potential generated on the surface of the stationary phase and its variation with pH. To comprehensively understand the adsorption process, adsorption isotherms for these solutes were individually acquired on the C18-WCX and reversed-phase C18 columns. The adsorption isotherms followed the Langmuir model for the uncharged solute and the electrostatically modified Langmuir model for charged solutes. The elution profiles for the single components were calculated from these isotherms using the equilibrium dispersion column model and were found to be in close agreement with the experimental elution profiles. To enable modelling of two-component cases involving charged solute(s), a competitive adsorption isotherm model based on electrostatic theory was derived. This model was later successfully used to calculate the elution profiles of two components for scenarios involving (a) a C18 Column: two charged solutes, (b) a C18 Column: one charged and one uncharged solute, and (c) a C18-WCX Column: two charged solutes. The strong alignment between the experimental and calculated elution profiles in all three scenarios validated the developed competitive adsorption model. 

This work introduces the Adsorption Energy Distribution (AED) calculation using competitive adsorption isotherm data, enabling investigation of the simultaneous AED of two components for the first time. The AED provides crucial insights by visualizing competitive adsorption processes, offering an alternative adsorption isotherm model without prior assuming adsorption heterogeneity, and assisting in model selection for more accurate retention mechanistic modeling. The competitive AED enhances our understanding of multicomponent interactions on stationary phases, which is crucial for understanding how analytes compete on the stationary phase surface and for selecting adsorption models for numerical optimization of preparative chromatography. Here, the two-component AED was tested on both synthetic and experimental data, and a very successful outcome was achieved.

In liquid chromatography (LC), adsorption heterogeneity arises from the distribution of adsorption sites on stationary phases with varying interaction energies, affecting retention and separation performance. This heterogeneity can cause peak tailing, reduced resolution, and unpredictable retention times in analytical chromatography, as well as broad, asymmetric elution profiles in preparative systems. Adsorption heterogeneity depends on the combined effects of the stationary phase, the mobile phase composition, the analyte properties, and the chromatographic conditions. Traditional adsorption isotherms often fail to fully describe these complex interactions because they assume uniform adsorption energies. The Adsorption Energy Distribution (AED) framework offers a powerful alternative by modelling adsorption as a sum of independent homogeneous sites, each with a specific energy, offering a realistic representation of heterogeneous adsorption. This review introduces the theoretical foundations of AED, including its mathematical formulation and computational approaches, and discusses its application in interpreting retention mechanisms in LC. AED analysis is illustrated through its use in both chiral and achiral separations, as well as its ability to explain peak tailing and surface heterogeneity. Practical considerations, such as the range of concentration data in the adsorption isotherm, the selection of a suitable kernel function, and the number of iterations and grid points in AED analysis, are discussed. Special emphasis is given on how to visualize and interpret the AED. This review aims to provide chromatographers with a comprehensive understanding of AED, emphasizing its practical value in characterizing the chromatographic system and elucidating retention mechanisms in liquid chromatography.