The collaborative effects between an applied orthogonal electrical field and the internal structure of polymer gels in gel electrophoresis is studied by using microscopic-based electrophoretic transport models that then are upscaled via the format of electro kinetics-hydrodynamics (EKHD). The interplay of the electrical field and internal gel morphology could impact the separation of biomolecules that, because of similar chemical properties, are usually difficult to separate. In this study, we focus on an irregular pore geometry of the polymer-gel structure by using an axially varying pore (i.e., an axially divergent section) and an orthogonal (to the main flow of solutes) applied electrical field. The microscopic-based conservation of species equation is formulated for the standard case of electrophoresis of charged particles within a geometrical domain, i.e., a pore, and upscaled to obtain macroscopic-based diffusion and mobility coefficients. These coefficients are then used in the calculation of the optimal time of separation to study the effect of the varying parameters of the pore structure under different values of the electrical field. The results are qualitatively consistent with those reported, in the literature, by using computational-based approaches as well as with experiments also reported in the literature, previously. The study shows the important collaborative effects between the applied electrical field and the internal geometry of the polymer gels that could lead to improving biomolecule separation in gel electrophoresis.

1. Introduction, Motivation, and Previous Efforts

The separation of biomacromolecules can be accomplished in various types of porous/fibrous media with electrophoresis, including polymer gels, microfluidic and field flow fractionation devices, and capillary electrophoresis, just to name a few [13]. Separation of these molecules has many potential applications, including the development of new pharmaceuticals, tissue scaffolding, and clinical diagnostic tests [4, 5]. In spite of excellent advances made among many of the electro-driven separation techniques, gel electrophoresis and its variations (2D gel electrophoresis, SDS-PAGE electrophoresis, temperature gradient gel electrophoresis, etc.) remains the widest used separation technique for proteomics, DNA/RNA fragments, and other macromolecules and it is a multimillion dollar business (Apogee Electrophoresis, 2014). Therefore, improvements in the separation efficiency will have a highly beneficial impact on the technology.

The use of an orthogonal field in gel electrophoresis can improve the separation of biomolecules. In particular, when the physicochemical properties of these biomolecules are similar and the separation can become very challenging, the presence of an orthogonal field is helpful. Sauer et al. [6] were the first contribution to explain the beneficial role of orthogonal fields found in experimental trends [7], from fundamental principles. Oyanader and Arce et al. [8] further identified the role of geometrical scales of the gel pore matrix in controlling the separation and, in addition, determining optimal time of separations. Pascal et al. [3] extended the analysis to electro-field flow fractionation in Couette-type separation devices.

In addition to the applied electrical field, another key aspect controlling the separation efficiency in electrophoresis is the internal structure, or morphology, of the gel media. Since gels are polymer matrices with pores of different sizes, shapes, and geometries and that only one of these is not enough to capture the separation behavior of the system [9, 10], it is important to gain a deeper understanding of the role that the gel pore geometry can play on the electrophoresis separation efficiency of macromolecules and their optimal time of separation. Early successful studies reported explained the potential impact of the geometry of the pores in the gel morphology, on electrophoresis separation, and included the use of a sudden change in the axial direction of the cross-sectional area of a pore domain of a gel medium [11]. Furthermore, a Monte Carlo simulation approach was used in either an isotropic or a nonisotropic gel morphology [12, 13] to determine the effective diffusivity of “point molecules” inside a model gel matrix. All of these studies show the crucial role that the pore-gel morphology could play in dictating the values of the effective transport parameters of biomolecules that, in turn, control their optimal time of separations.

As opposed to a sudden change in the cross-sectional area of flow in the gel-pore domain [11], Simhadri et al. [14] reported a computational analysis based on pores of rectangular geometry with an axially varying cross section. The current contribution examines the role of the same pore geometry, also with a rectangular diverging pore from an analytical point of view, i.e., a continuously axially varying cross section of the gel-pore domain (linear fashion) and its effect on the prediction of optimal separation times coupled with an orthogonal applied electrical field (constant). Through the use of the methodology known as “electrokinetic hydrodynamics,” or EKHD [15, 16], the relevant physics of the motion of the biomolecule inside the gel matrix is captured based on microscopic differential models of the gel-pore domain. Then, the model is upscaled to the macropore level by using the spatial-averaging method, a scaling technique, in conjunction with the fundamental principles of electrostatics and hydrodynamics. This approach is useful in obtaining the effective diffusivity and the effective mobility of the biomolecule (inside the gel matrix); then, these effective coefficients are used to compute optimal separation times as well as to assess the role of the axially varying morphology of the pore domain of the gel matrix on its separation performance. This can be accomplished by varying the various geometrical parameters of the diverging section, i.e., aspect ratio, angle of divergence, and axial position, and illustrating their impact on both the effective transport parameters as well as the optimal separation time as functions of the applied electrical field. One of the advantages of using the EKHD approach is that the universal or generalized dispersion coefficients (see Section 4) as well as the optimal separation time are explicit functions of the geometrical parameters controlling/affecting the morphology of the pore domain of the gel matrix and the applied orthogonal electrical field. This is an advantage with respect to computational simulations [14] where the relation of the effective transport coefficients, with geometrical parameters of the pore domain, is not explicit.

