General Solution of the Extended Plate Model including Diffusion, Slow Transfer Kinetics and Extra-column Effects for Isocratic Chromatographic Elution

The plate model proposed by Martin and Synge has been used for the characterization of columns up-to-date. In this approach, the column is divided into a large number (N) of identical theoretical plates. Mobile phase transference between plates takes place in infinitesimal steps with mixing of the solutions in the adjacent plates during the flow. The plate height is related to the band broadening that occurs in the mixing process due to the microscopic heterogeneities in the mobile phase flow. According to the original Martin and Synge model, solutes reach the equilibrium instantaneously in each theoretical plate, where dispersion is produced by: (i) convection or mixing of the mobile phase reaching a theoretical plate with that existing in that plate; and (ii) the equilibrium of the solute that is partitioned between mobile phase and stationary phase. In this work, a general method is proposed to solve the problem of chromatographic elution by means of an extended plate model assuming slow mass transfer, longitudinal diffusion in both mobile phase and stationary phase, and the extra column dispersion. The final equation was validated by comparing the results with those obtained through the numerical simulation of the solute migration using the finite differential approach. Experimental data were also used to check the validity of the derived equations.


Introduction
In 1941, Martin and Synge developed the principles of partition chromatography [1,2], and proposed the theoretical plate model to describe the elution in linear chromatography.These authors introduced the concept of "theoretical plate" in analogy to the description of fractional distillation [3].In 1952, both researchers were awarded with the Nobel Prize [4].Because of its simplicity, the plate model was useful to measure the chromatographic column performance through the proposal of the concept of number of column theoretical plates (N).Later, Craig [5] introduced the counter-current distribution and derived a mathematical treatment to explain the experimental results [6].In the plate model, the theoretical plate is represented as a column disk of very small thickness, where the equilibrium is reached instantaneously.In the Craig model, the elution is discontinuous and at each step the mobile phase inside a theoretical plate is replaced completely by the mobile phase in the preceding plate.In that time interval, the distribution equilibrium is assumed to be reached.This process results in a peak showing a binomial distribution [7], which can be approximated to a Gaussian with a variance equal to: where t R and t 0 are the retention time and dead time, respectively.According to this model, a non-retained compound would give rise to a peak with null width; hence, the peak broadening modeled by Equation ( 1) is due solely to the equilibrium between mobile phase and stationary phase.Instead, according to Martin and Synge, the elution is continuous and the mobile phase is transferred from one plate to the next in infinitesimal volumes, and the incoming phase is mixed with the theoretical plate receiver.Therefore, the proposed model implies an additional broadening due to this convective effect.Martin and Synge originally solved the plate model by the sequential study of the infinitesimal steps taking place during the elution process.Later, Said [8,9] performed a treatment based on the continuous flow of eluent through the column plates.More recently, we have studied the elution process by solving a system of differential equations consisting of mass balances that occur in each plate [10].For plate i: where k is the retention factor and n i the moles of solute in plate i at time t.The analytical solution of the system of differential equations yields the following peak profile versus time [7,10]: Equation ( 3) describes an Erlang distribution [11], which for sufficiently large values of N, tends to a normal distribution [7], with mean and variance equal to: In the Martin and Synge model, when a solute does not interact with the stationary phase, peak broadening is given by: The plate model theory described by Martin and Synge has been adopted as a basis for defining the efficiency and separation capability of chromatographic columns up-to-date.Despite its popularity, the plate model had only been solved for the simplest case that considers only the instantaneous equilibrium.By applying the Laplace transform, we managed to consider a slow mass transfer between mobile phase and stationary phase [10], proposing the following dispersion equation: where τ is a parameter related with the mass transfer due to the flow, and ν is the mass transfer constant between phases.
A more comprehensive study of the chromatographic elution process requires considering besides the equilibrium and slow mass transfer kinetics between mobile phase and stationary phase, the longitudinal diffusion in the mobile phase (and to a lesser extent, in the stationary phase), together with the eddy diffusion and heterogeneities in the flow [12,13].In addition, the dispersion due to extra-column sources should be included [14,15].The rate models [16] have been most widely used to consider some of these factors in both linear and non-linear elution [17], and for isocratic and gradient elution [18].However, its resolution can be difficult and requires several approximations [17,19], or the use of numerical methods [20][21][22].When an equilibrium with slow mass transfer and diffusion in the mobile phase is considered, the differential mass balance in an infinitesimal column segment can be expressed as: In the solid-film linear driving-force model, the rate of concentration change in the stationary phase is given by: Bq Bt " k f pq ˚´qq " k f pK c ´qq (9) and the differential equations can be solved as [18,23]: In Equations ( 8)-( 10), c and q are the concentration of the solute in the mobile phase and stationary phase, respectively; x is the axial coordinate along the column; u is the linear flow velocity; q* is the theoretical concentration of solute at equilibrium with c in the mobile phase; Φ is the volume phase ratio (Φ = V s /V m , where V s and V m are the volumes of stationary phase and mobile phase, respectively); D m is the solute diffusion coefficient in the mobile phase; k f is the mass transfer rate between phases; K is the partition constant; and L is the column length.
Band broadening is usually expressed in terms of column plate height (H): Accordingly, Equation ( 10) is developed as: where the first term has been added to take into account the eddy diffusion.A, B and C are constants that account for the contributions to band broadening from the eddy diffusion, longitudinal diffusion, and mass transfer resistance, respectively.This equation is known as "van Deemter equation" [24], and provides a basis to study the properties of chromatographic columns (column performance) that affect the peak shape, usually characterized by relating H to the linear mobile phase velocity (u), showing the different contributions to band broadening [25].
In another approach, proposed by Giddings and Eyring [26], the random migration of a single molecule is considered from a probabilistic point of view to describe the distribution function of the solute in the elution process [27], the elution profile being the probability density function of the retention time of the individual molecules.Recently, Chen [28] obtained identical results with this approach to those provided by the macroscopic models.
Other studies have described the band dispersion as a combination of independent processes, being the global variance the summation of the individual variances [29][30][31]: where σ 2 E , σ 2 d , σ 2 s and σ 2 ext are the variance associated to the eddy diffusion, longitudinal diffusion, slow mass transfer, and external factors, respectively.Therefore, each mechanism can be studied separately.Thus, for example, Knox et al. included the molecular diffusion in the stationary phase by considering it as a random walk process, where the theoretical plate height was given by [32]: Separations 2016, 3, 11 where the subscripts m and s account for the mobile phase and stationary phase, respectively, γ is the obstruction factor, D the diffusion coefficient, and the products γD quantify the effective diffusion coefficients.Following these approaches, several equations similar to Equation ( 12) have been developed [19,33,34] to explain the behavior of the theoretical plate height at varying mobile phase velocity.According to Knox, the most important contribution to dispersion under practical operating conditions comes from factors affecting the term A in Equation (12), the slow mass transfer within the porous particles being usually unimportant [35,36].On the other hand, Desmet et al. has reviewed the expression for term B [37], using a numerical procedure.A significant deviation was shown from the Knox type expression traditionally employed in high-performance liquid chromatography (HPLC).
In this work, a general approach is proposed to describe the peak broadening in liquid chromatography by applying an extended plate model that assumes slow mass transfer, the longitudinal diffusion in both mobile phase and stationary phase, and the extra-column contributions.The results have been validated by comparison with those obtained by simulating the elution using a finite difference approach.Experimental data were also used to check the validity of the derived equations.

