Abstract
An analogy relaxation model, able to describe the dynamics of Covid-19 pandemics, in terms of daily presented new cases, is obtained. The main objective here is the prediction of pandemic dynamics, as the forecast of herd immunity time, and provide quantitative criteria for some predictions. We propose an analogy based on the concept of fractal derivative and justify the appearance of fractal time. A Burr XII-shaped solution of a fractal-like equation was obtained. The data fitting shows that the obtained function is useful to describe the behavior of the pandemic in the studied countries. The temporal evolution of contagion over different countries and worldwide are presented and discussed. We also show that a deformed derivative model leads to the same analytical results.
PACS numbers 05.90.+m - Other topics in statistical physics, thermodynamics, and nonlinear dynamical systems, 11.10.Lm - Nonlinear or nonlocal theories and models, 02.90.+p - Other topics in mathematical methods in physics.
1. INTRODUCTION
The worldwide pandemic provoked by the SARS-CoV-2 coronavirus outbreak have attracted the attention of the scientific community due to, among other features, its fast spread. Its strong contamination capacity has created a fast growing population of people enduring COVID-19, its related disease, and a non small peak of mortality.
The temporal evolution of contagion over different countries and worldwide brings up a common dynamic characteristic, in particular, its fast rise to reach a maximum followed by a slow decrease (incidentally, very similar to other epidemic processes) suggesting some kind of relaxation process, which we try to deal with. In this case, it is in the contagion, or the number of infected cases, where the “relaxation process” occurs.
We are interested in studying the daily presented new cases, not the current cases by day. This must be noted to avoid confusion in the interpretation, i.e. we study not the cumulative number of infected patients reported in databases, but its derivative.
This relaxation process in this case is, for us, an scenario that, by analogy, will serve to model the dynamics of the pandemics. This is not an ordinary relaxation process. Due to the concurrence of many factors that make very complex its description but, just for that reason, its description must turn out to non classical description. So, we will consider that the dynamics of this pandemic is described by a “fractal” or internal time [1]. The network formed by the people in its daily activity forms a complex field of links very difficult, if not impossible, to describe. However, we can take a simplified model where all the nodes belong to a small world network, but the time of transmission from one node to other differs for each link. So, in order to study this process let us assume that spread occurs in “fractal time” or internal time [1, 2]. This is not a new tool in physics. In refs. [3-5] this concept has been successfully introduced and here, we keep in mind the possibility of a fractal-like kinetics [6], but generalizing as a nonlinear fractal-like kinetics.
Here we will follow to what we refer as a “relaxation-like” approach, to model the dynamics of the pandemic and that justify the fractal time. By analogy with relaxation, an anomalous relaxation, we build up a simple nonlinear equation with fractal-time. We also regain the analytical results using a deformed derivative approach, called conformable derivative (CD) [7]. In Ref. [8] one of the authors have shown intimate relation of this derivative with complex systems and nonadditive statistical mechanics. This was done without go into details of any kind of specific entropy definition.
Our article is outlined as follows. In Section 2, we present the relaxation model, to develop the relevant expressions to adjust data of COVID-19. In Section 3, we regain the relevant expressions by CD approach. In Section 4, we show the results and figures referring to the data fitting along with discussions. In section 5, we finally cast our general conclusions and possible paths for further investigations.
2. ANALOGY RELAXATION MODEL
Let us denote by F(t) the number of contagions up to time t.
A nonlinear relaxation model can be proposed here, based on a generalization of Brouers-Sotolongo fractal kinetic model (BSf) [3, 4, 9], represented here by the equation:
where τ is the relaxation time, q here is a real parameter and ς = tα is the fractal-time or internal time. Also, note that we do not impose any restrictive values to the parameters. There is no need to demand that the solution always converge. The equation to obtain Burr XII has to impose restrictions but this is not the case. In Burr XII the function was used as the probability distribution. But here the function describes a dynamic, which can be explosive, as will be shown for the curves of Brazil and Mexico. Therefore, if we consider infinite population, a peak will never be reached unless the circumstances change (treatments, vaccines, isolation, etc.). Our model does not impose finiteness of the solution.
Equation (1) has as a well known solution a function with the shape of Burr XII [10], with:
And so
The density (in a similar form of a probability density function (pdf), but here it is not a pdf) is, then,:
(for t0 = 0) which can be expressed as:
where the parameter are 
, 
, 
, a = α.
This is very similar, though not equal, to the function proposed by Tsallis [11] in an ad hoc way. Here, however, a physical representation by the method of analogy is proposed to describe the evolution of the pandemics. The possibility for a decay of the pandemic in a given region in this model requires the fulfillment of the condition
what expresses the property that
what means that the function has a local maximum. If this condition is not accomplished, the pandemic does not have a peak and, therefore, the number of cases increases forever in this model.
In this case, there is, apart from the change of propagation and development conditions, the possibility for a given country that does not satisfies condition (3), to reach “herd immunity”, i.e., when the number of contagions has reached about 60% of population, in which case we may calculate the time to reach such state using (3):
The parameter τ is considered here only as a fitting parameter. We will work with what T1000
3. CONFORMABLE DERIVATIVE APPROACH
Here we show that the analytical results of Section II can also be obtained by an alternative approach, that is, by conformable derivative (CD) approach.
The CD is defined as [7](9)
Note that the deformation is placed in the independent variable.
For differentiable functions, the CD can be written as
This can be obtained by a simple change of variable [7].
An important point to be noticed here is that the deformations affect different functional spaces, depending on the problem under consideration. For the conformable derivative [8, 12-15], the deformations are put in the independent variable, which can be a space coordinate, in the case of, e.g, mass position dependent problems, or even time or space-time variables, for temporal dependent parameter or relativistie problems.
Since we are dealing with a complex system, a search for a mathematical approach that could take into account some fraetality or hidden variables seems to be adequate. This thinking is also based in the fact that we do not have full information about the system under study. Here, we will show that a nonlinear relaxation equation for SARS-CoV-2 can be written with CD and the solution of this equation will give a powerful analytic function to make previsions, particularly for the pandemic’s behavior. Deformed derivatives, in the context of generalized statistical mechanics are present and connected [8]. There, the authors have also shown that the q-deformed derivative has also a dual derivative and a q-exponential related function [16].
Here, in the case under study, the deformation is considered for the solutions-space or dependent variable, that is, the number F(t) of contagions up to time t, since F has has even a fractal shape, similar to stock-market pattern. One should also consider that justification for the use of deformed derivatives finds its physical basis on the mapping into the fractal continuum [8, 17-19]. That is, one considers a mapping from a fractal coarse-grained (fractal porous) space, which is essentially discontinuous in the embedding Euclidean space, to a continuous one [12].
A nonlinear relaxation model can be proposed here, again based on a generalization of Brouers-Sotolongo fractal kinetic model (BSf) [3, 4, 9], but here represented by a nonlinear equation written in terms of CD:
where τ is the relaxation time and q and α here are real parameters. We do not impose any limit for the parameters, as explained in Section II.
Expliciting the time t, we have
Equation (11) has as a well known solution a function with the shape of Burr XII [10], with:
And so
The density (in a similar form of a pdf, but here it is not a pdf) is, then,:
(for t0 = 0) which can be expressed as:
where the parameter are 
, 
, 
, a = α.
All the results above, on this Section, are identical to those obtained in the Section II, which reinforces our results.
4. DATA FITTING
With eq. (5) let us fit the data of the epidemic worldwide. The data was extracted from Johns Hopkins University [20] and the website [21] to process the data for several countries.
We covered the infected cases taken at Jan 22 as day 1, up to May 19. The behavior of new infected cases by day is shown in figure 1. The fitting was made with gnuplot 5.2. As it seems, the pandemic shows some sort of “plateau”, so the present measures of prevention are not able to eliminate the infection propagation in a short term, but it can be seen that condition (6) is weakly fulfilled. On the other hand, eq. (3) can be used to predict the peak occurrence in different countries.
Worldwide infections from Jan, 22 to May 19 and fitting with eq. (5). The behavior fits well with A=1.35E-17, B=2.72E-24, a=13.1978, b=0.918342, P=7786246434. Condition (6) is satisfied.
In the particular case of Mexico the fitting is shown in figure 2. In this case condition (6) is not fulfilled. In terms of our model this means that the peak is not predictable within the present dynamics. Something similar occurs with Brazil, as shown in figure 3. The data for Brazil neither fulfill the condition (6). In this case there is neither the prevision of a peak and we can say that the data for Mexico and Brazil reveals a dynamics where the peak seems to be quite far if it exists.
But there are some illustrative cases where the peak is reached. Progression of the outbreak in Cuba and Iceland are shown in figure 4 and 5 respectively. Condition (6) is satisfied for both countries and we can see that the curve of infection rate descends at a good speed after past the peak.
Now let us take a look at United States data, shown in figure 6. The USA outbreak is characterized by a very fast growth until the peak and, then, very slow decay of the infection rate is evident. As discussed above, the outbreak will be controlled for almost infinite time in this dynamics.
There is also some intermediate cases as Spain and Italy, shown in figures 7 and 8. In this case the data exhibits the same behavior as in USA, a fast initial growth and a very slow decay after the peak. However, the outbreak is controlled in a finite amount of time.
Daily infections in Mexico and fitting with eq. (5) for A=0.0531581, B=0.0104261, a=3.65293, b=0.0555406, P=128932753, Thi = 675 days. Condition (6) is not satisfied.
Evolution of daily cases in Brazil and fitting with eq. (5) for A=0.00534297, B=0.0104326, a=4.63345, b=0.0651455, P=212559417, Thi = 199 days. Condition (6) is not satisfied.
Daily infections in Cuba. The peak was reached in an interval from the day 20 to day 40. The theoretical curve fits with data though with a poor correlation due to the dispersion of data. A=0.00181475, B=3.02E-08, a=3.87082, b= 9.02583, P 11320010. Condition(6) is satisfied.
Daily infections in Iceland, where the pandemic seems to have ceased. Here again, in spite of the relatively small correlation coefficient, the behavior of the pandemic in this country looks well described by eq. (5). A=9.84E-05, B=4.99E-09 a=5.05608, b 3.11973. P 311213. Condition (6) is satisfied.
Daily infections in USA, where the peak looks already surpassed. Here again, in spite of the relatively small correlation coefficient, the behavior of the pandemic in this country looks well described by eq. (5). A=2.06E-05, B=1.50E-11, a=7.18133, b=0.955859, P=331002651. Condition (6) is satisfied.
Daily infections in Spain. Here, the characteristics of the system reveal a large dispersion of data giving a poor correlation coefficient but the curve describes well the behavior. The peak developed between the 10th to 20th days. A=0.00244986, B=4.04E- 09, a=5.63775, b=1.39028, P=46754778. Condition (6) is satisfied.
Daily infections in Italy. A=0.000357966, B=6.84E-10, a 5.72952. 1. 1.32211. p=60461826. Condition (6) is satisfied.
In Table I we present the relevant fitting parameters, including herd immunity, Thi, for countries that have not reached the epidemic peak, Mexico and Brazil. We also include the population, P, of each country.
Relevant fitting parameters for each country.
As can be seen from fitting coefficients, the exponent b drives the behavior of infections in every country. Those countries that manage well the disease expansion have b values wide larger than 1. Countries with b values close to one, as Italy and Spain, have managed the pandemics but poorly and at high costs. The recovery in both countries will be long. The same is valid for USA, that manage poorly the outbreak and its struggling with an even longer recovery to normal life. Even worse scenario is taken place in Mexico and Brazil, with very low values of b. Those countries are experiencing a big outbreak where even can get herd immunity. This, however, implies very high values of infections and mortality for the near future.
But let us briefly comment about herd immunity. Those countries that have managed to stop the outbreak, even with relative high mortality as Spain and Italy, will not reach the heard immunity. As a matter of fact, Thi can not be calculated for those countries. Then, we can see countries like Brazil where, if the way of deal with the outbreak do not change, the herd immunity will be reached. Even when it seems desirable, the ability to reach the herd immunity brings with it a high payload. That is, for a country like Brazil the herd immunity would charge more than 100 million of infected people. That is, much the same as if a non small war devastates the country.
There is an alike scenario in Mexico, but the difference here is that the value for Thi is so high that SARS-CoV-2 could even turn into a seasonal virus, at least for some years. We can expect around the same mortality but scattered over a few years.
A special observation deserves USA, where Thi tends to infinity. Here we can expect a continuous infection rate for a very long time. The outbreak is controlled but not enough to eradicate the virus. Virus will not disappear in several years but maybe the healthcare system could manage it. The virus will get endemic, and immunity will never be reached. However the infections and mortality rate associated with it, can be, hypothetically, small if compared with Mexico and Brazil.
We can also compare the speed of the outbreak in different countries. As can be seen in table Table II we calculated the time to reach a rate of 1000 infected daily for some countries (T1000). However, it should be noticed that this time is not calculated from day 0, which is always January 22, but for the approximated day when the outbreak began in the correspondent country. By example, in Brazil there was no cases at January, 22 but the first cases were detected around March, 10. So both, data fitting and T1000, were calculated from March, 10.
It can be seen clearly that Italy, Spain and USA had a faster progression from first cases than Mexico and Brazil. This fact seems to be independent from the quarantine measures taken. That is, quarantine affects the posterior evolution of the SARS-CoV-2 spread but clearly not the early one. However there are other factors that can explain the virus behavior as demographics, social behavior, mobility patterns or even weather conditions. By example, population in first world is elder, and we could say more vulnerable, than in Mexico or Brazil which could make that even when virus spreads with the same rate, the reported C OVID-19 are more in former countries.
T1000 for some populated countries and worldwide.
5. CONCLUSIONS AND OUTLOOK FOR FURTHER INVESTIGATIONS
In this work, we presented a model built using the method of analogy, in this case with a nonlinear relaxation-like behavior. With this, a good fitting with the observed behavior of the daily number of cases with time is obtained. The explicit expressions obtained may be used as a tool to approximately forecast the development of the COVID-19 pandemic in different countries and worldwide. In principle, this model could be used as a help to elaborate or change actions. This model does not incorporate any particular property of this pandemic, so we think it could be used to study pandemics with different sources. With the collected data of the pandemics at early times, using this model, it can be predicted the possibility of a peak, indefinite growth, time for herd immunity, etc.
What seems to be clear from the COVID-19 data, the fitting and the values shown in the table Table I, is that SARS-CoV-2 is far from being controlled at world level. Even when some countries appear to control the outbreak, the virus is still a menace for its health system. Furthermore, in the nowadays interconnected world it is impossible for any country to keep closed borders and pay attention to what happens only inside. All isolation measures should be halted at some time and we can expect new outbreaks in countries like Spain or Italy even after the current one could be controlled. The only way to control the spread of SARS-CoV-2 seems to be the development of a vaccine that provides the so much desired herd immunity.
Indeed, the model made possible to make an approximate forecast of the time to reach the herd immunity. This may be useful in the design of actions and policies about the pandemic.
We have introduced the T1000, that gives information about the early infection behavior in populous countries.
With a physical justification, the CD leads to the appearance of a fractal time.
A possible improvement of this model is the formal inclusion of a formulation including the dual conformable derivative [15, 16]. This will be published elsewhere,
Data Availability
All data from: https://coronavirus.jhu.edu/data/new-cases. https://www.worldometers.info/.to.
https://www.medrxiv.org/content/early/2020/05/04/2020.04.24.20078154.full.pdf.
CONFLICT OF INTEREST
The authors declare that they have no conflict of interest.