This study shows the important interplay among model scaling, pore morphology of the gel matrix, charge of the biomolecule, and the applied orthogonal electrical field in assessing or improving electrophoresis separations. In addition, the optimal times of separation that are obtained are relatively simple mathematical functions, and therefore, they lead to a tunable model that can predict behaviors of the system performance for a wide range of variables (i.e., magnitude of the applied orthogonal electrical field, angle of divergence, aspect ratio, axial position, etc.). The effects of these geometrical parameters as well as the molecule charges and the magnitude of the applied orthogonal electrical field on the optimal separation time are illustrated. Furthermore, a qualitative agreement between the predictions reported here with those published by Simhadri et al. [14] is observed. Both predictions (those communicated here and those by Simhadri et al. [14]) are in agreement with experimental evidences reported by Thompson et al. [17].

This contribution is organized as follows. First, a detailed introduction and reference to previous work on the subject is given in Section 1. Next, an overview of the related literature and remarks where this contribution differentiates from others is presented in Section 2. The formulation of the so-called fluid problem [15, 16] is introduced in Section 3 and its subsections on electrostatics and hydrodynamic aspects. The “solute problem” equations are developed in Section 4, and the upscaling of the model is performed here leading to the effective diffusivity and effective mobility. Asymptotic solutions and validation are discussed in Section 5. Parametric illustrations and discussions of results are presented in Section 6. Section 7 is a discussion of the implications of the finding for the practical aspects of the problems; in particular, calculations of the optimal time of separation and the effect of the different parameters are included. Finally, concluding remarks are given in Section 8 of the study.

2. Overview of Relevant Literature

Biomacromolecules, in general, can often have similar physical or transport properties, i.e., electrophoretic mobility and/or diffusivity, causing the separation to become very challenging in certain cases [18]. In addition to the applied electrical field, for a given mixture of biomacromolecules, the internal gel structure, or morphology of the gel media, is one key property that can be exploited to enhance the separation efficiency of the material or device. Gels for electrophoresis are fibrous or porous matrices that can expose the biomacromolecules to a “differentiating” medium allowing the physical properties of the molecules, i.e., charge, size, shape, etc., to yield a unique motion for biomacromolecules. This unique motion is controlled by two different transport properties, i.e., effective diffusivity and effective mobility of the molecules inside the gel media. There are few experimental techniques that researchers have used to modify the morphology of the gel matrix. The traditional one is to change the concentration of the polymer cross-linker (see, for a review, Stellwagen [19]). A novel one is to use templating agents both within the polymer-chemical structure (molecular imprinting, see for a review Byrne and Salian [20]) and to manufacture gels by using molecular templates such as DNA and xanthan. For instance, Dharia et al. [21] and Rill et al. [22] proposed a modified architecture or morphology, by templating gels with DNA (or surfactant) templates that were removed from the gels once the polymerization took place, by electrophoresis action. The resulting material morphology, hypothesized as a “dual” porous structure [11, 13, 23, 24] displaying both small (of the order of the interfiber voids) and large (of the order of the templating agent) pores or voids, showed an improved efficiency in the electrophoretic separation of proteins. The most recent approach to modifying the morphology of the separation media is to add nanoparticles of different characteristics, i.e., nanocomposite gels, that could affect the separation efficiency (see, for example Thompson et al. [17, 23]). This subject has been reviewed by Simhadri et al. [25]; see also Texter [26].

Structurally, gel matrices can be viewed as consisting of pores of different geometries, sizes, shapes, etc. This situation leads to the fact that, realistically, only one idealized geometry, i.e., rectangular, cylindrical, or annular, for instance, is not enough to capture the internal architecture of the media [9, 10]. Instead, a collection of these may represent different limiting physical situations that can, collectively, aid to understanding the transport and separations of biomacromolecules. In addition, and based on the use of relatively simple geometrical models, Trinh et al. [1113] were able to show that the geometry of the pores or capillary channels could play a significant role in the determination of the effective transport properties of biomacromolecules.

The use of pores or capillary channels with axially varying cross sections is not uncommon in studies of transport phenomena. For example, the effect of an axially varying cross section on the transport of solute in a cylindrical capillary was studied by Horner [27] in order to predict the behavior of atherosclerosis in blood vessels. In this study, the lubrication approximation for the hydrodynamics in conjunction with the spatial averaging approach (Cwirko and Carbonell [28]; Paine et al. [29]; Slattery [30]; Whitaker [31]; Zanotti and Carbonell [3234]) was used in order to obtain both effective transport parameters as well as concentration profiles of the solute within the artery. Xuan et al. [35] further studied, experimentally, the effect of pressure on the movement of polystyrene particles in a converging diverging microchannel. The ratio of the particle velocity in the throat to that in a straight channel was studied in terms of pressure gradient, particle size, and particle trajectory.

Several contributions have been reported in the literature related to the electrokinetic transport in rectangular channels of varying cross section for applications in bioseparations and microfluidics. For instance, Ghosal [36] studied the transport of a polyelectrolyte through a nanopore with variable axial cross section and examined the effects of an (axially) applied electrical field. Through the use of electrostatics and hydrodynamics within a continuum mechanics framework, expressions for the electrophoretic speed of the polyelectrolyte, representative of a typical strand of DNA often used in gel electrophoresis, were determined. Four different modes of transport of the polyelectrolyte were found as a function of the dimensionless height of the channel. It was noted that the use of continuum hydrodynamics with water as the buffer is a reasonable approach that is valid for pores of sizes as low as tenths of a nanometer. In addition, Ghosal [37] also examined electroosmotic flow through a microfluidic channel of varying axial cross section and arbitrary zeta potential distribution. Supported by the use of Green’s function associated with the problem, a number of asymptotic cases for the electrohydrodynamic profile were studied as a function of the aspect ratio, , including cases for the long channel approximation (). Although many of these studies are more relevant for electrokinetic-based separations than for electrophoresis, one general conclusion is that the variable axial cross section increases mixing, worsening the resolution for an electro-based separation.

