Martin-Synge algorithm for the solution of equilibrium-dispersive model of liquid chromatography

An alternative method, called the Martin-Synge algorithm, is introduced to calculate numerical solutions of the equilibrium-dispersive (ED) model. The developed algorithm is based on the earlier work of Friday and Levan and on the continuous plate model of Martin and Synge. The column is divided evenly into a series of virtual vessels in which a simplified mass balance equation is solved accurately by the Runge-Kutta-Fehlberg method and the elution profile is given by the numerical solution for the last vessel. The dispersion of the compound during the elution process is controlled by adjusting the number of virtual vessels into which the column is divided. Solving the ED model under linear conditions with this method gives exactly the same profile as the analytical solution of the Martin-Synge plate model. The Martin-Synge method gives better results than the Rouchon method (1) when the isotherms involved are sigmoidal or anti-Langmuir; and, more importantly, (2) in the case of multi-component problems. Finally, the Martin-Synge method proves to be more robust and faster than the OCFE method that, until now, was considered to be one of the most robust and accurate algorithms. The developed algorithm was used for the calculation of the coefficients of the isotherm of butyl benzoate by the inverse method, using a simplex optimization algorithm.     

A systematic investigation of algorithm impact in preparative chromatography with experimental verifications

Computer-assisted optimization of chromatographic separations requires finding the numerical solution of the Equilibrium-Dispersive (ED) mass balance equation. Furthermore, the competitive adsorption isotherms needed for optimization are often estimated numerically using the inverse method that also solves the ED equations. This means that the accuracy of the estimated adsorption isotherm parameters explicitly depends on the numerical accuracy of the algorithm that is used to solve the ED equations. The fast and commonly used algorithm for this purpose, the Rouchon Finite Difference (RFD) algorithm, has often been reported not to be able to accurately solve the ED equations for all practical preparative experimental conditions, but its limitations has never been completely and systematically investigated. In this study, we thoroughly investigate three different algorithms used to solve the ED equations: the RFD algorithm, the Orthogonal Collocation on Finite Elements (OCFE) method and a Central Difference Method (CDM) algorithm, both for increased theoretical understanding and for real cases of industrial interest. We identified discrepancies between the conventional RFD algorithm and the more accurate OCFE and CDM algorithms for several conditions, such as low efficiency, increasing number of simulated components and components present at different concentrations. Given high enough efficiency, we experimentally demonstrate good prediction of experimental data of a quaternary separation problem using either algorithm, but better prediction using OCFE/CDM for a binary low efficiency separation problem or separations when the compounds have different efficiency. Our conclusion is to use the RFD algorithm with caution when such conditions are present and that the rule of thumb that the number of theoretical plates should be greater than 1000 for application of the RFD algorithm is underestimated in many cases.     

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.     

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.     

Theoretical study of the accuracy of the elution by characteristic points method for bi-langmuir isotherms

The bi-Langmuir equation has recently been proven essential to describe chiral chromatographic surfaces and we therefore investigated the accuracy of the elution by characteristic points method (ECP) for estimation of bi-Langmuir isotherm parameters. The ECP calculations was done on elution profiles generated by the equilibrium-dispersive model of chromatography for five different sets of bi-Langmuir parameters. The ECP method generates two different errors; (i) the error of the ECP calculated isotherm and (ii) the model error of the fitting to the ECP isotherm. Both errors decreased with increasing column efficiency. Moreover, the model error was strongly affected by the weight of the bi-Langmuir function fitted. For some bi-Langmuir compositions the error of the ECP calculated isotherm is too large even at high column efficiencies. Guidelines will be given on surface types to be avoided and on column efficiencies and loading factors required for adequate parameter estimations with ECP.     

Accurate and rapid estimation of adsorption isotherms in liquid chromatography using the inverse method on plateaus

The inverse method (IM) is an attractive approach for estimating adsorption isotherm parameters in liquid chromatography (LC), mainly due to its experimental simplicity and low sample consumption. This article presents a new experimental approach, the inverse method on plateaus (IMP), which uses elution profiles on concentration plateaus together with IM. This approach enabled us to obtain very accurate adsorption isotherms that agreed well with those estimated by means of frontal analysis over the entire concentration range under consideration. IMP is recommended when accurate adsorption isotherm estimates are required, and standard IM is insufficient.     

Numerical determination of the adsorption isotherms of tryptophan at different temperatures and mobile phase compositions

