{"title":"A spatial SIS point pattern model with movement.","authors":"Eka Suci Pramana Sari, Nanang Susyanto, Fajar Adi-Kusumo, Atina Husnaqilati, Fugo Takasu","doi":"10.3934/mbe.2026080","DOIUrl":"https://doi.org/10.3934/mbe.2026080","url":null,"abstract":"<p><p>Spatial epidemic dynamics are strongly affected by local interactions and by the movement of individuals. This study develops a susceptible-infected-susceptible (SIS) model based on point pattern dynamics with individual movement. The model extends the spatial SIS point pattern framework of Hamada and Takasu by allowing susceptible and infected points to move in continuous space. Disease transmission is represented by a distance-dependent infection kernel, and Gaussian, step-function, and exponential kernels are compared. Heuristic singlet and pair-density equations are derived to describe how infection, recovery, birth, death, and movement affect the spatial distribution of susceptible and infected individuals. Numerical simulations are used to illustrate the effects of kernel shape, infection range, and movement on singlet densities, pair correlations, and the temporal evolution of point patterns. The results show how spatial movement and distance-dependent transmission jointly shape infection spread in a recurrent-infection framework.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2208-2227"},"PeriodicalIF":2.6,"publicationDate":"2026-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664817","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Using active subspaces to explore discrepancies between global and local parameter sensitivities in a Lotka-Volterra system.","authors":"Huiyan Zou, Allison L Lewis","doi":"10.3934/mbe.2026079","DOIUrl":"https://doi.org/10.3934/mbe.2026079","url":null,"abstract":"<p><p>Global sensitivity metrics are essential tools for assessing parameter importance in complex models, particularly when precise information about parameter values is unavailable. In many cases, such metrics are used to provide parameter rankings that allow for necessary dimension reduction in moderate-to-high dimensional systems. However, globally derived sensitivity results may obscure localized variability in parameter sensitivities, resulting in misleading conclusions about parameter importance and ensuing consequences for subsequent tasks such as model calibration and surrogate model construction. In this study, we illustrated how discrepancies between globally and locally based sensitivity information may arise for an emerging sensitivity metric based on active subspace methodology, as well as for other commonly used sensitivity techniques. In response, we outlined a framework that exploits the active subspace to evaluate the stability of parameter sensitivities over the admissible parameter space. This analysis allows one to determine the subregions of the parameter space in which a globally derived sensitivity metric may be considered representative of local behavior. The proposed framework was illustrated on a collection of simple examples for ease of visualization, as well as a variation of the Lotka-Volterra model used to study interactions between competing tumor cell lines, for which we demonstrated how these issues may be exacerbated in higher dimensions. Our findings suggest that globally derived sensitivity information should be treated with caution, and that incorporating analysis on local subregions may improve robustness and accuracy in downstream modeling tasks.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2178-2207"},"PeriodicalIF":2.6,"publicationDate":"2026-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664797","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Mathematical analysis of a vector model using larvicides and <i>Wolbachia</i>-infected male mosquitoes: an optimal and impulsive control approach.","authors":"Wendpanga Birba, Boureima Sangaré, Abdoulaye Kaboré","doi":"10.3934/mbe.2026078","DOIUrl":"https://doi.org/10.3934/mbe.2026078","url":null,"abstract":"<p><p>In this paper, we develop and analyze a mathematical model for mosquito population control, incorporating the release of <i>Wolbachia</i>-infected male mosquitoes and the application of larvicides. Two complementary approaches are considered: An optimal control model with continuous-time control functions, and an impulsive control model accounting for periodic releases of <i>Wolbachia</i>-infected male mosquitoes and larvicide. The main objective is to minimize the population of fertilized females, which are responsible for the transmission of vector-borne diseases, while reducing the associated biological, social, and economic costs. Using the Pontryagin maximum principle, we characterize the optimal controls and demonstrate the existence of admissible solutions. Within the framework of the impulsive model, the existence and stability conditions for periodic solutions are rigorously derived, allowing the identification of critical thresholds for the larvicide concentration and the minimum density of <i>Wolbachia</i>-infected mosquitoes required to achieve a significant reduction in the target population. Numerical simulations illustrate the enhanced effectiveness of the combined intervention strategy, demonstrating both faster response times and a greater impact on population suppression. These findings emphasize the necessity of a scientifically grounded, optimized, and integrated approach in the design and implementation of vector control programs.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2132-2177"},"PeriodicalIF":2.6,"publicationDate":"2026-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664801","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An adaptive div-curl-grad least-squares finite element method for the Cahn-Hilliard equation.","authors":"Chad R Westphal","doi":"10.3934/mbe.2026077","DOIUrl":"https://doi.org/10.3934/mbe.2026077","url":null,"abstract":"<p><p>This paper develops a new least-squares finite element approach for the efficient approximation of solutions of the Cahn-Hilliard equation. We discretize in time with either implicit backward Euler or Crank-Nicolson and use a regularized Newton iteration for the nonlinearity. The innermost iteration is reformulated as a coupled first-order div-curl-grad system, which is efficiently solved by a least-squares minimization approach using standard Lagrange finite elements for each unknown. Theoretical findings include results on the optimal finite element convergence of the variational problem, convergence of the nonlinear iteration, error bounds for the time integration schemes, and discrete asymptotic mass conservation and energy dissipation in the least-squares sense. The basic iterative approach is expanded to include adaptivity in both time and space. Numerical examples are provided to confirm the theoretical results and to demonstrate effectiveness of the proposed method for a range of test problems.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2110-2131"},"PeriodicalIF":2.6,"publicationDate":"2026-07-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664730","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Learning contact policies for SEIR epidemics on networks: A mean-field game approach.","authors":"Weinan Wang","doi":"10.3934/mbe.2026076","DOIUrl":"https://doi.org/10.3934/mbe.2026076","url":null,"abstract":"<p><p>In this paper, we develop a mean-field game model for SEIR epidemics on heterogeneous contact networks, where individuals choose state-dependent contact effort to balance infection losses against the social and economic costs of isolation. The Nash equilibrium is characterized by a coupled Hamilton-Jacobi-Bellman/Kolmogorov system across degree classes. An important feature of the SEIR setting is the exposed compartment: the incubation period separates infection from infectiousness and changes incentives after infection occurs. In the baseline formulation, exposed agents optimally maintain full contact, while susceptible agents reduce contact according to an explicit best-response rule driven by infection pressure and the value gap. We also discuss extensions that yield nontrivial exposed precaution by introducing responsibility or compliance incentives. We establish existence of equilibrium via a fixed-point argument and prove the uniqueness under a suitable monotonicity condition. The analysis identifies a delay in the onset of precaution under longer incubation, which can lead to weaker behavioral responses and larger outbreaks. Numerical experiments illustrate how network degree and the cost exponent shape equilibrium policies and epidemic outcomes.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2083-2109"},"PeriodicalIF":2.6,"publicationDate":"2026-07-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664726","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Riccardo Sacco, Greta Chiaravalli, Giovanna Guidoboni, Anita Layton, Gal Antman, Alice Verticchio, Brent Siesky, Thomas A Ciulla, Alon Harris
{"title":"Modeling the relation between intraocular pressure and human trabecular meshwork fluid-mechanical properties.","authors":"Riccardo Sacco, Greta Chiaravalli, Giovanna Guidoboni, Anita Layton, Gal Antman, Alice Verticchio, Brent Siesky, Thomas A Ciulla, Alon Harris","doi":"10.3934/mbe.2026075","DOIUrl":"https://doi.org/10.3934/mbe.2026075","url":null,"abstract":"<p><p>Elevated introcular pressure (IOP) is a major risk factor for primary open angle glaucoma (POAG), the second leading cause of blindness worldwide. In this article, we investigate the relation between IOP increase and fluid-mechanical properties of the trabecular meshwork (TM), which is the main pathway of aqueous humor (AH) outflow in the human eye. TM is represented as an axisymmetric deformable biphasic porous medium (solid and fluid) with hydraulic permeability nonlinearly depending on the local AH pressure. Solving Darcy's law provides a TM hydraulic resistance (TMR) depending on the radial pressure drop, while TM stiffness is included multiplying the TMR by a shape function of TM Young's modulus constructed using postmortem measurements of human eye outflow facility. Simulations of AH flow based on an electric equivalent scheme of the eye indicate that TMR increases with TM stiffness, leading to a progressive increase of IOP until over the upper limit of the normal IOP range. The obtained results shed light on the impact of microscopic properties of the eye on its macroscopic function, and support the integration of experimental studies, mathematical modeling, and data analysis to develop an optimized patient-specific therapy to cure POAG.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2055-2082"},"PeriodicalIF":2.6,"publicationDate":"2026-07-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664799","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Chenyu Wu, Nuozhou Wang, Casey Garner, Kevin Leder, Shuzhong Zhang
{"title":"Novel optimization techniques for inferring heterogeneous population dynamics.","authors":"Chenyu Wu, Nuozhou Wang, Casey Garner, Kevin Leder, Shuzhong Zhang","doi":"10.3934/mbe.2026074","DOIUrl":"https://doi.org/10.3934/mbe.2026074","url":null,"abstract":"<p><p>In this paper, we introduce a new optimization algorithm that is well suited to solve parameter estimation problems that arise when inferring heterogeneous population dynamics. In these estimation problems, parameter estimation is complicated by the presence of two types of constraints: inequality constraints (e.g., non-negativity and boundedness of rates (so-called box-constraints)) and equality constraints that arise due to the need of the population fractions to sum to one. We call our new method cubic regularized Newton with affine scaling (CRNAS). In contrast to so-called first-order methods, which solely rely on the gradient of the objective function, our method utilizes the Hessian of the objective. As a result, it is able to focus on points that satisfy the second-order optimality conditions, as opposed to first-order methods that simply converge to critical points. This is an important feature in parameter estimation problems, where the objective function is often non-convex; as a result, there can be many critical points, which makes it nearly impossible to identify the global minimum. We use an affine scaling approach to handle a wide class of constraints, including equality constraints. We establish that CRNAS identifies a point that satisfies $ epsilon $-approximate second-order optimality conditions within $ O(epsilon^{-3/2}) $ iterations. Finally, we compare CRNAS with MATLAB's optimization solver <i>fmincon</i> on three different test problems. These test problems all feature mixtures of heterogeneous populations, a problem setting that CRNAS is particularly well-suited for. Our numerical simulations show that CRNAS has a favorable performance, thereby performing comparable, if not better than, <i>fmincon</i> in accuracy and computational cost for most of our examples.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"2018-2054"},"PeriodicalIF":2.6,"publicationDate":"2026-07-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664794","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
León A Valencia, Raúl Alejandro Morán-Vásquez, Duván H Cataño Salazar
{"title":"Modeling transmission intensity in SI epidemics via CIR and Jacobi processes: Asymptotic results and preliminary intervention strategies.","authors":"León A Valencia, Raúl Alejandro Morán-Vásquez, Duván H Cataño Salazar","doi":"10.3934/mbe.2026073","DOIUrl":"https://doi.org/10.3934/mbe.2026073","url":null,"abstract":"<p><p>This paper studies an SI epidemic model with stochastic transmission rates of the form $ (beta_t = varphi(t)P_t:tgeq0) $, where $ varphi(t) $ is a deterministic modulation function and $ P_t $ is a positive stochastic process. We show that the asymptotic behavior of the epidemic is determined by the integrated intensity process $ (H_t = int_0^t beta_s, ds:tgeq0) $. We consider two stochastic models for $ (P_t:tgeq0) $: the bounded Jacobi process and the Cox-Ingersoll-Ross (CIR) process. Both preserve positivity, but differ in the support of their sample paths. In the non-modulated regime $ (varphiequiv1) $, the CIR framework allows explicit expressions for Laplace transforms and probabilistic bounds associated with the integrated intensity process. Additionally, we present numerical simulations in two regimes: the non-modulated case $ (varphi(t) = 1) $ and the exponentially damped case $ (varphi(t) = e^{-alpha t}) $. The simulations show that the bounded and unbounded structures of the stochastic transmission processes produce different tail behaviors, particularly in high-volatility regimes.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"1995-2017"},"PeriodicalIF":2.6,"publicationDate":"2026-07-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664756","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Dynamic analysis of a SEAIRS model on co-evolution of human behavior changes in COVID-19 epidemic.","authors":"Xiaoyu Li, Nan Wang, Zhiming Li","doi":"10.3934/mbe.2026072","DOIUrl":"10.3934/mbe.2026072","url":null,"abstract":"<p><p>Considering that people usually take non-drug protective measures (making behavior changes), this paper proposes a new SEAIRS model by combining behavior changes during the spread of the COVID-19 epidemic. First, we prove the positivity and boundedness of the solution of the model. Furthermore, the existence and stability of the equilibrium points are discussed in three behavior scenarios $ (B_c = 0, B_c = 1, 0 < B_c < 1, $ where $ B_c $ represents the proportion of susceptible individuals ($ S $), exposed individuals ($ E $), and asymptomatic infected individuals ($ A $) who adopt preventive behavior changes). Finally, we illustrate the model's behaviour using COVID-19 data from Romania, and summarize some valuable results: (i) the spread of diseases can be effectively controlled by reducing the contact between susceptible, exposed, and infected individuals; (ii) increasing the proportion of behavioral changes can significantly reduce the number of infected individuals. The model suggests two main theoretical measures to promote behavior changes within the framework: one is to increase the return difference between the two behaviors by adjusting relevant parameters, and the other is to increase the value of the information function.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"1971-1994"},"PeriodicalIF":2.6,"publicationDate":"2026-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664751","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Paul Klockenkemper, Folashade Agusto, Wandi Ding, Gavin McGee, Rachel Leander
{"title":"Modeling the impact of host heterogeneity on the risk and dynamics of a West Nile virus epidemic.","authors":"Paul Klockenkemper, Folashade Agusto, Wandi Ding, Gavin McGee, Rachel Leander","doi":"10.3934/mbe.2026071","DOIUrl":"https://doi.org/10.3934/mbe.2026071","url":null,"abstract":"<p><p>West Nile virus (WNV) is a mosquito-borne arbovirus with significant ecological and public health implications. Its transmission cycle involves avian hosts and mosquito vectors. Many factors, including host diversity and mosquito population dynamics are known to shape epidemic patterns. In this study, we extend previous work by incorporating multiple host types, horizontal transmission among hosts, and mosquito feeding preference into a WNV model. Using a system of ordinary differential equations, we analyze the impact of variable host competence, host abundance, host community structure, and mosquito biting preference on epidemic metrics. We derive an expression for the basic reproduction number in multihost systems and qualify the impact of less competent hosts on its value. Numerical simulations explore the impact of structural variations in the model on epidemic metrics and assess the potential for less competent hosts and mosquito feeding preference to dilute the epidemic. Our findings elucidate the influence of host heterogeneity, horizontal transmission in hosts, and mosquito biting preference on the severity and dynamics of a WNV epidemic.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 7","pages":"1926-1970"},"PeriodicalIF":2.6,"publicationDate":"2026-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148664733","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}