Mathematical Modelling in Analytical Chromatography: Problems and Solutions

Based on an analysis of the results of original research performed in the Laboratory of Sorption Methods of the Vernadsky Institute of Geochemistry and Analytical Chemistry of the Russian Academy of Sciences within the project “Mathematical Chromatograph,” the review covers the aim and strategy of the imitation modeling of high-performance chromatography; associated problems of the theory of intermolecular interactions; classifications of polar stationary phases by their selectivity; descriptions of the kinetics and dynamics of sorption processes, choice of the composition of multicomponent mobile phases in HPLC and ion chromatography using the method of the dynamic map of a chromatographic system; and the development of alternating gradient modes using a mathematical experiment.

Chromatographic processes are characterized by high requirements to the composition of mobile and stationary phases and to the precision of implementation of complex modes of the separation of multicomponent mixtures. However, in practice, the choice of phases and conditions is often accidental with the parameters far from the optimum ones.
The state-of-the-art theory of sorption dynamics in combination with the results of the physicochemical description of selectivity and kinetics of sorption allow a researcher to make this choice conspicuous and more perfect. The software products developed on the basis of the mathematical modeling of chromatographic processes that take into account their specific features that allow a researcher to perform full-scale mathematical experiments are especially convenient for the search for optimal conditions and new separation modes. Such imitation modeling must implement high-level mathematical models using no adjustable parameters and, therefore, capable of a priori calculations.
For this purpose, we put forward a new concept of chromatographic analysis and developed a corresponding methodology, based on the maximum replacement of real experiments by simulations, performed using high-level mathematical models. Mathematical experiments require solution of the so-called direct modeling problem, which ensures the calculation of the result of a process based on specified input data. With an adequate model, the results of a priori calculations agree well with the results of test experiments. By varying the parameters of an instrument and solving the direct problem for each version, one can choose such values of parameters that satisfy restrictions of the method and the optimality criterion. This is the aim of using a mathematical model for the analysis optimization. The inverse problem of modeling, being an instrument for an analysis of the results of a real experiment, can be posed only at a successful solution of the direct problem. When performing chemical analyses, one should interpret an analytical signal as chemical information. This is a mathematically incorrect procedure, associated with the multiplication of facts based on the assumptions and intuition of an operator and associated with the risk of a false result. The use of the inverse problem of modeling, based on the theoretical knowledge of the features of the applied method, significantly increases the level of accuracy and reliability of the results of analysis. Thus, imitation modeling of chromatography allows a researcher to considerably reduce the volume of experimental work in developing analytical procedures, avoid mistakes in the identification of analytes, and reduce the systematic component of the experimental error.
The general aim of imitation modeling in high-performance chromatography can be subdivided into a number of problems for the development of mathematical models of methods of gas and ion chromatography and HPLC and the development of correspond-ing software products. This goal was considered within the initiative project no. 18-03-00382, supported by the Russian Foundation for Basic Research in 2018-2020 under the conditional name "Mathematical Chromatograph." Figure 1 presents a flowchart of correlations in the software for the imitation modeling of methods of analytical chromatography taking into account the following features: -the dynamic process in chromatography follows the laws of quasi-equilibrium thermodynamics and includes equilibrium and kinetic parameters; -equilibrium is determined by the nature of intermolecular interactions, in accordance with the participation of nonpolar, polar, and ionic interactions and hydrogen bonds; -kinetics is determined by the values of diffusion constants in the mobile phase (MP) and stationary phase (SP) and by the sizes of the substrate and the SP working area; -the input data of a model are parameters of the experimental system and the fundamental properties of substances forming the phases and analytes.