Single-component adsorption isotherm data were acquired by frontal analysis (FA) for tryptophan on a C(18)-Kromasil packed column, using acetonitrile-water solutions of various compositions (2.5, 5, and 7.5% ACN+1% acetic acid) and at five different temperatures between 25 and 65 degrees C. The adsorption isotherm model accounting best for these data is the bi-Moreau model, showing that two types of adsorption sites coexist on the surface and that strong adsorbate-adsorbate interactions take place. Large concentration band profiles of tryptophan were obtained for the three mobile phase compositions, at five different temperatures and the best values of the adsorption isotherm coefficients were determined by the inverse method (IM) of chromatography. The advantages and drawbacks of using the FA and the IM for determining the coefficients of the adsorption isotherm of tryptophan under the experimental conditions selected are discussed. The results of the FA and IM measurements are in good agreement. Both indicate that the retention time of tryptophan decreases rapidly with increasing acetonitrile concentration in the mobile phase as well as the saturation capacities of the two types of adsorption sites, with the highest values of the two saturation capacities being found for the lowest ACN content and the lowest temperature. The adsorption constant on the low-energy sites decreases with increasing acetonitrile content and temperature. In contrast, the adsorption constant on the high-energy sites increases with increasing ACN content of the mobile phase but decreases with increasing temperature. The solute-solute interaction parameters for the low and the high-energy adsorption sites increase rapidly with increasing ACN concentration in the mobile phase and with increasing temperature.     

Choice of Model for Estimation of Adsorption Isotherm Parameters in Gradient Elution Preparative Liquid Chromatography

The inverse method is a numerical method for fast estimation of adsorption isotherm parameters directly from a few overloaded elution profiles and it was recently extended to adsorption isotherm acquisition in gradient elution conditions. However, the inverse method in gradient elution is cumbersome due to the complex adsorption isotherm models found in gradient elution. In this case, physicochemically correct adsorption models have very long calculation times. The aim of this study is to investigate the possibility of using a less complex adsorption isotherm model, with fewer adjustable parameters, but with preserved/acceptable predictive abilities. We found that equal or better agreement between experimental and predicted elution profiles could be achieved with less complex models. By being able to select a model with fewer adjustable parameters, the calculation times can be reduced by at least a factor of 10.

     

Effect of the ionic strength of salts on retention and overloading behavior of ionizable compounds in reversed-phase liquid chromatography I. XTerra-C18

