Fast estimation of adsorption isotherm parameters in gradient elution preparative liquid chromatography. I: The single component case

•We adapted the inverse method for fast estimation of adsorption parameters under gradient runs.•The adapted inverse method was verified using both simulations and experiments.•For simulated data, almost exact Langmuir adsorption parameters could be determined.•Using experimental data, the inverse method was compared with the perturbation peak method.•The shapes of the bi-Langmuir adsorption isotherms were identical for the two methods but not the numerical parameters. The inverse method is a numeric method for fast estimation of adsorption isotherm parameters directly from overloaded elution profiles. However, it has previously only been used for isocratic experiments. Here we will extend the inverse method so it can be used for gradient elution too. This extended inverse method will make it possible to study the adsorption of substances whose retention factor vary strongly with the mobile-phase composition, like peptides and proteins, where the classic methods will fail. Our extended inverse method was verified using both simulations and real experiments. For simulated overloaded elution profiles we were able to determine almost exact Langmuir adsorption isotherm parameters with the new approach. From real experimental data, bi-Langmuir adsorption parameters were estimated using both the perturbation peak method and the extended inverse method. The shape of the acquired adsorption isotherms did match over the considered concentration range; however, the adsorption isotherm parameters found with the two methods were not the same. This is probably due to the fact that adsorption isotherm estimated with the inverse method is only a good approximation up to the highest eluted concentration in the used chromatograms. But this is not a serious drawback from a process point of view where the main objective is to make accurate predictions of elution profiles. The bi-Langmuir adsorption isotherm obtained with both methods could accurately predict the shape of overloaded elution profiles.