Pub Date : 2025-11-24DOI: 10.1016/j.matcom.2025.11.029
Xing Guo , Lianghao Ji , Shasha Yang , Rongjian Liu
This paper investigates the impulsive consensus control for nonlinear multiagent systems (MASs) with semi-Markov switching topologies (semi-MSTs) under deception attacks. A novel switching-triggered impulsive control scheme is proposed, which creatively uses topology switching as an event to drive impulsive control. This method decouples the execution of impulsive control from conventional time-triggered or even-triggered approaches by introducing topology switching as the primary determinant of control updates, thereby enabling full utilization of topological information. The network topology condition for achieving impulsive consensus of MASs with semi-MSTs under this scheme can be relaxed to only require that the union of all switching subtopologies contains a spanning tree. Furthermore, deception attacks occurring in communication channels are considered, which can cause incorrect state information transmission. The random variables describing whether deception attacks occur obey Bernoulli distribution. Sufficient conditions for realizing secure impulsive consensus control of MASs with semi-MSTs under deception attacks and the upper bound on the mean square error between the leader and the followers are given. Finally, an example is presented to validate the effectiveness of the main results.
{"title":"Switching-triggered control for multiagent systems with semi-MSTs under deception attacks","authors":"Xing Guo , Lianghao Ji , Shasha Yang , Rongjian Liu","doi":"10.1016/j.matcom.2025.11.029","DOIUrl":"10.1016/j.matcom.2025.11.029","url":null,"abstract":"<div><div>This paper investigates the impulsive consensus control for nonlinear multiagent systems (MASs) with semi-Markov switching topologies (semi-MSTs) under deception attacks. A novel switching-triggered impulsive control scheme is proposed, which creatively uses topology switching as an event to drive impulsive control. This method decouples the execution of impulsive control from conventional time-triggered or even-triggered approaches by introducing topology switching as the primary determinant of control updates, thereby enabling full utilization of topological information. The network topology condition for achieving impulsive consensus of MASs with semi-MSTs under this scheme can be relaxed to only require that the union of all switching subtopologies contains a spanning tree. Furthermore, deception attacks occurring in communication channels are considered, which can cause incorrect state information transmission. The random variables describing whether deception attacks occur obey Bernoulli distribution. Sufficient conditions for realizing secure impulsive consensus control of MASs with semi-MSTs under deception attacks and the upper bound on the mean square error between the leader and the followers are given. Finally, an example is presented to validate the effectiveness of the main results.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 69-81"},"PeriodicalIF":4.4,"publicationDate":"2025-11-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618618","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-24DOI: 10.1016/j.matcom.2025.11.026
Sayan Mandal, Pankaj Kumar Tiwari
In this study, we develop and analyze a deterministic prey–predator model where predators are generalist and follows modified Beverton–Holt-type growth dynamics due to additional foods, incorporating prey refuge. We also analyze system’s dynamics in the presence of seasonal and environmental fluctuations. Our key attention is on emphasizing the effects of density-dependent prey refuge and additional food availability on species coexistence and stability. Through theoretical analysis, we establish the feasibility of solutions under both autonomous and seasonal settings, identifying local stability criteria and the existence of positive periodic solutions. Our numerical results reveal that when there are no refuge and additional food, the system undergoes transcritical and supercritical Hopf bifurcations, leading to stable coexistence or population oscillations. However, the provision of prey refuge increases the number of coexistence equilibria, inducing bistability and, at higher levels, potential predator extinction. On variations of the levels of refuge and additional food, the system transitions from bistability to tristability, displaying complex dynamical shifts. However, the time variation of parameters significantly alter population stability, triggering periodic oscillations, chaotic regimes, and potential predator extinction under high-intensity of seasonal strengths. Sensitivity analysis confirms chaotic behavior under specific seasonal conditions, reinforcing the unpredictability of ecological dynamics. Notably, environmental noise can drive transitions between multiple equilibria, with moderate noise promoting coexistence and high noise leading to species extinction.
{"title":"Predator–prey interactions: How prey refuge, additional food, seasonality, and stochasticity shape ecological stability?","authors":"Sayan Mandal, Pankaj Kumar Tiwari","doi":"10.1016/j.matcom.2025.11.026","DOIUrl":"10.1016/j.matcom.2025.11.026","url":null,"abstract":"<div><div>In this study, we develop and analyze a deterministic prey–predator model where predators are generalist and follows modified Beverton–Holt-type growth dynamics due to additional foods, incorporating prey refuge. We also analyze system’s dynamics in the presence of seasonal and environmental fluctuations. Our key attention is on emphasizing the effects of density-dependent prey refuge and additional food availability on species coexistence and stability. Through theoretical analysis, we establish the feasibility of solutions under both autonomous and seasonal settings, identifying local stability criteria and the existence of positive periodic solutions. Our numerical results reveal that when there are no refuge and additional food, the system undergoes transcritical and supercritical Hopf bifurcations, leading to stable coexistence or population oscillations. However, the provision of prey refuge increases the number of coexistence equilibria, inducing bistability and, at higher levels, potential predator extinction. On variations of the levels of refuge and additional food, the system transitions from bistability to tristability, displaying complex dynamical shifts. However, the time variation of parameters significantly alter population stability, triggering periodic oscillations, chaotic regimes, and potential predator extinction under high-intensity of seasonal strengths. Sensitivity analysis confirms chaotic behavior under specific seasonal conditions, reinforcing the unpredictability of ecological dynamics. Notably, environmental noise can drive transitions between multiple equilibria, with moderate noise promoting coexistence and high noise leading to species extinction.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 121-148"},"PeriodicalIF":4.4,"publicationDate":"2025-11-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618620","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-22DOI: 10.1016/j.matcom.2025.11.027
Yi Zhang , Jin Song , Wenjie Zuo , Zhengdi Zhang
This paper aims to explore compound relaxation oscillations and underlying mechanisms in the dynamical system with periodic non-smoothness, focusing on the effect of non-smooth bifurcations on compound relaxation oscillations. Based on the Rayleigh–Duffing system with external excitation, a modified non-smooth dynamical system is developed by introducing a periodic term that represents discontinuous external influences, such as wave-induced forces in ship rolling dynamics. Various non-smooth bifurcation phenomena are systematically investigated, including non-smooth homoclinic bifurcation, C-bifurcation, persistence bifurcation, and non-smooth fold bifurcation. Five different oscillation modes are demonstrated through numerical simulations, and their mechanisms are revealed in combination with the slow–fast analysis. It is found that the non-smooth homoclinic bifurcation significantly alters the oscillation process and induces transitions between stable states. The C-bifurcation has less effect on the oscillation mode even though it changes the topology of limit cycles. Different types of boundary equilibrium bifurcations lead to substantial changes in the stability and structure of compound relaxation oscillations. In addition, two types of coexisting attractors are identified through the basin of attraction, indicating multistability that gives rise to different oscillation modes.
{"title":"Compound relaxation oscillations in a modified Rayleigh–Duffing system with periodic non-smoothness","authors":"Yi Zhang , Jin Song , Wenjie Zuo , Zhengdi Zhang","doi":"10.1016/j.matcom.2025.11.027","DOIUrl":"10.1016/j.matcom.2025.11.027","url":null,"abstract":"<div><div>This paper aims to explore compound relaxation oscillations and underlying mechanisms in the dynamical system with periodic non-smoothness, focusing on the effect of non-smooth bifurcations on compound relaxation oscillations. Based on the Rayleigh–Duffing system with external excitation, a modified non-smooth dynamical system is developed by introducing a periodic term that represents discontinuous external influences, such as wave-induced forces in ship rolling dynamics. Various non-smooth bifurcation phenomena are systematically investigated, including non-smooth homoclinic bifurcation, C-bifurcation, persistence bifurcation, and non-smooth fold bifurcation. Five different oscillation modes are demonstrated through numerical simulations, and their mechanisms are revealed in combination with the slow–fast analysis. It is found that the non-smooth homoclinic bifurcation significantly alters the oscillation process and induces transitions between stable states. The C-bifurcation has less effect on the oscillation mode even though it changes the topology of limit cycles. Different types of boundary equilibrium bifurcations lead to substantial changes in the stability and structure of compound relaxation oscillations. In addition, two types of coexisting attractors are identified through the basin of attraction, indicating multistability that gives rise to different oscillation modes.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 82-94"},"PeriodicalIF":4.4,"publicationDate":"2025-11-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618555","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This work introduces a ratio-dependent Holling–Tanner predator–prey model with the Allee effect in prey and then discretizes the introduced model through the Euler forward scheme. A brief discussion is held on the stability analysis for several fixed points in the discretized model. Several types of bifurcations, including codimension one and two bifurcations, are demonstrated in this study. Codimension-1 bifurcation, which covers Neimark–Sacker and flip bifurcations, and codimension-2 bifurcations, which include strong resonance 1:2, 1:3, and 1:4 at a positive fixed point. Various critical states under non-degeneracy conditions are computed using the critical normal form coefficient approach for each bifurcation. The model displays complex dynamical behaviours, like quasi-periodic orbits and chaotic sets. Additionally, the system’s chaos was managed by the development of control mechanisms, such as the OGY methodology. It has been established that bifurcation and chaos can be stabilized under certain circumstances. A thorough numerical simulation further supports our analytical findings, which include stability regions, bifurcation curves in 2D & 3D, phase plots, and the maximal Lyapunov exponent, etc.
{"title":"Multiple bifurcations and managing chaos: A discretized ratio-dependent Holling–Tanner predator–prey model with Allee effect in prey","authors":"Md. Jasim Uddin , Savita Boora , Sarker Md. Sohel Rana , Pradeep Malik","doi":"10.1016/j.matcom.2025.11.024","DOIUrl":"10.1016/j.matcom.2025.11.024","url":null,"abstract":"<div><div>This work introduces a ratio-dependent Holling–Tanner predator–prey model with the Allee effect in prey and then discretizes the introduced model through the Euler forward scheme. A brief discussion is held on the stability analysis for several fixed points in the discretized model. Several types of bifurcations, including codimension one and two bifurcations, are demonstrated in this study. Codimension-1 bifurcation, which covers Neimark–Sacker and flip bifurcations, and codimension-2 bifurcations, which include strong resonance 1:2, 1:3, and 1:4 at a positive fixed point. Various critical states under non-degeneracy conditions are computed using the critical normal form coefficient approach for each bifurcation. The model displays complex dynamical behaviours, like quasi-periodic orbits and chaotic sets. Additionally, the system’s chaos was managed by the development of control mechanisms, such as the OGY methodology. It has been established that bifurcation and chaos can be stabilized under certain circumstances. A thorough numerical simulation further supports our analytical findings, which include stability regions, bifurcation curves in 2D & 3D, phase plots, and the maximal Lyapunov exponent, etc.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 95-120"},"PeriodicalIF":4.4,"publicationDate":"2025-11-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618554","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-21DOI: 10.1016/j.matcom.2025.11.022
Mohammad Amini , Ramin Vatankhah , Mohammad Mehdi Arefi
Recent advancements in mathematical modeling have enhanced the analysis of cancer responses to treatments, particularly in the promising field of chemovirotherapy. The present research introduces a novel data-driven mathematical model of chemovirotherapy that comprehensively incorporates the immune response. This thorough consideration of the immune system enables a comparative analysis of chemovirotherapy with immunotherapy, specifically CD8+T cells, CD4+T cells, and IL-2 cytokine therapies. By estimating a Michaelis-Menten constant from empirical therapy data for virotherapy, the computational efficiency of the dynamical system is enhanced while maintaining high accuracy in capturing virotherapy dynamics. Parameters are estimated using the Unscented Kalman Filter based on data from human melanoma cell lines. A stability analysis investigates parameter-dependent equilibrium shifts of the model, revealing that treatments such as chemotherapy destabilize the system at any dosage, which may inform treatment scenarios. Numerical simulations conducted on the model demonstrate that the combination of chemotherapy and virotherapy yields superior outcomes, particularly in cases of high tumor burden and weakened immune systems. This study presents a comprehensive framework for comparing immunotherapy, chemotherapy, and virotherapy, thereby advancing cancer therapeutic modeling and facilitating the optimization of comparative treatment strategies.
{"title":"Enhancing chemovirotherapy through a data-driven model with detailed consideration of immune system response using Unscented Kalman Filter","authors":"Mohammad Amini , Ramin Vatankhah , Mohammad Mehdi Arefi","doi":"10.1016/j.matcom.2025.11.022","DOIUrl":"10.1016/j.matcom.2025.11.022","url":null,"abstract":"<div><div>Recent advancements in mathematical modeling have enhanced the analysis of cancer responses to treatments, particularly in the promising field of chemovirotherapy. The present research introduces a novel data-driven mathematical model of chemovirotherapy that comprehensively incorporates the immune response. This thorough consideration of the immune system enables a comparative analysis of chemovirotherapy with immunotherapy, specifically CD8<sup>+</sup>T cells, CD4<sup>+</sup>T cells, and IL-2 cytokine therapies. By estimating a Michaelis-Menten constant from empirical therapy data for virotherapy, the computational efficiency of the dynamical system is enhanced while maintaining high accuracy in capturing virotherapy dynamics. Parameters are estimated using the Unscented Kalman Filter based on data from human melanoma cell lines. A stability analysis investigates parameter-dependent equilibrium shifts of the model, revealing that treatments such as chemotherapy destabilize the system at any dosage, which may inform treatment scenarios. Numerical simulations conducted on the model demonstrate that the combination of chemotherapy and virotherapy yields superior outcomes, particularly in cases of high tumor burden and weakened immune systems. This study presents a comprehensive framework for comparing immunotherapy, chemotherapy, and virotherapy, thereby advancing cancer therapeutic modeling and facilitating the optimization of comparative treatment strategies.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 283-306"},"PeriodicalIF":4.4,"publicationDate":"2025-11-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145685416","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-19DOI: 10.1016/j.matcom.2025.11.021
Huiyan Zhang , Yu Huang , Ning Zhao , Kalidass Mathiyalagan , Peng Shi
This paper investigates the issue of adaptive neural non-fragile proportional and derivative (PD) feedback control for the singular systems with unknown nonlinear dynamics. First, considering the inaccuracy of controller implementation, the problem of non-fragile controller design is considered and solved by using a robust control strategy. Second, PD feedback control is established to transform the singular system into a normal system, which facilitates stability analysis of the system. Third, the adaptive proportional–derivative radial basis function neural network technique is used to approximate the unknown nonlinear function and resist its influence. Under this designed framework, the stability conditions of the closed-loop system are given by using the Lyapunov method. The designed methods of state feedback gains and observer-based gain matrices are presented, respectively. Last, three examples are employed to elucidate the feasibility of the developed control strategy.
{"title":"RBFNN-based adaptive control of singular systems via non-fragile proportional and derivative feedback method","authors":"Huiyan Zhang , Yu Huang , Ning Zhao , Kalidass Mathiyalagan , Peng Shi","doi":"10.1016/j.matcom.2025.11.021","DOIUrl":"10.1016/j.matcom.2025.11.021","url":null,"abstract":"<div><div>This paper investigates the issue of adaptive neural non-fragile proportional and derivative (PD) feedback control for the singular systems with unknown nonlinear dynamics. First, considering the inaccuracy of controller implementation, the problem of non-fragile controller design is considered and solved by using a robust control strategy. Second, PD feedback control is established to transform the singular system into a normal system, which facilitates stability analysis of the system. Third, the adaptive proportional–derivative radial basis function neural network technique is used to approximate the unknown nonlinear function and resist its influence. Under this designed framework, the stability conditions of the closed-loop system are given by using the Lyapunov method. The designed methods of state feedback gains and observer-based gain matrices are presented, respectively. Last, three examples are employed to elucidate the feasibility of the developed control strategy.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 51-68"},"PeriodicalIF":4.4,"publicationDate":"2025-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145571202","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-19DOI: 10.1016/j.matcom.2025.09.025
Yovan Singh , Bapan Ghosh , Suman Mondal
Time delays are integral to ecological processes. Population models incorporate time delays to account for the time required for maturation, gestation, dispersal, and many more. Time delay can induce various stability dynamics, including (i) stability invariance, (i) stability change, (iii) stability switching, (iv) instability invariance, and (v) instability switching. Even one of these dynamics can occur with multiple mechanisms based on the distribution of critical time delays. Generally, two or three types of dynamics are detected in many population models, but exhibiting all the above dynamics is not observed. In an ecological system, species form groups to improve their chances of survival. Taking inspiration from tuna’s forging behavior Cosner et al. (1999) developed the Cosner functional response. In this study, we propose a delayed predator–prey model with Cosner functional response. The non-delayed model can have up to four equilibria, two coexisting equilibria (anti-saddle and saddle), along with trivial and boundary equilibria. The stability of all equilibria is analyzed with time delay. Under certain parameter conditions, the boundary equilibrium remains globally stable for all delays. For increasing delay, the anti-saddle equilibrium may: (i) remain stable, (ii) undergo stability change (two possible scenarios), (iii) undergo stability switching, (iv) remain unstable (two possible scenarios), or (v) undergo instability switching. These seven stability scenarios are verified to exhibit, while an additional instability invariance scenario, where no critical delay exists, is analytically shown to be non-existent. Showing all these mentioned stability scenarios in a predator–prey model with a single delay is a novelty of this paper. If the anti-saddle equilibrium is stable in the absence of delay, then the degenerate case may occur, which implies the local stability between any two consecutive delay thresholds. Moreover, we have analytically proved that the degenerate case is not possible if the anti-saddle equilibrium is unstable in the absence of delay, which is a new observation in population dynamics. We have computed species survival basin for increasing delay. Our investigation reveals that increasing delay can change the shape and size of the basin, making delay beneficial or harmful for the species’ survival, depending on the initial populations of species. Finally, we have proposed an open question and outlined a couple of potential directions for future research.
时间延迟是生态过程不可或缺的一部分。种群模型包含了时间延迟,以解释成熟、孕育、扩散等所需的时间。时间延迟可以诱发各种稳定性动力学,包括(i)稳定性不变性,(i)稳定性变化,(iii)稳定性切换,(iv)不稳定性不变性和(v)不稳定性切换。甚至这些动态中的一种也可能发生在基于临界时间延迟分布的多种机制中。通常,在许多种群模型中检测到两种或三种类型的动态,但没有观察到表现出上述所有动态。在生态系统中,物种形成群体是为了提高生存的机会。Cosner et al.(1999)从金枪鱼的锻造行为中获得灵感,开发了Cosner功能反应。在本研究中,我们提出了一个具有Cosner功能响应的延迟捕食者-猎物模型。非延迟模型最多可以有四个平衡点,两个共存平衡点(反鞍态和鞍态),以及平凡平衡点和边界平衡点。用时滞分析了所有平衡点的稳定性。在一定的参数条件下,边界平衡对所有时滞保持全局稳定。对于增加的延迟,反鞍平衡可能:(i)保持稳定,(ii)经历稳定性变化(两种可能的情况),(iii)经历稳定性切换,(iv)保持不稳定(两种可能的情况),或(v)经历不稳定切换。这七个稳定性场景经过验证,而另一个不稳定不变性场景(不存在临界延迟)分析显示不存在。在具有单一延迟的捕食者-猎物模型中显示所有上述稳定性情景是本文的一个新颖之处。如果在没有延迟的情况下,反鞍平衡是稳定的,则可能出现退化情况,这意味着任意两个连续延迟阈值之间的局部稳定性。此外,我们还解析地证明了在没有时滞的情况下,如果反鞍平衡是不稳定的,就不可能出现退化情况,这是种群动力学中的一个新的观察结果。我们计算了增加延迟的物种生存盆地。我们的研究表明,延迟的增加可以改变盆地的形状和大小,使延迟对物种的生存有利或有害,这取决于物种的初始种群。最后,我们提出了一个开放性问题,并概述了未来研究的几个潜在方向。
{"title":"Delay-induced multiple stability scenarios, species coexistence, and predator extinction in an ecological system","authors":"Yovan Singh , Bapan Ghosh , Suman Mondal","doi":"10.1016/j.matcom.2025.09.025","DOIUrl":"10.1016/j.matcom.2025.09.025","url":null,"abstract":"<div><div>Time delays are integral to ecological processes. Population models incorporate time delays to account for the time required for maturation, gestation, dispersal, and many more. Time delay can induce various stability dynamics, including (i) stability invariance, (i) stability change, (iii) stability switching, (iv) instability invariance, and (v) instability switching. Even one of these dynamics can occur with multiple mechanisms based on the distribution of critical time delays. Generally, two or three types of dynamics are detected in many population models, but exhibiting all the above dynamics is not observed. In an ecological system, species form groups to improve their chances of survival. Taking inspiration from tuna’s forging behavior Cosner et al. (1999) developed the Cosner functional response. In this study, we propose a delayed predator–prey model with Cosner functional response. The non-delayed model can have up to four equilibria, two coexisting equilibria (anti-saddle and saddle), along with trivial and boundary equilibria. The stability of all equilibria is analyzed with time delay. Under certain parameter conditions, the boundary equilibrium remains globally stable for all delays. For increasing delay, the anti-saddle equilibrium may: (i) remain stable, (ii) undergo stability change (two possible scenarios), (iii) undergo stability switching, (iv) remain unstable (two possible scenarios), or (v) undergo instability switching. These seven stability scenarios are verified to exhibit, while an additional instability invariance scenario, where no critical delay exists, is analytically shown to be non-existent. Showing all these mentioned stability scenarios in a predator–prey model with a single delay is a novelty of this paper. If the anti-saddle equilibrium is stable in the absence of delay, then the degenerate case may occur, which implies the local stability between any two consecutive delay thresholds. Moreover, we have analytically proved that the degenerate case is not possible if the anti-saddle equilibrium is unstable in the absence of delay, which is a new observation in population dynamics. We have computed species survival basin for increasing delay. Our investigation reveals that increasing delay can change the shape and size of the basin, making delay beneficial or harmful for the species’ survival, depending on the initial populations of species. Finally, we have proposed an open question and outlined a couple of potential directions for future research.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 171-195"},"PeriodicalIF":4.4,"publicationDate":"2025-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618557","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This study investigates a fractional-order three-disk dynamo system incorporating time delay and viscous friction, enhancing its relevance to real-world phenomena. We analyze dynamics of the system with and without time delay, revealing richer behaviors in the delayed case. Through theoretical analysis, we investigate equilibrium points and their stability, identifying pitchfork and double-Hopf bifurcations that lead to complex dynamics, including three-dimensional torus structures. Numerical simulations validate these findings for both fractional and classical systems, highlighting the impact of fractional-order derivatives and time delays. A comparative analysis shows that the fractional-order system exhibits a broader stability region than its integer-order counterpart, underscoring the stabilizing role of fractional calculus. These results provide insights into modeling magnetic field dynamics in geophysical and astrophysical systems, with potential applications to geomagnetic reversals and stellar magnetic cycles.
{"title":"Stability and bifurcation of a time-delayed fractional three-disk system","authors":"Elham Ghafari , Reza Khoshsiar Ghaziani , Javad Alidousti , Khayyam Salehi","doi":"10.1016/j.matcom.2025.11.023","DOIUrl":"10.1016/j.matcom.2025.11.023","url":null,"abstract":"<div><div>This study investigates a fractional-order three-disk dynamo system incorporating time delay and viscous friction, enhancing its relevance to real-world phenomena. We analyze dynamics of the system with and without time delay, revealing richer behaviors in the delayed case. Through theoretical analysis, we investigate equilibrium points and their stability, identifying pitchfork and double-Hopf bifurcations that lead to complex dynamics, including three-dimensional torus structures. Numerical simulations validate these findings for both fractional and classical systems, highlighting the impact of fractional-order derivatives and time delays. A comparative analysis shows that the fractional-order system exhibits a broader stability region than its integer-order counterpart, underscoring the stabilizing role of fractional calculus. These results provide insights into modeling magnetic field dynamics in geophysical and astrophysical systems, with potential applications to geomagnetic reversals and stellar magnetic cycles.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 16-34"},"PeriodicalIF":4.4,"publicationDate":"2025-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145555217","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-19DOI: 10.1016/j.matcom.2025.11.019
Miao Chen , Shijie Zhao , Tianran Zhang , Xin Yu
Constrained multi-objective optimization problems (CMOPs) constitute a prevalent and ubiquitous class of optimization challenges that are frequently encountered across diverse field within science and engineering. To solve the complementary multi-objective optimization problem with narrow and disconnected feasible regions, dual-population two-archive evolutionary framework for constrained multi-objective optimization with constrained-archive solution phase-transition and auxiliary-population environment selection pause-termination (CAE_2SP) is proposed. The algorithm uses dual-population with different efficacy and two archives with different functions. To improve the problem of lower population diversity, constrained-archive solution phase-transition strategy is proposed. In this strategy, the diversity of solutions is emphasized in the early generation, so the non-dominated infeasible solutions generated by the evolution of main population are stored in the archive. In the late generation, the feasibility of solutions is taken into account, hence, constrained archive is transformed into storing non-dominated feasible solutions. In addition, this paper puts forward auxiliary-population environment selection pause-termination strategy. In this strategy, auxiliary population stop updating in the late generation and uses the optimal population information in the early generation to guide the evolution, to reduce the consumption of computing resources in the late generation and provide more computing resources for main population to help it search for potential feasible regions. The experimental results of nine comparison algorithms in three benchmark function suites demonstrate that CAE_2SP has superior performance in solving CMOPs compared with others. To validate the applicability of the proposed algorithm in solving practical problems, six real-world problems are employed for testing. The experimental results demonstrate that CAE_2SP exhibits competitive performance in addressing practical issues.
{"title":"Dual-population two-archive evolutionary framework for constrained multi-objective optimization","authors":"Miao Chen , Shijie Zhao , Tianran Zhang , Xin Yu","doi":"10.1016/j.matcom.2025.11.019","DOIUrl":"10.1016/j.matcom.2025.11.019","url":null,"abstract":"<div><div>Constrained multi-objective optimization problems (CMOPs) constitute a prevalent and ubiquitous class of optimization challenges that are frequently encountered across diverse field within science and engineering. To solve the complementary multi-objective optimization problem with narrow and disconnected feasible regions, dual-population two-archive evolutionary framework for constrained multi-objective optimization with constrained-archive solution phase-transition and auxiliary-population environment selection pause-termination (CAE_2SP) is proposed. The algorithm uses dual-population with different efficacy and two archives with different functions. To improve the problem of lower population diversity, constrained-archive solution phase-transition strategy is proposed. In this strategy, the diversity of solutions is emphasized in the early generation, so the non-dominated infeasible solutions generated by the evolution of main population are stored in the archive. In the late generation, the feasibility of solutions is taken into account, hence, constrained archive is transformed into storing non-dominated feasible solutions. In addition, this paper puts forward auxiliary-population environment selection pause-termination strategy. In this strategy, auxiliary population stop updating in the late generation and uses the optimal population information in the early generation to guide the evolution, to reduce the consumption of computing resources in the late generation and provide more computing resources for main population to help it search for potential feasible regions. The experimental results of nine comparison algorithms in three benchmark function suites demonstrate that CAE_2SP has superior performance in solving CMOPs compared with others. To validate the applicability of the proposed algorithm in solving practical problems, six real-world problems are employed for testing. The experimental results demonstrate that CAE_2SP exhibits competitive performance in addressing practical issues.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"243 ","pages":"Pages 196-220"},"PeriodicalIF":4.4,"publicationDate":"2025-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618556","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-19DOI: 10.1016/j.matcom.2025.11.006
Zhidong Guo
Paper (Shokrollahi et al. [Mathematics and Computers in Simulation 226 (2024) 172-183]), addresses the pricing of geometric Asian options with jumps using a method analogous to that employed in non-jump scenarios. In this comment, based on the probability distribution of path variables, we will point out that the main conclusion of the paper is incorrect. More importantly, we aim to show that the conclusions drawn from non-jump models do not directly carry over to jump models when dealing with path-dependent options.
论文(Shokrollahi et al. [Mathematics and Computers in Simulation 226(2024) 172-183]),使用类似于非跳跃情景的方法,解决了具有跳跃的几何亚洲期权的定价问题。在这篇评论中,我们将根据路径变量的概率分布,指出本文的主要结论是不正确的。更重要的是,我们的目标是表明,当处理路径依赖选项时,从非跳跃模型得出的结论并不直接适用于跳跃模型。
{"title":"Comments on “Pricing Asian options under the mixed fractional Brownian motion with jumps”.","authors":"Zhidong Guo","doi":"10.1016/j.matcom.2025.11.006","DOIUrl":"10.1016/j.matcom.2025.11.006","url":null,"abstract":"<div><div>Paper (Shokrollahi et al. [Mathematics and Computers in Simulation 226 (2024) 172-183]), addresses the pricing of geometric Asian options with jumps using a method analogous to that employed in non-jump scenarios. In this comment, based on the probability distribution of path variables, we will point out that the main conclusion of the paper is incorrect. More importantly, we aim to show that the conclusions drawn from non-jump models do not directly carry over to jump models when dealing with path-dependent options.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"242 ","pages":"Pages 121-124"},"PeriodicalIF":4.4,"publicationDate":"2025-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145580213","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}