General Description of the Chromatographic Elution Based on the Extended Plate Model
The next conditions will be considered to develop the approach: The solute crosses a fixed length L inside the column.(ii) The column is divided into N equal microscopic theoretical plates.(iii) The interaction between solute and stationary phase in each theoretical plate is determined by a linear isotherm: where rAs i and rBs i are the solute concentration in the mobile phase and stationary phase, respectively, at equilibrium inside the i theoretical plate, and b i and a i the corresponding number of moles.
In the theoretical plate theory, the column is divided into N microscopic and equal segments through which the mobile phase moves.The solute mass transfer between adjacent theoretical plates is influenced by the mobile phase flow, the mass transfer between mobile phase and stationary phase, and the longitudinal diffusion in both mobile phase and stationary phase.The mass balance between adjacent plates can be described by the differential equations given below.The transfer component due to the mobile phase flow is first quantified through [10]: where a i is the number of moles in the mobile phase inside the theoretical plate i, and τ the flow rate expressed as the number of theoretical plates per time unit: where F is the mobile phase flow rate.The mass transfer between mobile phase and stationary phase in each theoretical plate is described by the following equation [13]: Separations 2016, 3, 11 where ν is the mass-transfer rate constant and a i the moles of solute in the mobile phase when the equilibrium is reached.Considering that: and assuming a linear isotherm as in Equation ( 15), the next expression is derived: where s and m are the kinetic constants for the slow mass transfer of the solute from the stationary phase and mobile phase, respectively.Similarly, the change in the number of moles of solute in the stationary phase will be given by: Finally, the flow expressed as moles per unit area and time between two theoretical plates, due to the longitudinal diffusion, will be given by the Fick's law [7].For the mobile phase: where D m is the diffusion coefficient (length 2 ¨time -1 ), ∆c the difference in the solute concentration between two adjacent plates, and H the theoretical plate height; that is, the distance between the center of both plates.The change in the number of moles of solute in the mobile phase between two adjacent plates will be: where S is the mobile phase area, V m = SH, and ∆a is the difference in the number of moles between two adjacent plates.Therefore, taking into account the two theoretical plates adjacent to plate i, it results: for the mobile phase and for the stationary phase, where d m and d s are the diffusion transfer kinetic constants for the mobile phase and stationary phase, respectively: The change in the number of moles of solute due to all the described processes for plate i will be: db i dt " m a i ´s b i `ds pb i`1 ´bi q ´ds pb i ´bi´1 q (28) For the whole column, the following system of differential equations can be outlined: " ´τ pa i ´ai´1 q ´m a i `s b i `dm pa i`1 ´ai q ´dm pa i ´ai´1 q b p1q i " m a i ´s b i `ds pb i`1 ´bi q ´ds pb i ´bi´1 q . . .
where a  29) can be visualized as a kinetic system as shown in Figure 1, with the kinetic mass transfer constants τ, (τ + d m ), s, m, d m and d s (observe that the units for all these parameters are time ´1).These kinetics parameters are independent of the flow rate, except τ as shown in Equation (17).
For the whole column, the following system of differential equations can be outlined:  29) can be visualized as a kinetic system as shown in Figure 1, with the kinetic mass transfer constants τ, (τ + dm), s, m, dm and ds (observe that the units for all these parameters are time −1 ).These kinetics parameters are independent of the flow rate, except τ as shown in Equation ( 17).The solute monitored at the detector corresponds to the portion that flows out of the N plate to enter the detector in each dt.The peak profile at the detector is described by: Here we are assuming that the moles of solute that are detected at each time go from the last plate to the detector pushed by the flow.Therefore, the area of the function in Equation (30) will equal the moles of solute injected into the column.

