^{1}

^{2}

^{1}

^{*}

We show that the processes described by Avrami functions are self-similar. A comparative function characterizes a self-similar process by a certain Avrami exponent. We define the self-similar categories of some well-known biological processes. The method to determine the Avrami exponent by choosing the comparative function is demonstrated on the diffusion model of the growth of nuclei. We generalize the results.

The self-organizing behavior of the materials is a well-known and widely investigated topic in research [

The Avrami equation could serve as a mathematical model of different biological processes [

Our objective in this present work is to derive the universality of the Avrami- equation.

Self-similarity means the scaling possibility of all the measurable quantities of the studied object. Let us assume that r is the linear size of the scale and L is a scale-dependent property: L ( r ) . By lengthening the scale by a times, the property transforms as

L ( a r ) = k L ( r ) (1)

where the factor k depends only on a . The power type functional relationship is valid in the biology and the fractal theory:

L ( r ) = A r α (2)

which changes in the form of

L ( a r ) = A ( a r ) α = a α L ( r ) (3)

When Equation (2) is converted to the log-log scale, the result is a straight line with a gradient equal to the slope.

lg L ( r ) = lg A + α lg r (4)

The slope usually verifies the self-similarity of fractals.

Definition of the self-similarity in (1) means that the relationship is valid for any a and k independent from r . Hence k is the general function of a , expressing (1) in the form of

L ( a r ) = k ( a ) L ( r ) (5)

From this, we get a new and more general definition of self-similarity, namely

L : r ∈ I r ↦ L ( r ) ∈ R 1 (6)

the mapping is self-similar, if L ( a r ) may be expressed in a more general form where the scale-transformed function can be separated in the form of a product:

L ( a r ) = k ( a ) G ( r ) (7)

Solving the Equation (7), we suppose that a = 1 . Then,

L ( r ) = k ( 1 ) G ( r ) = c G ( r ) (8)

where c = k ( 1 ) is a constant. Consequently, L ( r ) and G ( r ) differ only in the fixed multiplicator. By the last property the Equation (7) can be rewritten:

L ( a r ) = k ( a ) k ( 1 ) L ( r ) = η ( a ) L ( r ) [ η ( a ) = k ( a ) k ( 1 ) ] (9)

Repeating the same procedure for the case of r = 1 , the following relationship is valid for any u , v variables:

L ( u v r ) = η ( u v ) L ( r ) = η ( u ) η ( v ) L ( r ) (10)

Therefore, the function η ( u v ) in Equation (10) is also a self-similar function:

η ( u v ) = η ( u ) η ( v ) (11)

The general solution of this function-equation is difficult. We may choose a = r , r = 1 for the common part of variables from the relationship Equation (9):

L ( r ) = k ( a ) k ( 1 ) L ( r ) = η ( r ) L ( 1 ) = C η ( r ) (12)

where C = L ( 1 ) is a constant. Consequently, when L ( r ) is a continuous function on this common set, then it is true also for η ( r ) . In this case, the solution of the function Equation (12) is known for positive r values [

η ( r ) = K r α (13)

where K and α are constants. Hence, we get from Equation (12):

L ( r ) = A r α (14)

where A is another constant. Consequently, the power-law principle shown in Equation (1), can be applied indeed in general case for the description of self- similarity.

The concept of self-similarity may also be extended from the structures to the processes [

f : t ∈ I = [ 0 , ∞ ) ↦ f ( t ) ∈ R (15)

function. Let us choose the time t expressed in the form of

t = τ t ′ (16)

with the aid of an earlier t ′ time. Our statement is that the process described by the function f is self-similar, if

f ( t ) = f ( τ t ′ ) = g ( τ ) h ( t ′ ) ∀ t , t ′ ∈ I = [ 0 , ∞ ) (17)

Since we have proved earlier that the solution of the equation of the above function has the form of

f ( t ) = a b t n , g ( τ ) = a τ n , h ( t ′ ) = b t ′ n (18)

thus, self-similarity can be written in a simpler form of:

f ( τ t ′ ) = f ( τ ) f ( t ′ ) ∀ t , t ′ ∈ I = [ 0 , ∞ ) (19)

Equation (19) means that the value of the physical quantity of the process can be established at a later point in time by the extension or shrinking of an earlier value, where the amount of extension or shrinking depends only on the ratio of the two time points. Therefore, the ratio of function values is the function of the ratio of time points. A solution of Equation (19) (Cauchy-type function-equa- tion) is [

f ( t ) = t α (20)

where α is an arbitrary number. Of course, in nature, not every process is self- similar. Those deserve attention, which can be transformed them into self-simi- lar processes. When f is not self-similar, we may apply a suitably chosen u transformation, when the composition of u ∘ f will be self-similar, namely:

u [ f ( τ t ′ ) ] = g ( τ ) h ( t ′ ) ∀ t , t ′ ∈ I = [ 0 , ∞ ) (21)

The u function is named: comparative function. It is proven [

A wide group of processes differ regarding mathematics; however, a process can be transformed into self-similar one with the help of a comparative function scarcely dissimilar in approximation capability.

We are going to show that the power type equation is the consequence of some simple and very general assumptions.

The Avrami equation [

x ( t ) = 1 − e − ( c t ) α (22)

where t is the time of the reference system, and α and c are constants. (The function is identical to the well-known Weibull function [

g ( t ) = ln 1 1 − x = ( c t ) α (23)

Moreover, it is transformed to a new time-scale by the

t ↦ υ = c t (24)

linear transformation. Following these changes, the Equation (22) Avrami-equ- ation takes its final form as:

f ( υ ) = ln 1 1 − x = ( υ ) α (25)

The function ln 1 1 − x is the comparative function for the Avrami processes,

(Avrami comparative function). Consequently, the Equation (24) linear transformation could be chosen when the Equation (25) Avrami comparative function is self-similar.

In the case when f is a continuously differentiable function, the function Equation (19) is identical to the following differential equation:

d f d υ = α f υ (26)

Substituting function Equation (25) into Equation (26), we get the differential equation:

d x d υ = α υ α − 1 ( 1 − x ) (27)

Transforming Equation (27) back to the reference time-scale, we obtain the differential equation describing the dynamic features of the process.

d x d t = c α ( c t ) α − 1 ( 1 − x ) (28)

It is proven [

A wide group of processes differs regarding mathematics. However, it can be transformed into a self-similar one with the help of a comparative function hardly dissimilar in approximation capability.

With our results, we could explain numerous important processes in biology. Let us describe first the growth of nuclei in a homogeneous thermodynamic system. For the sake of simplicity, consider the nuclei as a sphere, and the inside solid material has a constant number of particles by volume units n_{0}, as well as the concentration distribution of the particles in the neighborhood of the nuclei considered as having spherical symmetry. For the diffusion growth of these nuclei, the conservation of the particles could be formulated:

d N d t = I (29)

where N is the number of the particles in the nuclei, I is the particle current on the surface of the nuclei. If the nuclei radius is r and it is homogenous, then

N = n 0 4 π 3 r 3 (30)

moreover, the particle current

I = 4 π r r 2 j (31)

where j is the current density of particles. If j = const. then

r = j n 0 t (32)

The exponent, in this case, is trivially n = 1. Do not forget: this is a particular instance when both the concentration gradient and the thickness of the boundary of the nuclei are constant. Let us study when only the current density kept constant on the nuclei boundary. Then, based on (25):

n 0 4 π r 2 d r d t = I = 4 π r 2 D ∂ n ∂ r = all . (33)

where D is the diffusion constant. If the concentration distribution is stationary, then

n 0 4 π r 2 d r d t = I = 4 π D r [ n ∞ − n ( r ) ] (34)

where n ∞ is the particle concentration far from the nuclei. The particle concentration of the nuclei could be written [

n ( r ) = n plane ( 1 + 2 α n 0 k T r ) (35)

where n_{plane} is the concentration at a plane, α is the surface tension, and T is the temperature. If ( α / n 0 ) ≪ 1 , then

2 α n 0 k T r ≪ 1 (36)

Hence from the condition Equation (36) and Equation (35):

r ≈ t 1 2 (37)

so the exponent is n = 0.5. Contrarily, if ( α / n 0 ) ≫ 1 , then

r ≈ t 1 3 (38)

so the exponent is n = 1/3. Consequently, the volume of the nuclei depends on the following:

V grain ≈ t 3 n (39)

which is the Avrami exponent.

Let us calculate the Avrami comparative function for the same conditions. Denote the part of the old phase transforming into the new one by x. According to the Kolmogorov-Avrami assumptions [_{nucl}/dt) and the non-transformed part of the phase:

d x d t ≈ d V n u c l d t ( 1 − x ) (40)

Hence, by integration of (40), the V_{nucl}:

ln 1 1 − x ≈ V n u c l (41)

Using the k = const. factor, we get:

ln 1 1 − x = k t 3 n (42)

Consequently, the nuclei growth process by diffusion with the

ln 1 1 − x (43)

Avrami comparative function and on the

k 1 3 n t (44)

time-scale is a self-similar process and could be described by the Avrami exponent:

α = 3 n (45)

A summary of these approximations is shown in

To generalize the above, let us denote a cluster size with arbitrary shape by ξ. In this case, we can assume its volume is V n u c l ≈ ξ k . Hence the Avrami-expo- nent would be

α = β n

The actual value of β depends on the form of the cluster, e.g. in the needle

Character of the process | Avrami exponent (α) |
---|---|

Constant-current density of the diffusion | 3 |

Constant flux of the diffusion, with ( α / n 0 ) ≪ 1 | 3/2 |

Constant flux of the diffusion, with ( α / n 0 ) ≫ 1 | 1 |

form; it is β = 1; in the planar form, β = 2. The above considerations could be extended even for fractal dimensions: β = dfr, where dfr is the measured fractal dimension.

Various processes of class α = 1.0 are summarized in

The process character is the Avrami exponent α. Based on the extensive works of FW Cope [

The method could also be generalized on vectorial processes. Assuming that the processes are disjointed and the individual modes have the same comparative functions, then all the g i ( t ) comparative functions of the modes could be transformed into the form:

g i ( t ) = ( c i t ) α , ( i = 1 , 2 , ⋯ , N ) (47)

where c_{i} is the scale factor of the time-scale of the i-th mode and N is the number of modes. The example to study could be a bioprocess, chemical catabolism, exchanging materials with the actual neighborhood. There could be two reasons to change the number of particles in the system: a breakdown of a molecule or a molecular exchange at the surface of the system. We must consider not only one but numerous chemical components participating in this process. It could happen even in the simplest double-component system, not only because of the chemical potential of the breaking but also because the final component differs from its environmental values. Under these circumstances, the diffusion of both compounds starts and becomes coupled according to Onsager’s non-equilibrium [

d x 1 d t = J k 1 − λ x 1 d x 2 d t = J k 2 + λ x 1 (48)

where λ is the constant of the decomposition. On the other hand, for J k 1 and J k 2 , the Onsager equations determine particle currents:

J k 1 = − D 11 x 1 − D 12 x 2 J k 2 = − D 21 x 1 − D 22 x 2 (49)

where the matrix of D_{11}, D_{12}, D_{21} and D_{22} conduction constants is a positive defi- nite structure due to the second (entropy) law of thermodynamics. Furthermore,

Process | Descriptive function | Comparative function |
---|---|---|

Decomposition, breakdown | x = x ∞ e − c t | ln x ∞ x |

Infinite proliferation | x = x ∞ e c t | − ln x ∞ x |

Logistic process by Fermi distribution | x = 1 − x ∞ 1 + e − c t | ln 1 − x x ∞ + x − 1 |

Logistic process by Gompertz distribution | x = x 1 ∞ e − x 2 ∞ e − c t | ln x 2 ∞ ln x 1 ∞ x |

Demographic explosion | x = x ∞ π arccot ( c t ) | cot x x ∞ π |

Processes in solids [ | Avrami exponent (α) |
---|---|

Beginning of nucleation | 3 |

Nucleation with constant speed | 4 |

Nucleation with growing speed | <4 |

Staring nucleation and surface growth | 2 |

Diffusional growth from nuclei | 1.5 |

Diffusional growth from growing nuclei | 2.5 |

Diffusional growth in 1 or 2 dimensions | 1 |

Limited diffusional growth in 1 or 2 dimensions | 1.5 |

Biological processes [ | Avrami exponent (α) |
---|---|

Growth of bacteria | 4 |

K+ conductance decay in nerve axon | 1.9 |

K+ leakage from poisoned muscle | 1.7 |

Growth (weight) of rat | 2.0 |

Growth (length) of regular leg of salamander | 2.3 |

Growth (height) of sunflower plant | 1.9 |

Synthesis of chlorophyll in maize plant | 2.2 |

Muscle tension during tetanic contraction | 1.21 |

Myosin splitting of ATP | 1.24 |

Fresh green leaf IR phosphorescence decay | 1.0 |

Dried green leaf IR phosphorescence decay | 1.28 |

Melanin IR phosphorescence decay | 1.08 |

Cytochrome c IR phosphorescence decay | 1.16 |

the set of the equation is symmetrical (D_{12} = D_{21}) according to the Onsager’s theory, [

d d t [ x 1 x 2 ] + [ D 11 + λ ˙ D 12 D 21 − λ ˙ D 22 ] ⋅ [ x 1 x 2 ] = [ 0 0 ] (50)

describing the time variation of the particle currents. The modes correspond with the eigenvalues of the coefficient matrix. The time dependence of the modes:

e η t (51)

moreover, the characteristic equation to determine the eigenvalues:

det [ D 11 + λ − η D 12 D 21 − λ D 22 − η ] = 0 (52)

Hence, we get:

η 1 , 2 = − ( D 11 − D 22 + λ ) ± ( D 22 − D 11 − λ ) 2 + 4 D 2 12 − 4 D 12 λ 2 (53)

Both the modes of the vectorial process belong to the class α = 1.0, with the ln x i ∞ x i , ( i = 1 , 2 ) comparative function and with υ_{1} = η_{1}t and υ_{2} = η_{2}t time scale.

Description of the self-similar traveling waves is also possible in above format. Let us denote the variable for the traveling wave:

ξ : = v t + y (54)

which belongs to a wave traveling in the y-direction at velocity v. According to our previous results, the comparative function is:

g ( y , t ) = ( c ξ ) α (55)

The comparative functions are similar by an acceptable accuracy with appropriately chosen parameters. Let us choose the Avrami-like comparative function. By this, the traveling wave could be described by the following wave function:

x ( y , t ) = 1 − e − [ c ( v t + y ) ] α (56)

This wave is not a solution because the corresponding

d x d ξ = c α ( c ξ ) α − 1 ( 1 − x ) (57)

Avrami-equation leads to a linear equation with non-constant coefficients:

1 v ∂ x ∂ t + ∂ x ∂ y − 2 c α ( c ( v t + y ) ) α − 1 ( 1 − x ) = 0 (58)

The two-component autocatalytic reactions could be described by the following vector equation [

d x 1 d υ = x 1 ( 1 − x 2 ) d x 2 d υ = − x 2 ( 1 − x 1 ) (59)

Introducing the variable of the traveling wave by Equation (54) and regarding the process as self-similar by Equation (55), then we receive the following coupled vectorial wave-equations:

( ∂ ∂ t + v 1 ∂ ∂ y ) x 1 = x 1 ( 1 − x 2 ) ( ∂ ∂ t + v 2 ∂ ∂ y ) x 2 = − x 2 ( 1 − x 1 ) (60)

A solution of Equation (60) was made with bilinear transformation by Hirota [

x 1 = x 2 = 1 2 { 1 + tanh [ v 1 + v 2 2 ( y + ( v 1 − v 2 ) t ) ] } (61)

This equation is very similar to the various assumptions (e.g. [

In this paper, we show the possible transformation of the mathematical models of biological processes to the self-similar processes with the help of an appropriate comparative function. The mathematical expressions of the applied compa- rative functions are different, but by the suitable choice of their parameters, they approach each other well with little error. Consequently, in the practical solutions, to use one of them for the various tasks is satisfactory. The Avrami equation will describe the mathematical model of the processes with Avrami-like comparative function. We believe this is the reason for the unusual universality of the Avrami equation in biology. Consequently, the self-organizing conditions and the analogy of the self-organized processes could be a fruitful heuristic method in biological model calculations.

Szasz, O., Szigeti, G.P. and Szasz, A. (2017) On the Self- Similarity in Biological Processes. Open Journal of Biophysics, 7, 183-196. https://doi.org/10.4236/ojbiphy.2017.74014