Fuente 2024 model

\begin{align*}
  \frac{d Q\_active}{dt} &= \frac{k3 \cdot (1 - Q\_active)}{1 + {\frac{keq\_NPQ}{H\_lumen}}^{n\_NPQ}} \\
  & - Q\_active \cdot k4 \\ 
  \frac{d PQ}{dt} &= - 0.5 \cdot PQ \cdot RCII\_closed \cdot k1p - PQH\_2 \cdot RCII\_open \cdot k1m \\
  & + 0.5 \cdot PQH\_2 \cdot PSI\_ox \cdot k2p - PQ \cdot PSI\_red \cdot k2m \\
  & + PQH\_2 \cdot k\_X \\ 
  \frac{d PSI\_ox}{dt} &= \frac{L\_PSI \cdot osc\_light \cdot sigma\_PSI\_0 \cdot stoic\_PSI \cdot (stoic\_PSI - PSI\_ox)}{L\_PSI + osc\_light} \\
  & - PQH\_2 \cdot PSI\_ox \cdot k2p - PQ \cdot PSI\_red \cdot k2m \\ 
  \frac{d H\_lumen}{dt} &= \frac{bH}{N\_A \cdot V\_lumen} \cdot osc\_light \cdot sigma\_PSII \cdot stoic\_PSII \cdot (1 - RCII\_closed) \\
  & + \frac{bH}{N\_A \cdot V\_lumen} \cdot PQH\_2 \cdot PSI\_ox \cdot k2p - PQ \cdot PSI\_red \cdot k2m \\
  & - \frac{4.666666666666667 \cdot V\_stroma \cdot bH}{V\_lumen} \cdot k5 \cdot (ADP\_st - \frac{ATP\_st \cdot {\frac{H\_stroma}{H\_lumen}}^{4.666666666666667}}{cEqP}) \\
  & - k7 \cdot (H\_lumen - H\_stroma) \\ 
  \frac{d ATP\_st}{dt} &= k5 \cdot (ADP\_st - \frac{ATP\_st \cdot {\frac{H\_stroma}{H\_lumen}}^{4.666666666666667}}{cEqP}) \\
  & - ATP\_st \cdot k6
\end{align*}