Extra-Column Effects
The extra-column effects can be taken into account assuming that the solute moves through the external tubing similarly to through the column, but without interacting with the stationary phase.Thus, the tubes can be divided into M plates as shown in Figure 2. The solute monitored at the detector corresponds to the portion that flows out of the N plate to enter the detector in each dt.The peak profile at the detector is described by: Here we are assuming that the moles of solute that are detected at each time go from the last plate to the detector pushed by the flow.Therefore, the area of the function in Equation (30) will equal the moles of solute injected into the column.

Extra-Column Effects
The extra-column effects can be taken into account assuming that the solute moves through the external tubing similarly to through the column, but without interacting with the stationary phase.Thus, the tubes can be divided into M plates as shown in Figure 2.This involves a mass transfer process that can be described by the following system of equations.
(1) This involves a mass transfer process that can be described by the following system of equations.a p1q 1,ext " ´τext a 1,ext `dext pa 2,ext ´a1,ext q . . .a p1q i,ext " ´τext pa i,ext ´ai´1,ext q `dext pa i`1,ext ´ai,ext q ´dext pa i,ext ´ai´1,ext q . . .a p1q M,ext " ´τext pa M,ext ´aM´1,ext q ´dext pa M,ext ´aM´1,ext q (31) which would be added at the beginning of Equation System (29).The theoretical plate in the extra-column region has been considered to be different from the theoretical plate inside the column: where H ext is the extra-column plate height, t ext the extra-column time, and D the apparent longitudinal diffusion that can be described by the Taylor-Aris regime if the capillary diameter is large enough [38].
Regardless the longitudinal diffusion, the extra-column peak shape according to this approach is similar to that given by Equation (3). Figure 3 shows the fitting to Equation ( 3) of an extra-column peak obtained for a probe compound (sulfamerazine), which is fairly satisfactory.This involves a mass transfer process that can be described by the following system of equations. ( which would be added at the beginning of Equation System (29).The theoretical plate in the extracolumn region has been considered to be different from the theoretical plate inside the column: where Hext is the extra-column plate height, text the extra-column time, and D the apparent longitudinal diffusion that can be described by the Taylor-Aris regime if the capillary diameter is large enough [38].
Regardless the longitudinal diffusion, the extra-column peak shape according to this approach is similar to that given by Equation (3). Figure 3 shows the fitting to Equation ( 3) of an extra-column peak obtained for a probe compound (sulfamerazine), which is fairly satisfactory.3)) (solid line).Peak eluted using 1 mL/min flow rate and 20% acetonitrile.A.U. is arbitrary units.3)) (solid line).Peak eluted using 1 mL/min flow rate and 20% acetonitrile.A.U. is arbitrary units.

