Determination of the single component and competitive adsorption isotherms of the 1-indanol enantiomers by the inverse method

The inverse method of isotherm determination consists in calculating the numerical values of the coefficients of an isotherm model that give a set of chromatographic profiles in best possible agreement with the set of experimental profiles available. This method was applied to determine the adsorption isotherms of the 1-indanol enantiomers on a cellulose tribenzoate chiral stationary phase. Both single-component and competitive isotherms were determined by using no more than one or two overloaded band profiles. The isotherms determined from the overloaded band profiles agreed extremely well with the isotherms determined by frontal analysis. Several isotherm models were used and tested. The best-fit isotherm was selected by means of statistical evaluation of the results. The results show that the adsorption is best characterized with a model describing heterogeneous adsorption with bimodal adsorption energy distribution.     

Retention of ionizable compounds in high-performance liquid chromatography. 14. Acid-base pK values in acetonitrile-water mobile phases

Using competitive frontal analysis, the binary adsorption isotherms of the enantiomers of 1-phenyl-l-propanol were measured on a microbore column packed with a chiral stationary phase based on cellulose tribenzoate. These measurements were carried out using only the racemic mixture. The experimental data were fitted to four different isotherm models: Langmuir, BiLangmuir, Langmuir-Freundlich and Tóth. The BiLangmuir and the Langmuir-Freundlich models accounted best for the competitive adsorption data. An excellent agreement between the experimental and the calculated overloaded band profiles for various samples of racemic mixture was obtained when the equilibrium dispersive model of chromatography was used together with the BiLangmuir competitive isotherm. The isotherm parameters measured under competitive conditions were used to calculate the overloaded band profiles of large samples of the pure S- and R-enantiomers, too. A satisfactory agreement between the experimental and calculated band profiles was observed when using in the computation the corresponding single component BiLangmuir isotherm derived from the binary isotherm previously determined. Thus oniy data derived from the racemic mixture are required for computer optimization of the preparative chromatography separation of the enantiomers.     

Adsorption behavior and prediction of the band profiles of the enantiomers of 3-chloro-1-phenyl-1-propanol. Influence of the mass transfer kinetics

The single-component and competitive adsorption isotherms of the enantiomers of 3-chloro-1-phenyl-1-propanol were measured by frontal analysis. The stationary phase was a cellulose tribenzoate coated on silica, the mobile phase an n-hexane-ethyl acetate (95:5) solution. The adsorption data measured fitted well to the Langmuir isotherm model. The band profiles of single components and of their mixtures were calculated using the equilibrium-dispersive model. These profiles were found to match quite satisfactorily the experimental band profiles. However, the agreement between calculated and experimental band profiles was significantly improved when a more complex model taking into account the mass transfer kinetics was used. The mass transfer rate coefficients, k(f), for both single components were determined by using the transport-dispersive model of chromatography. The coefficients obtained were used to predict the band profiles of mixtures of the two enantiomers to good agreement.     

Determination of competitive isotherms of enantiomers by a hybrid inverse method using overloaded band profiles and the periodic state of the simulated moving-bed process

A procedure for determination of adsorption isotherms in simulated moving-bed (SMB) chromatography is presented. The parameters of a prescribed adsorption isotherm model and rate constants are derived using a hybrid inverse method, which incorporates overloaded band profiles of the racemic mixture and breakthrough data from a single frontal experiment. The latter are included to reduce the uncertainty on the estimated saturation capacity, due to the dilution of the chromatograms with respect to the injected concentrations. The adsorption isotherm model is coupled with an axially dispersed flow model with finite mass-transfer rate to describe the experimental band profiles. The numerical constants of the isotherm model are tuned so that the calculated and measured band profiles match as much as possible. The accuracy of the isotherm model is then checked against the cyclic steady state (CSS) of the target SMB process, which is readily and cheaply obtained experimentally on a single-column set-up. This experiment is as expensive and time consuming as just a few breakthrough experiments. If necessary, the isotherm parameters are adjusted by applying the inverse method to the experimental CSS concentration profile. The method is successfully applied to determine the adsorption isotherms of Tr?gers base enantiomers on Chiralpak AD/methanol system. The results indicate that the proposed inverse method offers a reliable and quick approach to determine the competitive adsorption isotherms for a specific SMB separation.     

Modeling of the separation of two enantiomers using a microbore column