Opposite to the observation made above, Ross [38] examined electrokinetic separations in a diverging microchannel and found, for a specific case, that the diverging channel could actually lead to a better separation than a straight channel could. Yariv and Dorfman [39] studied the electrophoretic transport in rectangular channels of periodically varying cross sections. Macrotransport theory, as described by Brenner and Edwards [1], was used to obtain expressions for the effective velocity and effective dispersion within the channel. The analysis yielded asymptotic results for the macroscopic transport parameters (i.e., effective mobility and effective diffusivity) for both small and large Peclet numbers. Xuan et al. [35] examined the electrophoretic motion and separation of particles within a converging-diverging microchannel fabricated with poly(dimethylsiloxane) chips. The dependence of the separation on electrical field strength, particle size, and morphology was studied by using the ratio of the velocity of the particle at the throat of the device to that in a straight microchannel as the main parameter. Ghosal [37] presented an analysis related to channels of varying cross sections; he examined that electroosmotic flow and sophisticated mathematical techniques were employed to obtain volume fluxes for various geometries and zeta potential distributions. In addition, Yariv and Dorfman [39] also computed expressions for dispersion and effective velocity but only for very large and small Peclet number values. None of these contributions, however, have reported any study of the role of the orthogonal electrical fields in conjunction with the characteristics of the pore or the separation domain nor have they shown any calculation of the optimal time of separation.

In this contribution, and for the particular case of gel-based electrophoretic applications, attention is given to the analysis of orthogonal (The word “orthogonal” here implies perpendicular to the axial (varying) direction of the capillary channel as it was used in Oyanader and Arce [8].) applied electric fields to potentially enhance the separation. As mentioned in Section 1, above, Sauer et al. [7] studied this effect on the transport of solute in a Poiseuille flow regime using a capillary channel or pore of rectangular geometry. Before the work reported by Sauer et al. [6], only ad hoc expressions were known for dispersion under the influence of an applied electrical field (Eringen and Maugin [40]; Hiemenz [41]; Hunter [42]) and the effects of an orthogonal electrical field were only known experimentally (Lee et al. [7]; see also Gajdos and Brenner [43], Heller et al. [44], and Hiemenz [41]). Sauer et al. [6] introduced a mathematically nontrivial and fundamental modeling approach as a tool for the analysis of the system under study. This effort was the first that systematically explained experimental trends (Lee et al. [7]) in enhancing separation efficiency when an orthogonal field is applied in electrophoretic separations. One limitation of the analysis is that the predicting formulae for the effective transport parameters are in terms of integrals that require numerical solution. Motivated by this limitation, and attempting to further apply and/or develop the methodology, Oyanader and Arce [8] have revisited this contribution and have derived simpler and mathematically friendlier expressions for the effective transport parameters that clearly show the dependence on the orthogonal electrical field and the scaling. More specifically, Oyanader and Arce [8] have shown that these expressions depend critically on the scaling of the geometrical parameters of the device. The present analysis is concerned with the effect of an applied orthogonal electrical field on a rectangular pore with (continuously) axially varying cross sections. To the best of our knowledge, this is the first analytical-based contribution in this subject. A computational-based approach was reported by Simhadri et al. [14].

In general, the effects of transverse or orthogonal applied electrical fields have not been extensively studied regarding the separation of biomacromolecules, except for the case of electro-FFF (see, for example, Pascal et al. [3]). However, Baldessari and Santiago [45] in a review related to electrophoresis in nanochannels stated that the presence of the thick electrical double layer (EDL) in nanochannels can produce transverse electrical fields. Because of this observation, the dispersion in nanochannels must then be a function of the applied orthogonal electrical field. Therefore, in addition to the previous work by Sauer et al. [6] and Oyanader and Arce [8], Baldessari and Santiago’s [45] work suggests that further examining the role of applied orthogonal electrical fields in the separation of biomacromolecules in various systems is important.

3. Problem Formulation: General Framework and Methodology

The system under study in this contribution is depicted in Figure 1. This figure is a general representation of the different scales presented in a gel-electrophoresis system. The transport associated with the motion of biomolecules in gel electrophoresis is a multiscale problem: level 1 (see Figure 1) describes the macroscopic scale in the system where in general, pores of different geometries, sizes, shapes, etc., are located at a smaller scale. Level 2 in the figure shows the pore-scale domain or microscopic scale where for example, an axially varying cross-section pore domain can be easily located in the gel matrix. In general, these domains could be oriented in many different ways and have different sizes and shapes. In this contribution, we are interested in a single type of pore orientation with potentially different shapes, all controlled by the so-called “shape factor” (). Level 2 is where microscopic transport phenomena based on continuum mechanics can be described. Level 3 is, therefore, a single-pore domain (representing all the pores of the gel matrix (A distribution of pore orientation can be used to analyze the effect of the microdomain on the macroscale level transport. However, this is beyond the scope of the present work.)) with the Cartesian coordinate system and other important dimensions anchored to the domain.

