The susceptible-infected-recovered/removed-vaccinated (SIRV) epidemics model is an important generalization of the SIR epidemics model as it accounts quantitatively for the effects of vaccination campaigns on the temporal evolution of epidemics outbreaks. Additional to the time-dependent infection ($a(t)$) and recovery ($\mu (t)$) rates, regulating the transitions between the compartments $S\to I$ and $I\to R$, respectively, the time-dependent vaccination rate $v(t)$ accounts for the transition between the compartments $S\to V$ of susceptible to vaccinated fractions. An accurate analytical approximation is derived for arbitrary and different temporal dependencies of the rates which is valid for all times after the start of the epidemics for which the cumulative fraction of new infections $J(t)\ll 1$. As vaccination campaigns automatically reduce the rate of new infections by transferring susceptible to vaccinated persons, the limit $J(t)\ll 1$ is even better fulfilled than in the SIR-epidemics model. The comparison of the analytical approximation for the temporal dependence of the rate of new infections $\jt(t)=a(t)S(t)I(t)$, the corresponding cumulative fraction $J(t)$, and $V(t)$, respectively, with the exact numerical solution of the SIRV-equations for different illustrative examples proves the accuracy of our approach. The considered illustrative examples include the cases of stationary ratios with a delayed start of vaccinations, and an oscillating ratio of recovery to infection rate with a delayed vaccination at constant rate. The proposed analytical approximation is self-regulating as the final analytical expression for the cumulative fraction $\jinf $ after infinite time allows us to check the validity of the original assumption $J(t)\le \jinf \ll 1$.

The dynamical equations of the susceptible-infected-recovered/removed (SIR) epidemics model play an important role to predict and/or analyze the temporal evolution of epidemics outbreaks. Crucial input quantities are the time-dependent infection ($a(t)$) and recovery ($\mu (t)$) rates regulating the transitions between the compartments $S\to I$ and $I\to R$, respectively. Accurate analytical approximations for the temporal dependence of the rate of new infections $\dot{J}(t)=a(t)S(t)I(t)$ and the corresponding cumulative fraction of new infections $J(t)=J(t_0)+\int_{t_0}^tdx\dot{J}(x)$ are available in the literature for either stationary infection and recovery rates and for a stationary value of the ratio $k(t)=\mu (t)/a(t)$. Here a new and original accurate analytical approximation is derived for general, arbitrary and different temporal dependencies of the infection and recovery rates which is valid for not too late times after the start of the infection when the cumulative fraction $J(t)\ll 1$ is much less than unity. The comparison of the analytical approximation with the exact numerical solution of the SIR-equations for different illustrative examples proves the accuracy of the analytical approach.

Monitored infection and vaccination rates during past past Corona waves are used to infer a posteriori two key parameter of the SIRV-epidemic model, namely the real time variation of the (i) ratio of recovery to infection rate and (ii) ratio of vaccination to infection rate. We demonstrate that using the classical SIR model the ratio between recovery and infection rates tends to overestimate the true ratio, that is of relevance in predicting the dynamics of an epidemics in the presence of vaccinations.

Monitored differential infection rates of past Corona waves are used to infer, a posteriori, the real time variation of the ratio of recovery to infection rate as key parameter of the SIR-epidemic model. From monitored Corona waves in five different countries it is found that this ratio exhibits a linear increase at early times below the first maximum of the differential infection rate before the ratios approach a nearly constant value close to unity at the time of the first maximum with small amplitude oscillations at later times. The observed time dependencies at early times and at times near the first maximum agree favorably well with the behavior of the calculated ratio for the Gaussian temporal evolution of the rate of new infections, although the predicted linear increase of the Gauss ratio at late times is not observed.

