Transmission dynamics of monkeypox virus with treatment and vaccination controls: a fractional order mathematical approach
Suvankar Majee, Soovoojeet Jana, Snehasis Barman, T. K. Kar
Abstract
Abstract Presently, monkeypox virus infection has spread worldwide in the ongoing outbreak that began in the UK. To study the transmission dynamics of monkeypox, we formulate here a seven-compartmental (five compartments for the human population and two compartments for animals or rodents) fractional-order mathematical model. The existence and uniqueness of the solution of the proposed fractional order model are examined here. The basic reproduction number for humans ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>h</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> ) and animals ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>a</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> ) are obtained through the next-generation matrix approach. Depending on the values of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>h</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>a</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> , we observed that the fractional order model has three equilibria, namely, monkeypox-free equilibrium, animal-free endemic equilibrium, and endemic equilibrium. Also, the stability of all equilibria is checked in this present article. We found that the model goes through transcritical bifurcation at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>a</mml:mi> </mml:mrow> </mml:msubsup> <mml:mo>=</mml:mo> <mml:mn>1</mml:mn> </mml:math> for any values of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>h</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> and at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>h</mml:mi> </mml:mrow> </mml:msubsup> <mml:mo>=</mml:mo> <mml:mn>1</mml:mn> </mml:math> for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>a</mml:mi> </mml:mrow> </mml:msubsup> <mml:mo><</mml:mo> <mml:mn>1</mml:mn> </mml:math> . Best of our knowledge, this is the first work where the fractional order optimal control for monkeypox is formulated and solved considering vaccination and treatment controls. Several feasible parameter values are used in the simulations to visualize and verify the findings, from which the results show that fractional order is more appropriate. Finally, parameters involved in the expression of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>h</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="fraktur">R</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>0</mml:mn> </mml:mrow> <mml:mrow> <mml:mi>a</mml:mi> </mml:mrow> </mml:msubsup> </mml:math> are scaled using the sensitivity index approach.