A microbore column packed with Chiralcel OB (cellulose tribenzoate coated silica) was used for the measurement of the single and competitive equilibrium-isotherm data of the 1-indanol enantiomers by frontal analysis. The amount of sample needed for the isotherm data acquisition was about 20 times less than that required with a conventional column. The data obtained were fitted to different single and competitive isotherm models. Both the single and the competitive data sets fitted best to the same Bilangmuir (BL) isotherm model with small differences in the numerical values of the parameters. The best fitted Bilangmuir single and competitive isotherm models were used to predict the overloaded experimental profiles of both pure enantiomers, of the racemic mixture, and of different enantiomeric mixtures. All the calculated profiles were in excellent agreement with the experimental ones. This agreement confirms that in many chiral separations, the competitive isotherms can be derived from data acquired from the mere racemic mixture with a sufficient accuracy for a correct prediction of the band profiles of all kinds of enantiomer mixtures, making possible the computer-assisted optimization of the experimental conditions.     

Numerical determination of competitive adsorption isotherm of mandelic acid enantiomers on cellulose-based chiral stationary phase

The use of inverse method for the determination of competitive adsorption isotherm of mandelic acid enantiomers on cellulose tris(3,5-diethylphenyl carbamate) stationary phase is proposed in this work. Non-dominated sorting genetic algorithm with jumping genes (NSGA-II-JG) was applied to acquire the isotherm parameters by minimizing the sum of square deviations of the model predictions from the measured elution profiles. Three different competitive isotherm models, i.e., Langmuir, biLangmuir and Tóth, combined with transport-dispersive chromatographic model were used in predicting the elution profiles. Orthogonal collocation on finite element (OCFE) method was applied to obtain the calculated elution profiles. Results indicate that biLangmuir isotherm and Tóth isotherm give remarkably similar equilibrium isotherms within the investigated liquid concentration range. Band profiles calculated from both isotherm models are in good agreement with the experimental data. The validity of the determined parameters was verified by comparing the model predictions with experimental elution profiles at various experimental conditions.     

Modeling of the adsorption behavior and the chromatographic band profiles of enantiomers : Behavior of methyl mandelate on immobilized cellulose

The adsorption isotherms of (−)- and (+)-methyl mandelate from a hexane-isopropanol (90:10) solution were measured on a chromatographic column packed with 4-methylcellulose tribenzoate coated on silica. These isotherms are accounted for by a bi-Langmuir isotherm model, the two Langmuir terms having widely different initial slopes and saturation capacities, but each term having the same saturation capacity for the two enantiomers. The competitive isotherms were also measured. They are in excellent agreement with the prediction of a competitive bi-Langmuir model based on the single-component isotherms. The individual band profiles are in agreement with the profiles calculated from these isotherms. Thus, a simplified competitive isotherm can be used to model a separation on a chiral stationary phase the recognition mechanism of which is not well identified and the adsorption behavior of which is certainly not ideal.     

Influence of the column hold-up time measurement accuracy on the prediction of chromatographic band profiles

The influence of the column hold-up time measurement accuracy on the determination of equilibrium isotherms by classical frontal analysis and the prediction of overloaded elution band profiles were investigated. The ideal model of chromatography in combination with a Langmuir isotherm was used. Breakthrough curves and overloaded elution profiles were computer generated with a known hold-up time value (true hold-up time). Then these data were evaluated the same way as it is done with experimental chromatographic data where the true hold-up time is unknown, i.e. to determine the equilibrium isotherm by the frontal analysis procedure, to fit the isotherm data to the Langmuir model and then to predict chromatographic band profiles using, e.g. the ideal model of chromatography. A comparison of overloaded elution profiles obtained with different deviations of the hold-up time from its true value shows that the effect of its measurement error is significant in preparative liquid chromatography because the isotherm is usually strongly nonlinear in this case.     

Thermodynamic analysis of the heterogenous binding sites of molecularly imprinted polymers

The thermodynamic interactions of two polymers, one Fmoc-L-Trp-imprinted (MIP), the other one an unimprinted reference (NIP), with the two Fmoc-tryptophan enantiomers were studied by frontal analysis, which allows accurate measurements of the adsorption isotherms. These isotherms were acquired at temperatures of 40, 50, 60, and 70 degrees C, for sample concentrations ranging between 0.005 and 40 mM. The mobile phase used was acetonitrile with one percent acetic acid as an organic modifier. Within the measured concentration ranges, the tri-Langmuir isotherm model accounts best for the isotherm data of both enantiomers on the MIP, the bi-Langmuir model for the isotherm data of Fmoc-L-Trp on the NIP. These isotherm models were selected using three independent processes: statistical tests on the results from regression of the isotherm data to different isotherm models; calculation of the affinity energy distribution from the raw isotherm data; comparison of the experimental and the calculated band profiles. The isotherm parameters obtained from these best selected isotherm models showed that the enantiomeric selectivity does not change significantly with temperature, while the affinity of the substrates for both the MIP and the NIP decrease considerably with increasing temperatures. These temperature effects on the binding performance of the MIP were clarified by considering the thermodynamic functions (i.e., the standard molar Gibbs free energy, the standard molar entropy of adsorption, and the standard molar enthalpy of adsorption) for each identified type of adsorption sites, derived from the Van't Hoff equation. This showed that the entropy of transfer of Fmoc-L-Trp from the mobile to the MIP stationary phase is the dominant driving force for the selective adsorption of Fmoc-L-Trp onto the enantioselective binding sites. This entropy does not change significantly with increasing temperatures from 40 to 70 degrees C.     