Peak Profile
The ordinary differential Equation System (29) is a first order system that can take the general form: where n is a vector containing the solute moles in both mobile phase and stationary phase in each column plate, t is the time, and Q is a square matrix (2N ˆ2N) containing the mass transfer constants.
The general solution for this system is given by: where m j are integers related to the multiplicity of the eigenvalues of matrix Q (m j takes a value unity if the multiplicity is one).The peak profile function only has a simple form when the elution follows the Martin and Synge ideal conditions.In this case, matrix Q has only one eigenvalue with N multiplicity and the peak equation is that shown in Equation (3).With other assumptions, the complexity of the peak function is extraordinarily increased [10].However, even assuming slow mass transfer, the peak equations tend to an almost Gaussian function (with very low asymmetry) if N is high enough.Therefore, given the high number of theoretical plates and the complexity reached by matrix Q for complex cases, the peak equation (Equation ( 35)) lacks of practical value, being the study through the calculation of moments preferable.

General Method
As explained in previous work [10], obtaining an explicit equation for the peak profile is highly difficult or unfeasible.Therefore, we characterized the peaks through the calculation of the associated moments, following the procedure described below: The equations describing the elution process (in this work, the sum of Equation Systems (31) and ( 29)) are first outlined.(ii) The system of differential equations is transformed to a system of algebraic equations in the Laplace space, which is a function of the variable r. (iii) The peak function is then obtained in the Laplace space (r-domain): F(r).(iv) To obtain the moments about the origin, the following property of the Laplace transformation is considered [10]: (v) For the moment of order n, the system of equations in the Laplace space describing the elution is derived n fold with regard to r.Then, it is solved to obtain the corresponding derivative after making r = 0 (see Equation (36)).(vi) The peak parameters are obtained from the moments about the origin.
Although the operations are simple and just consist in solving a linear equation system, the large number of implied equations makes the process rather complex.To follow it, it is interesting to understand the meaning of the derivatives at r = 0 in Equation (36).The derivative of order 0 is proportional to the amount of substance passing through the theoretical plate, the derivative of order 1 is proportional to the average time taken to cross the theoretical plate, and the derivative of order 2 is proportional to the average squared time to cross the theoretical plate.

Transformation of the System of Differential Equations into the Laplace Space
Obtaining the Laplace transform of the first-order Equation Systems ( 29) and ( 31) is rather simple.The Appendix of this work shows the whole process to calculate the first and second derivatives with regard to r.From Equation (30), the peak profile in the r-domain will be: where A N is the variable in the Laplace domain associated to the mobile phase moles in the last theoretical plate.Taking into account Equation (36) and considering that τ is a constant parameter: The peak position is given by the first moment about the origin: The variance is given by the second-order moment about the mean:

