# Does this dynamical system have another conserved quantity?

For the 3D system of ODEs:

$$\begin{eqnarray}\dot{x} &=& -\beta x y \\ \dot{y} &=& \beta x y + \hat{\beta} z y - \delta y \\ \dot{z} &=& -\hat{\beta} z y + \delta y, \end{eqnarray}$$

the quantities $$Q_1 = x+y+z$$ (assume $$Q_1=1$$) and $$Q_2 = (\delta/\hat{\beta}-z)(x)^{\hat{\beta}/\beta}$$ are conserved quantities (i.e., $$\dot{Q}_1 = \dot{Q}_2=0$$).

$$Q_2$$ can be found by first reducing the 3D system to 2D by using $$Q_1$$ (e.g. $$z = 1-x-y$$) and then computing $$dy/dx$$, which can be solved analytically.

My question is about the following generalization of the previous system to the following 6D system of ODEs:

$$\begin{eqnarray}\dot{x}_1 &=& -\beta x_1 y_2 \\ \dot{y}_1 &=& \beta x_1 y_2 + \hat{\beta} z_1 y_2 - \delta y_1 \\ \dot{z}_1 &=& -\hat{\beta} z_1 y_2 + \delta y_1, \end{eqnarray}$$

$$\begin{eqnarray}\dot{x}_2 &=& -\beta x_2 y_1 \\ \dot{y}_2 &=& \beta x_2 y_1 + \hat{\beta} z_2 y_1 - \delta y_2 \\ \dot{z}_2 &=& -\hat{\beta} z_2 y_1 + \delta y_2. \end{eqnarray}$$

It is easy to see that $$Q_{11} = x_1 + y_1 + z_1$$ and $$Q_{12} = x_2 + y_2 + z_2$$ are conserved quantities for the new system.

Are there any other conserved quantities in the new system? If there are, what are they?

Please assume that $$Q_{11} = Q_{12} = 1$$.

I have played with many analogous forms to $$Q_2$$ without success and I would like some help from the community to help me find the other conserved quantity(ies).

Even without assuming that $$Q_{11} = Q_{12} = 1$$, you can use $$y_2 = -\frac{1}{\beta} \frac{\dot{x}_1}{x_1}$$ (and the equivalent for $$y_1$$) to obtain the conserved quantity
$$Q_2 = \beta(y_1-y_2) + (\beta-\hat{\beta})(x_1-x_2) - (\delta + \hat{\beta} Q_{12}) \log x_2 + (\delta + \hat{\beta}Q_{11}) \log x_1 .$$
\begin{align} \dot{Q}_2 =& \beta (\dot{y}_1 - \dot{y}_2) + (\beta-\hat{\beta})(\dot{x}_1 - \dot{x}_2) - (\delta + \hat{\beta} Q_{12}) \frac{\dot{x}_2}{x_2} + (\delta + \hat{\beta} Q_{11}) \frac{\dot{x}_1}{x_1}\\ =& \beta\left(\beta x_1 y_2 + \hat{\beta}(Q_{11} - x_1 - y_1)y_2 - \delta y_1 - \beta x_2 y_1 - \hat{\beta}(Q_{12} - x_2 - y_2)y_1 + \delta y_2\right) \\ & + (\beta-\hat{\beta})(-\beta x_1 y_2 + \beta x_2 y_1) - (\delta + \hat{\beta} Q_{12})(-\beta y_1) + (\delta + \hat{\beta}Q_{11})(-\beta y_2)\\ =& 0. \end{align}