Showing 1–37 of 37 results for author: Silva, C J

Stability Analysis of Delayed COVID-19 Models

Authors: Mohamed A. Zaitri, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We analyze mathematical models for COVID-19 with discrete time delays and vaccination. Sufficient conditions for the local stability of the endemic and disease-free equilibrium points are proved for any positive time delay. The stability results are illustrated through numerical simulations performed in MATLAB. We analyze mathematical models for COVID-19 with discrete time delays and vaccination. Sufficient conditions for the local stability of the endemic and disease-free equilibrium points are proved for any positive time delay. The stability results are illustrated through numerical simulations performed in MATLAB. △ Less

Submitted 13 August, 2022; originally announced August 2022.

Comments: This is a preprint of a paper whose final and definite form is published Open Access in 'Axioms' at [https://doi.org/10.3390/axioms11080400]

MSC Class: 34C60; 92D30

Journal ref: Axioms 11 (2022), no. 8, Art. 400, 21 pp

A SIQRB delayed model for cholera and optimal control treatment

Authors: Ana P. Lemos-Paiao, Helmut Maurer, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We improve a recent mathematical model for cholera by adding a time delay that represents the time between the instant at which an individual becomes infected and the instant at which he begins to have symptoms of cholera disease. We prove that the delayed cholera model is biologically meaningful and analyze the local asymptotic stability of the equilibrium points for positive time delays. An opti… ▽ More We improve a recent mathematical model for cholera by adding a time delay that represents the time between the instant at which an individual becomes infected and the instant at which he begins to have symptoms of cholera disease. We prove that the delayed cholera model is biologically meaningful and analyze the local asymptotic stability of the equilibrium points for positive time delays. An optimal control problem is proposed and analyzed, where the goal is to obtain optimal treatment strategies, through quarantine, that minimize the number of infective individuals and the bacterial concentration, as well as treatment costs. Necessary optimality conditions are applied to the delayed optimal control problem, with a $L^1$ type cost functional. We show that the delayed cholera model fits better the cholera outbreak that occurred in the Department of Artibonite -- Haiti, from 1 November 2010 to 1 May 2011, than the non-delayed model. Considering the data of the cholera outbreak in Haiti, we solve numerically the delayed optimal control problem and propose solutions for the outbreak control and eradication. △ Less

Submitted 25 June, 2022; originally announced June 2022.

Comments: This is a preprint of a paper whose final and definite form is published Open Access in 'Math. Model. Nat. Phenom.' at [https://www.mmnp-journal.org]. Submitted to MMNP: Sept 23, 2021; Revised: March 15, 2022; Accepted: June 25, 2022

MSC Class: 34C60; 49K15; 92D30

Journal ref: Math. Model. Nat. Phenom. 17 (2022), Art. 25, 23 pp

Model-free based control of a HIV/AIDS prevention model

Authors: Loïc Michel, Cristiana J. Silva, Delfim F. M. Torres

Abstract: Controlling an epidemiological model is often performed using optimal control theory techniques for which the solution depends on the equations of the control system, objective functional and possible state and/or control constraints. In this paper, we propose a model-free control approach based on an algorithm that operates in 'real-time' and drives the state solution according to a direct feedba… ▽ More Controlling an epidemiological model is often performed using optimal control theory techniques for which the solution depends on the equations of the control system, objective functional and possible state and/or control constraints. In this paper, we propose a model-free control approach based on an algorithm that operates in 'real-time' and drives the state solution according to a direct feedback on the state solution that is aimed to be minimized, and without knowing explicitly the equations of the control system. We consider a concrete epidemic problem of minimizing the number of HIV infected individuals, through the preventive measure pre-exposure prophylaxis (PrEP) given to susceptible individuals. The solutions must satisfy control and mixed state-control constraints that represent the limitations on PrEP implementation. Our model-free based control algorithm allows to close the loop between the number of infected individuals with HIV and the supply of PrEP medication 'in real time', in such manner that the number of infected individuals is asymptotically reduced and the number of individuals under PrEP medication is below a fixed constant value. We prove the efficiency of our approach and compare the model-free control solutions with the ones obtained using a classical optimal control approach via Pontryagin maximum principle. The performed numerical simulations allow us to conclude that the model-free based control strategy highlights new and interesting performances compared with the classical optimal control approach. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: This is a preprint of a paper whose final and definite form is published open access by 'Mathematical Biosciences and Engineering' (ISSN: 1551-0018) at [https://www.aimspress.com/mbe/article/5926/special-articles]

Journal ref: Math. Biosci. Eng. 19 (2022), no. 1, 759--774

Fractional model of COVID-19 applied to Galicia, Spain and Portugal

Authors: Faical Ndairou, Ivan Area, Juan J. Nieto, Cristiana J. Silva, Delfim F. M. Torres

Abstract: A fractional compartmental mathematical model for the spread of the COVID-19 disease is proposed. Special focus has been done on the transmissibility of super-spreaders individuals. Numerical simulations are shown for data of Galicia, Spain, and Portugal. For each region, the order of the Caputo derivative takes a different value, that is not close to one, showing the relevance of considering frac… ▽ More A fractional compartmental mathematical model for the spread of the COVID-19 disease is proposed. Special focus has been done on the transmissibility of super-spreaders individuals. Numerical simulations are shown for data of Galicia, Spain, and Portugal. For each region, the order of the Caputo derivative takes a different value, that is not close to one, showing the relevance of considering fractional models. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: This is a preprint of a paper whose final and definite form is published by 'Chaos Solitons Fractals' (ISSN: 0960-0779). Paper Submitted 04/June/2020; Revised 22/July/2020 and 02/Aug/2020; Accepted 04/Jan/2021

MSC Class: 26A33; 34A08; 92D30

Journal ref: Chaos Solitons Fractals 144 (2021), Art. 110652, 7 pp

A New Compartmental Epidemiological Model for COVID-19 with a Case Study of Portugal

Authors: Ana P. Lemos-Paiao, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a compartmental mathematical model for the spread of the COVID-19 disease, showing its usefulness with respect to the pandemic in Portugal, from the first recorded case in the country till the end of the three states of emergency. New results include the compartmental model, described by a system of seven ordinary differential equations; proof of positivity and boundedness of solutions;… ▽ More We propose a compartmental mathematical model for the spread of the COVID-19 disease, showing its usefulness with respect to the pandemic in Portugal, from the first recorded case in the country till the end of the three states of emergency. New results include the compartmental model, described by a system of seven ordinary differential equations; proof of positivity and boundedness of solutions; investigation of equilibrium points and their stability analysis; computation of the basic reproduction number; and numerical simulations with official real data from the Portuguese health authorities. Besides completely new, the proposed model allows to describe quite well the spread of COVID-19 in Portugal, fitting simultaneously not only the number of active infected individuals but also the number of hospitalized individuals, respectively with a $L^2$ error of $9.2152e-04$ and $1.6136e-04$ with respect to the initial population. Such results are very important, from a practical point of view, and far from trivial from a mathematical perspective. Moreover, the obtained value for the basic reproduction number is in agreement with the one given by the Portuguese authorities at the end of the emergency states. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: This is a preprint of a paper whose final and definite form is published by 'Ecological Complexity' (ISSN: 1476-945X). Paper Submitted 05/June/2020; Revised 24/July/2020, 08/Sept/2020 and 02/Nov/2020; Accepted 17/Nov/2020

MSC Class: 34C60; 92D30

Journal ref: Ecological Complexity 44 (2020) Art. 100885, 8 pp

Optimal control of vaccination and plasma transfusion with potential usefulness for COVID-19

Authors: Juliana Couras, Ivan Area, Juan J. Nieto, Cristiana J. Silva, Delfim F. M. Torres

Abstract: The SEIR model is a compartmental model used to simulate the dynamics of an epidemic. In this chapter, we introduce two control functions in the compartmental SEIR model representing vaccination and plasma transfusion. Optimal control problems are proposed to study the effects of these two control measures, on the reduction of infected individuals and increase of recovered ones, with minimal costs… ▽ More The SEIR model is a compartmental model used to simulate the dynamics of an epidemic. In this chapter, we introduce two control functions in the compartmental SEIR model representing vaccination and plasma transfusion. Optimal control problems are proposed to study the effects of these two control measures, on the reduction of infected individuals and increase of recovered ones, with minimal costs. Up to our knowledge, the plasma transfusion treatment has never been considered as a control strategy for epidemics mitigation. The proposed vaccination and treatment strategies may have a real application in the challenging and hard problem of controlling the COVID-19 pandemic. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Comments: This is a preprint of a paper whose final and definite form is published in 'Analysis of Infectious Disease Problems (Covid-19) and Their Global Impact', Springer Nature Singapore Pte Ltd. Submitted 15/July/2020; revised 02/Oct/2020; accepted 08/Oct/2020

Complex network model for COVID-19: human behavior, pseudo-periodic solutions and multiple epidemic waves

Authors: Cristiana J. Silva, Guillaume Cantin, Carla Cruz, Rui Fonseca-Pinto, Rui Passadouro da Fonseca, Estevao Soares dos Santos, Delfim F. M. Torres

Abstract: We propose a mathematical model for the transmission dynamics of SARS-CoV-2 in a homogeneously mixing non constant population, and generalize it to a model where the parameters are given by piecewise constant functions. This allows us to model the human behavior and the impact of public health policies on the dynamics of the curve of active infected individuals during a COVID-19 epidemic outbreak.… ▽ More We propose a mathematical model for the transmission dynamics of SARS-CoV-2 in a homogeneously mixing non constant population, and generalize it to a model where the parameters are given by piecewise constant functions. This allows us to model the human behavior and the impact of public health policies on the dynamics of the curve of active infected individuals during a COVID-19 epidemic outbreak. After proving the existence and global asymptotic stability of the disease-free and endemic equilibrium points of the model with constant parameters, we consider a family of Cauchy problems, with piecewise constant parameters, and prove the existence of pseudo-oscillations between a neighborhood of the disease-free equilibrium and a neighborhood of the endemic equilibrium, in a biologically feasible region. In the context of the COVID-19 pandemic, this pseudo-periodic solutions are related to the emergence of epidemic waves. Then, to capture the impact of mobility in the dynamics of COVID-19 epidemics, we propose a complex network with six distinct regions based on COVID-19 real data from Portugal. We perform numerical simulations for the complex network model, where the objective is to determine a topology that minimizes the level of active infected individuals and the existence of topologies that are likely to worsen the level of infection. We claim that this methodology is a tool with enormous potential in the current pandemic context, and can be applied in the management of outbreaks (in regional terms) but also to manage the opening/closing of borders. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 23 pages, 10 figures, submitted 5-Oct-2020

MSC Class: 92D30 (Primary) 05C82; 34C60; 34D20 (Secondary)

Journal ref: J. Math. Anal. Appl. 514 (2022), no. 2, Art. 125171, 25pp

Optimal control of the COVID-19 pandemic: controlled sanitary deconfinement in Portugal

Authors: Cristiana J. Silva, Carla Cruz, Delfim F. M. Torres, Alberto P. Munuzuri, Alejandro Carballosa, Ivan Area, Juan J. Nieto, Rui Fonseca-Pinto, Rui Passadouro da Fonseca, Estevao Soares dos Santos, Wilson Abreu, Jorge Mira

Abstract: The COVID-19 pandemic has forced policy makers to decree urgent confinements to stop a rapid and massive contagion. However, after that stage, societies are being forced to find an equilibrium between the need to reduce contagion rates and the need to reopen their economies. The experience hitherto lived has provided data on the evolution of the pandemic, in particular the population dynamics as a… ▽ More The COVID-19 pandemic has forced policy makers to decree urgent confinements to stop a rapid and massive contagion. However, after that stage, societies are being forced to find an equilibrium between the need to reduce contagion rates and the need to reopen their economies. The experience hitherto lived has provided data on the evolution of the pandemic, in particular the population dynamics as a result of the public health measures enacted. This allows the formulation of forecasting mathematical models to anticipate the consequences of political decisions. Here we propose a model to do so and apply it to the case of Portugal. With a mathematical deterministic model, described by a system of ordinary differential equations, we fit the real evolution of COVID-19 in this country. After identification of the population readiness to follow social restrictions, by analyzing the social media, we incorporate this effect in a version of the model that allow us to check different scenarios. This is realized by considering a Monte Carlo discrete version of the previous model coupled via a complex network. Then, we apply optimal control theory to maximize the number of people returning to "normal life" and minimizing the number of active infected individuals with minimal economical costs while warranting a low level of hospitalizations. This work allows testing various scenarios of pandemic management (closure of sectors of the economy, partial/total compliance with protection measures by citizens, number of beds in intensive care units, etc.), ensuring the responsiveness of the health system, thus being a public health decision support tool. △ Less

Submitted 2 February, 2021; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: Paper whose final and definite form is published Open Access by 'Scientific Reports', ISSN 2045-2322, Springer Nature. Submitted 01/Sept/2020; Revised 16/Dec/2020; Accepted 27/Jan/2021; Published by [www.nature.com/scientificreports], see [https://doi.org/10.1038/s41598-021-83075-6]

Journal ref: Scientific Reports 11 (2021), Art. 3451, 15 pp

On SICA models for HIV transmission

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We revisit the SICA (Susceptible-Infectious-Chronic-AIDS) mathematical model for transmission dynamics of the human immunodeficiency virus (HIV) with varying population size in a homogeneously mixing population. We consider SICA models given by systems of ordinary differential equations and some generalizations given by systems with fractional and stochastic differential operators. Local and globa… ▽ More We revisit the SICA (Susceptible-Infectious-Chronic-AIDS) mathematical model for transmission dynamics of the human immunodeficiency virus (HIV) with varying population size in a homogeneously mixing population. We consider SICA models given by systems of ordinary differential equations and some generalizations given by systems with fractional and stochastic differential operators. Local and global stability results are proved for deterministic, fractional, and stochastic-type SICA models. Two case studies, in Cape Verde and Morocco, are investigated. △ Less

Submitted 24 April, 2020; originally announced April 2020.

Comments: This is a preprint of the following paper: C. J. Silva and D. F. M. Torres, On SICA models for HIV transmission, published in 'Mathematical Modelling and Analysis of Infectious Diseases', edited by K. Hattaf and H. Dutta, Springer Nature Switzerland AG. Submitted 28/Dec/2019; revised 21/Apr/2020; accepted 24/Apr/2020. arXiv admin note: text overlap with arXiv:1703.06446, arXiv:1612.00732, arXiv:1812.06965

MSC Class: 92D30

Journal ref: Studies in Systems, Decision and Control 302 (2020), 155--179

Optimal Control of Aquatic Diseases: A Case Study of Yemen's Cholera Outbreak

Authors: Ana P. Lemos-Paiao, Cristiana J. Silva, Delfim F. M. Torres, Ezio Venturino

Abstract: We propose a mathematical model for the transmission dynamics of some strains of the bacterium Vibrio cholerae, responsible for the cholera disease in humans. We prove that, when the basic reproduction number is equal to one, a transcritical bifurcation occurs for which the endemic equilibrium emanates from the disease-free point. A control function is introduced into the model, representing the d… ▽ More We propose a mathematical model for the transmission dynamics of some strains of the bacterium Vibrio cholerae, responsible for the cholera disease in humans. We prove that, when the basic reproduction number is equal to one, a transcritical bifurcation occurs for which the endemic equilibrium emanates from the disease-free point. A control function is introduced into the model, representing the distribution of chlorine water tablets for water purification. An optimal control problem is then proposed and analyzed, where the goal is to determine the fraction of susceptible individuals who should have access to chlorine water tablets in order to minimize the total number of new infections plus the total cost associated with the distribution of chlorine water tablets, over the considered period of time. Finally, we consider real data of the cholera outbreak in Yemen, from 27 April 2017 to 15 April 2018, choosing the values of the parameters of the uncontrolled model that fit the real data. Using our optimal control results, we show, numerically, that the distribution of chlorine water tablets could have stopped, in a fast way, the worst cholera outbreak that ever occurred in human history. Due to the critical situation of Yemen, we also simulate the case where only a small percentage of susceptible individuals has access to chlorine water tablets and obtain an optimal control solution that decreases, substantially, the maximum number of infective individuals attained at the outbreak. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: This is a preprint of a paper whose final and definite form is with 'J. Optim. Theory Appl.' (JOTA). Submitted 6-Apr-2019; Revised 12-Oct-2019 and 21-Jan-2020; Accepted 15-Apr-2020

MSC Class: 34C60; 49K15; 92D30

Journal ref: J. Optim. Theory Appl. 185 (2020), no. 3, 1008--1030

Numerical Optimal Control of HIV Transmission in Octave/MATLAB

Authors: Carlos Campos, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We provide easy and readable GNU Octave/MATLAB code for the simulation of mathematical models described by ordinary differential equations and for the solution of optimal control problems through Pontryagin's maximum principle. For that, we consider a normalized HIV/AIDS transmission dynamics model based on the one proposed in our recent contribution (Silva, C.J.; Torres, D.F.M. A SICA compartment… ▽ More We provide easy and readable GNU Octave/MATLAB code for the simulation of mathematical models described by ordinary differential equations and for the solution of optimal control problems through Pontryagin's maximum principle. For that, we consider a normalized HIV/AIDS transmission dynamics model based on the one proposed in our recent contribution (Silva, C.J.; Torres, D.F.M. A SICA compartmental model in epidemiology with application to HIV/AIDS in Cape Verde. Ecol. Complex. 2017, 30, 70--75), given by a system of four ordinary differential equations. An HIV initial value problem is solved numerically using the ode45 GNU Octave function and three standard methods implemented by us in Octave/MATLAB: Euler method and second-order and fourth-order Runge-Kutta methods. Afterwards, a control function is introduced into the normalized HIV model and an optimal control problem is formulated, where the goal is to find the optimal HIV prevention strategy that maximizes the fraction of uninfected HIV individuals with the least HIV new infections and cost associated with the control measures. The optimal control problem is characterized analytically using the Pontryagin Maximum Principle, and the extremals are computed numerically by implementing a forward-backward fourth-order Runge-Kutta method. Complete algorithms, for both uncontrolled initial value and optimal control problems, developed under the free GNU Octave software and compatible with MATLAB are provided along the article. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: Submitted 1-Oct-2019; Revised 7-Dec-2019; Accepted 17-Dec-2019; Published 19-Dec-2019 in Math. Comput. Appl., EISSN 2297-8747 [https://doi.org/10.3390/mca25010001]

MSC Class: 34K28; 49N90; 92D30

Journal ref: Math. Comput. Appl. 25 (2020), no. 1, 20 pp

Influence of the topology on the dynamics of a complex network of HIV/AIDS epidemic models

Authors: Guillaume Cantin, Cristiana J. Silva

Abstract: In this paper, we propose an original complex network model for an epidemic problem in an heterogeneous geographical area. The complex network is constructed by coupling nonidentical instances of a HIV/AIDS epidemiological model for which a disease-free equilibrium and an endemic equilibrium can coexist. After proving the existence of a positively invariant region for the solutions of the complex… ▽ More In this paper, we propose an original complex network model for an epidemic problem in an heterogeneous geographical area. The complex network is constructed by coupling nonidentical instances of a HIV/AIDS epidemiological model for which a disease-free equilibrium and an endemic equilibrium can coexist. After proving the existence of a positively invariant region for the solutions of the complex network problem, we investigate the effect of the coupling on the dynamics of the network, and establish the existence of a unique disease-free equilibrium for the whole network, which is globally asymptotically stable. We prove the existence of an optimal topology that minimizes the level of infected individuals, and apply the theoretical results to the case of the Cape Verde archipelago. △ Less

Submitted 12 November, 2019; originally announced November 2019.

MSC Class: 34A34; 34C60; 92B05

Journal ref: AIMS Mathematics, 2019, 4(4): 1145-1169

Stability of a fractional HIV/AIDS model

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a fractional order model for HIV/AIDS transmission. Local and uniform stability of the fractional order model is studied. The theoretical results are illustrated through numerical simulations. We propose a fractional order model for HIV/AIDS transmission. Local and uniform stability of the fractional order model is studied. The theoretical results are illustrated through numerical simulations. △ Less

Submitted 2 March, 2019; originally announced March 2019.

Comments: This is a preprint of a paper whose final and definite form is with 'Mathematics and Computers in Simulation', ISSN: 0378-4754. Submitted 21/Jan/2018; Revised 1/Feb/2019; Accepted for publication 1/March/2019. arXiv admin note: text overlap with arXiv:1703.06446, arXiv:1801.07059

MSC Class: 34C60; 34D23; 92D30

Journal ref: Math. Comput. Simul. 164 (2019), 180--190

A minimal HIV-AIDS infection model with general incidence rate and application to Morocco data

Authors: El Mehdi Lotfi, Marouane Mahrouf, Mehdi Maziane, Cristiana J. Silva, Delfim F. M. Torres, Noura Yousfi

Abstract: We study the global dynamics of a SICA infection model with general incidence rate. The proposed model is calibrated with cumulative cases of infection by HIV-AIDS in Morocco from 1986 to 2015. We first prove that our model is biologically and mathematically well-posed. Stability analysis of different steady states is performed and threshold parameters are identified where the model exhibits clear… ▽ More We study the global dynamics of a SICA infection model with general incidence rate. The proposed model is calibrated with cumulative cases of infection by HIV-AIDS in Morocco from 1986 to 2015. We first prove that our model is biologically and mathematically well-posed. Stability analysis of different steady states is performed and threshold parameters are identified where the model exhibits clearance of infection or maintenance of a chronic infection. Furthermore, we examine the robustness of the model to some parameter values by examining the sensitivity of the basic reproduction number. Finally, using numerical simulations with real data from Morocco, we show that the model predicts well such reality. △ Less

Submitted 15 December, 2018; originally announced December 2018.

Comments: This is a preprint of a paper whose final and definite form is with 'Statistics Opt. Inform. Comput.', Vol. 7, No 2 (2019). See [http://www.IAPress.org]. Submitted 16/Sept/2018; Revised 10 & 15/Dec/2018; Accepted 15/Dec/2018

MSC Class: 34D23; 92D30

Journal ref: Stat. Optim. Inf. Comput. 7 (2019), no. 3, 588--603

A cholera mathematical model with vaccination and the biggest outbreak of world's history

Authors: Ana P. Lemos-Paiao, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose and analyse a mathematical model for cholera considering vaccination. We show that the model is epidemiologically and mathematically well posed and prove the existence and uniqueness of disease-free and endemic equilibrium points. The basic reproduction number is determined and the local asymptotic stability of equilibria is studied. The biggest cholera outbreak of world's history began… ▽ More We propose and analyse a mathematical model for cholera considering vaccination. We show that the model is epidemiologically and mathematically well posed and prove the existence and uniqueness of disease-free and endemic equilibrium points. The basic reproduction number is determined and the local asymptotic stability of equilibria is studied. The biggest cholera outbreak of world's history began on 27th April 2017, in Yemen. Between 27th April 2017 and 15th April 2018 there were 2275 deaths due to this epidemic. A vaccination campaign began on 6th May 2018 and ended on 15th May 2018. We show that our model is able to describe well this outbreak. Moreover, we prove that the number of infected individuals would have been much lower provided the vaccination campaign had begun earlier. △ Less

Submitted 13 October, 2018; originally announced October 2018.

Comments: This is a preprint of a paper whose final and definite form is with 'AIMS Mathematics', available in open access from [http://www.aimspress.com/journal/Math]. Submitted 7-July-2018; Revised 14-Sept-2018; Accepted 12-Oct-2018. arXiv admin note: substantial text overlap with arXiv:1611.02195

MSC Class: 34C60; 92D30

Journal ref: AIMS Mathematics 3 (2018), no. 4, 448--463

A stochastic SICA epidemic model for HIV transmission

Authors: Jasmina Djordjevic, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a stochastic SICA epidemic model for HIV transmission, described by stochastic ordinary differential equations, and discuss its perturbation by environmental white noise. Existence and uniqueness of the global positive solution to the stochastic HIV system is proven, and conditions under which extinction and persistence in mean hold, are given. The theoretical results are illustrated vi… ▽ More We propose a stochastic SICA epidemic model for HIV transmission, described by stochastic ordinary differential equations, and discuss its perturbation by environmental white noise. Existence and uniqueness of the global positive solution to the stochastic HIV system is proven, and conditions under which extinction and persistence in mean hold, are given. The theoretical results are illustrated via numerical simulations. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: This is a preprint of a paper whose final and definite form is with 'Applied Mathematics Letters', ISSN 0893-9659. Submitted 22/Jan/2018; Revised 03/May/2018; Accepted for publication 03/May/2018

MSC Class: 34F05; 60H10; 92D30

Journal ref: Appl. Math. Lett. 84 (2018), 168--175

Uniform asymptotic stability of a fractional tuberculosis model

Authors: Weronika Wojtak, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a Caputo type fractional-order mathematical model for the transmission dynamics of tuberculosis (TB). Uniform asymptotic stability of the unique endemic equilibrium of the fractional-order TB model is proved, for any $α\in (0, 1)$. Numerical simulations for the stability of the endemic equilibrium are provided. We propose a Caputo type fractional-order mathematical model for the transmission dynamics of tuberculosis (TB). Uniform asymptotic stability of the unique endemic equilibrium of the fractional-order TB model is proved, for any $α\in (0, 1)$. Numerical simulations for the stability of the endemic equilibrium are provided. △ Less

Submitted 22 January, 2018; originally announced January 2018.

Comments: This is a preprint of a paper whose final and definite form is with Math. Model. Nat. Phenom., ISSN: 0973-5348, available at http://dx.doi.org/10.1051/mmnp/2018015 . Paper Submitted 12-08-2017; Revised 19-01-2018; Accepted 22-01-2018

MSC Class: 26A33; 34C60; 34D20; 92C60

Journal ref: Math. Model. Nat. Phenom. 13 (2018), no. 1, Art. 9, 10 pp

Mathematical modeling of Zika disease in pregnant women and newborns with microcephaly in Brazil

Authors: Faical Ndairou, Ivan Area, Juan J. Nieto, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a new mathematical model for the spread of Zika virus. Special attention is paid to the transmission of microcephaly. Numerical simulations show the accuracy of the model with respect to the Zika outbreak occurred in Brazil. We propose a new mathematical model for the spread of Zika virus. Special attention is paid to the transmission of microcephaly. Numerical simulations show the accuracy of the model with respect to the Zika outbreak occurred in Brazil. △ Less

Submitted 14 November, 2017; originally announced November 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Mathematical Methods in the Applied Sciences', ISSN 0170-4214. Submitted Aug 10, 2017; Revised Nov 13, 2017; accepted for publication Nov 14, 2017

MSC Class: 34D20; 92D30

Journal ref: Math. Methods Appl. Sci. 41 (2018), no. 18, 8929--8941

Optimal control of a delayed HIV model

Authors: Filipe Rodrigues, Cristiana J. Silva, Delfim F. M. Torres, Helmut Maurer

Abstract: We propose a model for the human immunodeficiency virus type 1 (HIV-1) infection with intracellular delay and prove the local asymptotical stability of the equilibrium points. Then we introduce a control function representing the efficiency of reverse transcriptase inhibitors and consider the pharmacological delay associated to the control. Finally, we propose and analyze an optimal control proble… ▽ More We propose a model for the human immunodeficiency virus type 1 (HIV-1) infection with intracellular delay and prove the local asymptotical stability of the equilibrium points. Then we introduce a control function representing the efficiency of reverse transcriptase inhibitors and consider the pharmacological delay associated to the control. Finally, we propose and analyze an optimal control problem with state and control delays. Through numerical simulations, extremal solutions are proposed for minimization of the virus concentration and treatment costs. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Discrete and Continuous Dynamical Systems - Series B' (DCDS-B), ISSN 1531-3492 (print), ISSN 1553-524X (online), available at [http://aimsciences.org/journals/home.jsp?journalID=2]. Paper Submitted 12-July-2016; Revised 26-Nov-2016 and 01-March-2017; Accepted 21-Aug-2017

MSC Class: 34C60; 49K15; 92D30

Journal ref: Discrete Contin. Dyn. Syst. Ser. B 23 (2018), no. 1, 443--458

Global Stability for a HIV/AIDS Model

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We investigate global stability properties of a HIV/AIDS population model with constant recruitment rate, mass action incidence, and variable population size. Existence and uniqueness results for disease-free and endemic equilibrium points are proved. Global stability of the equilibria is obtained through Lyapunov's direct method and LaSalle's invariance principle. We investigate global stability properties of a HIV/AIDS population model with constant recruitment rate, mass action incidence, and variable population size. Existence and uniqueness results for disease-free and endemic equilibrium points are proved. Global stability of the equilibria is obtained through Lyapunov's direct method and LaSalle's invariance principle. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat.', ISSN: 1303-5991. Submitted 08 Sept 2016; Article revised 15 Apr 2017; Article accepted for publication 19 Apr 2017. arXiv admin note: text overlap with arXiv:1703.06446

MSC Class: 34D23; 34C60; 92D30

Journal ref: Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. 67 (2018), no. 1, 93--101

Optimal Spraying in Biological Control of Pests

Authors: Cristiana J. Silva, Delfim F. M. Torres, Ezio Venturino

Abstract: We use optimal control theory with the purpose of finding the best spraying policy with the aim of at least to minimize and possibly to eradicate the number of parasites, i.e., the prey for the spiders living in an agroecosystems. Two different optimal control problems are posed and solved, and their implications discussed. We use optimal control theory with the purpose of finding the best spraying policy with the aim of at least to minimize and possibly to eradicate the number of parasites, i.e., the prey for the spiders living in an agroecosystems. Two different optimal control problems are posed and solved, and their implications discussed. △ Less

Submitted 31 May, 2017; v1 submitted 13 April, 2017; originally announced April 2017.

Comments: This is a preprint of a paper whose final and definite form is: Math. Model. Nat. Phenom. 12 (2017), no. 3, 51--64

MSC Class: 92D25; 49K05

Journal ref: Math. Model. Nat. Phenom. 12 (2017), no. 3, 51--64

Modeling and optimal control of HIV/AIDS prevention through PrEP

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: Pre-exposure prophylaxis (PrEP) consists in the use of an antiretroviral medication to prevent the acquisition of HIV infection by uninfected individuals and has recently demonstrated to be highly efficacious for HIV prevention. We propose a new epidemiological model for HIV/AIDS transmission including PrEP. Existence, uniqueness and global stability of the disease free and endemic equilibriums ar… ▽ More Pre-exposure prophylaxis (PrEP) consists in the use of an antiretroviral medication to prevent the acquisition of HIV infection by uninfected individuals and has recently demonstrated to be highly efficacious for HIV prevention. We propose a new epidemiological model for HIV/AIDS transmission including PrEP. Existence, uniqueness and global stability of the disease free and endemic equilibriums are proved. The model with no PrEP is calibrated with the cumulative cases of infection by HIV and AIDS reported in Cape Verde from 1987 to 2014, showing that it predicts well such reality. An optimal control problem with a mixed state control constraint is then proposed and analyzed, where the control function represents the PrEP strategy and the mixed constraint models the fact that, due to PrEP costs, epidemic context and program coverage, the number of individuals under PrEP is limited at each instant of time. The objective is to determine the PrEP strategy that satisfies the mixed state control constraint and minimizes the number of individuals with pre-AIDS HIV-infection as well as the costs associated with PrEP. The optimal control problem is studied analytically. Through numerical simulations, we demonstrate that PrEP reduces HIV transmission significantly. △ Less

Submitted 27 September, 2017; v1 submitted 19 March, 2017; originally announced March 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Discrete and Continuous Dynamical Systems -- Series S' (DCDS-S), ISSN 1937-1632 (print), ISSN 1937-1179 (online), available at [https://www.aimsciences.org/journals/home.jsp?journalID=15]. In this version, typos detected while reading the proofs were corrected

MSC Class: 34C60; 92D30; 34D23; 49K15

Journal ref: Discrete Contin. Dyn. Syst. Ser. S 11 (2018), no. 1, 119--141

Multiobjective optimization to a TB-HIV/AIDS coinfection optimal control problem

Authors: Roman Denysiuk, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We consider a recent coinfection model for Tuberculosis (TB), Human Immunodeficiency Virus (HIV) infection and Acquired Immunodeficiency Syndrome (AIDS) proposed in [Discrete Contin. Dyn. Syst. 35 (2015), no. 9, 4639--4663]. We introduce and analyze a multiobjective formulation of an optimal control problem, where the two conflicting objectives are: minimization of the number of HIV infected indiv… ▽ More We consider a recent coinfection model for Tuberculosis (TB), Human Immunodeficiency Virus (HIV) infection and Acquired Immunodeficiency Syndrome (AIDS) proposed in [Discrete Contin. Dyn. Syst. 35 (2015), no. 9, 4639--4663]. We introduce and analyze a multiobjective formulation of an optimal control problem, where the two conflicting objectives are: minimization of the number of HIV infected individuals with AIDS clinical symptoms and coinfected with AIDS and active TB; and costs related to prevention and treatment of HIV and/or TB measures. The proposed approach eliminates some limitations of previous works. The results of the numerical study provide comprehensive insights about the optimal treatment policies and the population dynamics resulting from their implementation. Some nonintuitive conclusions are drawn. Overall, the simulation results demonstrate the usefulness and validity of the proposed approach. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Computational and Applied Mathematics', ISSN 0101-8205 (print), ISSN 1807-0302 (electronic). Submitted 04-Feb-2016; revised 11-June-2016 and 02-Sept-2016; accepted for publication 15-March-2017

MSC Class: 90C29; 92C50

Journal ref: Comput. Appl. Math. 37 (2018), no. 2, 2112--2128

Ebola Model and Optimal Control with Vaccination Constraints

Authors: Ivan Area, Faical Ndairou, Juan J. Nieto, Cristiana J. Silva, Delfim F. M. Torres

Abstract: The Ebola virus disease is a severe viral haemorrhagic fever syndrome caused by Ebola virus. This disease is transmitted by direct contact with the body fluids of an infected person and objects contaminated with virus or infected animals, with a death rate close to 90% in humans. Recently, some mathematical models have been presented to analyse the spread of the 2014 Ebola outbreak in West Africa.… ▽ More The Ebola virus disease is a severe viral haemorrhagic fever syndrome caused by Ebola virus. This disease is transmitted by direct contact with the body fluids of an infected person and objects contaminated with virus or infected animals, with a death rate close to 90% in humans. Recently, some mathematical models have been presented to analyse the spread of the 2014 Ebola outbreak in West Africa. In this paper, we introduce vaccination of the susceptible population with the aim of controlling the spread of the disease and analyse two optimal control problems related with the transmission of Ebola disease with vaccination. Firstly, we consider the case where the total number of available vaccines in a fixed period of time is limited. Secondly, we analyse the situation where there is a limited supply of vaccines at each instant of time for a fixed interval of time. The optimal control problems have been solved analytically. Finally, we have performed a number of numerical simulations in order to compare the models with vaccination and the model without vaccination, which has recently been shown to fit the real data. Three vaccination scenarios have been considered for our numerical simulations, namely: unlimited supply of vaccines; limited total number of vaccines; and limited supply of vaccines at each instant of time. △ Less

Submitted 3 March, 2017; originally announced March 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Journal of Industrial and Management Optimization' (JIMO), ISSN 1547-5816 (print), ISSN 1553-166X (online). Submitted February 2016; revised November 2016; accepted for publication March 2017

MSC Class: 49J15; 92D30; 34C60; 49N90

Journal ref: J. Ind. Manag. Optim. 14 (2018), no. 2, 427--446

The effect of immigrant communities coming from higher incidence tuberculosis regions to a host country

Authors: Eugenio M. Rocha, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We introduce a new tuberculosis (TB) mathematical model, with $25$ state-space variables where $15$ are evolution disease states (EDSs), which generalises previous models and takes into account the (seasonal) flux of populations between a high incidence TB country (A) and a host country (B) with low TB incidence, where (B) is divided into a community (G) with high percentage of people from (A) plu… ▽ More We introduce a new tuberculosis (TB) mathematical model, with $25$ state-space variables where $15$ are evolution disease states (EDSs), which generalises previous models and takes into account the (seasonal) flux of populations between a high incidence TB country (A) and a host country (B) with low TB incidence, where (B) is divided into a community (G) with high percentage of people from (A) plus the rest of the population (C). Contrary to some beliefs, related to the fact that agglomerations of individuals increase proportionally to the disease spread, analysis of the model shows that the existence of semi-closed communities are beneficial for the TB control from a global viewpoint. The model and techniques proposed are applied to a case-study with concrete parameters, which model the situation of Angola (A) and Portugal (B), in order to show its relevance and meaningfulness. Simulations show that variations of the transmission coefficient on the origin country has a big influence on the number of infected (and infectious) individuals on the community and the host country. Moreover, there is an optimal ratio for the distribution of individuals in (C) versus (G), which minimizes the reproduction number $R_0$. Such value does not give the minimal total number of infected individuals in all (B), since such is attained when the community (G) is completely isolated (theoretical scenario). Sensitivity analysis and curve fitting on $R_0$ and on EDSs are pursuit in order to understand the TB effects in the global statistics, by measuring the variability of the relevant parameters. We also show that the TB transmission rate $β$ does not act linearly on $R_0$, as is common in compartment models where system feedback or group interactions do not occur. Further, we find the most important parameters for the increase of each EDS. △ Less

Submitted 31 January, 2017; originally announced January 2017.

Comments: This is a preprint of a paper whose final and definite form is with 'Ricerche di Matematica', ISSN 0035-5038 (Print) 1827-3491 (Online), available at [http://link.springer.com/journal/11587]. Paper Submitted 10-Feb-2016; Revised and Accepted 31-Jan-2017

MSC Class: 92D30

Journal ref: Ric. Mat. 67 (2018), no. 1, 89--112

A SICA compartmental model in epidemiology with application to HIV/AIDS in Cape Verde

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a new mathematical model for the transmission dynamics of the human immunodeficiency virus (HIV). Global stability of the unique endemic equilibrium is proved. Then, based on data provided by the "Progress Report on the AIDS response in Cape Verde 2015", we calibrate our model to the cumulative cases of infection by HIV and AIDS from 1987 to 2014 and we show that our model predicts well… ▽ More We propose a new mathematical model for the transmission dynamics of the human immunodeficiency virus (HIV). Global stability of the unique endemic equilibrium is proved. Then, based on data provided by the "Progress Report on the AIDS response in Cape Verde 2015", we calibrate our model to the cumulative cases of infection by HIV and AIDS from 1987 to 2014 and we show that our model predicts well such reality. Finally, a sensitivity analysis is done for the case study in Cape Verde. We conclude that the goal of the United Nations to end the AIDS epidemic by 2030 is a nontrivial task. △ Less

Submitted 2 December, 2016; originally announced December 2016.

Comments: This is a preprint of a paper whose final and definite form is with 'Ecological Complexity', ISSN 1476-945X, available at [http://dx.doi.org/10.1016/j.ecocom.2016.12.001]. Submitted 29/April/2016; Revised 21/Oct/2016; Accepted 02/Dec/2016

MSC Class: 34C60; 34D23; 92D30

Journal ref: Ecological Complexity 30 (2017), 70--75

An epidemic model for cholera with optimal control treatment

Authors: Ana P. Lemos-Paiao, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a mathematical model for cholera with treatment through quarantine. The model is shown to be both epidemiologically and mathematically well posed. In particular, we prove that all solutions of the model are positive and bounded; and that every solution with initial conditions in a certain meaningful set remains in that set for all time. The existence of unique disease-free and endemic e… ▽ More We propose a mathematical model for cholera with treatment through quarantine. The model is shown to be both epidemiologically and mathematically well posed. In particular, we prove that all solutions of the model are positive and bounded; and that every solution with initial conditions in a certain meaningful set remains in that set for all time. The existence of unique disease-free and endemic equilibrium points is proved and the basic reproduction number is computed. Then, we study the local asymptotic stability of these equilibrium points. An optimal control problem is proposed and analyzed, whose goal is to obtain a successful treatment through quarantine. We provide the optimal quarantine strategy for the minimization of the number of infectious individuals and bacteria concentration, as well as the costs associated with the quarantine. Finally, a numerical simulation of the cholera outbreak in the Department of Artibonite (Haiti), in 2010, is carried out, illustrating the usefulness of the model and its analysis. △ Less

Submitted 1 December, 2016; v1 submitted 7 November, 2016; originally announced November 2016.

Comments: This is a preprint of a paper whose final and definite form is with 'Journal of Computational and Applied Mathematics', ISSN 0377-0427, available at [http://dx.doi.org/10.1016/j.cam.2016.11.002]. Submitted 19-June-2016; Revised 14-Sept-2016; Accepted 04-Nov-2016. This version: some typos detected while reading the proofs were corrected

MSC Class: 34C60; 49K15; 92D30

Journal ref: J. Comput. Appl. Math. 318 (2017) 168--180

Stability and Optimal Control of a Delayed HIV Model

Authors: Diana Rocha, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose and investigate a delayed model that studies the relationship between HIV and the immune system during the natural course of infection and in the context of antiviral treatment regimes. Sufficient criteria for local asymptotic stability of the infected and viral free equilibria are given. An optimal control problem with time delays both in state variables (incubation delay) and control… ▽ More We propose and investigate a delayed model that studies the relationship between HIV and the immune system during the natural course of infection and in the context of antiviral treatment regimes. Sufficient criteria for local asymptotic stability of the infected and viral free equilibria are given. An optimal control problem with time delays both in state variables (incubation delay) and control (pharmacological delay) is then formulated and analyzed, where the objective consists to find the optimal treatment strategy that maximizes the number of uninfected $CD4^{+}$ T cells as well as CTL immune response cells, keeping the drug therapy as low as possible. △ Less

Submitted 24 September, 2016; originally announced September 2016.

Comments: This is a preprint of a paper whose final and definite form is with 'Mathematical Methods in the Applied Sciences', ISSN 0170-4214. Submitted 05-July-2016; accepted 24-Sept-2016

MSC Class: 34C60; 49K15; 92D30

Journal ref: Math. Methods Appl. Sci. 41 (2018), no. 6, 2251--2260

Optimal Control of a Tuberculosis Model with State and Control Delays

Authors: Cristiana J. Silva, Helmut Maurer, Delfim F. M. Torres

Abstract: We introduce delays in a tuberculosis (TB) model, representing the time delay on the diagnosis and commencement of treatment of individuals with active TB infection. The stability of the disease free and endemic equilibriums is investigated for any time delay. Corresponding optimal control problems, with time delays in both state and control variables, are formulated and studied. Although it is we… ▽ More We introduce delays in a tuberculosis (TB) model, representing the time delay on the diagnosis and commencement of treatment of individuals with active TB infection. The stability of the disease free and endemic equilibriums is investigated for any time delay. Corresponding optimal control problems, with time delays in both state and control variables, are formulated and studied. Although it is well-known that there is a delay between two to eight weeks between TB infection and reaction of body's immune system to tuberculin, delays for the active infected to be detected and treated, and delays on the treatment of persistent latent individuals due to clinical and patient reasons, which clearly justifies the introduction of time delays on state and control measures, our work seems to be the first to consider such time-delays for TB and apply time-delay optimal control to carry out the optimality analysis. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: This is a preprint of a paper whose final and definite form will be published in the journal 'Mathematical Biosciences and Engineering', ISSN 1547-1063 (print), ISSN 1551-0018 (online). Submitted 30/Oct/2015; revised 24/May/2016; accepted for publication 26/Jun/2016

MSC Class: 34D30; 92D30; 49M05; 93A30

Journal ref: Math. Biosci. Eng. 14 (2017), no. 1, 321--337

A TB-HIV/AIDS coinfection model and optimal control treatment

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose a population model for TB-HIV/AIDS coinfection transmission dynamics, which considers antiretroviral therapy for HIV infection and treatments for latent and active tuberculosis. The HIV-only and TB-only sub-models are analyzed separately, as well as the TB-HIV/AIDS full model. The respective basic reproduction numbers are computed, equilibria and stability are studied. Optimal control t… ▽ More We propose a population model for TB-HIV/AIDS coinfection transmission dynamics, which considers antiretroviral therapy for HIV infection and treatments for latent and active tuberculosis. The HIV-only and TB-only sub-models are analyzed separately, as well as the TB-HIV/AIDS full model. The respective basic reproduction numbers are computed, equilibria and stability are studied. Optimal control theory is applied to the TB-HIV/AIDS model and optimal treatment strategies for co-infected individuals with HIV and TB are derived. Numerical simulations to the optimal control problem show that non intuitive measures can lead to the reduction of the number of individuals with active TB and AIDS. △ Less

Submitted 14 January, 2015; originally announced January 2015.

Comments: This is a preprint of a paper whose final and definite form will be published in the journal 'Discrete and Continuous Dynamical Systems. Series A', ISSN 1078-0947 (Print), ISSN 1553-5231 (Online). Submitted 03/May/2014; revised 08/Nov/2014 and 23/Nov/2014; accepted 14/Jan/2015. arXiv admin note: text overlap with arXiv:1406.0877

MSC Class: 92D30; 93A30; 34D30; 49J15

Journal ref: Discrete Contin. Dyn. Syst. 35 (2015), no. 9, 4639--4663

Multiobjective approach to optimal control for a tuberculosis model

Authors: Roman Denysiuk, Cristiana J. Silva, Delfim F. M. Torres

Abstract: Mathematical modelling can help to explain the nature and dynamics of infection transmissions, as well as support a policy for implementing those strategies that are most likely to bring public health and economic benefits. The paper addresses the application of optimal control strategies in a tuberculosis model. The model consists of a system of ordinary differential equations, which considers re… ▽ More Mathematical modelling can help to explain the nature and dynamics of infection transmissions, as well as support a policy for implementing those strategies that are most likely to bring public health and economic benefits. The paper addresses the application of optimal control strategies in a tuberculosis model. The model consists of a system of ordinary differential equations, which considers reinfection and post-exposure interventions. We propose a multiobjective optimization approach to find optimal control strategies for the minimization of active infectious and persistent latent individuals, as well as the cost associated to the implementation of the control strategies. Optimal control strategies are investigated for different values of the model parameters. The obtained numerical results cover a whole range of the optimal control strategies, providing valuable information about the tuberculosis dynamics and showing the usefulness of the proposed approach. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: This is a preprint of a paper whose final and definite form will be published in Optimization Methods and Software, ISSN 1055-6788 (Print), 1029-4937 (Online). Submitted 09-May-2014; revised 19-Oct-2014; accepted 30-Nov-2014. arXiv admin note: text overlap with arXiv:1305.2145

MSC Class: 90C29; 49N90

Journal ref: Optim. Methods Softw. 30 (2015), no. 5, 893--910

Cost effectiveness analysis of optimal control measures for tuberculosis

Authors: Paula Rodrigues, Cristiana J. Silva, Delfim F. M. Torres

Abstract: We propose and analyse an optimal control problem where the control system is a mathematical model for tuberculosis that considers reinfection. The control functions represent the fraction of early latent and persistent latent individuals that are treated. Our aim is to study how these control measures should be implemented, for a certain time period, in order to reduce the number of active infect… ▽ More We propose and analyse an optimal control problem where the control system is a mathematical model for tuberculosis that considers reinfection. The control functions represent the fraction of early latent and persistent latent individuals that are treated. Our aim is to study how these control measures should be implemented, for a certain time period, in order to reduce the number of active infected individuals, while minimizing the interventions implementation costs. The optimal intervention is compared along different epidemiological scenarios, by varying the transmission coefficient. The impact of variation of the risk of reinfection, as a result of acquired immunity to a previous infection for treated individuals on the optimal controls and associated solutions, is analysed. A cost-effectiveness analysis is done, to compare the application of each one of the control measures, separately or in combination. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Comments: This is a preprint of a paper whose final and definite form is published in the Bulletin of Mathematical Biology, ISSN 0092-8240 (see http://dx.doi.org/10.1007/s11538-014-0028-6). Submitted 04/July/2014; Revised 07/Sept/2014; Accepted 11/Sept/2014

MSC Class: 92D30; 49M05

Journal ref: Bull. Math. Biol. 76 (2014), no. 10, 2627--2645

Optimal Control of Tuberculosis: A Review

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We review the optimal control of systems modeling the dynamics of tuberculosis. Time dependent control functions are introduced in the mathematical models, representing strategies for the improvement of the treatment and cure of active infectious and/or latent individuals. Optimal control theory allows then to find the optimal way to implement the strategies, minimizing the number of infectious an… ▽ More We review the optimal control of systems modeling the dynamics of tuberculosis. Time dependent control functions are introduced in the mathematical models, representing strategies for the improvement of the treatment and cure of active infectious and/or latent individuals. Optimal control theory allows then to find the optimal way to implement the strategies, minimizing the number of infectious and/or latent individuals and keeping the cost of implementation as low as possible. An optimal control problem is proposed and solved, illustrating the procedure. Simulations show an effective reduction in the number of infectious individuals. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: This is a preprint of a paper whose final and definite form will be published in the volume 'Mathematics of Planet Earth' that initiates the book series 'CIM Series in Mathematical Sciences' (CIM-MS) published by Springer. Submitted Aug 2013; Revised and Accepted June 2014

MSC Class: 92D30; 49M05

Journal ref: Dynamics, Games and Science, CIM Series in Mathematical Sciences 1 (2015), 701--722

Modeling TB-HIV syndemic and treatment

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: Tuberculosis (TB) and human immunodeficiency virus (HIV) can be considered a deadly human syndemic. In this article, we formulate a model for TB and HIV transmission dynamics. The model considers both TB and acquired immune deficiency syndrome (AIDS) treatment for individuals with only one of the infectious diseases or both. The basic reproduction number and equilibrium points are determined and s… ▽ More Tuberculosis (TB) and human immunodeficiency virus (HIV) can be considered a deadly human syndemic. In this article, we formulate a model for TB and HIV transmission dynamics. The model considers both TB and acquired immune deficiency syndrome (AIDS) treatment for individuals with only one of the infectious diseases or both. The basic reproduction number and equilibrium points are determined and stability is analyzed. Through simulations, we show that TB treatment for individuals with only TB infection reduces the number of individuals that become co-infected with TB and HIV/AIDS, and reduces the diseases (TB and AIDS) induced deaths. Analogously, the treatment of individuals with only AIDS also reduces the number of co-infected individuals. Further, TB-treatment for co-infected individuals in the active and latent stage of TB disease, implies a decrease of the number of individuals that passes from HIV-positive to AIDS. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: This is a preprint of a paper whose final and definite form is: Journal of Applied Mathematics (ISSN 1110-757X) 2014, Article ID 248407, http://dx.doi.org/10.1155/2014/248407

MSC Class: 34D30; 92D30; 93A30

Journal ref: J. Appl. Math. 2014 (2014), Art. ID 248407, 14 pp

An optimal control approach to malaria prevention via insecticide-treated nets

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: Malaria is a life threatening disease, entirely preventable and treatable, provided the currently recommended interventions are properly implemented. These interventions include vector control through the use of insecticide-treated nets (ITNs). However, ITN possession does not necessarily translate into use. Human behavior change interventions, including information, education, communication (IEC)… ▽ More Malaria is a life threatening disease, entirely preventable and treatable, provided the currently recommended interventions are properly implemented. These interventions include vector control through the use of insecticide-treated nets (ITNs). However, ITN possession does not necessarily translate into use. Human behavior change interventions, including information, education, communication (IEC) campaigns and post-distribution hang-up campaigns are strongly recommended. In this paper we consider a recent mathematical model for the effects of ITNs on the transmission dynamics of malaria infection, which takes into account the human behavior. We introduce in this model a supervision control, representing IEC campaigns for improving the ITN usage. We propose and solve an optimal control problem where the aim is to minimize the number of infectious humans while keeping the cost low. Numerical results are provided, which show the effectiveness of the optimal control interventions. △ Less

Submitted 9 June, 2013; originally announced June 2013.

Comments: This is a preprint of a paper whose final and definite form will appear in Conference Papers in Mathematics, Volume 2013, Article ID 658468, http://dx.doi.org/10.1155/2013/658468. Paper Submitted 10-May-2013; accepted after minor revision 09-June-2013

MSC Class: 92D30; 49M05

Journal ref: Conference Papers in Mathematics, vol. 2013, Article ID 658468, 8 pages, 2013

Optimal control for a tuberculosis model with reinfection and post-exposure interventions

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We apply optimal control theory to a tuberculosis model given by a system of ordinary differential equations. Optimal control strategies are proposed to minimize the cost of interventions, considering reinfection and post-exposure interventions. They depend on the parameters of the model and reduce effectively the number of active infectious and persistent latent individuals. The time that the opt… ▽ More We apply optimal control theory to a tuberculosis model given by a system of ordinary differential equations. Optimal control strategies are proposed to minimize the cost of interventions, considering reinfection and post-exposure interventions. They depend on the parameters of the model and reduce effectively the number of active infectious and persistent latent individuals. The time that the optimal controls are at the upper bound increase with the transmission coefficient. A general explicit expression for the basic reproduction number is obtained and its sensitivity with respect to the model parameters is discussed. Numerical results show the usefulness of the optimization strategies. △ Less

Submitted 9 May, 2013; originally announced May 2013.

Comments: This is a preprint of a paper whose final and definite form will appear in Mathematical Biosciences. Paper submitted 9-Dec-2011; revised 9-Jun-2012, 13-Jan-2013, and 7-May-2013; accepted for publication 9-May-2013. arXiv admin note: substantial text overlap with arXiv:1203.3255

MSC Class: 92D30; 49M05

Journal ref: Math. Biosci. 244 (2013), no. 2, 154--164

Optimal control strategies for tuberculosis treatment: a case study in Angola

Authors: Cristiana J. Silva, Delfim F. M. Torres

Abstract: We apply optimal control theory to a tuberculosis model given by a system of ordinary differential equations. Optimal control strategies are proposed to minimize the cost of interventions. Numerical simulations are given using data from Angola. We apply optimal control theory to a tuberculosis model given by a system of ordinary differential equations. Optimal control strategies are proposed to minimize the cost of interventions. Numerical simulations are given using data from Angola. △ Less

Submitted 14 March, 2012; originally announced March 2012.

Comments: This is a preprint of a paper whose final and definite form will appear in the international journal Numerical Algebra, Control and Optimization (NACO). Paper accepted for publication 15-March-2012

MSC Class: 92D30; 49M05

Journal ref: Numer. Algebra Control Optim. 2 (2012), no. 3, 601--617