Study of the adsorption behavior of the enantiomers of 1-phenyl-1-propanol on a cellulose-based chiral stationary phase

Using single-step frontal analysis, we measured single-component and competitive adsorption isotherm data for the two enantiomers of 1-phenyl-1-propanol (PP). These experimental data were fitted to several competitive bi-Langmuir models (with 8, 6, 5 and 4 parameters) and to the competitive Langmuir model. The latter model accounted well for the behavior of both PP enantiomers on Chiracel OB (cellulose tribenzoate coated on silica gel). The parameters obtained were used in numerical calculations to predict the band profiles of the two single components and of their mixtures under overloaded conditions. The equilibrium-dispersive model provides satisfactory results, with minor differences between the calculated and the experimental profiles. These differences became negligible when a more complex kinetic model was used, with a concentration-dependent rate coefficient.     

Adsorption of the enantiomers of 3-chloro-1-phenyl-propanol on silica-bonded chiral quinidine carbamate

The interactions of 3-chloro-1-phenyl-propanol with a quinidine carbamate-bonded chiral stationary phase under NPLC conditions were studied by measuring the adsorption isotherm data of its enantiomers by frontal analysis, modeling these data with a suitable isotherm model, and comparing the experimental overloaded elution band profiles with those calculated with this isotherm and the equilibrium dispersive model of liquid chromatography. The affinity energy distribution was calculated from the adsorption isotherm data. The results show that the surface of the adsorbent is heterogeneous and exhibits a bimodal adsorption energy distribution. This fact is interpreted in terms of the presence of two different types of adsorption sites on the stationary phase, nonselective and enantioselective sites. Albeit the bi-Langmuir isotherm model successfully accounts for the single-component data corresponding to both enantiomers, the competitive bi-Langmuir isotherm model does not allow an accurate prediction of the overloaded band profiles of the racemic mixture. Thermodynamic data are drawn for explanation. Some aspects of the retention mechanism are discussed in the light of the data obtained.     

Repeatability and reproducibility of high-concentration data in reversed-phase liquid chromatography. III. Isotherm reproducibility on Kromasil C18

Single component equilibrium isotherms of six compounds (aniline, caffeine, ethylbenzene, phenol, propranolol, and theophylline) were determined by the inverse method on 10 Kromasil-C18 columns, using water-methanol solutions as the mobile phase. This method offers an economic and fast isotherm determination by means of the overloaded band profiles of the compounds. Five out of the ten columns used in this test come from the same batch whilst the other five columns represent five additional batches. Statistical evaluation was used to assess the reproducibility of the isotherm parameters. We found that the column-to-column reproducibility of the isotherm parameters is of the same magnitude as the batch-to-batch reproducibility (with the exception of one outlier column). In most of the cases, the reproducibilities of the saturation capacities and that of the retention factors are excellent, they are typically between 1.2 and 3%, and very often below 2%. Within the limits of the experimental precision, these results agree with those obtained earlier, using a conventional method of isotherm determination.     

Repeatability and reproducibility of high concentration data in reversed-phase liquid chromatography. I. Overloaded band profiles on Kromasil-C18

Single-component adsorption-isotherm data were acquired by frontal analysis (FA) for six low-molecular-mass compounds (phenol, aniline, caffeine, theophylline, ethylbenzene and propranolol) on one Kromasil-C18 column, using water-methanol solutions (between 70:30 and 20:80, v/v) as the mobile phase. Propranolol data were also acquired using an acetate buffer (0.2 M) instead of water. The data were modeled for best agreement between calculated and experimental overloaded band profiles. The adsorption energy distribution was also derived and used for the selection of the best isotherm model. Widely different isotherm models were found to model best the data obtained for these compounds, convex upward (i.e. Langmuirian), convex downward (i.e. anti-Langmuirian), and S-shaped isotherms. Using the same sample size for all columns (loading factor, Lf approximately 10%), overloaded band profiles were recorded on four different columns packed with the same batch of Kromasil-C18 and five other columns packed with different batches of Kromasil-C18. These experimental band profiles were compared to the profile calculated from the isotherm measured by FA on the first column. The repeatability as well as the column-to-column and the batch-to-batch reproducibilities of the band profiles are better than 4%.     