MATHEMATICAL MODELS OF ANALYTICAL
CHROMATOGRAPHY Usually modeling in chromatography is based on abundant empirical data, obtained by the manufacturers of equipment using interpolation methods [1,2]. Such approaches are laborious, specific to the equipment, and are quite limited in the optimization and search for new modes. There are also semiempirical models, more justified, and to a certain extent capable of predicting the results [3][4][5]. However, because of a large number of empirical parameters and the inaccuracy of correlations, the search for new, more effective modes with their help is almost excluded. In contrast to them, within the project "Mathematical Chromatograph," we developed and constantly improve high-level models based on the fundamental theoretical description of chromatographic methods of analysis [6,7]. Such models open wide possibilities for performing mathematical experiments, also under very complex conditions that cannot be found in the random experimental search.
General dynamic properties of chromatographic processes. The description of a chromatographic process can be conditionally subdivided into two almost nonoverlapping fields: the thermodynamics and the   States of aggregation for phases and kinetic data macrokinetics of sorption. The relevance of the description of the macrokinetics of chromatography is in many respects associated with the successful solution of the problem of the a priori calculation of a sorption equilibrium. Is it reasonable to calculate peak shape, height, and width if we cannot predict its position in the chromatogram, which is determined by the equilibrium? The consideration of a dynamic problem of high-performance column chromatography, as a first approximation, gives a chromatogram equation in the form of an envelope of Gaussian curves constructed for the sample components (analytes): (1) where is the level of signal of the i-th component at the instant t in the output section of the column; is mass of the i-th component in the injected sample; is coefficient of detector sensitivity to the ith component; MP flow rate at the moment of elution of a component peak; and is the dispersion of the peak due to the kinetic properties of the MP, sample component, and structural parameters of the SP.
The position of the maximum of a chromatographic peak in the general case of variable conditions and a compressible MP is determined as from the integral equation: (2) where is the function obtained as a result of solution of Eq. (2); is dead time; L is column length; is linear velocity of the mobile phase in point x at the instant when the component front intersects it; function is determined by the boundary conditions from the eluotropic force of the MP and the temperature of the system. In particular, in the isocratic/isothermal mode, k = const, and retention time takes the usual form: The kinetic characteristics of systems of high-performance chromatography were theoretically well studied and are displayed in terms of the height equivalent to a theoretical plate (HETP) [8], which is described by the van Deemter equation [9]: HETP depends not only on hydrodynamic (d is sorbent mesh, is linear velocity of the MP) and kinetic parameters ( are analyte diffusion constants in the MP and SP respectively, r is working area of the SP), but also on equilibrium characteristics ( is retention factor related to the Henry constant of the analyte by the phase ratio). Calculations of HETP for each model possess their specific features, the account of which allows a researcher to avoid experiments on the determination of diffusion constants and other empirical parameters [7].
The value inverse to HETP has meaning of the specific efficiency of a column, measured in the number of theoretical plates per unit column length; therefore, the efficiency of the column (number of theoretical plates N) is determined by the integral (4) Modeling of gas chromatography. Version of gas adsorption chromatography. The instrument proposed in the works [10,11] and most completely described in the monographs [6,7] presents the new nonlinear molecular and statistical (NMS) theory of adsorption in the Henry region, based on the theory of generalized charges (TGC), which was developed based on the general principles of quantum statistics [12][13][14]. Unlike the semiempirical methods based on the rough postulate of the additivity atom-atomic potentials [15][16][17], in TGC another summation rule was derived [18][19][20], because of which an acceptable accuracy of the description of interatomic and intermolecular interactions without using adjustable parameters was attained. The expression for the potential of van der Waals interactions, coinciding with the Lennard-Jones potential in shape and requiring no empirical parameters, was derived from TGC. On this basis, a model of the critical state of substance, extremely important for the theory of first-order phase transitions, was proposed in [21,22]. Using TGC and the model of the critical state, a relationship between the parameters of intermolecular interactions and critical temperature was considered [23]. On an example of inert and simple gases, ab initio calculations were per-  formed, the results of which well agreed with the experimental data.
The application of the NMS theory to a priori calculations of the adsorption of several hundred of organic molecules (linear, branched, and cyclic hydrocarbons; oxygen-and halogen-containing compounds) on graphite gave a good agreement with the experimental data on the temperature dependence of the Henry constant for different adsorbates and adsorbents [6,7,11,14,24,25]. The precision with which this dependence is known strongly affects a possibility of the prediction of the components behavior under complex conditions of programmable temperature.
Block schemes of the GEOMOL and THENRY software, developed for calculations of the temperature dependence of the partition coefficient by the structural formula of a sample component of test are presented in Fig. 2. However, the resolution of component peaks in a mixture depends not only on their position, but also on their width. The correlations of a dynamic process of chromatography in the model of gas chromatography were complemented with estimates of diffusion coefficients and of HETP as a whole in the description of kinetics [25,26]. Therefore, the mathematical basis in modeling of molecular chromatography is provided by fundamental theoretical ideas of the adsorption and kinetic mechanisms, implemented in the chromatographic process. Based on the described mathematical fundamentals, the MOLCHROM software, ensuring the solution of the direct problem of modeling of gas adsorption and gasliquid chromatography with nonpolar phases (Fig. 3), was developed.
Version of gas-liquid chromatography. In the description of processes of gas-liquid chromatography, the poorly studied problems are mainly those associated with the imperfection of the characterization of selectivity of polar stationary phases. The characterization of the selectivity of an SP in gas chromatography faces the problem of the determination of desorption energy in the system with dispersion, electrostatic, and hydrogen bonds. For the very important value, the contribution of hydrogen bonds (H-bonds), there still has been no consensus about its nature, and, therefore, value. Based on the conclusions of the traditional quantum theory of chemical bonding about the lack of vacant orbitals for the formation of Hbonds, it was considered after Poling [27] that this bond is of electrostatic rather than quantum nature. Today this version is implemented in models forming the basis of semiempirical methods of calculations of intermolecular interactions (IMI). Contrary to the observed feature of the "classical" H-bond, the longwavelength shift of the IR spectrum of vibrations of the covalent bonds of hydride on the formation of a molecular complex [28], the assumption about the electrostatic nature of the hydrogen bond can only explain the effect of the "blue" shift of the specified spectral band, for which the special term "non-classical" or "improper" hydrogen bond was introduced [29]. The contradiction between the specific nature of the hydrogen bond and electrostatic forces for a particular system is not pronounced, because the reason for these types of interactions is common: the presence of a hydrogen atom and electronegative atoms; the strengths of interactions are also of equal level. Nevertheless, works on hydrogen bonds based on the statistical analysis of empirical data are known, in which electrostatic and covalent forces are independent components [30]. Because of basic falseness of the model considering only dispersion and electrostatic IMI, the convenience of the approaches implemented in computer programs is reduced by the random choice of the parameters of calculations, relating to the set of conformations, fixed angles, and the number of the formed H-bonds, which results in the unreliability of predictions, which, as a result, must be checked and corrected [31,32].
A fundamental difference between a classical electrostatic interaction of charges localized on atoms and a quantum-mechanical (electron) bond consists in the direct participation of electrons in the formation of the quantum-mechanical bond. The steady state of an atomic system is provided by a finite movement of a hydrogen electron in the region of such a bond; therefore, its kinetic energy cannot exceed the absolute value of the potential energy of the bond. Having reverted this relation, we will obtain that binding energy cannot be lower a certain value dictated by the quantum uncertainty. Therefore, the quantum nature of an H-bond manifests itself in a special effect: the entrance of a molecule in the region of a small solid angle with the vertex on a hydrogen atom of another molecule leads to a jump of the energy of IMI, which can be presented as a drop of a system in a narrow potential well. Therefore, in addition to the effects of polar nature, depending on dipole moments, hydrogen bonding, as dispersion forces, makes an independent contribution to IMI.
The potentials of the approach developed on the basis of TGC were shown not only in a priori calculations of adsorption and van der Waals interactions [18][19][20], but also in the description of covalent bonds [33][34][35]. Hydrides are obligatory participants of hydrogen bonds. The accuracy of the developed approach to the description of hydride characteristics was demonstrated in [35]. The derived simple expression for the length of bond in hydrides with a good accuracy corresponds to experimental values for a wide series of chemical elements: (5) where are the number of the period and atomic number of acceptor X and is Bohr's radius. According to the developed theory of H-bond, the hydride and hydrogen bonds are of common nature, and their characteristics for identical elements are correlated. The regularities of the H-bond were studied in [36], where expressions for its characteristics were derived and criteria of its existence were formulated. The theoretical graph (Fig. 4), constructed on the coordinates ( are the length and energy of an H-bond, respectively; Z is atomic number of the acceptor in the H-bond; and e is elementary charge), demonstrates the region of the existence of a hydrogen bond, the borders of which separate classes of substances by their ability to form hydrogen bonds. Based on quantitative estimates, explanations of the presence or absence of solubility in water are given for substances bearing electronegative atoms. It was shown that the energy of a hydrogen bond cannot be lower a certain value: Expressions for the limits of the minimum energy and maximum length of a hydrogen bond through the atomic number of an acceptor can be derived from the theory: (kJ/mol), The developed theory gives quantitative explanations for some "irregularities," for example, to the effect of the hydrophobicity of chloroalkanes, despite the ability of the hydrogen chloride to hydrogen bonding [37].
In [37,38], the new approach to the description of H-bond, as in the general model of IMI, was expanded to organic matter and to the interaction of vapor molecules with a liquid. The energy of IMI, according to the proposed model, contains three independent groups of values, describing nonpolar and polar forces and hydrogen bonds. Each of forces was presented by its particular molecular descriptor, generalized charge, dipole moment, and two numbers reflecting the ability of a molecule to be a donor or an acceptor of an H-bond, respectively. Using TGC, "structure-property" correlations were obtained for all components of molecular interaction energy. In particular, the contribution of an H-bond to the total energy of IMI as a product of the threshold value (4), being of quantum-mechanical nature, and probabilities of the regular arrangement of interacting molecules n H were described. The last value was determined as the probability of the entrance of a system of interacting molecules in a narrow potential well and was expressed through the parameters of the molecular structure (Fig. 5).  Therefore, new ideas about the H-bond, together with the conclusions of TGC, allowed us to express the energy of IMI of a general type, and, in particular, the energy of adsorption of polar analytes on polar phases, through structural descriptors.
Traditional empirical methods of the classification of stationary phases are known in chromatography [3,39]. The five-dimensional Rorschneider scheme (and similar procedures) is based on the anticipated independence of five types of energy for five reference samples. However, "pure" substances, which participate in only one type of interaction, do not exist. In addition, independent types of intermolecular energy are not numerous, i.e., less than five. From here follows the inaccuracy and redundancy of the scheme, because the space of Rorschneider parameters is of lower dimensionality, and the characteristics of reference substances summarize contributions of different types of energy and, therefore, cannot serve as basis vectors of the coordinate system. Other methods [40] utilize the one-dimensional scheme of a hydrophobic-hydrophilic equilibrium, in which stationary phases possess one estimation parameter, which is obviously insufficient because of the complex nature of intermolecular interactions.
The new method of three-parameter characterization (TPC), developed on the basis of the proposed model, is theoretically justified and in balance by the number of independent variables in the description of IMI energies [41]. In the TPC method, two characteristics of an SP, polarity ( is generalized charge, is squared dipole moment) and hydrophily (n H is probability of H-bond formation) are determined from three molecular properties (generalized charge, dipole moment, and the parameter of Hbonding), which can be conveniently plotted on a map, which is named the map of phase selectivity. The TPC method allows an analyst not only to perform a priori calculations of the characteristics of phase selectivity, but also to reveal new dependences of these characteristics on the specific features of the molecular structure of phases. In particular, a dependence of the characteristics of polarity and hydrophily of polyethyleneglycol (PEG) on the weight of the polymer molecule was found in [42]. The derived theoretical dependences of the polarity and hydrophily of PEG on molecular weight are presented in Fig. 6 as a parametrical dependence of hydrophily on polarity and well correlate with the experimental data plotted on the selectivity map.
The IMI model was also used in solving another important problem of analytical chromatography, the determination of the conditions and SP parameters of the fine chromatographic separation of the geometrical isomers of fatty acid methyl esters [43]. The inversion of the selectivity of different phases, observed in practice for transand cis-isomers, was explained by two sorption mechanisms, one of which was responsible for physical adsorption (and the cis/trans elution order), and the other, for the absorption of the adsorbate by the macromolecule of the polymer phase (and reverse retention order). The last mechanism is more efficient, but requires the obligatory formation of H-H n w Q =  bonds and a correspondence of the sizes of the macromolecule globule to the adsorbate length. The attained accuracy of the IMI model opens possibilities for the consideration of the direct modeling problem in gas chromatography: the determination of retention parameters by the specified structures of analytes and the stationary phase. Therefore, the model of gas-adsorption chromatography, developed for nonpolar substances and phases, is supplemented by the description of polar interactions, which ensures the development of a general model of gas-liquid chromatography. The adequacy of the model was checked at different stages of its development. The agreement between the results of calculations for the direct problem and experimental data for gas chromatography is so satisfactory that there appears a possibility for the creation of a simulator program of a stateof-the-art gas chromatograph with wide possibilities in the selection of a carrier gas, a column, a sorbent, test sample components, temperature, and gasdynamic modes.
Modeling of liquid chromatography. In HPLC, problems of the description of competitive sorption, multicomponent eluents, solvation, dissociation, etc., peculiar to ion chromatography, are added to the problems of modeling of gas chromatography. The phenomenological model of HPLC [6,44], summarizing present-day retention theories [3][4][5], in combination with TGC allowed us to describe real systems almost without empirical parameters rather well, except for such transferable parameters as sorption capacity and standard energy of the adsorbent. Figure 7 presents an example of a comparison of the results of a priori calculations of the retention factor of naphthalene in the system of reversed-phase HPLC depending on the percentage of modifier in the MP with the experimental data [45].   1 and 2) with H-bond atoms A 1 (acceptor) and D 2 (donor) and independent motion of molecules relative to the center of a dispersion interaction (shown by arrows) that must form for the system of two molecules fall in a potential well. The concept of a narrow potential well illustrates the graph for potential energy (shown in the bottom-right corner). The specified model was complemented with an equation for the energy of IMI, derived for polar systems in gas chromatography, for the description of various versions of solvation, affecting retention in HPLC [46]. Because of its extreme complexity, this problem requires further investigation.
An equation for the eluotropic force of a multicomponent MP in liquid chromatography for the case of effective eluents was derived in [47]; it is similar to the equation for an eluent in ion chromatography [48,49]. Eluotropic force ϕ is a unified characteristic of a multicomponent mixture, similar to the concentration of a single-component MP. It obeys the following equation, common for ion and liquid chromatography: Modeling of ion chromatography. The principles and possibilities of an approach to modeling of ion chromatography (IC), based on the achievements of the theory of the dynamics of sorption, adapted to the specific features of ion chromatography, were described in [48,49]. As a result of computer modeling, the IONCHROM© software (Fig. 8) was developed [50].
Using the IONCHROM software, problems of multicolumn IC with isocratic and gradient elution were considered [51,52]. For the determination of the parameters of a new chromatographic system, necessary for modeling, its dynamic map is built, by which the conditions of the required separation in different pH and concentration ranges are calculated, after which three or four experiments with the chosen eluents are sufficient to correct the parameters. The model refined with the found parameters is used to predict the behavior of the system and, in particular, to optimize the analysis [53].

