The Curie-von Schweidler law relates to relaxation current in dielectric when a step DC voltage is applied and is given by , where and the power (exponent) i.e. n is called relaxation constant or decay constant, where     . We note that n is non-integer. This relaxation law is taken as universal law, at least for dielectric relaxations. Whereas we are used to Debye type of relaxation i.e. exponential decay law given by or where is the relaxation time constant while denotes the relaxation rate of the process with . The radioactive decay is example of ideal Debye law where the exponential decay is governed by “one-lumped” decay constant i.e. . The Curie-von Schweidler behavior has been observed in many instances, since late 19th Century, such as those shown in dielectric studies and experiments  -  .
This power law relaxation of the non-Debye type i.e. has been interpreted as a many-body problem but can also be formulated as an infinite number of independent relaxing bodies meaning infinite number of time constants or relaxation rates varying from near zero to infinity   . The observations of power law relaxation are also made in the experiments and studies with super-capacitors      . These studies also indicate the fractional calculus is used as constituent expression to describe super-capacitors. The use of empirical power law i.e. Curie-von Schweidler Law of relaxation of current to a step input of voltage to get constituent relation with fractional derivative was proposed in  . Apart from relaxation of current decay in dielectrics and super-capacitors, the power law type or non-Debye relaxation is observed in visco-elastic experiments strain relaxation in     .
In this paper, we are giving the derivation of the distribution of relaxation rates ( ) particularly for Curie-von Schweidler law and we observe the distribution nature as Zipf’s distribution  -  . We try to reason out as to why this distribution of relaxation rates takes Zipfian nature. We also show that Curie-von Schweidler law has time varying rate of relaxation. This paper will not deal with the mathematics of Zipfian distribution (or power law distribution) like probability density function, cumulative probability density function, and the conditions of finding finite mean, variance or standard deviation for power law distribution. This paper describes finding the distribution function of relaxation rates (or histogram) by formulating Laplace integral, and show that the distribution thus obtained is a Zipf’s power law.
We extend this mathematical approach to get the distribution function for time constants ( ). We observe that time constants are also distributed as Zipf’s power law; but this observation points to a contrary physical interpretation derived from this obtained power law distribution for the relaxation rate ( ) distributions. Thus we can conclude that this method developed by Laplace integral approach is restricted to get only distribution of relaxation rates i.e. and not to get the distribution of time constants i.e. . Though we are discussing especially Curie-von Schweidler law, yet we will tabulate relaxation rate distributions obtained for some other relaxation functions which are obtained via this Laplace integral method.
We shall demonstrate the formation of fractional derivative in the expression relating current and voltage considering the relaxation rates as Zifian distribution; and thus forming a scale dependent power law for relaxation rates as the scale varies from zero to infinity. Though by experiments one cannot make histogram directly for the rates of relaxation for any non-Debye processes, yet this mathematical procedure that we develop helps in extracting this information from the observations relaxation function. This is new treatment, and much more research is required, across various dynamic processes.
2. Obtaining Fractional Derivative Directly from Curie-Von Schweidler Law for Capacitor
Practically on applying a step input voltage Volts at to a capacitor which is initially uncharged; we get a power-law decay of current given by empirical Curie-von Schweidler as  . That we write in following way as indicated by experimental studies  -  :
The parameter is constant. This is from observation and the evaluation of order of power-law function is       . Let the capacitor be excited by a step input of VBB Volts, i.e. written as , where is unit step function at time . The Laplace transform of step input is . Then taking Laplace transform of above power-law decay current (1), we obtain
 . Then using the formula