.3. Solution for the Mean and Variance of the Peak Profile
The derivatives of the Laplace function associated to the last theoretical plate have been obtained using the approach described above, as shown in the Appendix.The derived equations are: By substituting Equations ( 41)-(43) in the corresponding equations (Equations ( 38)-( 40)), the different moments are obtained: which agrees with the characteristics of the Laplace transform and the fact that, along the chromatographic elution (where τ > d m ), all solute molecules cross the whole column.
The mean elution time is obtained from Equation (39): where t ext is the extra-column time, t col = t 0 (1 + k) is the column time, and t R the experimental retention time.For this calculation, Equations ( 17) and ( 32) have been considered.
To calculate the variance, we will start from Equation (40), and operate considering Equations ( 41)-(43): Taking into account Equations ( 17), ( 26), ( 32) and ( 33), the following results: This equation agrees with Equation ( 5) when t ext = 0, D m = 0, D s = 0 and ν = 8.The first term is the extra-column variance.Considering that in the extra-column capillary tubes, the axial dispersion should be controlled by the Taylor-Aris regime [38]: where r ext is the radius of the extra-column tubing.Finally, the extra-column effect is given by: The column variance produced by a solute travelling in a unique path will be described by the second term in Equation ( 47), which could also include a term for considering the Taylor-Aris dispersion in the column: However, owing to the small interstitial diameter within the column and the smaller linear velocity, this term would only be significant at high velocities, and will not be considered.

Reagents and Columns
Five sulfonamides were used as probe compounds: sulfamerazine, sulfamethizole, sulfachloropyridazine, sulfisoxazole and sulfadimethoxine, purchased from Sigma (St. Louis, MO, USA), which were dissolved in a small amount of acetonitrile and diluted with water up to reach the concentrations of the stock and injected solutions (100 and 10 µg/mL, respectively), and stored in the darkness at 4 ˝C.
The mobile phase was prepared with 20% (v/v) HPLC-grade acetonitrile from Scharlau (Barcelona, Spain) and buffered at approximately pH 3 using 0.01 M sodium dihydrogen phosphate from Panreac (Barcelona, Spain) and HCl from Scharlau.Nanopure water obtained with a Barnstead purification system from Sybron (Boston, MA, USA) was used throughout.The mobile phase and the solutions of the probe compounds were filtered through 0.45 µm Nylon membranes with a diameter of 47 mm (Magna type) and 17 mm (Cameo type), respectively, from Osmonics (Herentals, Belgium).
Chromatographic analysis was carried out with a Spherisorb (150 mm ˆ4.6 mm ID) column with 5 µm particles from Scharlau.The column was protected with a pre-column of similar characteristics (30 mm ˆ4.0 mm).

Apparatus and Measurement of Peak Parameters
The HPLC system was equipped with the following modules from Agilent (Waldbronn, Germany): quaternary pump, autosampler, thermostated column compartment and UV-visible detector (Series 1200 and 1260).The signals were monitored at 254 nm.The flow rate was 1 mL/min and the dead time, 1.18 min.Dead time and extra-column contributions were measured by injecting a solution containing 20 µg/mL KBr from Acros Organics (Fair Lawn, NJ, USA).
The chromatographic system was controlled by a ChemStation from Agilent (B.04.03).The peak mean time and variance were calculated from the peak moments measured by numerical integration of experimental peaks using the Simpson rule.Data treatment was implemented in QB64.
The pH was measured with a potentiometer from Crison (Model MicropH 2002, Barcelona, Spain) and a combined glass electrode containing Ag/AgCl reference electrodes with 3.0 M KCl aqueous solution as salt bridge from Orion (Model 8102, Barcelona, Spain).The electrode was calibrated with aqueous buffers, and the pH of the mobile phases measured after the addition of the organic solvent.