The influence of the salt concentration (potassium chloride) on the retention and overloading behavior of the propranolol cation (R'-NH2+ -R) on an XTerra-C18 column, in a methanol:water solution, was investigated. The adsorption isotherm data were first determined by frontal analysis (FA) for a mobile phase without salt (25% methanol, v/v). It was shown that the adsorption energy distribution calculated from these raw adsorption data is bimodal and that the isotherm model that best accounts for these data is the bi-Moreau model. Assuming that the addition of a salt into the mobile phase changes the numerical values of the parameters of the isotherm model, not its mathematical form, we used the inverse method (IM) of chromatography to determine the isotherm with seven salt concentrations in the mobile phase (40% methanol, v/v; 0, 0.002, 0.005, 0.01, 0.05, 0.1 and 0.2 M). The saturation capacities of the model increase, q(s,1) by a factor two and q(s,2) by a factor four, with increasing salt concentration in the range studied while the adsorption constant b1 increases four times and b2 decreases four times. Adsorbate-adsorbate interactions vanish in the presence of salt, consistent with results obtained previously on a C18-Kromasil column. Finally, besides the ionic strength of the solution, the size, valence, and nature of the salt ions affect the thermodynamic as well as the mass transfer kinetics of the adsorption mechanism of propranolol on the XTerra column.     

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.     

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.     

Simultaneous determination of two‐component isotherm parameters and lumped mass transfer coefficients in RPLC with the 0‐1 model‐inverse method

The 0‐1 model‐inverse method of nonequilibrium nonlinear chromatography was developed to simultaneously determine the isotherm parameters and the lumped mass transfer coefficients of the two‐component systems in RPLC. By comparing the simulated elution curves with experimental curves with regard to profiles and areas, the suitable isotherm parameters and the lumped mass transfer coefficients were obtained with the 0‐1 model‐inverse method. With a solute cell unit width of cm, the average errors of the peak areas were 0.178% for one component and ?0.40% for two components, and the numerical diffusions of the 0‐1 model for the contribution to band broadening may be negligible. In addition, the results showed that the lumped mass transfer coefficients decrease as the solute concentration increases. The 0‐1 model‐inverse method has not only the advantages of high calculation speed (less than 10 min for one‐component systems or approximately 3 h for two‐component systems using an ordinary computer) and high accuracy in simultaneously obtaining thermodynamic parameters and kinetic parameters of two‐component systems, but this method also possesses the potential to optimally design and control the time‐variant preparative chromatographic system due to the thermodynamic state recursion and the Lagrangian‐Eulerian presentation of the 0‐1 model.     

Numerical determination of the competitive isotherm of enantiomers

A numerical method was developed and used to determine adsorption isotherms in chromatography. The numerical parameters of an isotherm model were derived from the recorded band profiles of the racemic mixture of the 1-phenyl-1-propanol enantiomers, by means of a nonlinear least-squares method. We used the equilibrium-dispersive model of chromatography with several isotherm models. The numerical constants of the isotherm models were tuned so that the calculated and the measured band profiles match as much as possible. We show that this numerical inverse method can be applied even without the knowledge of the individual band profile of the pure enantiomers. The isotherms determined from the--usually unresolved--overloaded band profiles matched extremely well the isotherms determined by frontal analysis. Several isotherm models were used and tested--such as Langmuir, biLangmuir, Tóth, Langmuir-Freundlich. The best-fit isotherm was selected by means of statistical evaluation of the results.     

Modeling of overloaded gradient elution of nociceptin/orphanin FQ in reversed-phase liquid chromatography

The Reversed-phase (RP) gradient elution chromatography of nociceptin/orphanin FQ (N/OFQ), a neuropeptide with many biological effects, has been modeled under linear and non-linear conditions. In order to do this, the chromatographic behavior has been studied under both linear and nonliner conditions under isocratic mode at different mobile phase compositions--ranging from 16 to 19% (v/v) acetonitrile (ACN) in aqueous trifluoracetic acid (TFA) 0.1% (v/v)-on a C-8 column. Although the range of mobile phase compositions investigated was quite narrow, the retention factor of this relatively small polypeptide (N/OFQ is a heptadecapeptide) has been found to change by more than 400%. In these conditions, gradient operation resulted thus to be the optimum approach for non-linear elution. As the available amount of N/OFQ was extremely reduced (only a few milligrams), the adsorption isotherms of the peptide, at the different mobile phase compositions examined, have been measured through the so-called inverse method (IM) on a 5 cm long column. The adsorption data at different mobile phase compositions have been fitted to several models of adsorption. The dependence of the isotherm parameters on the mobile phase composition was modeled by using the linear solvent strength (LSS) model and a generalized Langmuir isotherm that includes the mobile phase composition dependence. The overloaded gradient separation of N/OFQ has been modeled by numerically solving the equilibrium-dispersive (ED) model of chromatography under a selected gradient elution mode, on the basis of the previously determined generalized Langmuir isotherm. The agreement between theoretical calculations and experimental overloaded band profiles appeared reasonably accurate.     

Effect of the flow rate on the measurement of adsorption data by dynamic frontal analysis

The adsorption data of propyl benzoate were acquired by frontal analysis (FA) on a Symmetry-C18 column, using a mixture of methanol (65%, v/v) and water as the mobile phase, at three different flow rates, 0.5, 1.0 and 2.0 mL/min. The exact flow rates Fv were measured by collecting the mobile phase in volumetric glasses (deltaFv / Fv < or = 0.2%). The extra-column volumes and the column hold-up volume were accurately measured at each flow rate by tracer injections. The detailed effect of the flow rate on the value of the amount adsorbed was investigated. The best isotherm model accounting for the adsorption data was the same BET isotherm model at all three flow rates. Only slight differences (always less than 5%) were found between the three different sets of isotherm parameters (saturation capacity, q(s), equilibrium constant on the adsorbent, b(s) and equilibrium constant on successive layers of propyl benzoate, bL). The reproducibility of the same isotherm parameters measured by the inverse method (IM) is less satisfactory, leading to R.S.D.s of up to 10%. A flow rate increase is systematically accompanied by a slight increase of the amount adsorbed. This phenomenon is consistent with the influence of the pressure on the equilibrium constant of adsorption due to the difference between the partial molar volumes of the solute and the adsorbate. The larger average pressure along the column that is required to achieve a larger flow rate causes a larger amount of solute to be adsorbed on the column at equilibrium. This result comforts the high sensitivity and versatility of the FA method for isotherm determination under any kind of situation.     

Numerical approximation of a two-dimensional model of non-isothermal reactive liquid chromatography involving slow rates of adsorption–desorption kinetics

Abstract

A two-dimensional general rate model of non-isothermal reactive column chromatography is formulated considering homogenous and heterogeneous reaction rates, slow rates of adsorption–desorption kinetics, and enthalpies of adsorption and reaction. The model is expressed by a system of six nonlinear partial differential equations (PDEs) coupled with algebraic expressions for the adsorption and reaction rates. The nonlinearity of adsorption isotherm and reaction term hinders the derivation of analytical solutions. For that reason, a flux-limiting high-resolution finite volume scheme is suggested to numerically approximate the model equations. The effects of several kinetic and thermodynamic parameters are rigorously analyzed on the reactant conversion and components separation.     

Theoretical study of the accuracy of the pulse method,frontal analysis,and frontal analysis by characteristic points for the determination of single component adsorption isotherms

The adsorption isotherms of selected compounds are our main source of information on the mechanisms of adsorption processes. Thus, the selection of the methods used to determine adsorption isotherm data and to evaluate the errors made is critical. Three chromatographic methods were evaluated, frontal analysis (FA), frontal analysis by characteristic point (FACP), and the pulse or perturbation method (PM), and their accuracies were compared. Using the equilibrium-dispersive (ED) model of chromatography, breakthrough curves of single components were generated corresponding to three different adsorption isotherm models: the Langmuir, the bi-Langmuir, and the Moreau isotherms. For each breakthrough curve, the best conventional procedures of each method (FA, FACP, PM) were used to calculate the corresponding data point, using typical values of the parameters of each isotherm model, for four different values of the column efficiency (N=500$N=500$, 1000, 2000, and 10,000). Then, the data points were fitted to each isotherm model and the corresponding isotherm parameters were compared to those of the initial isotherm model. When isotherm data are derived with a chromatographic method, they may suffer from two types of errors: (1) the errors made in deriving the experimental data points from the chromatographic records; (2) the errors made in selecting an incorrect isotherm model and fitting to it the experimental data. Both errors decrease significantly with increasing column efficiency with FA and FACP, but not with PM.     

Automated instrumental method for on-line fraction analysis and peak deconvolution in gradient multicomponent overloaded high performance liquid chromatography

An automated method for deconvolution of overloaded band profiles in gradient elution is described. The instrumental set-up consists of a pseudo-bidimensional HPLC system, where overloaded band profiles generated in the first direction are sampled on the second one. The method, previously employed under isocratic conditions, has been now extended to gradient elution, where one has to face the problem of band compression (and possibly band interference) during gradient, which decreases the time-window available along the first direction for sampling the overloaded profiles. The effect of the gradient steepness on the problem of defining a minimum number of sampling points to reconstruct single component bands from overloaded profiles is investigated. This approach is particularly useful in the framework of using inverse method (IM) to determine adsorption isotherm parameters for preparative gradient elution. In fact, it allows for the gathering of the information necessary to run IM calculations with minimum efford, also in cases where the individual component forming the mixture have different UV spectra.     

General theory of indirect detection in chromatography

Indirect detection is used in liquid chromatography to detect sample components that more or less lack detectable properties. A detectable component, the probe, is added to the mobile phase and the sample containing the components is injected. Every solute, even those not normally detectable, will give rise to detectable peaks due to the presence of the probe. The theory currently used to describe this situation is different for different adsorption isotherm models and also restrictions have to be imposed on the adsorption isotherm parameters. We present a general theory that describes both the retention times and areas of the detected peaks without making any a priori adsorption isotherm assumptions. Our general theory is compared with the current theory and validated by computer simulations. An alternative quantification method is suggested based on the measurement of the relative areas of the sample components without the requirement of having standard curves of known solutions of the sample components. Using this approach, the new theory is able to predict the solute concentrations in cases when the current theory did not, i.e., when the saturation capacities of the sample components and the probe differed from each other. In addition, a new application is proposed: the estimation of adsorption isotherm parameters.     

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

Single-component adsorption isotherm data of l-tryptophan on a C(18)-bonded silica column were acquired by frontal analysis (FA), with aqueous mobile phases containing 5% of acetonitrile at five different temperatures between 23 and 62 degrees C. The non-linear fitting of these data provided the bi-Moreau model for all temperatures as the best isotherm model. The inverse method (IM) was used to derive the parameters at these temperatures from the parameters of the 25 degrees C isotherm. The adsorption constants and the saturation capacities of the low and high-energy sites decreases by increasing the temperature, while the adsorbate-adsorbate parameters of both sites increase. An excellent agreement was found between the experimental and calculated overloaded band profiles at all the temperatures used. The breakthrough curves obtained and the overloaded band profiles obtained were found to have different shapes according to the range of concentration studied and the temperatures. At low concentration 0.05-0.5 g/L the breakthrough curves and the overloaded band profiles have a front shock and diffuse rear, which indicates langmuirian behavior, but at intermediate 1-2 g/L and high concentration 8 g/L they start to have diffuse fronts and shocks at the rear or more than one shock at the rear which indicates non-langmuirian behavior. At 23 degrees C the isotherm has another langmuirian part, which appears at high concentration. The behavior of the breakthrough curves is explained by the shape of the isotherm in which all of the isotherms have a langmuirian part (the isotherm is concave upward) and an antilangmuirian part (the isotherm is concave downward). The temperature affected the breakthrough curves by decreasing the time of the appearance of the fronts for all concentration ranges studied, and by decreasing the time difference between the highest concentration and lowest concentration of the fronts, especially the low concentration range at 0.5 g/L. The fronts of the breakthrough curves at high concentration seems to be the most affected by temperature.