Estimation of single solute adsorption isotherms applying the nonlinear frequency response method using non-optimal frequencies

Abstract

In order to estimate single solute adsorption isotherms, the nonlinear frequency response (FR) of a chromatographic column is analyzed experimentally and evaluated using the concept of higher order frequency response functions (FRFs) based on the Volterra series and generalized Fourier transform. In this case study, it has been investigated the adsorption of ethyl benzoate on octadecyl silica from a mixture of methanol and water (60:40) as a solvent. Experiments are performed using a standard gradient HPLC unit. For estimation of adsorption isotherms by the nonlinear FR method the column inlet concentration is changed in a nearly sine waveform around several steady-state concentrations. Using this method the first three local derivatives of a single solute adsorption isotherm are estimated from the low frequency asymptotes of the corresponding functions, i.e. the phase and first order derivative of the FRFs. For an accurate estimation of isotherm coefficients periodical experiments should be preformed for frequencies below a certain critical frequency. This is the frequency needed for approaching the low frequency asymptotic behaviour of the corresponding functions close enough, so that errors due to the non-feasibility of experiments with zero frequency can be neglected. Unfortunately, depending on the properties of the system, it can happen (as for the system investigated here) that experiments for the critical frequency would be too long and cannot be realized. In order to study the loss of accuracy of the nonlinear FR method, when it is applied for non-optimal frequencies, experiments are performed for frequencies approximately one order of magnitude higher than the critical frequency required to evaluate the FRF phases. The obtained isotherm model coefficients are compared with the ones estimated using conventional frontal analysis as a reference method. The isotherms determined by two methods are similar, however a closer look reveals that peaks predicted under overloading conditions differ.