Determination of single component isotherms and affinity energy distribution by chromatography

Adsorption isotherm data were acquired by frontal analysis (FA) and large sample-size band profiles were recorded for phenol and caffeine. For both compounds, the isotherm data fit well to the Langmuir, Toth, and Bi-Langmuir models of adsorption. The Langmuir model must be dismissed because it does not predict accurately the overloaded band profiles. However, profiles calculated using the unimodal Toth and the bimodal Bi-Langmuir models are indistinguishable. The expectation-maximization procedure was used to calculate directly the affinity energy distribution (AED) from the raw FA data points. For both compounds, the AED converges to a bimodal distribution at high numbers of iterations. This result, which shows the high sensitivity of the EM method, suggest that the Bi-Langmuir model makes better physical sense than the Toth model. This model also permits a detailed investigation of the properties of active sites, a feature often evoked in chromatography but so far rarely the topic of a quantitative investigation.     

Determination of competitive adsorption isotherm parameters of pindolol enantiomers on alpha1-acid glycoprotein chiral stationary phase

In this paper, inverse method (IM) was used to determine the binary competitive adsorption isotherm of pindolol enantiomers by a least-square fitting of the proposed model to the experimentally measured elution curves of racemic pindolol. The isotherm parameters were determined by minimizing the least-square error using an adaptation of genetic algorithm, non-dominated sorting genetic algorithm with jumping genes (NSGA-II-JG). An equilibrium dispersive (ED) model combined with bi-Langmuir isotherm was used in predicting the elution profiles. The determined parameters show good agreement with the experimental profiles at various experimental conditions such as sample volume, concentration and flow rates of the racemic mixture. Robustness and validity of the isotherm parameters were also verified by frontal analyses at various step inputs. Results from both the pulse tests and the frontal analysis indicate that adsorption isotherm derived from the inverse method is quite reliable. This method requires relatively less number of experiments to be performed and therefore, lower experimental costs confirming that inverse method is an attractive alternative approach of experimental technique in determining the competitive adsorption isotherm for binary systems.     

Features of the adsorption of naproxen enantiomers on weak chiral anion-exchangers in nonlinear chromatography

The retention mechanism of the enantiomers of naproxen on a Pirkle-type chiral stationary phase (CSP) was studied. This CSP is made of a porous silica grafted with quinidine carbamate. It can interact with the weak organic electrolyte naproxen either by adsorbing it or by ion-exchange. Using frontal chromatography, we explored the adsorption equilibrium under such experimental conditions that naproxen dissociates or cannot dissociate. Under conditions preventing ionic dissociation, the adsorption isotherms were measured, the adsorption energy distributions determined, and the chromatographic profiles calculated. Three different types of the adsorption sites were found for both enantiomers. The density and the binding energy of these sites depend on the nature of the organic modifier. Different solute species, anions, neutral molecules, solvent-ion associates, and solute dimers can coexist in solution, giving rise to different forms of adsorption. This study showed the unexpected occurrence of secondary steps in the breakthrough profiles of S-naproxen in the adsorption mode at high concentrations. Being enantioselective, this phenomenon was assumed to result from the association of solute molecules involving a chiral selector moiety. A multisite Langmuir adsorption model was used to calculate band profiles. Although this model accounts excellently for the experimental adsorption isotherms, it does not explain all the features of the breakthrough profiles. A comparison between the calculated and experimental profiles allowed useful conclusions concerning the effects of the adsorbate-adsorbate and adsorbate-solvent interactions on the adsorption mechanism.     

New theoretical expressions for the five adsorption type isotherms classified by BET based on statistical physics treatment

New theoretical expressions to model the five adsorption isotherm types have been established. Using the grand canonical ensemble in statistical physics, we give an analytical expression to each of five physical adsorption isotherm types classified by Brunauer, Emett, and Teller, often called BET isotherms. The establishment of these expressions is based on statistical physics and theoretical considerations. This method allowed estimation of all the mathematical parameters in the models. The physicochemical parameters intervening in the adsorption process that the models present could be deduced directly from the experimental adsorption isotherms by numerical simulation. We determine the adequate model for each type of isotherm, which fixes by direct numerical simulation the monolayer, multilayer, or condensation character. New equations are discussed and results obtained are verified for experimental data from the literature. The new theoretical expressions that we have proposed, based on statistical physics treatment, are rather powerful to better understand and interpret the various five physical adsorption type isotherms at a microscopic level.     