Comparison of the Proposed Model with a Simulation Approach
The proposed approach was checked by comparison with the profile obtained by simulation of the solute migration using the finite difference approach [18], dividing the process inside a theoretical plate in 100 time intervals.This was repeated for all theoretical plates up to achieve full elution from the column.The variance values calculated with the simulation approach and predicted with Equation (46) are compared in Figure 4 for hypothetical compounds exhibiting different retention times and assuming several values of column and extra-column theoretical plates, flow rate, mass transfer, and longitudinal diffusion.The high agreement between both approaches in different situations proves the validity of Equations ( 46) and (47).and the dead time, 1.18 min.Dead time and extra-column contributions were measured by injecting a solution containing 20 μg/mL KBr from Acros Organics (Fair Lawn, NJ, USA).
The chromatographic system was controlled by a ChemStation from Agilent (B.04.03).The peak mean time and variance were calculated from the peak moments measured by numerical integration of experimental peaks using the Simpson rule.Data treatment was implemented in QB64.
The pH was measured with a potentiometer from Crison (Model MicropH 2002, Barcelona, Spain) and a combined glass electrode containing Ag/AgCl reference electrodes with 3.0 M KCl aqueous solution as salt bridge from Orion (Model 8102, Barcelona, Spain).The electrode was calibrated with aqueous buffers, and the pH of the mobile phases measured after the addition of the organic solvent.

Comparison of the Proposed Model with a Simulation Approach
The proposed approach was checked by comparison with the profile obtained by simulation of the solute migration using the finite difference approach [18], dividing the process inside a theoretical plate in 100 time intervals.This was repeated for all theoretical plates up to achieve full elution from the column.The variance values calculated with the simulation approach and predicted with Equation (46) are compared in Figure 4 for hypothetical compounds exhibiting different retention times and assuming several values of column and extra-column theoretical plates, flow rate, mass transfer, and longitudinal diffusion.The high agreement between both approaches in different situations proves the validity of Equations ( 46) and (47).Figure 5 depicts simulated peaks along with the corresponding Gaussian peaks, obtained using the peak profile parameters (mean time and variance) according to Equations ( 45) and (46).As observed, the agreement is highly satisfactory and the peaks were highly symmetrical, independently of the variation in mass transfer or diffusion (asymmetry values, measured as B/A at 10% peak high, were in the 1.04-1.09range).Figure 5 depicts simulated peaks along with the corresponding Gaussian peaks, obtained using the peak profile parameters (mean time and variance) according to Equations ( 45) and (46).As observed, the agreement is highly satisfactory and the peaks were highly symmetrical, independently of the variation in mass transfer or diffusion (asymmetry values, measured as B/A at 10% peak high, were in the 1.04-1.09range).

Comparison of Predicted Variance with Experimental Data
It is not possible to know the linear velocities needed in Equation (47), because the effective lengths are unknown.A transformation to volume and flow units is thus convenient, according to the following relationship: The final equation, obtained by combination of Equations ( 49) and (50), would be: Figure 6 shows the behavior of parameters AV, BV and CV according to Equations ( 53)-(55).As observed in the figure, the parameters show the behavior predicted in the equations.Thus, parameters AV and CV are linear against (1 + k) 2 and k(1 + k), respectively, while parameter BV gives rise to a third-order polynomial function against k.The quality of the fittings should be emphasized,

Comparison of Predicted Variance with Experimental Data
It is not possible to know the linear velocities needed in Equation (47), because the effective lengths are unknown.A transformation to volume and flow units is thus convenient, according to the following relationship: The final equation, obtained by combination of Equations ( 49) and (50), would be: where Figure 6 shows the behavior of parameters A V , B V and C V according to Equations ( 53)-(55).As observed in the figure, the parameters show the behavior predicted in the equations.Thus, parameters A V and C V are linear against (1 + k) 2 and k(1 + k), respectively, while parameter B V gives rise to a third-order polynomial function against k.The quality of the fittings should be emphasized, especially for the C V term, considering that there are other more complex contributions to the band broadening that have not been taken into account [38].In addition, the diffusion constants and mass transfer may vary with the eluted solute.Finally, a similar linear behavior of parameter C V with (1 + k) 2 could be found if the dispersion followed the Taylor-Aris regime.
Separations 2016, 3, 11 13 of 23 especially for the CV term, considering that there are other more complex contributions to the band broadening that have not been taken into account [38].In addition, the diffusion constants and mass transfer may vary with the eluted solute.Finally, a similar linear behavior of parameter CV with (1 + k) 2 could be found if the dispersion followed the Taylor-Aris regime.