The temporal evolution of pandemics described by the susceptible-infectious-recovered (SIR)-compartment model is sensitively determined by the time dependence of the infection (a(t)) and recovery (u(t)) rates regulating the transitions from the susceptible to the infected and from the infected to the recovered compartment, respectively. Here approximated SIR-solutions for different time dependencies of the infection and recovery rates are derived which are based on the adiabatic approximation assuming time-dependent ratios k(t)=u(t)/a(t) varying slowly in comparison to the typical time characteristics of the pandemic wave. For such slow variations the available analytical approximations from the KSSIR-model, valid for a stationary value of the ratio k, are used to insert a-posteriori the adopted time-dependent ratio of the two rates. Instead of investigating endless different combinations of the time dependencies of the two rates a(t) and u(t) a suitably parameterized reduced time dependence of the ratio k(tau) is adopted. Together with the definition of the reduced time this parameterized ratio k(tau) allows us to cover a great variety of different time dependencies of the infection and recovery rates. The agreement between the solutions from the adiabatic approximation in its four different studied variants and the exact numerical solutions of the SIR-equations is remarkably good providing strong confidence in the accuracy of the proposed adiabatic approximation.

Based on hospital capacities, facts from past experience with the coronavirus disease 2019 (COVID-19) virus and the number of dark infections during the second wave (DII = 2D2), a reasonable limiting value of 140/D2 for the 7-day incidence per 100,000 persons (MSDIHT) and a second wave herd immunization threshold fraction value of 0.26 in Germany were calculated. If the MSDIHT is held below this limiting value, the German hospital system can cope with the number of new seriously infected persons without any triage decisions. On the basis of the SIRV epidemics model, the classical threshold values for herd immunization were calculated for 18 countries. For these countries, the dates regarding when herd immunization against the second COVID-19 wave will be reached were estimated

With the now available vaccination against Covid-19 it is quantitatively explored how vaccination campaigns influence the mathematical modeling of epidemics. The standard susceptible-infectious-recovered/removed (SIR) epidemic model is extended to the fourth compartment V of vaccinated persons and the vaccination rate v(t) that regulates the relation between susceptible and vaccinated persons. The vaccination rate v(t) competes with the infection (a(t)) and recovery (\mu(t)) rates in determining the time evolution of epidemics. In order for a pandemic outburst with rising rates of new infections it is required that k+b<1-2\eta, where k=\mu_0/a_0 and b=v_0/a_0 denote the initial ratios of the three rates, respectively, and \eta << 1 is the initial fraction of infected persons. Exact analytical inverse solutions t(Q) for all relevant quantities Q=[S,I,R,V] of the resulting SIRV-model in terms of Lambert functions are derived for the semi-time case with time-independent ratios k and b between the recovery and vaccination rates to the infection rate, respectively. These inverse solutions can be approximated with high accuracy yielding the explicit time-dependences Q(t) by inverting the Lambert functions. The values of the three parameters k, b and \eta completely determine the reduced time evolution the SIRV-quantities Q(\tau). The influence of vaccinations on the total cumulative number and the maximum rate of new infections in different countries is calculated by comparing with monitored real time Covid-19 data. The reduction in the final cumulative fraction of infected persons and in the maximum daily rate of new infections is quantitatively determined by using the actual pandemic parameters in different countries. Moreover, a new criterion is developed that decides on the occurrence of future Covid-19 waves in these countries. Apart from Israel this can happen in all countries considered.

We start out by deriving simple analytic expressions for all measurable amounts of cases and fatalities during a pandemic evolution exhibiting multiple waves, described by the semi-time SIR model. The approximant shares all relevant features with the exact solution, including time and position of the peak of daily new infections, as well as the asymptotic behaviors at small and large times. We derive exact analytic expressions for the early doubling time, late half decay time, and a half-early peak law, characterizing the dynamical evolution. We show, in particular, how the asymmetry of the first epidemic wave and its exponential tails are affected by the initial conditions; a feature that has no analogue in the all-time SIR model. We apply the approach to available data from different continents. Our analysis reveals that the immunity is very strongly increasing during the 2nd wave, while it was still at a very moderate level of a few percent in several countries at the end of the first wave. The wave-specific SIR parameters describing the infection and recovery rates we find to behave in a similar fashion, while their ratio k was decreasing only by a about 5% for most countries. Still, an apparently moderate change of k can have significant consequences for the relevant numbers like the final amount of infected or deceased population. As we show, the probability for an additional wave is however low in several countries due to the fraction of immune inhabitants at the end of the 2nd wave, irrespective the currently ongoing vaccination efforts. We compare with alternate approaches.