for generalization of factorial i.e.    , we get the following expressions
We get Transfer function  of capacitor as following expression
This Expression (3) i.e. is admittance expression in complex frequency s-domain of a capacitor. From here we write impedance expression for capacitor as following
From the obtained Expressions (3) (4) i.e. and by taking inverse Laplace transform by using the identity i.e. fractional derivative operation   , we get the constituent relation for capacity as following
This fractional derivative Expression (5) gives a new capacitor theory  and we utilize this above Formula (5) to find characteristics of super-capacitors, like the variation of n with the current excitation, and the efficiency of the energy discharged to the energy stored      . Classically the expression of capacitor is i.e. with integer “one-whole” order classical derivative. Therefore Curie-von Schweidler law gives a different approach based on fractional calculus   . In experimental observations we find that capacitor has fractional order impedance  -      ,  . The impedance (is obtained by writing Laplace variable s in (4) as , i.e. considering steady state analysis   . This fractional impedance observed in  -       , has implication in dissipation  theory of di-electrics that we will not cover here. We state that while classical capacitor unit is in Farad, the the “fractional capacity” is in units of  .
This section gives us the understanding that this Curie-von Schweidle law i.e. the empirical law gives a relation of voltage and current of capacitor by using fractional derivative. In this paper we will show how we get the same relation (i.e. via use of fractional derivative) by considering Zipfian distribution of relaxation rates ( ) that we get for Curie-von Schweidler relaxation law.
3. About the Zipf’s Power Law Distribution and Probable Hypothesis for Its Mechanism
The Zipf’s law is widely referred in linguistic studies, economics studies, population studies  -  . We use this for a dielectric relaxation law (i.e. Curie-von Schweidler law), which is observed as , since late 19th century. We derived histogram of relaxation rates for relaxation function and show that it follows Zipf’s power law. We try to give possible reasons as to why Zipfian distribution is observed for the distribution of relaxation rates. The histogram function of Zipf’s law is a power law type. In this section we assume that relaxation rates l’s follow a Zipfian histogram say . This we will derive in subsequent section. The l’s are relaxation rates of infinite number of relaxing bodies, simultaneously relaxing as per Debye law i.e. .
Zipf’s Law is an empirical law formulated using mathematical statistics that refers to the fact that many types of data studied in the physical and social sciences can be approximated with a Zipfian distribution. This distribution is one of a family of related discrete power law probability distributions  -  . This power law distribution help to describe phenomena where large events are rare, but small ones are quite common. For example, there are few large earthquakes but many small ones. There are a few mega-cities, but many small towns. There are few words, such as “and” and “the” that occur very frequently, but many which occur rarely.
The emergence of a complex language is one of the fundamental events of human evolution, and several remarkable features suggest the presence of fundamental principles of organization. These principles seem to be common to all languages. The best known is the so-called Zipf’s law, which states that the frequency of a word decays as a (universal) power law of its rank. The possible origins of this law have been controversial, and its meaningfulness is still an open question. One of the early hypotheses of Zipf of a principle of least effort for explaining the law is shown to be sound   . But still the exact mechanism how the Zipf’s distribution manifests is debated.
Many of the things that we measure have a typical size or “scale”. We ask ourselves why the relaxation rates cannot be arranged as simple “normal distribution”. Like while we plot the height of person in X-axis and the percentage of occurrence of that particular height in Y-axis, we get a “normal distribution” peaked around mean height with a spread both ways, that is a histogram. We find that ratio of maximum height and minimum height of a person is finite (or relatively low value). For example as per Guinness book of records tallest person was having height 272 cm and shortest person was having the height of 57 cm, making this ratio 4.8. This ratio is relatively low value. We see the most adults are about 170 cm tall-there is some variations around this figure notably depending on sex, but we never measure persons having height of 1 cm or 1000 cm.
But not all things we measure are peaked around a typical value. Some may over a very large dynamic range, sometimes many orders of magnitude. For example the ratio of population of largest town to population of smallest town is about 250,000. The histogram if plotted for X-axis with population of cities and Y-axis with percentage of cities having that population; the distribution will not show the “normal-distribution”. The histogram of cities & population is highly “right-skewed”, meaning that while the bulk of distribution occurs for fairly small sizes―i.e. most cities have small population-there is small number of cities with population much higher than a said typical value, producing the long tail to the right of histogram. This “right skewed” form is qualitatively quite different from histogram of person’s height. That is because we know that there is large dynamic range from smallest to largest city sizes, we can immediately infer that there can only be a small number of very large cities. The histogram of this sort is like a function i.e. . The distribution of this nature is called Zip’s power law distribution.
The same we observe when relaxation rates call them having a large (ideally infinite) spreads follow Zipfian distribution, we call that and will show that the histogram follows the function i.e. . Thus one reason that this non-Debye relaxation (explained in subsequent sections) of Curie-von Schweidler Law ( ) in dielectric is having infinite spread of relaxation rates of l’s-thus forming a Zipfian power law.
4. Zipfian Power Law Distribution as a Result of Connected Exponential Processes―A Postulate
Having discussed the formation of a histogram as power law type, when there is very large dynamic spreads amongst the relaxation rates of a complex relaxing process we move to a probable postulate of explaining this process via exponentially distributed processes. A much more common distribution than power law is the exponential distribution. In this complex relaxation mechanism i.e. that we are discussing we consider infinite number of bodies relaxing simultaneously, in different time scales ( ). We consider that a complex relaxation mechanism and a quantity say survival time of a relaxing body, has exponential distribution of probability . This means that a probability for a body having very large survival time (age) is very low; and vice-versa. Then indicates the fraction of survival numbers of bodies between survival time and . Now suppose that the real quantity that we are interested is not but other quantity , say the relaxation rate of discharge which is exponentially related to ; thus . That implies the surviving bodies with very large time of survival (age) have a very low rate of relaxation. This also states that . Then if probability distribution of is ; then we have (effect of conservation of probability  ). The negative sign indicates opposite movement, as is increased from to , then is decreased from to . This means that number of discharging units having relaxation rates between and is equal to number of surviving bodies having survival time between and . Thus we write following steps:
The above discussion in steps (6) gives a power law distribution for relaxation rates l’s where there is combination of exponential processes. Thus we expect that in our complex relaxation process governed by Curie-von Scweidler Law which is having infinite number of simultaneously discharging bodies will have a power law distribution for relaxation rates as a histogram . This we will derive subsequently. We proceed with this explanation and hypothesis. This could be one explanation in physical sense, in line with exponential distribution in the Boltzmann distribution of energies in statistical mechanics.
5. Complex Relaxation of Non-Debye Type Composing with Several Exponentially Relaxing Functions of Debye Type
We call the Curie-von Shweidler relaxation law as complex process, of non-Debye type. Where a Debye type relaxation is a decay function given by function of exponential type as . We mention that Curie-von Shweidler relaxation law is not exponential relaxation process of Debye type.
In this section we formulate the method to extract the histogram of the relaxation rates call it , for a complex non-Debye relaxation process , which we assume to be composed of several Debye type exponential relaxation functions , with varying from zero to infinity. The complex decay may be expressed as following with several rate constants with weights , where is having units in sec−1 i.e. “per second”, and is equal to inverse of time constant i.e. . We write following composite relaxation expression as sum of several “discrete” relaxations of Debye type i.e.
The coefficients ’s in (7) can be positive or negative that we will elucidate in later section. In the continuum limit, we may write the above discrete expression as following integral equation
Where the function i.e. is the distribution-function of the rate of the relaxation ( ) of the process, or we may call it as histogram of relaxation rates. The function can be positive or negative that we will elucidate in later section.
While for the case with discrete set of relaxation rates i.e. the rate distribution function would be having discrete delta functions ( ) at points ; which we write like following expression
From above formulation (9) we infer that if we have only one single Debye relaxation i.e. having only one rate constant say i.e. then . This is verified in the following expression
In above derivation (10) we used the property of delta function    i.e. .
6. Extraction of Rate Distribution Function by Formulating Laplace Integral
In this section we formulate Laplace integral of the complex decay given by Curie-von Shweidler relaxation law and then getting by inverse Laplace transform of time domain response i.e. we get relaxation rate distribution function i.e. . Conventionally we are used to get inverse Laplace transform of a frequency domain function to time domain function; we note here we will be inverting a time domain function i.e. .
The Laplace transform of a function in time domain is defined as following integral transform relation    , i.e. called Laplace integral
This Expression (11) is standard integral transform of a function from a time domain (t) to a complex frequency domain i.e. ; where real part is significant in the transient response and the imaginary part of the frequency corresponds to “steady-state” response; in classical “Control Science”  . Here is “inverse Laplace transform” of , and we write and .
We have in earlier Section 8 derived . Compare this with defined Laplace integral expression as follows
Both expressions in (12) above are Laplace transform expressions, (or Laplace integrals). The first expression is transforming the function from domain to “complex” t time domain; while the second one is transforming from t domain to “complex” s frequency domain. Thus both expressions are Laplace integral expressions with change of variable and symbol. Therefore we can say is inverse Laplace Transform of in the first expression, i.e. ; as is inverse Laplace of in the second expression, i.e. .
Therefore in order to get the rate distribution-function from the decay curve (or relaxation-function ), we need to perform inverse Laplace Transform of the time function . The definition of inverse Laplace Transform is described as following integral expressions
In the above Expression (13) is real number larger than , where being such that has some form of singularity on the real line but is analytic in the complex plane to the right of that line, i.e. for ,    . Thus in this formulation we treat time variable as complex quantity say in the Expression (13) of inverse Laplace Transform i.e.
. Though we cannot explain presently physical
meaning of concept an “imaginary time” in the expression of complex time quantity , yet mathematically there is no restriction in assuming time to be complex number. We thus proceed in mathematical sense to invert a function in complex time variable, by techniques of Laplace inversion.
7. Extracting the Rate Distribution Function by Performing Inverse Laplace Transform on Relaxation Function of Time Variable
Table 1 lists several types’ relaxation functions and its inverse Laplace describing the rate distribution function; mostly got from standard Laplace transform tables  . The integral representations of , shown in Table 1 i.e. for entries 12 to 16 is got via Berberan-Santos method   . The entry 12 is for stretched exponential decay function and entry 13 is Becquerel’s compressed hyperbolic radioactive decay function; the entry 15 and 16 is for Mittag-Leffler function and the entry 14 is general power law relaxation. These integral representations of are difficult to solve but are easy to plot via use of numerical integration techniques.
We have observed in the previous section that a Debye relaxation of has rate distribution as i.e. it is given by a delta function at point (10). This we verify with known Laplace relation i.e.  , where . Also we have ; thus we can write . From here we can write with change of variable for the inverse Laplace of this time domain function in domain we get as , i.e. the rate distribution function.
If there is no decay then say ; the rate distribution function is delta function at origin i.e. .
If the relaxation function is of say ; then we have
; giving . From this observation we say that for our earlier Expressions (7) and (9) i.e. ,
Table 1. Several relaxation functions and corresponding rate distribution functions.
the coefficients ’s can have negative values as well for some type of relaxation function.
For example, if is a relaxation function that initially grows to a maximum value and then starts falling as time increases, it has rate distribution function as , a oscillatory one. Thus in this case the distribution function i.e. can take positive as well as negative values (8).
One interesting observation is for a relaxation function the relaxation function is -a “uniform distribution”, for . All these are listed in Table 1.
The inverse Laplace transformation is usually carried out by contour integration. But the very modern technique of Berberan-Santos   method is the analytical Laplace inversion without the usual contour integration. We describe this now briefly.
Our aim is to evaluate Laplace inverse which is given as Laplace inversion (13) integral expression i.e.
Here we describe Berberan-Santos method formulas for evaluation of the Laplace inversion without going for contour integration. First is change of variable i.e. from “real time variable” to “complex time variable” as ; with . Here the real part i.e. is constant as a vertical line calls it a constant. The formulas are following  
Consider a very simple case of decay function and convert to complex time by putting as   . We know from standard Laplace pair that is . Thus, for we should get via inverse Laplace the rate distribution functions as . The application of the Berberan-Santros formula   with yields the following steps
Here we say that has integral representation as
Particularly for , we have . The condition enables us to choose we get following integral representation for   which is also rate distribution function is following
8. Derivation of Rate Distribution Function for Curie-Von Schweidler Law
For the Curie-von Schweidler relaxation of type function i.e. then rate distribution function is . With the known Laplace pair
i.e. , we can write the following steps
Therefore above discussion suggests that for a power law type relaxation, i.e. Curie-von Schweidler law i.e. , the relaxation rates l’s are having a power law distribution of type i.e. , , , . This is Zipf’s power law with .
For dielectric relaxation as observed that in Curie-von Schweidler relaxation , the rate relaxation distribution function has exponent in power in the range i.e. . Considering graph of as histogram, we infer that for Curie-von Schweidler relaxation function i.e. there are very large number of relaxations with small i.e. large number of slower decay takes place, compared to fewer faster decay rates-and the histogram is highly right skewed with long tail.
From the above discussion (18) and using our Laplace integral (8) i.e. we write for Curie-von Schweidler relaxation function the following
The above Expression (19) is integral representation of the shows weighted averaging of infinite Debye relaxations i.e. with weight applied for all from zero to infinity.
Using Berberan-Santos method   we get the Laplace inversion of . In reality of decay functions, we can take ; in complex time variable i.e. as the decay function is not expected to have singularity at time . Choosing in we have , that is
where we used
. The real part of the complex function is . Now using the Berberan-Snatos formula   , we
get following steps
From Laplace transform tables  (and Table 1) we have
i.e. from inverse Laplace transformed of , therefore we write following representation
From (21) we write the following
Now by changing variable to , to , and rearrange above
expression (22) to get .
Considering now as we write another integral representation of as follows
This means that if we chose basic relaxation function as , then Curie-von Schweidler relaxation is weighted sum of all ’s with weights , as is varied from zero to infinity.
9. Zipf’s Distribution for Relaxation Time Constants for Curie-Von Schweidler Law―A Contradiction
Now converting to , we assume the Distribution of time-constants call it , is the Zipf’s power law distribution. However direct taking of reciprocal of obtained inversion of i.e. the rate distribution function got via Laplace inversion of is not possible. This we demonstrate in this section.
As we have formulated Laplace integral (8) i.e. , just
by replacing , we will get
that is . This is not Laplace integral. Now we do the following steps, for and obtained
We write the two representations of as following integrals
Thus we have , as we have . Now we verify the above obtained result in the subsequent discussion.
By the logic that we had constructed which is Laplace integral (8); we will similarly get the integral which is not a direct Laplace Transform formula. Following steps will convert this expression into the Laplace Transform formula, and from there we will extract :
Proceeding further we obtain following result
Now we take different approach to verify the above (26) (27) obtained expression for . Let us have set of relaxation functions with various time constants ranging from 0 to infinity that is , comprising of infinite number of functions, in continuum in . The relaxation function varies from very-very quick decay (when ) to very-very slow decay curve (when ). We construct a weighted decay function as . This shows that we are multiplying by weight the decay function . We are assuming Zipf’s type distribution of , in form of , meaning the lowest time constant i.e. fastest decay occurs more frequent than slow decay i.e. large time constant. The time constant parameter let vary from 0 to infinity and construct the following integral , i.e.
The integral of (28) i.e. gives notion of weighted average of
infinite relaxation functions. We do the substitution i.e. , i.e. and
in the above integral to get following steps
By using the definition of the Gamma function   in integral form i.e. , we write the above integral as following
Putting in above we get integral representation of the power law and we represent this by time constant distribution function in following expressions
Earlier in (19) we have obtained ; where we called
rate distribution function as . Now from above (31) weighted average logic we get . We note that these two are not
reciprocal of each other.
Therefore we can conclude that Curie-von Schweidler law ( ) relates to weighted averaging of several classical Debye relaxations (of type ) over several time constants from zero to infinity, that is having Zipf’s power-law with time constant distribution as . What does it say for , that is also a right-skewed distribution, where the lower time constants (faster decay) appear more than larger time constant (slower decay). This is contradiction to what we inferred for .
10. Demonstration of Contradiction of Obtained Zipf’s Power Law Distribution for Relaxation Rates and Time Constant Distributions
This contradiction we demonstrate now. For we got ; we expect for ; let us see what happens in following steps
Though mathematically we can get integral representation for any relaxation function as but physically it will be contradictory to Laplace integral i.e. . Hence we will deal with the relaxation rate distribution function that we extracted as from via our devised method of Laplace inversion.
11. Experimental Validation of Range of Relaxation Exponent in Curie-Von Schweidler Law
The Curie-von Schweidler empirical law of power law relaxation, i.e. states that . This is validated via experiments on dielectric relaxations. A 100V step input applied to a completely discharged capacitor of 0.47 μF having metalized paper dielectric, and the current decay is recorded with time. The graphs of log-log plot i.e. vs. show a straight line of average slope −0.86  -  . This experiment indicates a Curie-von Schweidler law, with , having . This gives . The exponent n is in the range of in several di-electric relaxation experiments  -  . The experiments with super-capacitors   , show range as . A very low value of exponent n is found in relaxation of Laponite studies averagely  . Thus in case of Laponite studies we have relaxation rate distribution function as . In this Laponite study  though the exponent n was obtained on “self-discharge” curves with various charging time history-showing memory effect, the expression obtained for self-discharge decay of voltage assumes fractional capacity-that in turn assumes Curie-von Schweidler law as current relaxation function.
12. Time Variant Relaxation Rate for Non-Debye Relaxation & Curie-Von Schweidler Law
Any decay function is written as a general formulation in following way
where is the time ( ) dependent rate coefficient. When the relaxation is pure exponential, one has as constant say described as , expressed in following steps
Thus we get a Debye relaxation for a system having constant rate of relaxation. To extract that is time dependent rate coefficient we have to follow the following steps by taking logarithm of (33) and then differentiating both sides
We use the above rule for Mittag-Leffler relaxation function i.e. ,
this is defined as    . For
the is extracted as in following steps
For Curie-von Schweilder relaxation law , we have time dependent rate relaxation rate as
Therefore we have two observations that Curie-von Schweidler relaxation has time constant distributed as Zipfian power law , while the relaxation rate constant is variable in time as a function . Thus implying the relaxation starts with very-very fast relaxation at a very-very high rate ( ) and as the time goes the rate constant decreases indicating slow rate of current decay. This is a case of “equivalent” single body relaxation where the rate is varying with time; whereas the multi-body relaxation gives simultaneous relaxations with rates distributed as Zipf’s power law.
13. Scale Dependence Relaxation Rates Give Capacitors Charging Current as per Curie-Von Schweidler Law
Let a uncharged capacitor be connected to a voltage source Volts, at time ; obviously this capacitor will get charged to the battery voltage. Let this capacitor is uncharged at , thus there is no charge held by it, therefore the voltage across the capacitor is zero at , and the circuit current is . The voltage balance equation assuming be the total resistance of the circuit (including internal resistance of Capacitor) at is the following
where is the charging current flowing into the capacitor. The above integral Equation (38) is differentiated and is put as following, for
The RHS of above Expression (39) first order system indicates “forcing
function” which is . The forcing function is in this case.
If we take i.e. a constant, then considering and , i.e. “unit-step function”, we have for RHS of the above Equation (39) as following
Thus substituting the above (40) into (39) we get following equation for for the following
We see that forcing function of above first order Equation (41) is by delta function . The solution to the above equation gives Debye relaxation function i.e.
This solution is the “impulse response” of the circuit equation. The relaxation current of the above system (41) follows Debye’s relaxation, with one relaxation rate (also termed as Debye law). The rate distribution function is ; that we discussed in previous sections. With this as Green’s function call it i.e. solution of differential equation with “unit” impulse excitation ( ) or say Homogeneous solution, i.e.
Now we find if the input is step function at time , call it , where for and for ; then we get relaxation function for current as convolution integral, i.e. depicted as in following steps
We saw in earlier sections the relaxation rates ( ) distribution, for a Curie-von Schweidler relaxation law, i.e. is ; for relaxations in dielectrics. This is histogram of rates. It says that the relaxation of current is with several relaxation rates, which are distributed as discussed in Zipf’s law fashion with right-skewed-histogram. Thus if we represent the equivalent relaxation rate say with as scale of relaxation where the scale varies from zero to infinity; we will not be incorrect in assuming this. That is as we slide from a low scale to high scale the equivalent relaxation rate will be different at different scales of relaxation. If the index parameter i.e. then we have single rate constant system given by always at all scales of relaxation i.e. , and with solution as , i.e. Debye relaxation function.
We thus modify the capacitor discharge current equation, with i.e. with one relaxation rate at any scale of relaxation (λ), i.e. to following i.e. variable at any scale of relaxation rate(λ)
The initial condition is given as for . The above equation is having a free “scale” parameter varying from zero to infinity. The solution of the above is . We call this as “impulse response function” at a particular scale , i.e. we call it
The above Expression (46) actually is valid for all scale varying from zero to infinity. Thus on integrating this “impulse response function” on the free variable ( ) from 0 to , we get the function of time and that is called “impulse response” or the Green’s function as depicted in following derivation
To get above Expression (47), we substitute in ,
that makes following changes:
Then by using definition of Gamma function i.e. , and its property the following steps are followed to get the desired
By changing q to n we get integral representation of as following
For q = 1 case we have scale invariance λ thus , where at all scales. For this case “impulse response” or Green’s function is same as “impulse response function” i.e. .
We find that for a system where the equivalent relaxation rate is ; similar to a distribution function that we obtained as gives relaxation current as . We write the two currents expressions obtained as following for
Therefore we infer that the Curie-von Schweidler relaxation current for dielectric excited by a step voltage that follows the relation has distribution function a power law or Zipfian distribution, with scale dependent relaxation rate described as .
Here we mention that in this study, we have derived mathematically several integral representations of Curie-von Schweidler relaxation function , (19), (23), (24), (50). Those we list as
and . These are mathemati-
cally equivalent representations of the function , yet using for physical explanations one needs care. However, these are formulations for few definite integrals giving the same result as .
14. Appearance of Fractional Derivative―In the System Having Zipfian Power Law Distribution in Relaxation Rates, where the Equivalent Relaxation Rate Is Scale Dependent
The delta-function for excitation as shown in above section gives homogeneous system with solution as i.e. described as following
Now let the system described above be excited by a signal proportional to , a derivative of voltage excitation function ; so we write this as following
Note that if , that is unit-step-function at time then , we recover the above homogeneous differential Equation (52). Then the response to this new excitation function is convolution of
Green’s function obtained i.e. above (52), with the forcing
function i.e. now . We write the following steps to get for a forcing function
Multiplying and dividing the above Expression (54) with and using the definition of fractional integral   that is
we get the following derivation
In above derivation (56) we have used the relation and the composition rule i.e.   . This derivation implies the appearance of fractional derivative for cases where several relaxation rates (ideally infinite of them) define a relaxation process; which are having a scale dependence behavior, i.e. with histogram distributed as Zipf’s power law i.e. , and the relaxation is by Curie-von Schweidler law i.e. . Thus we have current through a system (having a complex relaxation process with several rate distributed as power law excited by a voltage as fractional derivative of it, i.e. .
Let this system with scale dependent relaxation rates with for i.e.
be excited by a source which is a delta function say ; at . This means ; where is unit step function at . With this excitation the relaxation current would be fractional integral of the input excitation that is from as depicted in above derivation (56) i.e. . We have fractional integration of delta
function   as ; and using this formula we get
This was what was derived in (49) above as impulse response (where ) i.e. .
If the excitation source of (57) is a step function as at ; meaning where the unit step function is ; then the relaxation current is fractional integration of order ; that is . Using the formula for
fractional integration of a constant i.e.   we have;
the relaxation current as following
The empirical law that is Curie-von Schweidler law, which is a type of non-Debye relaxation, (that is also considered to be universal law of dielectric relaxation of current, since late 19th century), states when dielectric is stressed with a constant voltage, gives relaxation current as . We derived the rate distribution-function (the histogram function) for current relaxation as Zipf’s power law distribution. The histogram function we found out to be of a function of type . We infer the Curie-von Schweidler relaxation is simultaneous multi-body relaxations which has a distribution for relaxation rates i.e. right-skewed. That is with large number of relaxations with lower value of rate (slow rates) followed with long tail of small number of relaxations with faster relaxation rates, relaxing simultaneously. We noted that the possibility of having Zipfian distribution arises due to very-very large ratio of maximum to minimum spread in the relaxation rates l’s, and possibility of connected exponential distribution of many body simultaneous relaxations. The method we obtained for getting rate distributions of relaxation rates via formation of Laplace integral. When this method is extended for finding distribution of time constants though mathematically correct yet gave contrary physical interpretation. Thus we carried out the entire discussion with rate distribution functions i.e. and not the time constant distribution function i.e. . We also showed that Curie-von Schweidler law gives constituent of current and voltage of capacitor via use of fractional derivative, i.e. , unlike classical capacitor relation i.e. . We verified by using obtained by Zipf’s distribution as power law for Curie-von Schwidler current relaxation law, assuming the scale dependence equivalent relaxation rate in the classical charging equation of capacitor with scale of relaxation varying from zero to infinity, i.e. . We also related the Curie-von Schweidler relaxation law which gives a time varying rate i.e. , indicating that the relaxation starts with very-very high rate, and becomes slower and slower with elapse of time. The paper gives a possible foundation for further studies in obtaining the rate relaxation distribution functions for other non-Debye type relaxation functions, and new type of explanation regarding reasons of Zipfian distributions.
I sincerely acknowledge support received from Prof Sujata Tarafdar, Prof Tapati Dutta, Research Scholars Smoasri Hazra, Moutushi Dutta Choudhury, Tania Basu, Simantini Mazumdar (Dept. of Physics, CMPRC Jadavpur University) students Adreja Mondol Rivu Gupta of St Xaviers Univ. Kolkata to have put to use concept of fractional calculus in relaxation experiments in visco-elastic studies, cooling law dynamics, impedance spectroscopy, studies in crack formation under electric stress. I also acknowledge Dr N C Pramanik Scientist CMET Thrissur, Prof Subhojit Ghosh and Research Scholar Manoranjan Kumar of Dept. of EE NIT Raipur for using fractional calculus in super-capacitor characterization experiments.