We introduce the mathematical model of electromigration of electrolytes in free solution together with free software Simul, version 5, designed for simulation of electrophoresis. The mathematical model is based on principles of mass conservation, acid-base equilibria, and electroneutrality. It accounts for any number of multivalent electrolytes or ampholytes and yields a complete picture about dynamics of electromigration and diffusion in the separation channel. Additionally, the model accounts for the influence of ionic strength on ionic mobilities and electrolyte activities. The typical use of Simul is: inspection of system peaks (zones), stacking and preconcentrating analytes, resonance phenomena, and optimization of separation conditions, in either CZE, ITP, or IEF.

We are introducing a computer implementation of the mathematical model of zone electrophoresis (CZE) described in Stedry, M., Jaros, M., Hruska, V., Gas, B., Electrophoresis 2004, 25, 3071-3079 program PeakMaster. The computer model calculates eigenmobilities, which are the eigenvalues of the matrix tied to the linearized continuity equations, and which are responsible for the presence of system eigenzones (system zones, system peaks). The model also calculates other parameters of the background electrolyte (BGE)-pH, conductivity, buffer capacity, ionic strength, etc., and parameters of the separated analytes--effective mobility, transfer ratio, molar conductivity detection response, and relative velocity slope. This allows the assessment of the indirect detection, conductivity detection and peak broadening (peak distortion) due to electromigration dispersion. The computer model requires the input of the BGE composition, the list of analytes to be separated, and the system instrumental configuration. The output parameters of the model are directly comparable with experiments; the model also simulates electropherograms in a user-friendly way. We demonstrate a successful application of PeakMaster for inspection of BGEs having no stationary injection zone.

A mathematical model of capillary zone electrophoresis (CZE) based on the conception of eigenmobilities, which are the eigenvalues of a matrix M tied to the linearized governing equations is presented. The model considers CZE systems, where constituents, either analytes or components of the background electrolyte (BGE), are weak electrolytes--acids, bases, or ampholytes. There is no restriction on the number of components nor on the valence of the constituents nor on pH of the BGE. An electrophoretic system with N constituents has N eigenmobilities. In most BGEs one or two eigenmobilities are very close to zero so their corresponding eigenzones move very slowly. However, there are BGEs where no eigenmobility is close to zero. The mathematical model further provides: the transfer ratio, the molar conductivity detection response, and the relative velocity slope. This allows the assessment of the indirect detection, conductivity detection and peak broadening (distortion) due to electromigration dispersion. Also, we present a spectral decomposition of the matrix M to matrices allowing the assessment of the amplitudes of system eigenpeaks (system peaks). Our model predicted the existence of BGEs having no stationary injection zone (or water zone, gap, peak, dip). A common practice of using the injection zone as a marker of the electroosmotic flow must fail in such electrolytes.

Simul 5 Complex is a one-dimensional dynamic simulation software designed for electrophoresis, and it is based on a numerical solution of the governing equations, which include electromigration, diffusion and acid-base equilibria. A new mathematical model has been derived and implemented that extends the simulation capabilities of the program by complexation equilibria. The simulation can be set up with any number of constituents (analytes), which are complexed by one complex-forming agent (ligand). The complexation stoichiometry is 1:1, which is typical for systems containing cyclodextrins as the ligand. Both the analytes and the ligand can have multiple dissociation states. Simul 5 Complex with the complexation mode runs under Windows and can be freely downloaded from our web page http://natur.cuni.cz/gas. The article has two separate parts. Here, the mathematical model is derived and tested by simulating the published results obtained by several methods used for the determination of complexation equilibrium constants: affinity capillary electrophoresis, vacancy affinity capillary electrophoresis, Hummel-Dreyer method, vacancy peak method, frontal analysis, and frontal analysis continuous capillary electrophoresis. In the second part of the paper, the agreement of the simulated and the experimental data is shown and discussed.

The Kohlrausch regulating function (KRF) is a conservation law (conservation function), which is held in electrophoresis and which enables calculation of the so-called adjusted concentrations of constituents. The KRF is not the only conservation function and, depending on the complexity of the electrophoretic system, other conservation laws may be obeyed having a broader range of applicability. The conservation laws are tightly related to system eigenmobilities and system zones (system peaks). In principle, no system eigenmobility is exactly zero, but in most practical cases at least one system's eigenmobility is close to zero. The existence of the close-to-zero eigenmobility inherently points to the existence of a conservation function and a system zone which is stationary. The stationary system zone is called injection zone, stagnant zone, water peak, or solvent dip. Electrophoretic (electromigration) systems can be divided into two types: (i) conservation systems, in which the absolute value of at least one system eigenmobility is close to zero and where at least one conservation law is obeyed and (ii) nonconservation systems, where no system eigenmobility is close to zero and no conservation law is obeyed. The paper reviews work dealing with conservation functions in electromigration, derives some "historical" conservation functions in a new way, derives several conservation functions for systems of multivalent electrolytes, and discusses electrophoretic systems that have nonconservation behavior. In some typical instances, the conservation functions are simulated by means of a dynamic simulation tool and depicted graphically.

The complete mathematical model of electromigration in systems with complexation agents introduced in the Part I of this article (V. Hruška et al., Eletrophoresis, 2012, 33, this issue), which was implemented into our simulation program Simul 5, was verified experimentally. Three different chiral selector (CS) systems differing in the type of the CS, the magnitude of the complexation constants as well as in the experimental conditions were selected for verification. The experiments and simulations were performed at various concentrations of the CSs in order to discuss the influence of the concentration of the CS on the separation. The simulated and experimental electropherograms show very good agreement in the position, shape and amplitude of the analyte peaks. The new Simul 5 Complex offers a deep insight into electrophoretical separations that take place in systems containing complexing agents, for example into enantiomer separations. Using Simul 5 Complex we were able to predict and explain the significant electromigration dispersion of analyte peaks. It was clarified that the electromigration dispersion in these systems results directly from complexation. The new Simul 5 Complex was also shown to be a useful and powerful tool for the prediction of the results of enantioseparations.

We extended the linearized model of electromigration, which is used by PeakMaster, by calculation of nonlinear dispersion and diffusion of zones. The model results in the continuity equation for the shape function ϕ(x,t) of the zone: ϕ(t) = -(v(0) + v(EMD) ϕ)ϕ(x) + δϕ(xx) that contains linear (v(0)) and nonlinear migration (v(EMD)), diffusion (δ), and subscripts x and t stand for partial derivatives. It is valid for both analyte and system zones, and we present equations how to calculate characteristic zone parameters. We solved the continuity equation by Hopf-Cole transformation and applied it for two different initial conditions-the Dirac function resulting in the Haarhoff-van der Linde (HVL) function and the rectangular pulse function, which resulted in a function that we denote as the HVLR function. The nonlinear model was implemented in PeakMaster 5.3, which uses the HVLR function to predict the electropherogram for a given background electrolyte and a composition of the sample. HVLR function also enables to draw electropherograms with significantly wide injection zones, which was not possible before. The nonlinear model was tested by a comparison with a simulation by Simul 5, which solves the complete nonlinear model of electromigration numerically.

No abstract