OPTIMIZATION AND DYNAMIC MAP OF A CHROMATOGRAPHIC SYSTEM
An approach to the description of a multicomponent MP was used to optimize HPLC analysis by an economics criterion [54]. Based on example of a chromatographic experiment from a published source, reference parameters of an HPLC system used in there were found and the optimum composition of the MP  The choice of chromatographic phases is not limited by the attainment of the required selectivity. Also research into the kinetics of sorption are performed for optimizing stationary phases by their resolving capacity and the rate of separation. To compare a new SP, developed for analytical chromatography, a new criterion, characterizing the efficiency and rapidity of the SP, was proposed [55]: (8) where is the diffusion constant of a component in the MP, and are the duration and efficiency of the determination of the component, and is the mesh of column packing. The meaning of the characteristic of the SP by efficiency and performance is given the name "sorbent performance coefficient:" the higher the value of the criterion, the higher the product of efficiency on the throughput of analysis by the time of longitudinal diffusion along the sorbent grain size.
The choice of the best conditions for the separation of mixtures taking into account the selectivity and efficiency of the chromatographic system is the central problem of optimization in chromatography. In the published data, a set of optimality criterions aimed at the achievement of the best separation within a certain time of analysis, the maximum capacity of analysis, etc., was developed [3,56]. An important criterion of optimization is the concept of peak capacity, introduced in the Kaiser works [57], as the number (it may be non-integer) of additional peaks that can be placed in chromatogram ranges free from peaks. The search for a minimum of peak capacity corresponds to the aim of increasing the performance of chromatography. However, peak capacity is not a very convenient optimization criterion, because precisely zero peak capacity is almost unattainable, and the result with nonzero peak capacity is degenerate and requires additional criteria. A more precise criterion is the duration of a single analysis, which is minimized in the fulfillment of requirements for the degree of component separation [48,51]. The effective use of this criterion requires the availability of a mathematical model of the chromatographic method, leading to the creation of a dynamic map of a chromatographic system.
The system of equations of a model of liquid chromatography can be presented in a graphical form, as a dynamic map of a chromatographic system (DMCS), which is convenient as a characteristic of the chromatographic behavior of the components of interest and an instrument for the search for the optimum separation modes.   The abscissa in the dynamic map is here the eluotropic force of MP is calculated by Eq. (7). The logarithm of relative retention is plotted on the axis of ordinates (axis y) if the characteristic of the benchmark component is taken equal to zero: (9) For each sample component, in the DMCS there is a corridor, constructed around the medium line described by the equation: (10) where unit values of the constant and charge are considered for the BC. The borders of the corridor are related to the characteristics of the chromatogram, i.e., corrected retention time and peak half-width , by the equation (11) where is the half-width of -th corridor, is the degree of peak resolution (the required degree of separation of a component with the log , . The zones of corridor crossing in the DMCS correspond to peaks in the chromatogram that are not separated to a necessary degree. As follows from Eq. (10), analytes with different charges can in all cases be separated by selecting a corresponding mode in the dynamic map, because their corridors are at different angles and, starting from certain point, do not intersect. We can also note that, on different sides from the point of intersection of medium lines, there are areas of conditions with different elution orders of the components. The medium lines of two equally charged ions of a sample are parallel, so that their elution order cannot be changed, and only band width depends on the eluent strength.
A DMCS for a mixture of anions in a model system with characteristics of a DIONEX (United States) ion chromatograph of is shown in Fig. 9 [58].
Gradient elution is the convenient method of analysis of mixtures with components strongly differing on retention. Simple models of retention, ensuring calculations of unicomponent eluents and linear gradients, were developed in not numerous works on modeling of gradient chromatography [4,5]. However, the level of these models is obviously insufficient for the consideration of complex gradient modes; the degree of the separation of sample components attained in the calculations is taken into account poorly, which strongly depreciates the role of modeling in the development of procedures and in analysis (correspondingly, in solving of direct and inverse problems). The main equa-i n Fig. 9. Dynamic map of a chromatographic system for the separation anion mixtures by the criterion R = 1: (1) F -, (2) HCOO -, (3) Cl -, (4) , (5) Br -, (6) , , (8) PO 3 F 2-, (9) at pH 10.82 and the following characteristics: separation column 180 × 3 mm; sorbent with the structure of the surface-grafted AS-14a ionite, mesh 5 ± 0.5 μm, thickness of spherical shell 140 nm, capacity 0.038 mEq/mL; ASRS-ULTRA2-mm suppressor electrodialyzer; flow rate of carbonate eluent 0.5 mL/min. [59]. Vertical lines correspond to two isocratic modes: I1, to the experiment presented in the DIONEX catalog; I2, to the optimum eluent under the specified conditions. The shaded rectangle (in the region forbidden for the isocratic mode with the criterion R = 1) corresponds to the borders of the weak gradient "lens" mode shown in Fig. 10.
, (8) PO 3 F 2-, (9) ] on an ICS-5000 system with column AS-14a. The system in the isocratic mode I2 (see caption to Fig. 9) can separate the mixture within as short as 9 min. with the elution order of the components presented above. SO − tions for the equilibrium and kinetic characteristics of gradient IC were given in [7,53]; examples of solutions of direct and inverse problems for systems of DIONEX ion chromatographs using both isocratic and gradient elution were presented. Further research [59,60] showed that the accuracy of the description of gradient chromatography can be improved. The level of refinement of the model of gradient chromatography in the IONCHROM software is rather high for the consideration of alternating gradient modes [61,62]. Figure 10 presents an example of a chromatogram calculated for the alternating "lens" mode. leading to the complete separation of a mixture of anions within threefold shorter time than in the best of isocratic mode in the same system.
Work [62] indicated a possibility of the creation of a universal algorithm for solving an extremely complex problem of the optimization of the gradient mode of ion chromatography and HPLC on an example of the so-called "weak" gradient mode.
All important components and logical relations of the developed Mathematical Chromatograph software were characterized in the presentation [63], and the main references, some developments and software are given on the project website https://www.grompiac.com. * * * The Mathematical Chromatograph project includes mathematical models of analytical methods of gas, liquid, and ion chromatography. A theoretical basis for the development of high-level models, requiring no adjustable parameters and allowing mathematical experiments in wide limits, is created. The authors believe that a combination of the presented abundant arsenal of developments in the field of modeling of methods of analytical chromatography in a uniform program has no fundamental hindrances and will require rather short time.

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit ttp://creativecommons. org/licenses/by/4.0/.