Variance Components
From Equation (47), we can see that in the elution process the sum of variances established in Equation ( 13) is fulfilled and taken as an initial premise to obtain the usual equations to calculate the theoretical plate.According to this point of view, the variance components come from the following sources: (i) Equilibrium between the mobile phase and stationary phase according to a Craig process: (ii) Convection or flow through the column: This variance is exclusively due to the flow, and can be obtained considering that it follows Equation ( 6), which describes the variance of a solute that does not interact with the stationary phase.Assuming that this broadening is only experienced by molecules found in the mobile phase: where a* and b* are the moles of solute at equilibrium in the mobile phase and stationary phase, respectively.Finally, the ideal variance for a column without slow mass transfer or diffusion will be: (iii) Slow mass transfer between the phases: (iv) Mobile phase diffusion: (vi) Extra-column dispersion given by Equation (49) that includes convection, diffusion and dispersion due to the Taylor-Aris regime.
It is interesting to note that the components of band broadening due to mass transfer and diffusion do not depend on the theoretical plate height, but on the residence time inside the column.

Equation Describing the Theoretical Plate Height
The total variance is usually expressed in terms of the theoretical plate height, which includes all column band broadening factors.However, in the Martin and Synge elution model, the theoretical plate is a fixed parameter.Therefore, transforming Equation (47) to the usual terms, band broadening should be expressed as: Equation ( 63) gives an apparent value for the theoretical plate height (in the plate model H = L/N is a system constant).Therefore, the calculated theoretical plate (H app ) agrees with the real value only in the absence of longitudinal diffusion, multipath diffusion and slow mass transfer dispersion.In these conditions, as can be seen in Equation (10), the rate models would predict a null value for the band broadening.Thus, the extended plate approach and the rate models only agree when N tends to infinite.Equation (63) can be rewritten as: which agrees with the van Deemter equation except that the term A is associated to the dispersion due to convection.However, it can also be related to the particle diameter (d p ): (H = n H ˆdp , being n H the number of particles inside a theoretical plate).Besides, other sources of dispersion independent of the flow (as the multipath dispersion) should be included in this term.On the other hand, the term C in Equation ( 63) seems also different from that in Equation (10).Taking into account that ν = s + m and k = m/s, the mass transfer term in Equation ( 63) may be expressed with regard to the retention factor as: where each term assumes a constant kinetic parameter.

Conclusions
The original Martin and Synge approach has been extended to describe longitudinal diffusion in both mobile phase and stationary phase, slow mass-transfer and extra-column dispersion.To our knowledge, this is the first time the plate model has been solved including longitudinal diffusion in both phases, together with the slow mass transfer.Moreover, the extra-column dispersion has been included as a part of the dispersion model.
A systematic approach to obtain the moments of the elution profile applying the extended plate model has been developed.The final equation agrees with the van Deemter description for longitudinal diffusion in the mobile phase, and with the Knox equation, which also includes the diffusion in the stationary phase.The slow mass transfer results in a similar term to that in the van Deemter approach, but including a k/(1 + k) factor that fits better to the experimental data.The equations, expressed as a function of the theoretical plate height, provide an explanation for the three terms in the van Deemter equation.The term A is associated to the dispersion produced by the flow due to the mixing process inside each plate, and would only be null for N values tending to infinite.Other factors contributing to term A, such as the multipath dispersion, should be considered separately.
Simulated experiments showed that neither the slow mass transfer nor the longitudinal diffusion yield significant peak asymmetry.The plate model is assimilated to a kinetic mass transfer process coupling the various processes that take place both extra-column and inside the column.The proposed equation demonstrates the assumption of the sum of variances in the case of linear elution.Therefore, the plate model could be used to study separately each factor affecting the dispersion, adding in a further step each contribution to the final variance.The approach can be applied to study the factors affecting the elution process individually for educational purposes, only requiring knowledge on the solving of algebraic equations.