In level 3 of Figure 1, the pore’s geometry is described by Cartesian coordinates with length and height , with a variable axial cross section (i.e., a diverging pore domain), the degree to which is controlled by the angle . Since the cross section is not uniform, the fluid motion is in both the axial and orthogonal directions; however, only the flow in the axial direction is relevant for the spatial averaging approach used in this study (Cwirko and Carbonell [28]; Paine et al. [29]; Slattery [30]; Whitaker [31]; Zanotti and Carbonell [3234]). An electrical field is applied in the orthogonal direction (Oyanader and Arce [8]), denoted by the potential, , on both the upper and lower pore walls. It is also shown in Figure 1 that the walls of the capillary can be mathematically described by the equation of a line, , as follows: where is the slope of the line and is the -intercept. Equation (1) above can be expressed in dimensionless form as follows: where , , and , which is the aspect ratio or “shape factor” of the pore.

The key assumptions for conducting this analysis are as follows: (i) the buffer solution is Newtonian; (ii) the flow is incompressible; (iii) the system is under steady state; (iv) the capillary walls vary slowly (axially) so that the lubrication approximation can be applied to the system (Probstein [4]); (v) no significant charges are present on the walls (i.e., this is an electrophoretic process) and; (vi) there is long channel approximation (The value of will be parametrically varied within the range of the approximation.) ().

The assumptions above will be employed to develop microscopic governing equations for the “fluid problem,” i.e., the electrostatics and hydrodynamics of the buffer solution flow, through a systematic approach called “electrokinetic hydrodynamics,” or EKHD, that links the two domains in the system under study (i.e., “solute” and “fluid” or buffer domains) ([15, 16]). EKHD emphasizes scaling and involves a series of steps which can be used to obtain upscaled model equations for the electrostatics and hydrodynamics of the system (fluid domain) and by utilizing the spatial averaging technique to connect them with upscaled microscopic models such as the species continuity equation (solute domain) (Cwirko and Carbonell [28]; Paine et al. [29]; Slattery [30]; Whitaker [31]; Zanotti and Carbonell [3234]). From this method, constitutive equations for the effective transport coefficients (i.e., effective mobility and effective diffusivity) can be analytically obtained in order to predict the macroscopic or pore-level behavior of the system (A material or gel scale is also possible (see Simhadri [46]); however, in this contribution, we focus only on the pore-level scale. The material scale is equivalent to the “pellet scale” in catalytic systems (see Arce et al. [47]).). These coefficients are a function of fundamental transport parameters, such as the molecular diffusivity, molecular charges, and species electrophoretic mobility. Once these two parameters, i.e., the effective mobility and effective diffusivity, are computed, the optimal separation time for, for example, the two ideal biomacromolecules, or, in general, for a mixture can be determined. This allows for the optimization of the system with respect to the gel morphology shown in Figure 1 (see Simhadri et al. [48]).

In summary, by following the EKHD approach, once the electrostatics and hydrodynamic (formulation and scaling) aspects have been taken care of, and in order to aid in extracting information from the (solute/analyte) species continuity equation (Bird et al. [49]), without actually solving it, the method of spatial (area) averaging is used (Cwirko and Carbonell [28]; Paine et al. [29]; Slattery [30]; Whitaker [31]; Zanotti and Carbonell [3234]). In particular, this methodology states that any variable of interest, , can be decomposed into a simple addition of two expressions, i.e., the cross-sectional area average, , and a “deviation” function from such an average value, . Equation (3), below, mathematically expresses this definition while equation (4) further characterizes the mathematical operation needed to obtain the area average of the variable under study: where and where is any variable of interest such as concentration, velocity, electrical potential, or temperature. It has to be noted that the term is only a function of the axial coordinate “” and remains, in general, as a function of all spatial coordinates of the system, i.e., “”, “,” and “.

Because the governing equations that describe the fluid domain and solute domain are coupled sequentially, the area-averaging approach on the solute domain can be used to reduce the number of independent variables through an upscaling analysis. In particular, the use of the methodology avoids completely solving the species continuity equation; instead, the approach yields “effective” coefficients that help to predict certain macropore-level performance information (such as the optimal time of separation) as a function of the different conditions of the operating parameters (electrical field strength, divergent angle, pore-domain shape, etc.). In the sections below, the analysis of the governing equation for the electrical field as well as the coupling between this variable and the species continuity equation are studied. Finally, the use of the spatial area-averaging approach also helps in determining effective transport parameters in “closed-form” mathematical equations that are handy in predicting the behavior of the system performance under a variety of parameter values. This type of analytical equations, for example, cannot be accomplished by other methods such as the Monte Carlo simulation where an algorithmic approach must be used every time that a system parameter is changed (Trinh et al. [12]).

4. Electrostatics

The gel-pore domain depicted in Figure 1 (see level 3) is undergoing electrophoresis; such a domain has an applied electrical field that is governed by the principle of conservation of charges. The pore domain is represented by using Cartesian coordinates. Under the general assumptions typically used in gel electrophoresis [6, 49], the applied electrical field potential for the system shown in Figure 1 is mathematically described by the Laplacian of such an electrical potential. where “” is the applied electrical field potential, “” is the axial coordinate, and “” is the orthogonal coordinate to the channel axial direction, . By neglecting entrance effects as suggested by the analysis described in Oyanader and Arce [8], the orthogonal surface potential, , dominates the variations of the electrical fields, , and such a field may be assumed as , explicitly. In addition, Oyanader and Arce [8] showed that for cases where the geometrical ratio satisfies , equation (5) simplifies for all practical cases to the following (The simplification is typical of an order of magnitude analysis in Laplacian equations describing heat and mass transfer processes as well (see Whitaker [50]).):