The earlier analytical analysis (part A) of the Susceptible-Infectious-Recovered (SIR) epidemics model for a constant ratio k of infection to recovery rates is extended here to the semi-time case which is particularly appropriate for modeling the temporal evolution of later (than the first) pandemic waves when a greater population fraction from the first wave has been infected. In the semi-time case the SIR model does not describe the quantities in the past; instead they only hold for times later than the initial time t=0 of the newly occurring wave. Simple exact and approximative expressions are derived for the final and maximum values of the infected, susceptible and revovered/removed population fractions as well the daily rate and cumulative number of new infections. It is demonstrated that two types of temporal evolution of the daily rate of new infections j(tau) occur depending on the values of k and the initial value of the infected fraction I(0)=eta: in the decay case for k > 1-2 eta the daily rate monotonically decreases at all positive times from its initial maximum value j(0)=eta (1-eta). Alternatively, in the peak case for k<1-2 eta the daily rate attains a maximum at a finite positive time. By comparing the approximated analytical solutions for j(tau) and J(tau) with the exact ones obtained by numerical integration, it is shown that the analytical approximations are accurate within at most only 2.5 percent. It is found that the initial fraction of infected persons sensitively influences the late time dependence of the epidemics, the maximum daily rate and its peak time. Such dependencies do not exist in the earlier investigated all-time case.

We revisit the Susceptible-Infectious-Recovered/Removed (SIR) model which is one of the simplest compartmental models. Many epidemological models are derivatives of this basic form. While an analytic solution to the SIR model is known in parametric form for the case of a time-independent infection rate, we derive an analytic solution for the more general case of a time-dependent infection rate, that is not limited to a certain range of parameter values. Our approach allows us to derive several exact analytic results characterizing all quantities, and moreover explicit, non-parametric, and accurate analytic approximants for the solution of the SIR model for time-independent infection rates. We relate all parameters of the SIR model to a measurable, usually reported quantity, namely the cumulated number of infected population and its first and second derivatives at an initial time t=0, where data is assumed to be available. We address the question on how well the differential rate of infections is captured by the Gauss model (GM). To this end we calculate the peak height, width, and position of the bell-shaped rate analytically. We find that the SIR is captured by the GM within a range of times, which we discuss in detail. We prove that the SIR model exhibits an asymptotic behavior at large times that is different from the logistic model, while the difference between the two models still decreases with increasing reproduction factor. This part A of our work treats the original SIR model to hold at all times, while this assumption will be released in part B. Releasing this assumption allows to formulate initial conditions incompatible with the original SIR model.

The Gauss model for the time evolution of the first corona pandemic wave allows to draw conclusions on the dark number of infections, the amount of heard immunization, the used maximum capacity of breathing apparati and the effectiveness of various non-pharmaceutical interventions in different countries. In Germany, Switzerland and Sweden the dark numbers are 7.4 +/- 6.1, 11.1 +/- 8.5 and 25 +/- 25, respectively. Our method of estimating dark numbers from modeling both, infection and death rates simultaneously spares these countries the laborious, time-consuming and costly medical testing for antibodies of large portions of the population. In Germany the total number of infected persons, including the dark number of infections by the first wave is estimated to be 1.06 +/- 0.60 million, corresponding to 1.28 +/- 0.72 percent of the German population. We work out direct implications from these predictions for managing the 2nd and further corona waves.

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the presence of both steady-state and transient shear and uniaxial elongational flows, supplemented by a variable link tension coefficient. Here, we provide the transient and stationary solutions of the tumbling-snake model under biaxial elongation both analytically, for small and large elongation rates, and via Brownian dynamics simulations, for the case of planar elongational flow over a wide range of rates, times, and the model parameters. We show that both the steady-state and transient first planar viscosity predictions are similar to their uniaxial counterparts, in accord with recent experimental data. The second planar viscosity seems to behave in all aspects similarly to the shear viscosity, if shear rate is replaced by elongation rate.

