Nauman Ahmed, Muhammad Jawaz, Ali Raza, Shahah Almutairi
{"title":"Conformable pertussis dynamics with waning immunity and natural boosting: stability regions and a structure-preserving NSFD scheme","authors":"Nauman Ahmed, Muhammad Jawaz, Ali Raza, Shahah Almutairi","doi":"10.1140/epjp/s13360-026-08365-2","DOIUrl":null,"url":null,"abstract":"<div><p>Pertussis recurrence is influenced by the interaction of vaccine-induced immunity, waning post-infection immunity, natural immune boosting, and secondary susceptibility. A six-compartment model is formulated using the difference–quotient representation of the conformable relaxation equation, with primary susceptible <i>S</i>, vaccinated <i>V</i>, infected <i>I</i>, recovered <i>R</i>, waning immunity <i>W</i>, and secondary-susceptible <i>F</i> populations. For differentiable state variables, the conformable operator is equivalent to an ordinary derivative under the nonlinear time transformation <span>\\(u=t^\\alpha /\\alpha\\)</span>. Hence, the conformable order changes the timing of transient dynamics without altering the equilibrium equations or the reproduction threshold. Existence and uniqueness, non-negativity, boundedness, disease-free and endemic equilibria, the control reproduction number, and local stability are established. Two-parameter stability regions for the disease-free regime <span>\\(\\mathscr {R}_{c}<1\\)</span> and the endemic regime <span>\\(\\mathscr {R}_{c}>1\\)</span> are obtained numerically from the spectral abscissa of the corresponding Jacobian. A nonstandard finite difference (NSFD) scheme is then constructed using the DQR/ESDDFD denominator obtained from the conformable relaxation structure. The discrete scheme preserves non-negativity, while boundedness on finite computational intervals is established by mathematical induction. Numerical simulations for several conformable orders support the analytical threshold behavior: below the threshold, infection approaches the disease-free state, whereas above the threshold, the solution approaches a positive endemic state. The conformable order changes the transient time scale but not the limiting equilibrium coordinates. The analysis therefore provides a consistent link between conformable pertussis dynamics and structure-preserving numerical approximation.\\query{Please check the edit made in the article title.}</p></div>","PeriodicalId":792,"journal":{"name":"The European Physical Journal Plus","volume":"141 10","pages":""},"PeriodicalIF":2.9000,"publicationDate":"2026-10-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"The European Physical Journal Plus","FirstCategoryId":"4","ListUrlMain":"https://link.springer.com/article/10.1140/epjp/s13360-026-08365-2","RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"PHYSICS, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
Pertussis recurrence is influenced by the interaction of vaccine-induced immunity, waning post-infection immunity, natural immune boosting, and secondary susceptibility. A six-compartment model is formulated using the difference–quotient representation of the conformable relaxation equation, with primary susceptible S, vaccinated V, infected I, recovered R, waning immunity W, and secondary-susceptible F populations. For differentiable state variables, the conformable operator is equivalent to an ordinary derivative under the nonlinear time transformation \(u=t^\alpha /\alpha\). Hence, the conformable order changes the timing of transient dynamics without altering the equilibrium equations or the reproduction threshold. Existence and uniqueness, non-negativity, boundedness, disease-free and endemic equilibria, the control reproduction number, and local stability are established. Two-parameter stability regions for the disease-free regime \(\mathscr {R}_{c}<1\) and the endemic regime \(\mathscr {R}_{c}>1\) are obtained numerically from the spectral abscissa of the corresponding Jacobian. A nonstandard finite difference (NSFD) scheme is then constructed using the DQR/ESDDFD denominator obtained from the conformable relaxation structure. The discrete scheme preserves non-negativity, while boundedness on finite computational intervals is established by mathematical induction. Numerical simulations for several conformable orders support the analytical threshold behavior: below the threshold, infection approaches the disease-free state, whereas above the threshold, the solution approaches a positive endemic state. The conformable order changes the transient time scale but not the limiting equilibrium coordinates. The analysis therefore provides a consistent link between conformable pertussis dynamics and structure-preserving numerical approximation.\query{Please check the edit made in the article title.}
百日咳复发受疫苗诱导免疫、感染后免疫减弱、自然免疫增强和继发易感性的相互作用影响。采用符合松弛方程的差商表示建立了一个六室模型,其中主要易感人群S,接种疫苗的V,感染的I,恢复的R,减弱的免疫力W和次要易感人群F。对于可微状态变量,合形算子等价于非线性时间变换\(u=t^\alpha /\alpha\)下的普通导数。因此,符合的顺序改变了瞬态动力学的时间,而不改变平衡方程或繁殖阈值。建立了系统的存在唯一性、非负性、有界性、无病和地方性平衡、控制繁殖数和局部稳定性。从对应的雅可比矩阵的谱横坐标数值计算得到了无病状态\(\mathscr {R}_{c}<1\)和地方病状态\(\mathscr {R}_{c}>1\)的双参数稳定区域。然后利用从适形松弛结构得到的DQR/ESDDFD分母构造非标准有限差分格式(NSFD)。离散格式保持了非负性,并通过数学归纳法建立了有限计算区间上的有界性。对几个符合序的数值模拟支持分析阈值行为:低于阈值,感染接近无病状态,而高于阈值,解决接近正地方性状态。顺序改变了暂态时间尺度,但没有改变极限平衡坐标。因此,分析提供了一致性百日咳动力学和结构保持数值近似之间的联系。\query{Please check the edit made in the article title.}
期刊介绍:
The aims of this peer-reviewed online journal are to distribute and archive all relevant material required to document, assess, validate and reconstruct in detail the body of knowledge in the physical and related sciences.
The scope of EPJ Plus encompasses a broad landscape of fields and disciplines in the physical and related sciences - such as covered by the topical EPJ journals and with the explicit addition of geophysics, astrophysics, general relativity and cosmology, mathematical and quantum physics, classical and fluid mechanics, accelerator and medical physics, as well as physics techniques applied to any other topics, including energy, environment and cultural heritage.