Thus, only the orthogonal field direction is important for such a case. Boundary conditions for the orthogonal direction are as follows: where is the value of the applied electrical potential at the top surface of the pore domain.

In order to study the influence of the electrical field on the system, a solution of equation (6) is needed. Furthermore, this solution must be averaged in the cross section of the capillary to obtain the averaged applied electrical field. This information is needed in the solute transport equation discussed in Section 4, below. Equations (6) & (7) produce the following analytical solution: where , , and are dimensionless variables defined previously. Equation (8) shows that the point, i.e., the “local” applied electrical field potential, , is a function of both , the dimensionless axial coordinate, and , the dimensionless orthogonal coordinate, due to the axially varying cross section shown in Figure 1. By applying the area-averaging definition (see equation (4)), the area-averaged applied electrical field, , is simply

The (spatial) area-averaged applied electrical field, , can be related to the point or local potential, , by the decomposition technique (Gray [2]) to produce where the variable is the spatial deviation field defined as

In order to predict the effective transport parameters of the system, equations (8), (9), (10), and (11) are useful in solving the species continuity equation as shown in Section 5.

5. Hydrodynamics

The next aspect in studying the convective-diffusive transport of solute/analyte in capillary channels is the hydrodynamics of the buffer solution. As shown in Figure 1, the flow regime under study is laminar- and pressure-driven. Due to the presence of the variable axial cross section, the lubrication approximation becomes useful to obtain an analytical velocity profile. Four requirements, or assumptions, are needed in order to utilize such an approximation. They include that the pressure is a function of the axial coordinate, , only, the velocity component in the axial direction, , is a function of only, and the convective terms in the Navier-Stokes equations are negligible. The first two conditions are a consequence of the long channel approximation, and the last two conditions are a consequence of low Reynolds number flow, as it has been clearly reported in the literature [51, 52]. The following three equations, i.e., the continuity equation and the Navier-Stokes equations for the velocity components in the - and -directions, will be simplified through the lubrication approximation to produce a differential equation for the velocity in the -direction as a function of the orthogonal coordinate, . where and are velocity in the - and -directions, respectively, is the density of the fluid, is the viscosity, and is the pressure. Once the above equations have been identified (Here, the equations are the usual ones for the case of viscous flows since the electrokinetic-based forces are negligible for the buffer solution for the case of electrophoretic transport (Arce et al. [53]).), the lubrication approximation can be used and requires defining the following dimensionless variables [51, 52]: where , , and are characteristic values for the pressure, velocity in the -direction, and velocity in the -direction, respectively. After equations (15), (16), & (17) are substituted into equations (12), (13), and (14), and the assumptions required for the lubrication approximation (stated previously) are employed, the final differential equation for the velocity field in the -direction as a function of is shown below. (As previously noted, the effect of the velocity in the -direction is not relevant for the purposes of this analysis.).

Equation (18) has the usual no-slip boundary conditions at each wall, expressed as follows:

Equations (18) and (19) produce the following solution for the velocity profile for the system depicted in Figure 1:

Equation (20) is derived under the assumptions that there are no charged walls, the fluid is Newtonian, the flow is incompressible, and the system is under steady-state condition. Next, the pressure gradient as a function of the axial coordinate, , needs to be determined. This can be accomplished by employing the definition of the volumetric flow rate (Bird [49]):

By using equations (20) and (21), the pressure gradient as a function of can be found in terms of the volumetric flow rate as shown below:

Substituting equation (22) into equation (20) produces the following dimensionless velocity profile for the system shown in Figure 1: where and is the velocity evaluated at and .

From equation (23) and by using the spatial area averaging method, as used in the electrostatics section, one can compute the area-averaged velocity for the flow to be

The area-averaged hydrodynamic velocity indicated in equation (24) is also a function of the axial coordinate, , as opposed to a constant value in the straight pore domain. Similarly, the deviation equation for the hydrodynamic velocity field is given by

The above equations will be useful along with the electrostatics equations from Section 3, to study the convective-diffusive transport of solutes in gel-pore domains by performing the spatial averaging, which will be outlined in the next section, below.

6. Convective-Diffusive Molar Transport Equation: Formulation and Analysis

Once the expressions for the area-averaged and deviation of the electrical potential and the hydrodynamic velocity profile of the buffer solution have been computed, as suggested by the EKHD approach [4, 52], it is necessary to link the solute/analyte transport aspects with both the electrostatics (The electrostatics is important in this analysis since the solutes/analytes are charged particles [4]. In particular, it is useful in the computation of the electromigration term of the species continuity equation.) and the hydrodynamics of the buffer solution. This will be accomplished through the use of the spatial averaging approach as applied to the molar species continuity equation as in Sauer et al. [6]. This approach yields expressions for the effective mobility and effective diffusivity which are required for the determination of the optimal separation time, . The steps to obtain these effective parameters are the focus of this section, and the optimal separation time is discussed in detail in Section 6.

