Abstract
Adopting an early doubling time of three days for the rate of new infections with the omicron mutant the temporal evolution of the omicron wave in different countries is predicted. The predictions are based on the susceptible-infectious-recovered/removed (SIR) epidemic compartment model with a constant stationary ratio k = μ(t)/a(t) between the infection (a(t)) and recovery (μ(t)) rate. The fixed early doubling time then uniquely relates the initial infection rate a0 to the ratio k, which therefore determines the full temporal evolution of the omicron waves. For each country three scenarios (optimistic, pessimistic, intermediate) are considered and the resulting pandemic parameters are calculated. These include the total number of infected persons, the maximum rate of new infections, the peak time and the maximum 7-day incidence per 100000 persons. Among the considered European countries Denmark has the smallest omicron peak time and the recently observed saturation of the 7-day incidence value at 2478 is in excellent agreement with the prediction in the optimistic scenario. For Germany we predict peak times of the omicron wave ranging from 32 to 38 and 45 days after the start of the omicron wave in the optimistic, intermediate and pessimistic scenario, respectively, with corresponding maximum SDI values of 7090, 13263 and 28911, respectively. Adopting Jan 1st, 2022 as the starting date our predictions implies that the maximum of the omicron wave is reached between Feb 1 and Feb 15, 2022. Rather similar values are predicted for Switzerland. Due to an order of magnitude smaller omicron hospitalization rate, due to the high percentage of vaccinated and boostered population, the German health system can cope with maximum omicron SDI value of 2800 which is about a factor 2.5 smaller than the maximum omicron SDI value 7090 in the optimistic case. By either reducing the duration of intensive care during this period of maximum, and/or by making use of the nonuniform spread of the omicron wave across Germany, it seems that the German health system can barely cope with the omicron wave avoiding triage decisions. The reduced omicron hospitalization rate also causes significantly smaller mortality rates compared to the earlier mutants in Germany. In the optimistic scenario one obtains for the total number of fatalities 7445 and for the maximum death rate 418 per day which are about one order of magnitude smaller than the beta fatality rate and total number.
I. INTRODUCTION
After being exposed to several Covid-19 outbursts the recently identified omicron mutant threatens many societies wordwide.1,2 Not many details are known sofar about its infection characteristics3,4 apart from alarming hints (1) that it is spreading at least four times quicker than the beta mutant with a short doubling time of t2 = 3 days, and (2) that the existing vaccines, taylored to prevent infections from the earlier alpha, beta, gamma and delta mutants, are less efficient against the action of the omicron mutant especially without the current booster campaigns.5–8 The alpha (α), beta (β), gamma (γ) and delta (δ) mutants have caused the first four Covid-19 waves, respectively. Positively, the omicron mutant seems to lead to on average milder symptons and thus to smaller hospitalization fractions compared to the earlier mutants.
Even with so little details known today it is of high interest to explore quantitatively the future time evolution of the omicron mutant under realistic scenarios of currently taken non-pharmaceutical interventions (NPIs). Of particular interest are reliable estimates of the maximum and total percentage of infected persons from this mutant in order to compare with the available medical capacities in different countries. In the following we provide these estimates by modeling the time evolution of the omicron wave with the susceptible-infectious-recovered/removed (SIR) epidemic compartment model.9
As in our earlier analysis10,11 – hereafter referred to as KSSIR-model – we adopt a constant stationary ratio k = μ(t)/a(t) =const. between the infection (a(t)) and recovery (μ(t)) rate regulating the transition from susceptible to infected persons and infected to recovered/removed persons in the semi-time case, respectively. As it is so far unclear whether earlier vaccinated persons are not infected by the omicron mutant, we adopt the worst case scenario here and treat the vaccinated persons as fully susceptible to the omicron mutant. However, when calculating hospitalization and mortality rates below we will account for the influence of boosted (with vaccines) persons.
As proven by Table 1 and Figs. 1 and 2 of ref.10 the KSSIR-model predicted the temporal evolution of the second wave in several countries convincingly good including the maximum rate 
of new infections and the total cumulative number (J∞) of infections as well as the initial and final second wave time dependence and the time of maximum. For the considered countries the maximum deviation in the total number of infected persons is at most 13 percent off from the later recorded values. An outstanding property of the KSSIR-model is that basically only one parameter, the ratio k of recovery and infection rates, fully determines the wave evolution in reduced time τ, whereas the influence of the initial fraction η of infected people at the onset of the modeled mutant at time t0 is only minor especially for values of η much smaller than unity. Here the reduced time
can be calculated for any arbitrary but given real time dependence of the infection rate a(t).
Second wave parameters 
in days−1, kβ, initial fraction ηβ, and the inferred second doubling time 
in days. For the omicron mutant in all countries we adopt 
to calculate the ratio of the the two doubling times 
.
(a) J∞, (b) tmax − t0, and (c) 
as a function of the only parameter k for different values of the initial fraction of infected persons η in the case of an early 3-day doubling time.
Location of the considered values of a0 and k for the omicron mutant in the investigated 12 countries with the adopted 3-day early doubling times. The symbols β represent the values for the earlier beta mutant. The solid black line is Eq. (23).
Adopting a constant infection rate a(t) = a0 is a good approximation for rapidly evolving mutant waves and not only for its initial phases, so that in this case the simple relation τ = a0(t − t0) holds between the reduced and the real time. Moreover, by determining the then two decisive parameters k and a0 from the early monitored real time evolution then allows us the accurate determination of all relevant quantities of the considered outburst. The two parameters k and a0 differ among different societies depending besides specific virus mutant properties also on the NPIs taken, the quality and ability of the health care system, and the discipline of the people in keeping distances, wearing masks and following quarantine measures. As the latter are mainly unchanged during different mutant actions it makes sense to relate the k and a0 parameters of the omicron mutant to those of the earlier beta mutant as we will adopt below.
II. RESULTS FROM THE SIR-MODEL
In terms of the reduced time (1) the KSSIR model equations read
obeying the sum constraint
at all times. In Eqs. (2)-(3) S, I and R denote the fractions of susceptible, infected and recovered/removed persons in a population, respectively, subject to the semi-time initial conditions
The rate of new infections and its cotresponding cumulative number are given by j(τ) = S(τ)I(τ) and 
, respectively, whereas 
and J(t) = J(τ).
A. Exact results
In terms of J the exact solution of the KSSIR model in the semi-time case is given by11
with ϵ = − ln(1 − η). The remaining SIR quantities are given by J(τ) as S(τ) = 1 − J(τ), I(τ) = J(τ)+kϵ+k ln[1 − J(τ)] and R(τ) = − k[ϵ + ln(1 − J(τ))]. Differentiating Eq. (5) with respect to τ readily yields for the rate of new cases
As shown before11 without the explicit inversion of the solution (5) to J(τ) one obtains for the final cumulative fraction of infected persons
with α = − (1 − η)k−1e−1/k, and for the maximum rate of new infections
occurring at
in terms of the principal (W0) and non-principal (W−1) solution of Lambert’s equation12, the well-known and documented Lambert functions. We emphasize that for small values of η ≪ 1 the results (7) and (9) are basically independent of the value of η and only determined by the parameter k. The first Eq. (8) implies
B. Approximate results
Very accurate approximations have been obtained12 for
so that Eq. (10) provides the approximation
Note that for small k < 1/8 the exact expressions (7), (9), and (10) are still useful as the approximation gives values slightly larger than unity for J∞.
The occurrence of the maximum rate of new infections (8) at positive values of the reduced peak time τmax > 0 requires values of k < 1 −2η. In this case11 the reduced peak time is well approximated by
with c0 = η(1 − η) and c1 = 1 − k − 2η,
The reduced time dependence of the rate of new infections is well approximated as
with d1 = J ∞ − (1 − k).
For a stationary infection rate a0 the corresponding real peak time is given by
Likewise, the early asymptotic reduced time behavior is well approximated12 by
corresponding to the early asymptotic real time behavior
In the considered case of a stationary infection rate Eq. (18) reduces to
implying for the early doubling time defined by 
that
where we inserted k = μ0/a0 in the case of stationary infection and recovery rates. Equation (20) will be used in the following two sections in two different ways.
The maximum 7-day incidence value per 105 persons is calculated by integrating
It is only slightly smaller than the estimate SDI ≃ 7 × 105 
from the maximum rate.
The late at times after the maximum half-decay time is given by
III. CONSEQUENCES OF EARLY 2-DAY DOUBLING TIME
For the omicron mutant the early doubling time of t2,omicron = 3 days has been reported3,4 in South Africa, Great Britain and Denmark. Adopting this value for all countries considered then provides according to Eq. (20) for the omicron mutant the relation
throughout. Using this relation in all results of the last section to eliminate a0 we find that all quantities of interest are solely determined by the parameter k. Particularly for the peak time (16) we obtain
whereas the maximum rate of new infections
Likewise the real time dependence of the rate of new infections with Eq. (14) is given by
In Fig. 1 we display the resulting dependence of J∞, 
and 
as a function of the parameter komicron ∈ [0, 1]. It can be seen that J∞ and 
decrease with increasing values of k almost independent of the initial fraction η of infected persons, except at very large values of komicron close to unity. Obviously, for comparatively small values of the total number of infected persons J∞ and the maximum rate of newly infected persons 
large values of the ratio k are required. Alternatively, the reduced time of maximum τmax decreases with increasing values of k as long as k is much smaller than 1 − 2η.
IV. OMICRON FORECAST IN INDIVIDUAL COUNTRIES
With the earlier inferred parameter values10 
and kβ for the second wave caused by the β-mutant we calculate the second wave doubling time 
for the countries listed in Table I. Adopting for all countries the short omicron doubling time 
we obtain the ratio of the doubling times also listed in the last column of Table I. With Eq. (20) we then infer for this ratio
yielding readily the relation
For each country we then consider 3 possible omicron scenarios:
the optimistic case with komicron = kβ so that the increase in the ratio r is solely due to the increase in the stationary infection rate
As noted earlier the larger the value of komicron the smaller the total cumulative number of infections J∞ and the maximum rate of new infections jmax will be. This justifies the classification of this case as optimistic.
the pessimistic scenario with
so that the increase in the ratio r is solely due to the decrease in the ratio k
Clearly, with these small values of komicron the resulting total cumulative number of infections J∞ and the maximum rate of new infections jmax will be highest, justifying the classification of this case as pessimistic. In four countries (ITA, FRA, RUS, USA) the resulting komicron,pess is negative which cannot be. In these cases we use komicron,pess = 0 and
.the intermediate case with
where half of the increase in the ratio r stems from the increase in the stationary infection rate. Then as a consequence
In Tables II, III and IV we calculate the forecast for the omicron mutant for these three scenarios, respectively. Figure 2 visualizes the relationship between a0 and k and the location of the three regimes for the 12 countries, and Fig. 3 shows the time dependence of 
and cumulative fraction J(t) of infected persons for all 12 countries.
Forecast of the omicron mutant for the optimistic case, i.e., 
, and initial fraction η = ηβ from Tab. I for this table. Columns list the final cumulative fraction J∞ of infected persons, the maximum (dimensionless) rate jmax of new infections, the cumulative fraction J0 of infected persons at peak time, the reduced peak time τmax, the peak time tmax − t0 in days, and the SDI, the maximum 7-day incidence per 105 persons. Country names are abbreviated by their α3 codes.
Forecast of the omicron mutant for the pessimistic case, i.e., 
and k = komicron,pess.
Forecast of the omicron mutant for the intermediate case, i.e., 
and k = komicron,inter.
Time dependence of the daily rate of newly infected persons, 
, as well as the cumulative fraction of infected persons, J (t), for all 12 countries.
It is obvious from these three tables that in European countries, apart from Russia with limited data reliability, Denmark has the shortest peak time of the omicron wave ranging from 16 to 22 and 27 days after the start of the omicron wave in the optimistic, intermediate and pessimistic scenario, respectively. The corresponding predicted maximum 7-day incidence values per 105 persons (SDI) are 2424, 4462 and 7148, respectively. Presently on January 10, 2022 the well-monitored data of Denmark13 indicate that the SDI has saturated at its maximum value at 2478 which is in excellent agreement with our predicted value in the optimistic case. Although preliminary this outstanding agreement is definitely encouraging and argues in favour of the optimistic scenario
Regarding Germany we predict peak times of the omicron wave ranging from 32 to 38 and 45 days after the start of the omicron wave in the optimistic, intermediate and pessimistic scenario, respectively, with corresponding maximum SDI values of 7090, 13263 and 28911, respectively. Adopting Jan 1st, 2022 as the starting date our predictions implies that the maximum of the omicron wave is reached between Feb 1 and Feb 15, 2022. In the optimistic case the total cumulative number of omicron infections will be 0.180 but can go up high to 0.812 in the pessimistic case. The late half decay times are 3.2 to 3.5 and 5.2 days in the optimistic, intermediate and pessimistic case, respectively.
Rather similar values are predicted for Switzerland. Here the peak times of the omicron wave range from 30 to 36 and 42 days after the start and the corresponding maximum SDI values are 8148, 15060 and 29259, respectively. Here with the same starting date the maximum of the omicron wave is reached between Jan 31 and Feb 13, 2022. Here, in the optimistic case the total cumulative number of omicron infections will be 0.208 but can go up high to 0.824 in the pessimistic case. The late half decay times are 3.2 to 3.6 and 5.3 days in the optimistic, intermediate and pessimistic case, respectively.
V. MEDICAL CONSEQUENCES FOR GERMANY
A. Tolerable maximum 7-day incidence value
We have argued earlier14 that the German health system can cope with maximum SDI values of 280/(hm) without any triage decisions, where m in months denotes the average duration of intensive care with access to breathing apparati for seriously infected persons, and h in units of percent indicates the percentage of people seriously infected needing access to breathing apparati in hospitals. For the earlier α and β mutants the value of h = 1 has been reasonable. Fortunately, for the omicron variant from studies in South Africa15 and Great Britain16 substantial 70–90 percent reductions as compared to earlier mutants in Covid-19 hospitalization have been reported. This strong reduction is predominantly caused by the high percentage of persons with boostered vaccination.
We therefore adopt here the value of homicron = 0.1, i.e. only one out of 1000 new infections with the omicron mutant needs to be hospitalised. Consequently, the German health system can cope with maximum omicron SDI value of 2800/m which is about a factor 2.5 smaller than the maximum omicron SDI value 7090 in the optimistic case. By either (1) reducing the duration of intensive care during this period of maximum to m = 0.5, and/or (2) by making use of the nonuniform spread of the omicron wave across Germany, appearing first in the northern states and considerably later in the southern and eastern states, combined with mutual help in hospital capacities, it seems that the German health system can cope with the omicron wave avoiding triage decisions.
B. Fatality rates and total number of fatalities
As before10,14 we assume that every second hospitalised person eventually dies from the omicron virus so that the omicron mortality rate is fomicron = 0.5homicron = 5 × 10−4 which is one order of magnitude smaller than the mortality rates of the earlier mutants. Consequently, the total fatality rate is given by D∞ = fomicronJ∞N, where N = 82.7 million denotes the German population. Likewise the maximum death rate is 
. In the optimistic scenario one obtains D∞ = 7445 and dmax = 418 per day which are about one order of magnitude smaller than the beta fatality rate and total number of fatalities of the second wave10. The main reason for these comparatively small numbers is the order of magnitude smaller hospitalization rate of the omicron mutant compared to the earlier more deadly mutants.
However, in the less likely pessimistic scenario the fatality numbers increase by a factor 4.5 to to D∞ = 33576 and dmax = 1708 which are about half of the fatality values of the second wave.
VI. SUMMARY AND CONCLUSIONS
Adopting an early doubling time of three days for the rate of new infections with the omicron mutant the temporal evolution of the omicron wave in different countries is predicted. The predictions are based on the susceptible-infectious-recovered/removed (SIR) epidemic compartment model with a constant stationary ratio k = μ(t)/a(t) between the infection (a(t)) and recovery (μ(t)) rate. The fixed early doubling time then uniquely relates the initial infection rate a0 to the ratio k, which therefore determines the full temporal evolution of the omicron waves.
As all considered countries have been exposed to earlier waves of the Covid-19 virus we relate the parameters a0 and k to those of the well-studied second wave. In the optimistic case we assume that the decrease in the early doubling time of the omicron mutant as compared to the beta mutant is solely due to a corresponding increase in the initial infection rate a0 whereas the ratio k is the same as for the beta mutant. In the pessimistic case we assume that the decrease in the early doubling times is fully caused by a corresponding decrease in the ratio k whereas the initial infection rate is the same as for the beta mutant. In the intermediate scenario half of the decrease in the early doubling time is assigned to a corresponding increase of the initial infection rate and a corresponding decrease of the ratio k. For 12 countries these three scenarios (optimistic, pessimistic, intermediate) are considered and the resulting pandemic parameters are calculated. These include the total number of infected persons, the maximum rate of new infections, the peak time and the maximum 7-day incidence per 100000 persons.
Among the considered European countries Denmark has the smallest omicron peak time and the recently observed saturation of the 7-day incidence value at 2478 is in excellent agreement with our prediction in the optimistic case. For Germany we predict peak times of the omicron wave ranging from 32 to 38 and 45 days after the start of the omicron wave in the optimistic, intermediate and pessimistic scenario, respectively, with corresponding maximum SDI values of 7090, 13263 and 28911, respectively. Adopting Jan 1st, 2022 as the starting date our predictions implies that the maximum of the omicron wave is reached between Feb 1 and Feb 15, 2022. In the optimistic case the total cumulative number of omicron infections will be 0.180 but can go up high to 0.812 in the pessimistic case. The late half decay times are 3.2 to 3.5 and 5.2 days in the optimistic, intermediate and pessimistic case, respectively.
Rather similar values are predicted for Switzerland. Here the peak times of the omicron wave range from 30 to 36 and 42 days after the start and the corresponding maximum SDI values are 8148, 15060 and 29259, respectively. Here with the same starting date the maximum of the omicron wave is reached between Jan 31 and Feb 13, 2022. Here, in the optimistic case the total cumulative number of omicron infections will be 0.208 but can go up high to 0.824 in the pessimistic case. The late half decay times are 3.2 to 3.6 and 5.3 days in the optimistic, intermediate and pessimistic case, respectively.
Adopting an order of magnitude smaller omicron hospitalization rate thanks to the high percentage of vaccinated and boostered population we conclude that the German health system can cope with maximum omicron SDI value of 2800 which is about a factor 2.5 smaller than the maximum omicron SDI value 7090 in the optimistic case. By either reducing the duration of intensive care during this period of maximum, and/or by making use of the nonuniform spread of the omicron wave across Germany, it seems that the German health system can barely cope with the omicron wave avoiding triage decisions.
The reduced omicron hospitalization rate also causes significantly smaller mortality rates compared to the earlier mutants in Germany. In the optimistic scenario one obtains for the total number of fatalities D∞ = 7445 and for the maximum death rate dmax = 418 per day which are about one order of magnitude smaller than the beta fatality rate and total number. In the less likely pessimistic scenario these numbers increase by a factor 4.5.
Data Availability
All data produced in the present work are contained in the manuscript