Band splitting in overloaded isocratic elution chromatography II. New competitive adsorption isotherms

The equations of two new binary competitive isotherms models are derived. The first of these models assumes that the isotherms of the two pure, single compounds have distinct monolayer capacities. Its derivation is based on kinetic arguments. The ideal adsorbed solution (IAS) framework was applied to derive the second model that is a thermodynamically consistent competitive isotherm. This second model predicts the competitive adsorption isotherm behavior of a mixture of two compounds that have single-component adsorption behavior following a BET and/or a Langmuir isotherms. Both models apply well to the binary adsorption of ethylbenzoate and 4-tert.-butylphenol on a Kromasil-C18 column (with methanol-water, 62:38, v/v, as the mobile phase). The best single-solute adsorption isotherms of these two compounds are the liquid-solid extended multilayer BET and the Langmuir isotherms, respectively. The kinetic and thermodynamic new competitive models were compared, regarding the accuracy of their prediction of the elution band profiles of mixtures of these two compounds. A better agreement between experimental and calculated profiles was observed with the kinetic model. The IAS model failed because the behavior of the ethylbenzoate/4-tert.-butylphenol adsorbed phase mixture is probably non-ideal. The most striking result is the qualitative prediction by these models of the peak splitting of 4-tert.-butylphenol during its elution in presence of ethylbenzoate.     

Effect of the mobile phase composition on the adsorption behavior of tryptophan in reversed-phase liquid chromatography

Single-component adsorption isotherm data of l-tryptophan on a C18-bonded silica column were acquired by frontal analysis (FA), with aqueous mobile phases containing 2.5, 5, and 7.5% of acetonitrile (ACN) or 7, 10, 15, and 20% of methanol (MeOH). Most of these isotherms have two inflection points and three different parts. The low and the high concentration parts exhibit langmuirian behavior. The intermediate part exhibits anti-langmuirian behavior. The inflection points shift toward higher concentrations with increasing mobile phase concentration in ACN or MeOH, which causes the differences in the isotherm profiles. The nature of the organic modifier and its concentration affect only the isotherm profile and the numerical values of its parameters, not the nature of the best model, which is the bi-Moreau model in all cases. The isotherm profiles depend on the experimental conditions because they affect the intensity of the adsorbate-adsorbate interactions. Overloaded band profiles of tryptophan were recorded with the seven mobile phase compositions. They were used to determine the best values of the isotherm coefficients by the inverse method (IM) of chromatography. There is an excellent agreement between the values of these parameters obtained by FA and by IM. Increasing the concentration of either ACN or MeOH in the mobile phase causes a slight decrease in the saturation capacities of the low and the high energy sites, and in the adsorption constant of the low energy sites. The adsorption constant of the high energy sites increases with increasing concentration of either solvent or is little affected. The adsorbate-adsorbate interaction constants of both low and high energy sites increase for both solvents. Saturation capacities of the high energy sites are higher for ACN than for MeOH.     

Fitting adsorption isotherms to the distribution data determined using packed micro-columns for high-performance liquid chromatography

Knowing the adsorption isotherms of the components of a mixture on the chromatographic system used to separate them is necessary for a better understanding of the separation process and for the optimization of the production rate and costs in preparative high-performance liquid chromatography (HPLC). Currently, adsorption isotherms are usually measured by frontal analysis, using conventional analytical columns. Unfortunately, this approach requires relatively large quantities of pure compounds, and hence is expensive, especially in the case of pure enantiomers. In this work, we investigated the possible use of packed micro-bore and capillary HPLC columns for the determination of adsorption isotherms of benzophenone, o-cresol and phenol in reversed-phase systems and of the enantiomers of mandelic acid on a Teicoplanin chiral stationary phase. We found a reasonable agreement between the isotherm coefficients of the model compounds determined on micro-columns and on conventional analytical columns packed with the same material. Both frontal analysis and perturbation techniques could be used for this determination. The consumption of pure compounds needed to determine the isotherms decreases proportionally to the second power of the decrease in the column inner diameter, i.e. 10 times for a micro-bore column (1 mm I.D.) and 100 times for capillary columns (0.32 mm I.D.) with respect to 3.3 mm I.D. conventional columns.