The first step in the analysis of the system involves the identification of the species continuity equation on a molar basis as well as the wall boundary conditions associated with such an equation. From Bird et al. [49], the following equation is written for the system under study (see Figure 1, level 3): where the notation, , is used to denote the cation valence in a given buffer solution, is the molar concentration of the species under study, is the axial mass average velocity of the fluid, is the dilute solution diffusion coefficient of the species, is the electrophoretic mobility of the species, and is the applied electrical field potential. Equation (26) has the following boundary conditions at the walls of the channel, implying that the walls are impermeable (Since the spatial area-averaging method is being used, only boundary conditions in the orthogonal direction are specified.):

By applying the definition of area-averaging to equation (26), the following differential expression is obtained (for details, please see Sauer et al. [6]):

Equation (29) requires the information related to both the deviation function of the velocity profile, , and the deviation of the applied electrical field, . These functions were derived in Sections 3 and 4 of this contribution, and they are given by equations (11) and (25), respectively. In addition, equation (29) requires a solution that yields as a function of the spatial area-averaged concentration variable, . This solution is obtained by using a closure approach such as the one proposed by Cwirko and Carbonell [28] and Paine [29]. Their approach seeks a solution of the following form: where and are two algebraic functions of the basic parameters of the system. In the system under study, however, and will be implicit functions of the axial coordinate, , also. In terms of nondimensional variables, and are computed from the differential model and boundary conditions to be where

Equation (30) is frequently referred to as the “closure” equation of the spatial area-averaging procedure, and in its derivation the assumption of a long capillary, i.e., large value of the geometrical ratio, is used [6, 31]. This is not a restriction in this analysis since the electrical field was computed for a case where (Oyanader and Arce [8]). In particular, the nondimensional applied field () in the orthogonal direction of the flow is the most relevant electrical field parameter of the equation in addition to geometrical dimensions. Based on the spatial/area-average approach [31], the general expressions for the effective transport parameters are given by

Furthermore, now equations (33) & (34) can be used to obtain the dimensionless effective parameters, and , defined as follows: where is the Peclet number of the system based on the averaged (axial) velocity as identified in equation (24). The function identified in equation (35) is referred to sometimes as the nondimensional axial dispersion of the capillary and has the advantage of being independent of the Peclet number. In short, it is convenient to view as the “universal dispersion” function for every value of the applied orthogonal electrical field since it includes all possible values of the Peclet number as a scaling factor (Mathematically, is a self-similar function with Pe, the self-similar parameter of the system. Sometimes, is also term as the generalized dispersion coefficient.). Please note that an increase in implies a decrease in and vice versa. Both and reflect the combined effects of the hydrodynamics of the buffer and the orthogonal electrical field applied to the system, under the assumptions stated, on the macrotransport parameters of the system. Both the characteristics of the electrical and hydrodynamic fields are included in these equations. In addition, geometrical dimensions and geometrical scaling for the electrical field (long channel approximation) and hydrodynamics (lubrication approximation) have been incorporated in this formulation.

Finally, by replacing the expressions of the variables identified in equations (31) & (32) and completing the spatial averaging indicated, the following analytical functions of the dimensionless effective parameters, and , are found:

It can be observed that the above expression for is an explicit function of the line, , representing the walls of the diverging capillary as well as the cation valence, , and the dimensionless applied orthogonal electric field, , whereas the expression is not a function of . This observation is useful to study the roles that these parameters play on the universal dispersion and effective velocity. This analysis is performed below.

7. Asymptotic Solutions: Validation of the Results for the Effective Transport Coefficients

Based on the equation for the effective dispersion coefficient, see Equation (37) (As and approach zero, the limit of is -2/105, the Taylor-Aris result for purely Poiseuille flow (see Probstein [4]).), the following relation can be written: where

Alternatively, Equation (40) can be written as

It is clear that for the case of , the relation reduces to the equation published by Oyanader and Arce [8]. In short, the effective coefficient reproduces the asymptotic value for the strait pore geometry. In addition, equation (38) (Please note that as approaches zero, the limit of is unity.) can be written in the exponential form as the following:

By comparison, Equation (42) with Equation (22b) from Oyanader and Arce [49] leads to the conclusion that they are exactly the same. Again, the straight-pore geometry result is reproduced for the asymptotic case of . In conclusion, the fact that the limiting case for the straight-pore domain () coincides exactly with the results reported in the previous work indicates a validating case for the general situation at hand in this contribution, i.e., .

8. Parametric Analysis of the Effective Transport Coefficients and Discussion

This section includes graphical illustrations for the analytical results found in the previous section of this contribution. These illustrations show the dependence of the effective parameters, , the universal dispersion, and , the nondimensional effective mobility, on the geometrical and physicochemical properties of the system. Specifically, in this section, the direct dependence of these two effective transport coefficients with key parameters dictating the geometry of the pore (i.e. , the aspect ratio, , the angle of divergence, and , the dimensionless axial position), as well as physicochemical parameters such as the molecular charges (), and the applied nondimensional orthogonal electrical field, , will de illustrated. However, the importance of the role that these parameters could play on determining the optimal separation time of biomacromolecules will be discussed in detail in Section 7.

Figures 2(a) and 2(b) show the influence of the molecular charges on the effective mobility coefficient and on the universal dispersion, with respect to the applied orthogonal field, , respectively. In the plots, the values for have been obtained for the values of the parameters within the range selected previously by Sauer [6]. Figure 2(a) shows the variation of the effective transport coefficient,, with respect to , the dimensionless applied orthogonal electrical field, for various values of the valence, . It can be observed that as the magnitude of the electrical field increases, the effective velocity asymptotes to a value of zero, whereas, when tends to a value of zero (i.e., no applied orthogonal electrical field), the effective velocity approaches a value of one. These results imply that as the magnitude of the electrical field increases, the molecules within the system will approach the velocities of the two pore walls (i.e., zero). The degree to which the particles approach the walls can also be deduced from Figure 2(b) as shown by the effect of the valence.