The Gauss model for the time evolution of the first corona pandemic wave rendered useful in the estimation of peak times, amount of required equipment, and the forecasting of fade out times. At the same time it is probably the simplest analytically tractable model that allows to quantitatively forecast the time evolution of infections and fatalities during a pandemic wave. In light of the various descriptors such as doubling times and reproduction factors currently in use to judge about lock-downs and other measures that aim to prevent spreading of the virus, we hereby provide both exact, and simple approximate relationships between the two relevant parameters of the Gauss model (peak time and width), and the transient behavior of two versions of doubling times, and three variants of reproduction factors including basic reproduction numbers.

We propose a Gauss model (GM), a map from time to the bell-shaped Gauss function to model the deaths per day and country, as a quick and simple model to make predictions on the coronavirus epidemic. Justified by the sigmoidal nature of a pandemic, i.e. initial exponential spread to eventual saturation, we apply the GM to existing data, as of April 2, 2020, from 25 countries during first corona pandemic wave and study the model's predictions. We find that logarithmic daily fatalities caused by Covid-19 are well described by a quadratic function in time. By fitting the data to second order polynomials from a statistical chi2-fit with 95\% confidence, we are able to obtain the characteristic parameters of the GM, i.e. a width, peak height and time of peak, for each country separately, with which we extrapolate to future times to make predictions. We provide evidence that this supposedly oversimplifying model might still have predictive power and use it to forecast the further course of the fatalities caused by Covid-19 per country, including peak number of deaths per day, date of peak, and duration within most deaths occur. While our main goal is to present the general idea of the simple modeling process using GMs, we also describe possible estimates for the number of required respiratory machines and the duration left until the number of infected will be significantly reduced.

We study the deterministic dynamics of N point particles moving at constant speed in a 2D table made of two polygonal urns connected by an active rectangular channel, which applies a feedback-control on the particles, inverting the horizontal component of their velocities, when their number in the channel exceeds a fixed threshold. Such a bounce--back mechanism is non-dissipative: it preserves volumes in phase space. An additional passive channel closes the billiard table forming a circuit in which a stationary current may flow. Under specific constraints on the geometry and on the initial conditions, the large N limit allows nonequilibrium phase transitions between homogeneous and inhomogeneous phases. The role of ergodicity in making a probabilistic theory applicable is discussed both for rational and irrational urns. The theoretical predictions are compared with the numerical simulation results. Connections with the dynamics of feedback-controlled biological systems are highlighted.

The susceptible-infected-recovered-vaccinated-deceased (SIRVD) epidemic compartment model extends the SIR model to include the effects of vaccination campaigns and time-dependent fatality rates on epidemic outbreaks. It encompasses the SIR, SIRV, SIRD, and SI models as special cases, with individual time-dependent rates governing transitions between different fractions. We investigate a special class of exact solutions and accurate analytical approximations for the SIRVD and SIRD compartment models. While the SIRVD and SIRD equations pose complex integro-differential equations for the rate of new infections and the fractions as a function of time, a simpler approach considers determining equations for the sum of ratios for given variations. This approach enables us to derive fully exact analytical solutions for the SIRVD and SIRD models. For nonlinear models with a high-dimensional parameter space, such as the SIRVD and SIRD models, analytical solutions, exact or accurately approximative, are of high importance and interest: not only as suitable benchmarks for numerical codes, but especially as they allow us to understand the critical behavior of epidemic outbursts as well as the decisive role of certain parameters. In the second part of our study, we apply a recently developed analytical approximation for the SIR and SIRV models to the more general SIRVD model. This approximation offers accurate analytical expressions for epidemic quantities, such as the rate of new infections and the fraction of infected persons, particularly when the cumulative fraction of infections is small. The distinction between recovered and deceased individuals in the SIRVD model affects the calculation of the death rate, which is proportional to the infected fraction in the SIRVD/SIRD cases but often proportional to the rate of new infections in many SIR models using a posteriori approach. We demonstrate that the temporal dependence of the infected fraction and the rate of new infections differs when considering the effects of vaccinations and when the real-time dependence of fatality and recovery rates diverge. These differences are highlighted for stationary ratios and gradually decreasing fatality rates. \Revised{The case of stationary ratios allows one to construct a new powerful diagnostics method to extract analytically all SIRVD model parameters from measured COVID-19 data of a completed pandemic wave.