A Nested Multi-Scale Model for COVID-19 Viral Infection

On the reduced SEI system with Nh>0 treated as constant, the paper establishes positivity/boundedness and global stability of the disease-free equilibrium E0 for R0<1 via the Castillo–Chávez framework, consistent with the candidate’s result. It also derives the endemic equilibrium E1 and sets up a Routh–Hurwitz test, simplifying C1 to μ(μ+π+γ1)(μ+γ2+dNh)(R0−1)>0, but leaves the key inequality A1B1−C1>0 as an explicit condition rather than proving it; elsewhere it asserts LAS of E1 “whenever R0>1,” creating a logical gap. The candidate solution fills this gap by proving A1B1−C1>0 in closed form, thereby completing the local stability proof for all R0>1, and provides a clean Lyapunov/LaSalle proof of global stability of E0 for R0<1. The paper also contains several typographical/notation errors (e.g., ω vs Λ in F(S,0), a stray α in the characteristic polynomial, and a self-referential formula for β*), but its main conclusions (DFE GAS for R0<1 and forward bifurcation at R0=1) remain correct. Overall, the model solution is complete and rigorous; the paper’s argument is correct where complete, but it is incomplete on the unconditional LAS of E1 for R0>1 and needs corrections to typos and derivations. Key places in the PDF: positivity/boundedness and Ω construction, R0 and equilibria, E0 global stability via Castillo–Chávez, E1 stability setup and condition, and forward bifurcation analysis with a<0,b>0.