In Figure 2(b), the universal dispersion coefficient is plotted as a function of for different values of valence, , and selected values of , , and . It can be seen that a maximum in dispersion is reached when the applied electrical field approaches zero, and the dispersion approaches a value of zero as the magnitude of the applied electrical field increases. As the magnitude of the charges increases (negatively), the shape of the bell curve of the dispersion gets narrower implying that a given value of the function can be obtained for smaller values of the applied orthogonal filed, .

Figures 3(a) and 3(b) represent the impact on the universal dispersion, , of two geometrical parameters, the nondimensional axial position within the pore domain, , and the shape factor () of such a pore domain. Clearly, positions closer to the pore entrance () show less values in the magnitude of the parameter () for values of , , and . This behavior shows a decrease in the nondimensional effective diffusion coefficient, . Figure 3(b) shows the parametric effect of the shape factor, , for the same value of the divergence angle and molecular charges as before and at . As the values of the shape factor decreases (longer pore domains), the values of () increases implying an increase in the nondimensional effective diffusion coefficient, , related to the pore domain.

Figures 4 illustrates the impact of the divergence angle, , on the universal dispersion coefficient, , with respect to the nondimensional applied orthogonal field, , parametrically, for different values of the divergence angle, , that have been selected; i.e., , , , and . All the values of the divergence angle have been selected such as less than to comply with the restriction of the lubrication approximation. Values of the shape factor are kept at and at a position inside the pore domain with a value of the molecular charge or valence of . The general trends are that (absolute) values of the dispersion coefficient decreases when the values of the divergent angle also decrease. This implies that the effective dispersion decreases as the pore tends to be more strait in their shape. The results are consistent (qualitatively) with those published by Simhadri [14]. This type of results is a clear indication of the dramatic role played by the geometry of the domain on the values of the effective transport coefficients that control the motion of the macromolecules inside the pores of the gel matrix.

9. Practical Implications of the Results: Optimal Time of Separation

The impact on the practical aspects of the results, obtained in Section 4 and illustrated in Section 6, deserves illustration to gain a better understanding of the role played by the irregular shape of the pore geometry on the electrophoresis; and an attempt to illustrate some of their useful implications for the separation of biomacromolecules is included in this section. The focus of this discussion is on the potential effects of the orthogonal electrical field and the (axially varying) morphology of the capillary, on the time required for optimal separation of two solutes with given physicochemical properties. In this regard, one of the most difficult examples of separation is the case of two species of approximately equal molecular diffusivities. This specific case is an ideal candidate for applying orthogonal electric fields in order to take advantage of other physicochemical properties (i.e., electrophoretic mobility, , and cation valence, ) or the morphology of the gel. Sauer and collaborators [6] and Oyanader and Arce [8] calculated the time required to achieve an “optimal” separation (resolution of two) to illustrate practical applications of applied electric fields. Following this approach, analytical expressions for effective transport parameters in order to estimate the optimal time of separation of two hypothetical species based on various parameters, including the angle of divergence, , aspect ratio, , and dimensionless axial position, , that characterize the pore geometry under analysis, have been computed. In addition, since the equations predicting the effective parameters show a direct dependence on the orthogonal electrical field, , the effect of this parameter on the optimal time, , for the separation of biomacromolecules is determined.

Applying the definition of the dimensionless orthogonal electric field (please see equation (32)), , to two hypothetical species with very similar molecular diffusivity values, but different electrophoretic mobilities, the following expression can be utilized to relate the effects of between the two species: where the term has been defined as the electrophoretic mobility ratio. Equation (43) indicates that for a given value of applied electric field to the system, species B is affected by a dimensionless electric field proportional to the electric field affecting species A, , by a factor of the mobility ratio . Another important parameter for analysis of the optimal separation time variation with is the cation valence ratio. For two species, similarly to equation (43), we have defined the following valence ratio: .

The time of separation of two species, for a resolution of 2, has been defined by the following expression (see Sauer et al. [6]): where all the terms have been previously identified and/or defined (see equations (33) and (34)) as functions of the dimensionless orthogonal field, , and the subscripts indicate the species for which the term must be computed, either using or . Please note that and are functions of the valence of each cation species, so that and may be specified to compute the optimal separation time,. Furthermore, the combination of equations (43) and (44) can be used to analyze the separation process under electric field stress.

Figures 5(a)5(d) show the (optimal) separation time, , as a function of the nondimensional orthogonal field, , taking species A as the reference. In particular, Figure 5(a) presents as a function of for different values of the valence ratio, , and Figure 5(b) shows the same variables for different values of the angle of divergence, , for fixed values of the other parameters (i.e., mobility ratio, valence ratio, aspect ratio, and dimensionless axial position). Figures 5(c) & 5(d) show the same function for different and dimensionless axial position values,, respectively. It can be seen from Figure 5(a) that a minimum optimal separation time, , exists for values of . This minimum is further influenced by the value of the valence ratio, of the two species, A and B. It is evident that as increases, the time for separation decreases, due to the increase in the difference in susceptibility of the two species to the applied orthogonal electrical field. Therefore, if the two species have vastly different charges, or valences, the model predictions indicate that they will be easier to separate under the action of an orthogonal electrical field.