"
da i {dt, b p1q i " db i {dt, and a 1 ... a N and b 1 ... b N are the moles of solute in the mobile phase and stationary phase associated to each theoretical plate, respectively, which change with time; a M,ext is the solute that is entering the column from the M plate of the extra-column tubing; and m and s are the mass-transfer constants between mobile phase and stationary phase, and vice versa.The Equation System ( and a1 ... aN and b1 ... bN are the moles of solute in the mobile phase and stationary phase associated to each theoretical plate, respectively, which change with time; aM,ext is the solute that is entering the column from the M plate of the extra-column tubing; and m and s are the mass-transfer constants between mobile phase and stationary phase, and vice versa.The Equation System (

Figure 1 .
Figure 1.Kinetic representation of the elution inside the column as a combination of mass transfer processes from plate to plate (see Glossary for symbol explanation).

Figure 1 .
Figure 1.Kinetic representation of the elution inside the column as a combination of mass transfer processes from plate to plate (see Glossary for symbol explanation).

Figure 2 .
Figure 2. Kinetic representation of the extra-column mass transfer (see Glossary for symbol explanation).

Figure 2 .
Figure 2.Kinetic representation of the extra-column mass transfer (see Glossary for symbol explanation).

Figure 2 .
Figure 2. Kinetic representation of the extra-column mass transfer (see Glossary for symbol explanation).

Figure 3 .
Figure 3. Extra-column experimental peak for sulfamerazine (symbols), compared with that predicted according to the plate model (Equation (3)) (solid line).Peak eluted using 1 mL/min flow rate and 20% acetonitrile.A.U. is arbitrary units.

Figure 3 .
Figure 3. Extra-column experimental peak for sulfamerazine (symbols), compared with that predicted according to the plate model (Equation (3)) (solid line).Peak eluted using 1 mL/min flow rate and 20% acetonitrile.A.U. is arbitrary units.

Figure 6 .
Figure 6.Behavior of parameters AV, BV and CV calculated according to Equations (53)-(55).The plotted values were obtained by the non-linear fitting to Equation (52) of the experimental variances for the five sulfonamides eluted at several flow rates.

Figure 6 .
Figure 6.Behavior of parameters A V , B V and C V calculated according to Equations (53)-(55).The plotted values were obtained by the non-linear fitting to Equation (52) of the experimental variances for the five sulfonamides eluted at several flow rates.
the slow mass transfer τ mass transfer kinetic parameter related to the flow in the column τ ext mass transfer kinetic parameter related to the flow in the extra-column tubing ν mass-transfer rate constant between both phases (ν = m + s) A, B, C van Deemter equation constants A V , B V , C V constants in the variance-flow rate equation, where the variance is expressed in volume units a 0 initial moles of solute in the mobile phase associated to the first theoretical plate a i moles of solute in the mobile phase associated to the i theoretical plate a i moles of solute in the mobile phase at equilibrium in the i theoretical plate rAs i solute concentration in the mobile phase at equilibrium in the i theoretical plate b i moles of solute in the stationary phase associated to the i theoretical plate b i moles of solute in the stationary phase at equilibrium with the mobile phase in the i theoretical plate rBs i solute concentration in the stationary phase at equilibrium associated to the i theoretical plate c solute concentration in the mobile phase as a function of t and x D m , D s solute longitudinal diffusion coefficient in the mobile phase and stationary phase d m , d s diffusion transfer kinetic constants in the mobile phase and stationary phase d ext diffusion transfer kinetic constant in the mobile phase in the extra-column tubing change in the total moles of solute per time unit q solute concentration in the stationary phase as a function of t and x q* solute concentration in the stationary phase at equilibrium with the solute in the mobile phase at concentration c r Laplace variable s kinetic constant for the solute slow mass transfer from the stationary phase to the mobile phase t time coordinate t 0 column dead time t col solute retention time associated to the column t ext solute retention time associated to the extra-column tubing t R solute retention time or peak mean time (time at the peak maximum for symmetrical peaks) u linear flow velocity V m mobile phase volume associated to a theoretical plate V s stationary phase volume associated to a theoretical plate x axial coordinate γ m , γ s obstruction to the diffusion transfer factors in the mobile phase and stationary phase