## Primary growth models included in biogrowth

### Modified Gompertz model under static conditions

Zwietering et al. (1990) proposed a reparameterization of the Gompertz model with more meaningful parameters parameters. This model predicts the population size \(N(t)\) for storage time \(t\) as a sigmoid using the following algebraic equation

\[ \log_{10} N(t) = \log_{10} N_0 + C \left( \exp \left( -\exp \left( 2.71 \frac{\mu}{C}(\lambda-t)+1 \right) \right) \right) \]

where \(N_0\) is the initial population size, \(\mu\) is the maximum growth rate, \(\lambda\) is the duration of the lag phase and \(C\) is the difference between the initial population size and the maximum population size.

### Logistic growth model

The logistic growth model can be parameterized by the following equation (Zwietering et al. 1990)

\[ \log_{10} N(t) = \log_{10} N_0 + \frac{C}{1 + \exp{ \left(\frac{4 \mu}{C}(\lambda - t)+2 \right) } } \]

where \(N_0\) is the initial population size, \(\mu\) is the maximum growth rate, \(\lambda\) is the duration of the lag phase and \(C\) is the difference between the initial population size and the maximum population size.

### Richards growth model

The Richards growth model can be parameterized by the following equation (Zwietering et al. 1990)

\[ \log_{10} N(t) = \log_{10} N_0 + C \left[1+\nu \cdot \exp{ \left(1 + \nu + \frac{\mu}{A}(1+\nu)^{1+1/\nu} (\lambda - t) \right)} \right]^{-1/\nu} \]

where \(N_0\) is the initial population size, \(\mu\) is the maximum specific growth rate, \(\lambda\) is the duration of the lag phase and \(C\) is the difference between the initial population size and the maximum population size, and \(\nu\) defines the sharpness of the transition between growth phases.

### Baranyi model under dynamic conditions

Baranyi and Roberts (1994) proposed a system of two differential equations to describe microbial growth:

\[ \frac{dN}{dt} = \frac{Q}{1+Q}\mu'\left(1 - \frac{N}{N_{max}} \right)N \]

\[ \frac{dQ}{dt}=\mu' \space Q \]

Note that the maximum specific growth rate is written as \(\mu'\). The reason for this is that **biogrowth** makes calculations in log10 scale for the population size. Therefore, for consistency with the equations for static conditions, we use a different notation in these equations. Both parameters are related by the identity \(\mu' = \mu \cdot \log(10)\).

In the Baranyi model, the deviations with respect to exponential growth are justified based on two hypotheses. It introduces the variable \(Q(t)\), which represents a theoretical substance that must be produced before the microorganism can enter the exponential growth phase. Hence, its initial value (\(Q_0\)) defines the lag phase duration (under static conditions) as \(\lambda = \frac{\log (1+1/Q_{0})}{\mu}\). On the other hand, the stationary growth phase is defined by the logistic term \((1-N/N_{max})\), which reduces the growth rate as the microorganisms reach the maximum count.

Note that the original paper by Baranyi and Roberts included an exponent, \(m\), in the term defining the stationary growth phase. However, that term is usually set to 1 by convention in predictive microbiology and, consequently, has been omitted from **biogrowth**. Also, in its original paper the specific growth rate of \(Q(t)\) was defined by a different parameter (\(\nu\)). However, because this variable does not correspond to any known substance, it is a convention in the field to set \(\nu = \mu\).

### Baranyi model under static conditions

Oksuz and Buzrul (2020) calculated the solution of the Baranyi model for static conditions given by the following equation

\[ \log_{10} N = \log_{10} N_{max} + \log_{10}{ \frac{1 + \exp{ \left( \ln (10) \mu (t-\lambda) \right)} - \exp{- \ln (10) \mu \lambda}} {\exp \left( \ln (10) \mu (t-\lambda) \right)- \exp{ \left( - \ln (10) \mu \lambda \right) + 10^{\log_{10} N_{max} - \log_{10} N_0}} } } \]

where \(N_0\) is the initial population size, \(\mu\) is the maximum specific growth rate and \(N_{max}\) is the maximum growth rate, and \(\lambda\) is the lag phase.

### Trilinear model under static conditions

Buchanan et al. (1997) proposed a trilinear model as a more simple approach to describe the growth of microbial populations. This model is defined by the piece-wise equation.

The lag phase is defined considering that, as long as \(t < \lambda\), there is no growth (i.e. \(N = N_0\)).

\[ \log_{10} N(t) = \log_{10} N_0; t \leq \lambda \] The exponential phase is described considering that during this phase, the specific growth rate is constant, with slope \(\mu_{max}\).

\[ \log_{10} N(t) = \log_{10} N_0 + \mu(t-\lambda); t\in(\lambda,t_{max}) \]

Finally, the stationary phase is modeled considering that once \(N\) reaches \(N_{max}\), it remains constant.

\[ \log_{10} N(t) = \log_{10} N_{max}; t \geq t_{max} \]

where \(t_{max}\), defined is the time required for the population size to reach \(N_{max}\) ($t_{max} = ( *{10} N*{max} - _{10} N_0 )/+ $).

### Gathering primary models metadata directly from biogrowth

The **biogrowth** package includes the `primary_model_data()`

function, which provides information about the primary models included in the package. It takes just one argument (`model_name`

). By default, this argument is `NULL`

, and the function returns the identifiers of the available models.

If a model identifier is passed, it returns a list with mode meta-information.

It includes the full reference

```
meta_info$ref
#> [1] "Buchanan, R. L., Whiting, R. C., and Damert, W. C. (1997). When is simple good enough: A comparison of the Gompertz, Baranyi, and three-phase linear models for fitting bacterial growth curves. Food Microbiology, 14(4), 313-326. https://doi.org/10.1006/fmic.1997.0125"
```

or the identifiers of the model parameters