As can be observed from Figure 5(b), optimal separation times decrease with decreasing values of the angle of divergence of the pore, . In particular, for , the optimal separation time is much lower than for or (Although the lubrication approximation limits the validity of the approach to small angles (less than 20°), computational simulations show that the most dramatic changes for the system’s behavior occur within the range of 0°-20° (Simhadri et al. [14]).). This result implies that a straight pore will create a better separation than a pore with a varying axial cross section will, since there will be less mixing of the two species. It is evident from Figure 5(c) that as the aspect ratio, , increases, the optimal separation time decreases. A large value of implies that the pore is shorter than it is long, whereas a small value of means that the pore is long. Therefore, from Figure 4(c), a shorter pore will lead to a lower optimal separation time.

Lastly, Figure 5(d) illustrates the effect of the dimensionless axial position within the pore on the optimal separation time. It can be seen from this figure that for low values of (i.e., near the entrance of the pore), the optimal separation time decreases, and for high values of (i.e., near the exit of the pore), the optimal separation time increases. This is due (most likely) to the fact that near the entrance of the pore (low values of ), the species do not have enough time to mix and therefore will be significantly easier to separate with an applied orthogonal electrical field, than if they were at a position near the exit of the pore. From the results discussed in this section and within the assumptions of the analysis, it can be concluded that a straight, short capillary will produce a lower optimal separation time. This result is indicative of the potentially important interplay of scaling, capillary morphology, and the direction of the applied field in order to improve separation. Moreover, an optimal operation requires the consideration (most likely) of all of them together as opposed to the “in isolation” fashion customarily used. The prediction of the optimal time of separation is consistent with those reported by Simhadri et al. [14] and with experimental evidences reported by Thompson et al. [17]. A complete scaling analysis and a quantitative comparison are outside the scope of this contribution; however, it will be the subject matter of the future communication.

10. Summary and Concluding Remarks

For the case of gel electrophoresis of charged solutes in gels modeled as a pore bundle, with uniform pore of an irregular shape, the results reported here clearly show the explicit dependence of the effective diffusivity and the optimal separation time (within the assumptions of the analysis) on the orthogonal field as well as the geometrical parameters of the capillary, including the aspect ratio, the axial position, and the angle of divergence. The case of electrophoresis separation of two solutes under an applied, orthogonal, and constant field is analyzed for a diverging pore cavity. These types of geometrical voids are hypothesized as possibly being housed in a nanocomposite gel [17, 53]. Expressions for the effective system parameters have been derived by using the EKHD method (Tijaro-Rojas et al. [17]) without actually solving for the concentration profiles. Taking advantage of this methodology, constitutive equations of the effective, or alternatively “macrotransport” transport parameters, i.e., the effective diffusivity or, even more general, the universal dispersion parameter, and the effective velocity, are obtained as a function of fundamental variables such as the magnitude of the electrical field, cation valence, angle of divergence, aspect ratio, and axial position (This type of functionality is, in general, tied to the characteristics of the electrical and hydrodynamic fields and our claims are by no means valid outside the specific assumptions and characteristics used in this analysis. The variation of the effective transport parameters with the location is consistent with analysis conducted for other geometries (see Nagarkar et al. [54]).). These functions are useful in identifying conditions for minimum values of the optimal separation time for two solutes. The potential benefits of the predicting equations for the design of optimal conditions for the electrophoresis separations of biomacromolecules are clearly shown. For a gel with diverging pore geometry, and to the best of our knowledge, these analytical results are not available in the current literature and it is the hope of the authors that they contribute to the better understanding of the fundamental role of morphology and the orthogonal field in electrophoretic separations with relevance for gels, with a modified morphology, microcapillary electrophoresis, and other similar cases.

Finally, and since the EKHD method uses concepts from continuum mechanics approaches, the results and observations of this contribution (most likely) will be also valid for the microchannel or microcapillary scales in separation devices such as microfluidics. However, continuum mechanics approaches have been validated for the nanoscale in contributions reported in the literature. For example, Conlisk [55] conducted an analysis comparing the validity of the Debye-Hückel approximation in microchannels (i.e., channel height is greater than electrical double-layer thickness) and nanochannels (i.e., channel height is on the order of the electrical double-layer thickness). The results using the continuum approach compare well with experimental data for rectangular channels. On the other hand, it was reported that the continuum approach is not valid for nanochannels with a height of 6 nm or less. However, for most separations conducted in nanochannels, the height is at least 10 nm; therefore, the continuum assumptions should be valid. We believe that the results presented here could be a very useful first approximation for understanding the role of the morphology of the separation media beyond the microscopic scale and, possible, give a qualitative idea at behaviors of the media also valid for smaller scales.


Jennifer A. Pascal’s present address is Department of Chemical and Biomolecular Engineering, University of Connecticut, USA. Preliminary versions of the results were presented at the American Electrophoresis Society (AES) Annual Meetings.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors are indebted to helpful discussions with Dr. Sharon G. Sauer of the Rose-Hulman Institute of Technology, Terre-Hautte, IN, and Dr. Y. Liu of the Department of Mathematics, TTU. Jennifer Pascal is grateful for a Diversity Fellowship from the Office of Research and Graduate Studies at Tennessee Technological University.