Arxiv:2103.08399V2 [Math.OC] 9 Apr 2021 Abstract Eiru Foeaos .Kt 4 Losuidasize-Stru a Studied Also [4] Kato N

arXiv:2103.08399v2 [math.OC] 9 Apr 2021 Abstract eiru foeaos .Kt 4 losuidasize-stru a studied also [4] Kato N. as operators, model a of population linear semigroup size-structured of evolving semigroup various an involves discussed considers analysis [3] Webb their and earliest. h the structure extensive are an have [1,2] on) Lotka so and and maturity, and characteristics length statistical diameter, or physiological shows which Words ones. Key existing g the complement a and minimizes new which are model article population our given principle the variational opti Ekeland’s to the controller Using obtain sets. We environm different conditions. three spatial rat optimality fertility a derive takes we in which system, problem position control optimal a the qualitatively occupy study to individuals semigroups which of method in the apply We ronment. nsrcue ouainmdl scmae ounstructure to compared as developed. models population di structured in or in an stages developmental of different characteristics through life progression and and complexi behaviour the optimal enhances Fur the stages that developmental pattern. or behavioral size complex of a terms have w usually natural population the a in processes evolutionary long-term Modeling Introduction 1 Classification Subject AMS operators. of Semigroup principle, nlsso iuieSz-tutrdPplto oe and Model Population Size-Structured Diffusive of Analysis ahmtclmdl fpplto yaisincorporating dynamics population of models Mathematical hswr drse h pia it oto rbe o i for problem control birth optimal the addresses work This : ∗ orsodn uhrseal:sba.ikgalcm abbas@iitmand [email protected], : email author’s Corresponding iesrcue ouainwt iuin pia it co birth Optimal diffusion, with population structured Size : 32;4J5 71;65M25 97M10; 47J35; 93C20; : ninIsiueo ehooyMandi, Technology of Institute Indian aad(.. 705 India 175005, - (H.P.) Kamand ao ua,Se Abbas Syed Kumar, Manoj colo ai Sciences, Basic of School it Control Birth 1 fidvdas o xml,sz a eae mass, age, be can size example, for individuals, of h xsec,aduiuns fotmlbirth optimal of uniqueness and existence, the , iiulsol rdal le ihmaturation with alter gradually should dividual sacnrlvral.Wt h epo adjoint of help the With variable. control a as e yaia ytmi natatsae Using space. abtract an in system dynamical a dnnieroeaos ehdo semigroups of Method operators. nonlinear nd rdi hlegn akbcueidvdasin individuals because task challenging a is orld eetszs u omr opeiisinvolved complexities more to Due sizes. fferent trdpplto oe ihdffso term. diffusion with model population ctured iesrcue ouainmdl ihspatial with models population size-structured iso vltoaymdl.As,i sevident is it Also, models. evolutionary of ties ouainmdl,eitn oesaeless are models existing models, population d aiycniin yfiigtebrhrt on rate birth the fixing by conditions mality hroe tutrn aua ouain in populations natural structuring thermore, soy h g-tutrdmdl ySharpe by models age-structured The istory. vncs ucinli hw.Otoe of Outcomes shown. is functional cost iven nlz iesrcue ouainmodel population size-structured a analyze n.Wt netpplto nmn,we mind, in population insect With ent. iesrcue(iei otnosvariable continuous a is (size size-structure ∗ vsv pce nasailenvi- spatial a in species nvasive to,Eeadsvariational Ekeland’s ntrol, i.ac.in Optimal Control in a size-structured population model is the process of forcing a population through a controller to obtain a certain behaviour on it [5, 6]. The objective in these types of control problems is to obtain a controller that is optimal in the sense that it takes minimum cost to reduce vermin population or to increase total harvest. So, the task in these control problems is to either minimize a cost functional or maximize total harvest. For size-structured population models, these problems can be categorized either as optimal harvesting problems (optimal harvesting of natural or farmed populations such as fish or plants) or optimal control problems (optimal control of vermin or pest population). There is enough literature available on the optimal control problems which take fertility rate as a control variable [5,7–17] and also on optimal harvesting problems [18–21]. Brokate [22] studied Pontryagin’s principle for a general problem of optimal cotrol for an age structured population model. Also recent work on optimal behavioural strategies in structured populations is discussed in [23]. Motivated by the work of N. Kato [4], in this paper, we establish important estimates on mild solutions in terms of fertility rates and initial population distributions. The main focus of the N. Kato’s work is to show the existence and uniqueness of mild solution, but in this work our main focus is to study the optimal birth control problem. To study optimal birth control problem, we need some estimates on mild solution and on the solution of adjoint equation. Rong et al. [7] studied the least cost size and least cost deviation problem for a nonlinear size-structured population model which takes the fertility rate of vermin population as a control variable. Ze-Rong et al. in [8] studied the optimal birth control problem without the diffusion term. The technique used to find the optimality conditions is quite different from the technique used by [7] and [8]. The technique used by [7] and [8] are tangent-normal cone technique, but in our case we construct three sets and by fixing birth rate on each set, we obtain the optimality conditions. Using Ekeland’s variational [24] principle and under some some assumptions, we show the existence and uniqueness of optimal controller as well. This paper is divided into five sections. In section 2, we formulated the abstract form of our model with some assumptions. In section 3, we discuss some preliminary results which we will require in other sections. In section 4, we discuss existence and uniqueness results and also derived some estimates on mild solution. The last section is devoted to the study of optimal birth control problem.

2 Model Formulation

Let us consider an evolving size-structured population living in a habitat Ω Rn with smooth boundary ⊂ ∂Ω. The spatial movement of individuals is controlled by diﬀusion process. Let p(s,t,x) be the population density of individuals which depends on size variable s [0,sf ], time t [0, T ] and spatial variable x Ω. ∈ ∈ ∈ We assume that size of individuals increases in same manner and also individuals do not move out of domain Ω through ∂Ω. Integrating population density w.r.t. size and spatial variable at any time t in the size interval (s ,s ) and in a subset Ω Ω will give the total population 1 2 1 ⊂ s2 P (t)= p(s,t,x)dxds. Zs1 ZΩ1 Let γ(s,t) be size and time dependent growth rate of individuals (growth rate here means a change in size with respect to time), µ(s,t,x), β(s,t,x) are mortality and reproduction rates respectively of size s individuals at time t in position x. Let f(s,t,x) and C(t,x) respectively be the number of s-size and zero size individuals coming from outside the domain Ω and let r(s,t,x) be ratio of female individuals. For our convenience, we will use the following notations in our article ΩT = (0, T ) Ω,Ωs = (0,sf ) Ω, × ×

2 ΩTs = (0,sf ) (0, T ) Ω and ΣTs = (0,sf ) (0, T ) ∂Ω. Now, size structured population model with × × × × diﬀusion can be formulated as

∂p(s,t,x) ∂ ∂t + ∂s (γ(s,t)p(s,t,x)) = k∆p(s,t,x) µ(s,t,x)p(s,t,x)+ f(s,t,x) in ΩTs sf − γ(0,t)p(0,t,x)= C(t,x)+ 0 r(s,t,x)β(s,t,x)p(s,t,x)ds in ΩT (2.1)  ∂p  ∂ν (s,t,x)=0 in ΣTs,R p(s, 0,x)= p0(s,x) in Ωs. We assume β(s,t,x) to be zero for those individuals who are not participating in the reproduction process. ∂ Constant k is a diffusion coefficient and ∂ν denotes the derivative in the direction of outward normal to the boundary ∂Ω. The homogeneous Neumann boundary condition shows that there is no flux of the population through the boundary of Ω and p0 is the initial population density. Also if we take the Dirichlet boundary condition i.e. p(s,t,x)=0onΣTs would mean individual will die as soon as they reached the boundary, which does not make sense for most of the species. So, the Neumann boundary condition is an appropriate choice. Let ∞ = u L (ΩTs) φl(s,t,x) u(s,t,x) φm(s,t,x) a.e in ΩTs (2.2) U { ∈ | ≤ ≤ } ∞ β where φl and φm are non negative functions lies in L (ΩTs). Let p (s,t,x) be the solution of (2.1), then our problem can be formulated as

T sf 1 minimize J(β)= pβ(s,t,x) ρ(β(s,t,x)2) dsdtdx (2.3) − 2 Z0 Z0 ZΩ subject to β ( [8] have also considered similar cost functional but without diﬀusion term). Here, ρ is ∈U a positive constant which is the weight factor of the cost to execute the control. J(β) can be considered as an energy function and our task is to show the existence of birth function β which minimizes J(β). It is clear from the energy function that if the cost of implementing control is high, then the fertility rate of invasive species will be lower. So, we want to minimize J(β) given in (2.3) subject to (2.1). Some authors ( [5–7]) have also taken energy function as

T sf 1 J(β)= pβ(s,t,x)+ ρ(β(s,t,x)2) dsdtdx. (2.4) 2 Z0 Z0 ZΩ Same analysis can be applied to this energy function with some minor modiﬁcations.

2.1 Abstract Formulation This model without spatial diﬀusion can be analyzed without converting it into abstract form, but to handle spatial diﬀusion we will convert this model into abstract form [4] in an appropriate Banach space. Let A be realization of Laplacian in Lq(Ω), q in (1, ) with the boundary conditions of Neumann type i.e. ∞ ∂v D(A)= v W 2,q(Ω) = 0 a.e on ∂Ω ∈ | ∂ν Aψ = k∆ψ for ψ D(A). ∈

3 It is well known that A will generate an analytic semigroup S(t),t 0 and there exist a positive constant ωt q ≥ C and ω R, such that S(t)ψ q Ce ψ q for ψ L (Ω). ∈ k kL (Ω)≤ k kL (Ω) ∈ For ψ Lq(Ω), let us deﬁne the operators M and B which are bounded linear by ∈ [M(s,t)ψ](x)= µ(s,t,x)ψ(x), [B(s,t)ψ](x)= β(s,t,x)ψ(x)

also [C(t)](x)= C(t,x), [f(s,t)](x)= f(s,t,x) and [p0(s)](x)= p0(s,x). Note that M ∞ µ ∞ , B ∞ β ∞ . (2.5) k kL (STf ;L(Z)) ≤ k kL (ΩTs ) k kL (STf ;L(Z)) ≤ k kL (ΩTs ) Above inequalities will be useful while proving existence and uniqueness of optimal control. Let Z be a Banach space with norm . In abstract form our model reduces to kk ∂p ∂ (s,t)+ (γ(s,t)p(s,t)) = [A M(s,t)]p(s,t)+ f(s,t), (s,t) ST = [0,sf ] [0, T ] ∂t ∂s − ∈ f × γ(0,t)p(0,t)= C(t)+ sf r(s,t)B(s,t)p(s,t)ds, t [0, T ] (2.6)  0 ∈ p(s, 0) = p (s), s [0,sf ]. 0 ∈ R Our main task here is to analyze p(s,t) the Z valued function, which describes the distribution of population at time t with respect to size s. We need the following assumptions to do qualitative analysis of our model:

(A1) γ: ST [0, ) is continuously diﬀerentiable w.r.t both s and t. Also γ(s,t) is positive on ST i.e. f 7→ ∞ f size always increases after crossing initial size upto ﬁnal size sf . We assume that the growth rate is non-negative and for each t [0, T ], we consider the following cases: ∈ (a) γ(0,t) > 0 and γ(sf ,t) > 0 (b) γ(0,t) > 0 and γ(sf ,t) = 0 (c) γ(0,t)=0 and γ(sf ,t) > 0 (d) γ(0,t)=0 and γ(sf ,t) = 0.

(A2) A is the inﬁnitesimal generator of a C semigroup S(t) t 0 in Z and S(t)ψ Z C˜ ψ Z , where 0 { | ≥ } k k ≤ k k C˜ is a generic constant.

∞ (A3) M,B L (ST ; (Z)), where (Z) is the space of all bounded linear operators in Z. ∈ f L L 1 1 1 (A4) f L (ST ; Z),C L (0, T ; Z) and p L (0,sf ; Z). ∈ f ∈ 0 ∈ (A5) 0 < r(s,t,x) < 1 (s,t,x) ΩT . ∀ ∈ s We can extend the function γ(s,t) on R [0, T ] as follows: × Assume that the growth of individuals with size less than zero is same as size zero individuals. And the growth of individuals with size greater than sf is same as size sf individuals. Due to assumption (A1), there exist a unique solution of the initial value problem

d s(t)= γ(s(t),t), s(t )= s where s R. (2.7) dt 0 0 0 ∈ Remark 2.1. The unique solution of initial value problem (2.7) is called characteristic curve of the system (2.1).

4 Let us define s(t)= ψ(t; t0,s0), z0(t)= ψ(t; 0, 0), and z1(t)= ψ(t; T,sf ). Now, we will define the initial and final time using the idea of [4] Case 1: For the case (A1)-(a), we have z (t) > 0 for t> 0 and z (t)

τ0, if s z0(t) τ0(t,s)= ≤ (2.8) (0, if s > z0(t).

Similarly for (s,t) ST satisfying z (t) s, there exists a unique τ [0, T ] such that ψ(t; τ ,sf )= s and ∈ f 1 ≤ 1 ∈ 1 the final time is defined by τ1, if s z1(t) τ1(t,s)= ≥ (2.9) (T, if s < z1(t). Similarly, we can define initial and final time for the other cases.

3 Preliminaries

The definition of lower semicontinuity and other conditions for a function to achieve infimum are given in [5]. Definition: Let Z be a topological space and z Z, then the function 0 ∈ f: Z R , is said to be lower semi-continuous function at z if 7→ ∪ {−∞ ∞} 0

lim inf f(z) f(z0). z→z0 ≥ Propostion 3.1. ( [5]) Let Z be a Banach space which is also reﬂexive and f: Z R , be a 7→ ∪ {−∞ ∞} lower semi-continuous convex function. If K is a bounded, closed and convex subset of Z, then f attains its inﬁmum on K.

Corollary 3.2. ( [5]) Let Z be a reﬂexive Banach space and f: Z R , be a lower semi- 7→ ∪ {−∞ ∞} continuous convex function such that lim f(z)=+ , then f attains its inﬁmum on Z. kzk→z0 ∞ The Ekeland variational principle is very useful to prove the existence of a minimum point in the absence of standard compactness condition [24]. The principle can be stated as follows:

Theorem 3.3. ( [24]) Let (Z, d) be a complete metric space and f: Z ( , ] be a lower semicontinuous 7→ −∞ ∞ function, bounded from below and nonidentically + . Also let ǫ> 0 and z Z be such that ∞ ∈ f(z) inf f(z) z Z + ǫ. ≤ { | ∈ }

Then for any λ> 0 there exist zǫ Z such that ∈

f(zǫ) f(z), ≤

d(zǫ, z) λ), ≤ −1 f(zǫ)

5 Now, we will state two important inequalities called Gronwall inequalities [5].

Lemma 3.4. ( [5]) Let f: [a, b] R (a, b R, a < b) be a continuous function, ϕ L∞(a, b) and 7→ ∈ ∈ λ L1(a, b), λ(t) 0 a.e. t (a, b). If ∈ ≥ ∈ t f(t) ϕ(t)+ λ(s)f(s)ds, ≤ Za for each t [a, b], then ∈ t t f(t) ϕ(t)+ ϕ(s)λ(s)exp λ(τ)dτ ds, ≤ Za Zs for each t [a, b]. ∈ Lemma 3.5. ( [5]) If in addition to Lemma 3.4, let ϕ(t)= C, where C is constant, then

t f(t) Cexp λ(s)ds . ≤ Za 4 Existence, Uniqueness and Estimates

In this section, our task is to obtain important estimates on mild solution. These estimates will be important while studying optimal birth control problem. Firstly, we will give key steps to derive mild solution. Similar form of mild solution is also given by N. Kato [4], but without the ratio r. The key steps to derive mild solution are given for the convenience of readers. Suppose p(s,t) satisﬁes (2.1) in strict sense and let

ζ u(ζ; t,s) = exp ∂sg(ψ(η; t,s), η))dη p(s(ζ),ζ), where ψ(ζ; t,s)= s(ζ). Zτ0 Diﬀerentiating u(ζ; t,s) with respect to ζ, we get

d ζ (u(ζ; t,s)) = [A M(ψ(ζ; t,s),ζ)]u(ζ; t,s) +exp ∂sg(ψ(η; t,s), η))dη f(s(ζ),ζ). (4.1) dζ − Zτ0 Using the variation of constant formula and after substituting ζ = t, we obtain

t p(s,t)= Q(τ )S(t τ )p(0,τ )+ Q(η)S(t η)[ M(s(η), η)p(s(η), η)) + f(s(η), η)]dη, 0 − 0 0 − − Zτ0 where t Q(τ ) = exp ∂sg(ψ(η; t,s), η)dη . 0 − Zτ0 For the case (A1)-(a) and (A1)-(b), γ(0,t) > 0, so p(0,t) is deﬁned by

1 sf p(0,t)= C(t)+ r(s,t)B(s,t)p(s,t)ds for t (0, T ). (4.2) γ(0,t) ∈ Z0

6 Also, for the case (A1)-(a) and (A1)-(b), we have

τ (t,s)= τ for s (0, z (t)) and τ (t,s)=0 for s (z (t),sf ). 0 0 ∈ 0 0 ∈ 0 Moreover, in this case for a.e t (0, T ), p(s,t) is given by ∈ Q(τ )S(t τ )p(0,τ )+ t Q(η)S(t η)[ M(s(η), η)p(s(η), η)) + f(s(η), η)]dη ,s (0, z (t)) 0 0 0 τ0 0 p(s,t)= − t − − ∈ (Q(0)S(t)p0 + Q(η)S(t η)[ M(s(η), η)p(s(η), η)) + f(s(η), η)]dη ,s (z0(t),sf ). 0 R− − ∈

Note that p0 = p0(ψ(0; t,sR )) and for the case (A1)-(c) and (A1)-(d), there is no reproduction rate, so for these cases, we have

t p(s,t)= Q(0)S(t)p + Q(η)S(t η)[ M(s(η), η)p(s(η), η)) + f(s(η), η)]dη for (s,t) ST . 0 − − ∈ f Z0 We deﬁne the above solution with (4.2) as mild solution.

∞ Remark 4.1. Suppose µ, β L (ΩTs) and assume that µ, β are non negative almost everywhere in ΩTs. 1 ∈ 1 Let f L (ΩTs),f(s,t,x) 0 a.e. (s,t,x) ΩTs and C L (ΩT ),C(t,x) 0 a.e. (t,x) ΩT and also ∈1 ≥ ∈ ∈ ≥ ∈ p0 L (Ωs),p0(s,x) 0 a.e. (s,x) Ωs. If the semigroup S(t) t 0 satisﬁes S(t)Z+ Z+, we have ∈ ∞ ≥ 1 ∈ { | ≥ } ⊂ existence of p L (0, T ; L (0,sf ; Z )) and p(s,t) Z a.e. (s,t) STf , where Z is a positive cone with ∈ + ∈ + ∈ + vertex 0.

k Let Cψ(STf ; Z) be the collection of those functions which are Z-valued and k times continuously dif- ferentiable along the characteristic curve ψ. Suppose Dψp(s,t) is the derivative along the characteristic ψ [p(ψ(t + h; t,s),t + h) p(s,t)] Dψp(s,t) := lim − . h→0 h Then, it is easy to see that d D p(ψ(η; t,s), η)= p(ψ(η; t,s), η). ψ dη ∞ 1 Theorem 4.2. For the Case (A1)-(a) and Case (A1)-(b), any mild solution p L (0, T ; L (0,sf ; Z)) ∈ of (2.1) is continuously diﬀerentiable along almost every characteristic ψ and satisﬁes

Dψp(s,t)= Ap(s,t) ∂sγ(s,t)p(s,t) M(s,t)p(s,t)+ f(s,t) a.e (s,t) ST − − ∈ f γ(0,t)p(0,t)= C(t)+ sf r(s,t)B(s,t)p(s,t)ds, a.e t (0, T ) (4.3)  0 ∈ p(s, 0= p (s), a.e s (0,sf ) 0 R ∈ where p(s, 0) and p(0,t) are given by the following limits:

p(s, 0) = lim p(ψ(η; 0,s), η) a.e s (0,sf ), η→0+ ∈ p(0,t) = lim p(ψ(η; t, 0), η) a.e t (0, T ), η→t ∈ and in case (A1)-(c) and case (A1)-(d)

Dψp(s,t)= Ap(s,t) ∂sγ(s,t)p(s,t) M(s,t)p(s,t)+ f(s,t) a.e (s,t) ST − − ∈ f (4.4) p(s, 0) = p (s), a.e s (0,sf ). ( 0 ∈ 7 ∞ 1 1 Conversely, if p L (0, T ; L (0,sf ; Z)) C (STf ; Z) satisﬁes (4.3) and (4.4), then p is a mild solution ∈ ∩ ψ of (2.6).

Proof. The proof of necessary part is easy and follows just by using the fact that A generates an analytic ∞ 1 1 semigroup. Conversely, suppose that p L (0, T ; L (0,sf ; Z)) C (STf ; Z) and satisﬁes (4.3). So, ∈ ∩ ψ Dψp(s,t)= Ap(s,t) ∂sγ(s,t)p(s,t) M(s,t)p(s,t)+ f(s,t) a.e (s,t) ST − − ∈ f because d D p(ψ(η; t,s), η)= p(ψ(η; t,s), η). ψ dη Diﬀerentiating p with respect to η, we obtain d p(ψ(η; t,s), η)= Ap(ψ(η; t,s), η) ∂sγ(ψ(η; t,s), η)p(ψ(η; t,s), η) dη − M(ψ(η; t,s), η)p(ψ(η; t,s), η)+ f(ψ(η; t,s), η) − p(ψ(τ0; t,s),τ0)= p(0,τ0).

Therefore, by variation of constants formula, the solution starting from initial time τ0 with η = t is given by t Q(τ )S(t τ )p(0,τ )+ Q(η)S(t η)[ M(s(η), η)p(s(η), η)) + f(s(η), η)]dη, 0 − 0 0 − − Zτ0 t where Q(τ0) = exp ∂sg(ψ(η; t,s), η)dη . Also p(ψ(0; t,s), 0) = p0(ψ(0; t,s)), so p is a mild solution − τ0 to (2.6). R

One can observe that s = ζ = τ0(t,s), which implies s = ψ(t; ζ, 0) and using (2.7), we have ds t = γ(0,ζ)exp ∂sg(ψ(η; ζ, 0), η)dη . dζ − − Zζ Moreover s = φ = ψ(ζ; t,s), which implies s = ψ(t; ζ, φ) and again using (2.7), we have ds t = exp ∂ γ(ψ(η; t,s), η)dη . dφ s Zζ ∞ 1 Suppose that the assumptions (A1)-(A5) holds, and suppose that p ,p L (0, T ; L (0,sf ; Z)) be mild 1 2 ∈ solutions to (2.1). Then for t [0, T ] ∈ z0(t) p ( ,t) p ( ,t) 1 Q(τ )S(t τ )(p (0,τ ) p (0,τ )) Z ds k 1 · − 2 · kL (0,sf ;Z)≤ k 0 − 0 1 0 − 2 0 k Z0 z0(t) t sf t + Q(η)S(t η)M(s(η), η)(p p ) Z dηds + Q(η)S(t η)M(s(η), η)(p p ) Z dηds. k − 1 − 2 k k − 1 − 2 k Z0 Zτ0 Zz0(t) Z0 Using the above transformation s = ζ to ﬁrst integral on R.H.S and s = φ to second and third integral on R.H.S,we will get t p ( ,t) p ( ,t) 1 C˜ p ( , η) p ( , η) 1 dη, k 1 · − 2 · kL (0,sf ;Z)≤ k 1 · − 2 · kL (0,sf ;Z) Z0 8 where C˜ is generic constant which depends on B ∞ , M ∞ , T and the bound of k kL (STf ;L(Z)) k kL (STf ;L(Z)) semigroup S(t). By using Gronwall’s inequality, we get

p ( ,t) p ( ,t) 1 = 0, k 1 · − 2 · kL (0,sf ;Z) which shows that the mild solution is unique. ∞ 1 Theorem 4.3. Suppose that the assumptions (A1)-(A5) hold and p ,p L (0, T ; L (0,sf ; Z)) be mild 1 2 ∈ solutions corresponding to the initial data p01 and p02 respectively. Then

p ( ,t) p ( ,t) 1 C˜ p p 1 , k 1 · − 2 · kL (0,sf ;Z)≤ k 01 − 02kL (0,sf ;Z) where C˜ is a generic constant which depends on B ∞ , M ∞ , T and the bound of the k kL (STf ;L(Z)) k kL (STf ;L(Z)) semigroup S(t),t 0 . { ≥ } ∞ 1 Proof. For t [0, T ] and p ,p L (0, T ; L (0,sf ; Z)), we have ∈ 1 2 ∈ z0(t) 1 p1( ,t) p2( ,t) L (0,sf ;Z) Q(τ0)S(t τ0)(p1(0,τ0) p2(0,τ0)) Z ds k · − · k ≤ 0 k − − k Z sf + Q(0)S(t)(p p ) Z ds k 01 − 02 k Zz0(t) z0(t) t + Q(η)S(t η)M(s(η), η)(p p ) Z dηds k − 1 − 2 k Z0 Zτ0 sf t + Q(η)S(t η)M(s(η), η)(p p ) Z dηds. k − 1 − 2 k Zz0(t) Z0 Let us denote z0(t) I = Q(τ )S(t τ )(p (0,τ ) p (0,τ )) Z ds, 1 k 0 − 0 1 0 − 2 0 k Z0 sf I = Q(0)S(t)(p p ) Z ds, 2 k 01 − 02 k Zz0(t) z0(t) t I = Q(η)S(t η)M(s(η), η)(p p ) Z dηds, 3 k − 1 − 2 k Z0 Zτ0 sf t I = Q(η)S(t η)M(s(η), η)(p p ) Z dηds. 4 k − 1 − 2 k Zz0(t) Z0 Now, again using the transformation s = ζ = τ0(t,s) on I1, s = φ = ψ(0; t,s) on I2 and s = φ = ψ(ζ; t,s) on I3 and I4, we get t p ( ,t) p ( ,t) 1 C˜ p p 1 +C˜ p ( ,ζ) p ( ,ζ) 1 , k 1 · − 2 · kL (0,sf ;Z)≤ k 01 − 02kL (0,sf ;Z) k 1 · − 2 · kL (0,sf ;Z) dζ Z0 where C˜ is generic constant which depends on B ∞ , M ∞ , T and the bound of the k kL (STf ;L(Z)) k kL (STf ;L(Z)) semigroup S(t). By using the Gronwall’s inequality, we obtain

p ( ,t) p ( ,t) 1 C˜ p p 1 . k 1 · − 2 · kL (0,sf ;Z)≤ k 01 − 02kL (0,sf ;Z)

9 In a similar manner, we can state the following result:

β β ∞ 1 Theorem 4.4. Suppose that the assumptions (A1)-(A5) hold and p 1 ,p 2 L (0, T ; L (0,sf ; Z)) be ∈ solutions corresponding to the fertility rates β1 and β2 respectively. Then

β1 β2 p ( ,t) p ( ,t) C˜exp T B1 B2 ∞ , (4.5) 1 L (STf ;L(Z)) · − · L (0,sf ;Z) ≤ k − k where C˜ is a positive constant which depends on M ∞ , T and the bound of the semigroup k kL (STf ;L(Z)) S(t),t 0 . { ≥ } Moreover, if C˜ C1 B ,B and for some positive constant C , then kB −B k ∞ L 1 2 1 ≤ 1 2 L (STf ; (Z)) ∀ ∈U

β1 β2 p ( ,t) p ( ,t) M1 B1 B2 ∞ , (4.6) 1 L (STf ;L(Z)) · − · L (0,sf ;Z) ≤ k − k

where M1 is a positive constant. β β ∞ 1 Proof. For t [0, T ] and p 1 ,p 2 L (0, T ; L (0,sf ; Z)), we have ∈ ∈

z0(t) β1 β2 β1 β2 p ( ,t) p ( ,t) 1 Q(τ )S(t τ )(p (0,τ ) p (0,τ )) Z ds k · − · kL (0,sf ;Z) ≤ k 0 − 0 0 − 0 k Z0 z0(t) t β1 β2 + Q(η)S(t η)M(s(η), η)(p p ) Z dηds k − − k Z0 Zτ0 sf t β1 β2 + Q(η)S(t η)M(s(η), η)(p p ) Z dηds. k − − k Zz0(t) Z0 Let us denote z0(t) β1 β2 I = Q(τ )S(t τ )(p (0,τ ) p (0,τ )) Z ds 1 k 0 − 0 0 − 0 k Z0 z0(t) t β1 β2 I = Q(η)S(t η)M(s(η), η)(p p ) Z dηds 2 k − − k Z0 Zτ0 sf t β1 β2 I = Q(η)S(t η)M(s(η), η)(p p ) Z dηds 3 k − − k Zz0(t) Z0 In the integral I1, substituting s = ζ = τ0(t,s), we get

t sf sf I S(t ζ) B (s(η), η)pβ1 (s(η), η)dη B (s(η), η)pβ2 (s(η), η)dη dζ. 1 ≤ − 1 − 2 Z0 Z0 Z0 Z

β2 β1 In [4], it has been shown that if β2 β1, then p p . So, here we prove the other case. Without loss ≤ ≤ of generality assume that β β , then 2 ≥ 1 β pβ1 β pβ2 β pβ1 + β pβ2 β pβ2 β pβ1 = (β β )(pβ1 pβ2 ). 1 − 2 ≤ 1 2 − 1 − 2 1 − 2 − Therefore, t β1 β2 I1 C˜ B1 B2 ∞ p ( ,ζ) p ( ,ζ) dζ. L (STf ;L(Z)) 1 ≤ k − k 0 · − · L (0,sf ;Z) Z

10 Similarly, using s = φ = ψ(ζ; t,s), we can always ﬁnd C˜ such that

t β1 β2 I2 + I3 C˜ B1 B2 ∞ p ( ,ζ) p ( ,ζ) dζ. L (STf ;L(Z)) 1 ≤ k − k 0 · − · L (0,sf ;Z) Z

Therefore t β1 β2 β1 β2 p ( ,t) p ( ,t) C + C˜ B1 B2 ∞ p ( ,ζ) p ( ,ζ) dζ, (4.7) 1 L (STf ;L(Z)) 1 · − · L (0,sf ;Z) ≤ k − k 0 · − · L (0,sf ;Z) Z where C is a positive constant and C˜ is a generic constant which depends on M ∞ , T and the k kL (STf ;L(Z)) bound of the semigroup S(t),t 0 . Now using the Gronwall inequality, we get { ≥ }

β1 β2 p ( ,t) p ( ,t) Cexp CT˜ B1 B2 ∞ . 1 L (STf ;L(Z)) · − · L (0,sf ;Z) ≤ k − k

Since C is a positive constant, we can ﬁnd another constant C (just for convenience) such that

C C B B ∞ . ≤ k 1 − 2kL (STf ;L(Z)) Using the estimate on C˜ and Gronwall’s lemma, we get

β1 β2 C1T p ( ,t) p ( ,t) C B1 B2 ∞ e . 1 L (STf ;L(Z)) · − · L (0,sf ;Z) ≤ k − k

C1T Taking M1 = Ce , we obtain our desired result. Remark 4.5. Applying the same procedure as we followed to prove the uniqueness of mild solution to (2.1), we can prove that the bound of mild solution depends on vital rates, inﬂow of zero and s-size individuals and the initial population distribution. So, if we assume inﬂow of zero and s-size individuals as constant and also vital rates (mortality, fertility and growth rate) as constant, then mild solution will be global mild solution in time variable. Now, let us consider the following adjoint system to (4.3) and (4.4): In case (A1)-(a) and (A1)-(c),

∗ ∗ ∗ ∗ Dψφ(s,t)= A φ(s,t)+ M (s,t)φ(s,t) + [1 B (s,t)r(s,t)φ(0,t)] + f (s,t) a.e (s,t) STf − − ∈ φ(sf ,t) = limh φ(ψ(t h; t,sf ),t h) = 0 a.e t (0, T )  →+0 − − ∈ φ(s, T ) = limh φ(ψ(T h; T,s), T h) = 0, a.e s (0,sf ) →+0 − − ∈ (4.8)  and for case (A1)-(b) and (A1)-(d)

∗ ∗ ∗ ∗ Dψφ(s,t)= A φ(s,t)+ M (s,t)φ(s,t) + [1 B (s,t)r(s,t)φ(0,t)] + f (s,t) a.e (s,t) STf − − ∈ φ(sf ,t) = limh φ(ψ(t h; t,sf ),t h) = 0 a.e t (0, T ), ( →+0 − − ∈ (4.9) ∞ ∗ 1 ∗ ∗ where φ L (STf ; D(A )) C (STf ; Z ) is unknown. In adjoint system, A is the adjoint operator of ∈ ∩ ψ A and will generate the adjoint semigroup S∗(t) t 0 in the dual space Z∗ of Z, also B∗ and M ∗ are { | ≥ } adjoint operators for B and M respectively. M ∗ and B∗ are deﬁned by

[M ∗(s,t)φ](x)= µ(s,t,x)φ(x), [B∗(s,t)φ](x)= β(s,t,x)φ(x).

11 Note that ∗ ∗ M L∞(S ;L(Z∗)) µ L∞(Ω ), B L∞(S ;L(Z∗)) β L∞(Ω ). (4.10) k k Tf ≤ k k Ts k k Tf ≤ k k Ts

Now, our task is to convert the data φ(sf ,t) = φ(s, T ) = 0 into the data φ(0,t) = φ(s, 0) = 0 to relate system (4.8) and (4.3). So, the natural choice of transformation will be φ(s,t) = φ(sf s, T t). Our − − next job is to ﬁnd the characteristic curve along which φ satisfy an ordinary diﬀerential equation. Let γ(s,t)= γ(sf s, T t) and ψ(t; t ,s ) be the solution to the system − − 0 0 ds(t) = γ(s(t),t), s(t )= s , t [0, T ]. dt 0 0 ∈

By following the idea of [4], one can easily show that D φ(s,t)= Dψφ(sf s, T t). Then applying the ψ − − − same procedure as used in Theorem 4.2, we can show the existence of solution to the dual system (4.8) and (4.9). Also the uniqueness of the solution to the adjoint system follows from the following estimates

T ∗ ∗ ∗ φ(s,t) Z∗ C B ( ,ζ)r( ,ζ)φ(0,ζ) r( ,ζ)φ(0,ζ) f ( ,ζ) ∞ ∗ dζ. k k ≤ k · · − · − · kL (0,sf ;Z ) Z0

Following the same steps as followed in theorem 4.4, we will get the following estimates for dual problem:

1 Theorem 4.6. Suppose that the assumptions (A1)-(A5) hold and φβ , φβ2 be solutions to dual problem corresponding to fertility rates β1, β2 respectively, then

β1 β2 ∗ ∗ ∗ φ (s,t) φ (s,t) Z M2 B1 B2 L∞(S ;L(Z∗)). (4.11) k − k ≤ k − k Tf Remark 4.7. Using (4.10) estimate in theorem 4.6 can be written as

β β φ 1 (0,t,x) φ 2 (0,t,x) M β β ∞ t (0, T ), x Ω. | − |≤ 2k 1 − 2kL (ΩTs ) ∀ ∈ ∈ Let us deﬁne the concept of weak solution to (2.1) as follows 1 In case (A1)-(a) and case (A1)-(b), a function p L (STf ; Z) is called a weak solution to (2.1) in the ∈ following sense

T sf p(s,t,x)( Dψφ(s,t,x) k∆φ(s,t,x)+ µ(s,t,x)φ(s,t,x) r(s,t,x)β(s,t,x)φ(0,t,x))dxdsdt − − − Z0 Z0 ZΩ

sf T T sf = p0(s,x)φ(s, 0,x)dxds + C(t,x)φ(0,t,x)dxdt + f(s,t,x)φ(s,t,x)dxdsdt Z0 ZΩ Z0 ZΩ Z0 Z0 ZΩ where φ is any absolutely continuous function along almost every characteristic line and in case (A1)-(a) satisﬁes ∞ φ L (ΩTs) ∈ ∞ Dψφ + k∆φ µφ + βφ(0, , ) L (ΩTs)  − · · ∈ (4.12) φ(sf ,t) = 0 a.e t (0, T )  ∈ φ(s, T ) = 0 a.e s (0,sf ),  ∈   12 and in case (A1)-(b) satisﬁes

∞ φ L (ΩTs) ∈ ∞ Dψφ + k∆φ µφ + βφ(0, , ) L (ΩTs) (4.13)  − · · ∈ φ(s, T ) = 0 a.e s (0,sf ). ∈  Denote set of all φ satisfying (4.12) by Φ1 and satisfying (4.13) by Φ2. Since φ is absolutely continuous along almost every characteristic line, φ(sf ,t) and φ(s, T ) should be understood as

φ(sf ,t) = lim φ(ψ(t h; t,sf ),t h) and φ(s, T ) = lim φ(ψ(T h; T,s), T h). h→+0 − − h→+0 − −

1 Also in case (A1)-(c) and case (A1)-(d), a function p L (STf ; Z) is called a weak solution to (2.1) in the ∈ following sense

T sf p(s,t,x)( Dψφ(s,t,x) k∆φ(s,t,x)+ µ(s,t,x)φ(s,t,x))dxdsdt − − Z0 Z0 ZΩ

sf T sf = p0(s,x)φ(s, 0,x)dxds + f(s,t,x)φ(s,t,x)dxdsdt Z0 ZΩ Z0 Z0 ZΩ where φ is any absolutely continuous function along almost every characteristic line and lies in Φ1 for case (A1)-(c) and in Φ2 for case (A1)-(d). Now, question arises, what guarantees us the existence of weak solutions. For this we have the following result:

∞ 1 Theorem 4.8. The mild solution p L (0, T ; L (0,sf ; Z)) of (2.1) is also a weak solution of (2.1). ∈ ∞ 1 Proof. Let us assume that p L (0, T ; L (0,sf ; Z)) be a mild solution of (2.1). Also assume that φ Φ ∈ ∈ 1 in case (A1)-(a) and φ Φ in case (A1)-(b). Also let us deﬁne the integral ∈ 2 T −h sf 1 I = (p(ψ(t + h; t,s),t + h, x) p(ψ(t; t,s),t,x)) φ(s,t,x)dxdsdt. (4.14) h − Z0 Z0 ZΩ Then as h 0, I will converge to → T sf Dψp(s,t,x)φ(s,t,x)dxdsdt Z0 Z0 ZΩ

T sf = (∆p(s,t,x) ∂sγ(s,t)p(s,t,x) µ(s,t,x)p(s,t,x)) φ(s,t,x)dxdsdt − − Z0 Z0 ZΩ T sf + f(s,t,x)φ(s,t,x)dxdsdt (4.15) Z0 Z0 ZΩ T sf = (∆φ(s,t,x) ∂sγ(s,t)φ(s,t,x) µ(s,t,x)φ(s,t,x)) p(s,t,x)dxdsdt − − Z0 Z0 ZΩ T sf + f(s,t,x)φ(s,t,x)dxdsdt (4.16) Z0 Z0 ZΩ 13 Now, using the change of variables (s,t) = (y,ζ), where (s,t) is related to (y,ζ) by the relation (s,t) = (ψ(ζ h; ζ,y),ζ h). Then I can be written as − − 1 T ψ(ζ;ζ−h,sf ) ζ I = p(y,ζ,x)φ(ψ(ζ h; ζ,y),ζ h, x)exp ∂sγ(ψ(η; ζ,y), η)dη dxdydζ h − − Zh Zψ(ζ;ζ−h,0) ZΩ Zζ−h 1 T −h ψ(ζ;ζ−h,sf ) = p (φ(ψ(ζ h; ζ,y),ζ h, x) φ(y,ζ,x)) h − − − Zh Zψ(ζ;ζ−h,0) ZΩ ζ exp ∂sγ(ψ(η; ζ,y), η)dη dxdydζ × Zζ−h 1 T −h ψ(ζ;ζ−h,sf ) + pφ(ψ(ζ h; ζ,y),ζ h, x) h − − Zh Zψ(ζ;ζ−h,0) ZΩ ζ exp ∂sγ(ψ(η; ζ,y), η)dη 1 dxdydζ × − Zζ−h 1 T ψ(ζ;ζ−h,sf ) ζ + pφ(ψ(ζ h; ζ,y),ζ h, x)exp ∂sγ(ψ(η; ζ,y), η)dη dxdydζ h − − ZT −h Zψ(ζ;ζ−h,0) ZΩ Zζ−h 1 T −h ψ(ζ;ζ−h,sf ) p(y,ζ,x)φ(y,ζ,x)dxdydζ − h Zh Zψ(ζ;ζ−h,0) ZΩ 1 h sf p(y,ζ,x)φ(y,ζ,x)dxdydζ − h Z0 Z0 ZΩ Now, combining ﬁrst two integrals and taking other integrals separately, I will converge to

T sf p(s,t,x) ( Dψφ(s,t,x) ∂sγ(s,t)φ(s,t,x)) dxdsdt − − Z0 Z0 ZΩ T sf + γ(0,ζ)p(0,ζ,x)φ(0,ζ,x)dxdζ + p0(s,x)φ(s, 0,x)dxds (4.17) Z0 ZΩ Z0 ZΩ as h 0. Now, from (4.16) and (4.17) it is clear that p is a weak solution. Similar arguments will work → in case (A1)-(c) and case (A1)-(d).

5 Optimal Control Problem

In this section, we need to characterize tangent cones and normal cones. In [6], the following characterization is given for tangent cone and normal cone. The Tangent cone u( ) to the convex set (2.2) has the following characterization: v u( ) iﬀ almost T U 1 ∈T U everywhere on ΩTs v (s,t,x) 0 if u(s,t,x)= φl(s,t,x) 1 ≥ v (s,t,x) 0 if u(s,t,x)= φm(s,t,x). 1 ≤ The Normal cone u( ) to the convex set (2.2) has the following characterization: v u( ) iﬀ N U 1 ∈ N U almost everywhere on ΩTs v (s,t,x) 0 if u(s,t,x)= φm(s,t,x) 1 ≥

14 v1(s,t,x)=0 if φl(s,t,x) < u(s,t,x) < φm(s,t,x)

v (s,t,x) 0 if u(s,t,x)= φl(s,t,x). 1 ≤ Now, let us consider the following optimal birth control problem

T sf 1 Minimize J(β) = Minimize pβ(s,t,x) ρ(β(s,t,x))2 dxdsdt, (5.1) − 2 Z0 Z0 ZΩ where pβ(s,t,x) is the solution of (2.1) corresponding to the birth function β. The term 1 ρ(β)2 in the − 2 objective functional implies that if cost of birth control is higher, the fertility rate will be lower which also matches with our intuition. In this section we will assume that the assumptions (A1)-(A5) hold and also the assumptions of (Theorem 4.4) hold. Definition 5.1. A pair (β,pβ) is said be optimal pair for the birth control problem (5.1) if the following conditions are satisfied (1) β . ∈U (2) pβ solves the system(2.1). (3) (β,pβ) minimizes the functional J(β). Let us define the functional Ψ in as follows U T sf pβ(s,t,x) 1 ρ(β(s,t,x))2 dxdsdt, β Ψ(β)= 0 0 Ω − 2 ∈U (5.2) (+ , otherwise. R ∞R R Our first task is to prove that Ψ is lower semicontinuous. Lemma 5.2. The functional Ψ is lower semicontinuous.

1 Proof. Let βn be a sequence in L (ΩTs) converges to β. Without loss of generality assume that βn βn ∈ . Since p is bounded (Theorem4.4) and the bound depends upon βn, we can ﬁnd a subsequence U βn β βn(s,t,x) β(s,t,x) and p (s,t,x) p (s,t,x) in ΩTs. So, we have → →

βn 1 2 β 1 2 p (s,t,x)+ ρ(βn(s,t,x)) p (s,t,x)+ ρ(β(s,t,x)) . 2 → 2 Now, using the Lebesgue dominated convergence theorem, we can conclude that

T sf βn 1 2 p (s,t,x) ρ(βn(s,t,x)) dxdsdt − 2 Z0 Z0 ZΩ converges to T sf 1 pβ(s,t,x) ρ(β(s,t,x))2 dxdsdt − 2 Z0 Z0 ZΩ Therefore, by Fatou’s lemma,

T sf βn 1 2 lim inf p (s,t,x) ρ(βn(s,t,x)) dxdsdt n →∞ − 2 Z0 Z0 ZΩ T sf 1 pβ(s,t,x) ρ(β(s,t,x))2 dxdsdt, ≥ − 2 Z0 Z0 ZΩ 15 which is lim inf Ψ(βn) Ψ(β). n→∞ ≥ Therefore Ψ is lower semicontinuous.

β∗ β∗ Theorem 5.3. Let (β∗,p ) be optimal pair for the least cost size problem and φ be the solution to the adjoint system (4.8) in case (A1)-(a), (A1)-(c) and adjoint system (4.9) in case (A1)-(b), (A1)-(d) respectively. Then for f ∗(s,t,x)= c, where c is a positive constant − r(s,t,x)φβ∗ (0,t,x)pβ∗ (s,t,x) β (s,t,x)= F ∗ − cρ 1 ∞ F : L (ΩTs) L (ΩTs) is deﬁned by 7→

φl, if h(s,t,x) < φl

(F h)(s,t,x)= h(s,t,x), if φl h(s,t,x) φm (5.3)  ≤ ≤ φm, if h(s,t,x) > φm.

1 β∗+ǫδ∗ β∗  ∞ Proof. Let zǫ = [p p ], where δ = β β , then zǫ converges to z L (STf ; Z) and z satisﬁes ǫ − ∗ − ∗ ∈

Dψz(s,t)= Az(s,t) ∂sγ(s,t)z(s,t) M(s,t)z(s,t) a.e (s,t) ST − − ∈ f sf sf β∗ γ(0,t)z(0,t)= r(s,t)B∗(s,t)z(s,t)ds + δ∗p (s,t)ds, a.e t (0, T ) (5.4)  si 0 ∈ z(s, 0 = 0, a.e s (0,sf ), R ∈ R ∗ β∗ ∗ β∗ ∗ β∗ ∞ where δ p (s,t)= δ (s,t, )p (s,t, ) and also δ p L (STf ; Z). · · ∈ Because z(s,t) is a weak solution of (5.4), we have

T sf z(s,t,x)( Dψ φ(s,t,x) k∆φ(s,t,x)+ µ(s,t,x)φ(s,t,x) r(s,t,x)β(s,t,x)φ(0,t,x))dxdsdt − − − Z0 Z0 ZΩ T sf β∗ = δ∗(s,t,x)p (s,t,x)φ(s,t,x)dxdsdt. Z0 Z0 ZΩ Because φ(s,t) satisﬁes

∗ ∗ ∗ ∗ Dψφ(s,t)= A φ(s,t) M (s,t)φ(s,t)+ r(s,t)B (s,t)φ(0,t)+ f (s,t) 1, − − − we have

T sf T sf z(s,t,x)(f ∗(s,t,x) 1)dxdsdt = δ r∗(s,t,x)pβ∗ (s,t,x)φ(0,t,x)dxdsdt. (5.5) − ∗ Z0 Z0 ZΩ Z0 Z0 ZΩ β∗ Because (β∗,p ) is optimal pair, we have 1 lim [Ψ(β∗ + ǫδ∗) Ψ(β∗)] 0, ǫ→+0 ǫ − ≥ which implies T sf β∗+ǫδ∗ β∗ 2 p p ρǫ(δ∗) lim − ρβ∗δ∗ dxdsdt 0, ǫ→+0 ǫ − 2 − ≥ Z0 Z0 ZΩ 16 which further implies T sf (z(s,t,x) ρβ δ ) dxdsdt 0. − ∗ ∗ ≥ Z0 Z0 ZΩ For f ∗(s,t,x)= c + 1, where c is a nonnegative constant not equal to 1 and using (5.5), we get − T sf r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) (β β) ρβ + dxdsdt 0, ∗ − ∗ c ≥ Z0 Z0 ZΩ which implies

T sf r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) ρ(β β) β + dxdsdt 0. ∗ − ∗ cρ ≥ Z0 Z0 ZΩ Let r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) = (s,t,x) ΩTs < φl . V1 ∈ | − cρ Choose β = φl on . So, in this case V1 r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) ρ(β φl) β + dxdsdt = 0, ∗ − ∗ cρ ZV1 because ρ is a positive constant, we have β∗ = φl.

Similarly, let r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) = (s,t,x) ΩTs > φm V2 ∈ | − cρ and choose β = φm on . Then, we have V2 r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) ρ(β φm) β + dxdsdt = 0. ∗ − ∗ cρ ZV2

Since ρ is a positive constant, we have β∗ = φm. If r(s,t,x)pβ∗ (s,t,x)φβ∗ (0,t,x) = (s,t,x) ΩTs φl φm , V3 ∈ | ≤− cρ ≤ then on 3 V r(s,t,x)φβ∗ (0,t,x)pβ∗ (s,t,x) β = . ∗ − cρ

Theorem 5.4. Suppose that the assumptions (A1) (A5) hold and assumptions of theorem 5.3 hold with − M M + M M 1 4 2 3 < 1, cρ then the optimal control problem (5.1) has a unique optimal control β¯ ΩT where M and M are ∈ s 3 4 supremum of p and φ respectively. | | | |

17 Proof. We know that the functional Ψ deﬁned in (5.2) is lower semicontinuous, so using Ekeland’s variational principle, for given ǫ> 0, there exists βǫ such that ∈U

Ψ(βǫ) inf Ψ(β)+ ǫ (5.6) ≤ β∈U

Ψ(βǫ) inf Ψ(β)+ √ǫ β βǫ L1(Ω ) . (5.7) ≤ β∈U{ k − k Ts } Thus the perturbed functional

Ψǫ(β)=Ψ(β)+ √ǫ β βǫ 1 (5.8) k − kL (ΩTs ) attain its inﬁmum at βǫ. Therefore,

1 ′ lim Ψǫ(βǫ + ǫ δ∗) Ψǫ(βǫ) 0, ǫ′→0 ǫ′ − ≥ which implies 1 ′ ′ lim Ψ(βǫ + ǫ δ∗)+ √ǫ ǫ δ∗ Ψ(βǫ 0. ǫ′→0 ǫ′ k k− ≥ Following the same steps as in theorem 5.3, we get

T sf r(s,t,x)pβǫ (s,t,x)φβǫ (0,t,x) T sf ρδ βǫ + dxdsdt + √ǫ δ dxdsdt 0. (5.9) ∗ cρ | ∗| ≥ Z0 Z0 ZΩ Z0 Z0 ZΩ

We know that v lies in tangent cone u( ) iﬀ almost everywhere on ΩT 1 T U s

v (s,t,x) 0 if u(s,t,x)= φl(s,t,x) 1 ≥

v (s,t,x) 0 if u(s,t,x)= φm(s,t,x). 1 ≤ Thus, δ∗ u( ) and also δ∗ = β βǫ depends on β which lies in , which means (5.9) holds for any ∈ T U − U ∞ δ u( ). Therefore using the structure of normal cones, there exists θ(s,t,x) L (ΩT ), θ 1 ∗ ∈ T U ∈ s k k∞≤ such that r(s,t,x)pβǫ (s,t,x)φβǫ (0,t,x) ρβǫ + + √ǫθ β ( ). c ∈N ǫ U Therefore r(s,t,x)pβǫ (s,t,x)φβǫ (0,t,x) √ǫθ βǫ(s,t,x)= F + . − cρ ρ Now, we show the uniqueness of optimal control. Deﬁne

∞ G: L (ΩT ) U ⊂ s 7→ U by r(s,t,x)pβ(s,t,x)φβ(0,t,x) √ǫθ (Gβ)(s,t,x)= F + . − cρ ρ

18 Then for any (s,t,x) ΩT , we have ∈ s (Gβ )(s,t,x) (Gβ )(s,t,x) | 1 − 2 | r(s,t,x)pβ1 (s,t,x)φβ1 (0,t,x) √ǫθ r(s,t,x)pβ2 (s,t,x)φβ2 (0,t,x) √ǫθ = F + F + − cρ ρ − − cρ ρ r(s,t,x) pβ1 (s,t,x)φβ1 (0,t,x) pβ2 (s,t,x)φβ2 (0,t,x) ≤ cρ − r(s,t,x) pβ1 (s,t,x) pβ2 (s,t,x) φβ1 (0,t,x) + pβ2 (s,t,x ) φβ1 (0,t,x) φβ2 (0,t,x) ≤ cρ | − || | | || − |

1 (M M + M M ) β β ∞ , ≤ cρ 1 4 2 3 k 1 − 2kL (ΩTs ) because 1 (M M + M M ) < 1. cρ 1 4 2 3 Hence G is a contraction mapping, so have a unique ﬁxed point β¯ . Optimality condition in theorem ∈ U 5.3 gives us the uniqueness of the optimal control. Now, we prove the existence of optimal controller. Let

Ψ(β¯) = inf Ψ(β) : β . (5.10) { ∈ U} Computing the norm, we obtain

Gβǫ βǫ ∞ k − kL (ΩTs ) r(s,t,x)pβǫ (s,t,x)φβǫ (0,t,x) √ǫθ r(s,t,x)pβǫ (s,t,x)φβǫ (0,t,x) = F + F − cρ ρ − − cρ √ǫ . ≤ ρ Therefore

β¯ βǫ ∞ = Gβ¯ βǫ ∞ k − kL (ΩTs ) k − kL (ΩTs ) Gβ¯ Gβǫ ∞ + Gβǫ βǫ ∞ ≤ k − kL (ΩTs ) k − kL (ΩTs ) 1 √ǫ (M1M4 + M2M3) β¯ βǫ ∞ + . ≤ cρ k − kL (ΩTs ) ρ Which further gives 1 1 −1 β¯ βǫ ∞ 1 (M1M4 + M2M3) √ǫ, k − kL (ΩTs )≤ ρ − cρ ∞ + which implies βǫ β¯ in L (ΩT ) as ǫ 0 . Hence by (5.6), we have → s → Ψ(β¯)= inf Ψ(β). β∈U

19 References

[1] Sharpe, F., Lotka, A.J. A problem in age-distribution. Philos. Mag., 6:435-438, (1911).

[2] Lotka, A.J. The stability of the age distribution. Proc. Natl. Acad. Sci., 8:330-345, (1922).

[3] Webb,G. F. Population models structured by age, size, and spatial position. In Structured population models in biology and epidemiology, volume 1936 of Lecture Notes in Math., pages 1-49. Springer, Berlin, (2008).

[4] Kato, N. Abstract linear partial diﬀerential equations related to size-structured population models with diﬀusion. J. Math. Anal. Appl., 436(2):890–910, (2016).

[5] Anita, S. Analysis and control of age-dependent population dynamics, volume 11 of Mathematical Modelling: Theory and Applications. Kluwer Academic Publishers, Dordrecht, (2000).

[6] Barbu, V., Iannelli, M. Optimal control of population dynamics. J. Optim. Theory Appl., 102(1):1–14, (1999).

[7] Liu, R., Liu, G. Optimal birth control problems for a nonlinear vermin population model with size- structure. J. Math. Anal. Appl., 449(1):265–291, (2017).

[8] He, Z. R., Liu, Y. An optimal birth control problem for a dynamical population model with size- structure. Nonlinear Anal. Real World Appl., 13(3):1369–1378, (2012).

[9] Chan, W. L., Zhu, G. B.. Optimal birth control of population dynamics. J. Math. Anal. Appl., 144(2):532–552, (1989).

[10] Chan, W. L., Guo, B. Z. Optimal birth control of population dynamics. II. Problems with free ﬁnal time, phase constraints, and mini-max costs. J. Math. Anal. Appl., 146(2):523–539, (1990).

[11] Luo, Z. Optimal birth control for an age-dependent competition system of n species. J. Syst. Sci.Complex., 20(3), 403-415., (2007).

[12] Chan, W. L., Zhu, G. B. Pareto Optimal Birth Control of Age-Dependent Populations. IFAC Proc. Vol., 22(4), 147-151., (1989).

[13] Luo, Z. Optimal birth control for competition system of three species with age-structure. J. Appl. Math. Comput., 24(1), 49-64., (2007).

[14] Luo, Z., He, Z. R., Li, W. T. Optimal birth control for predator–prey system of three species with age-structure. Appl. Math. Comput., 155(3), 665-685., (2004).

[15] Zhao, C., Wang, M. S., Zhao, P., Wan, A. H. Optimal birth control for periodic age-dependent population dynamics. In Proceedings of 2004 International Conference on Machine Learning and Cybernetics (IEEE Cat. No. 04EX826)., (Vol. 1, pp. 383-386)., (2004).

[16] He, Z. R., Wang, M. S., Ma, Z. E. Optimal birth control problems for nonlinear age-structured population dynamics. Discrete Cont. Dyn. Sys. ser. B., 4(3), 589., (2004).

20 [17] Huang, Y., Zhao, Y. Optimal Birth Control of A Nonlinear Population Diﬀusion with External Constraint. Control Theory Appl., 05,534–544., (1994).

[18] Kato, N. Optimal harvesting for linear size-structured population models with diﬀusion. J. Nonlinear Convex Anal., 18(7):1335–1348, (2017).

[19] Kato, N. Optimal harvesting for nonlinear size-structured population dynamics. J. Math. Anal. Appl., 342(2), 1388-1398, (2008).

[20] Iannelli, M., Marinoschi, G. Harvesting control for an age-structured population in a multilayered habitat. J. Optim. Theory Appl., 142(1):107–124, (2009).

[21] Getz, W. M. Harvesting discrete nonlinear age and stage structured populations. J. Optim. Theory Appl., 57(1):69–83, (1988).

[22] Brokate, M. Pontryagin’s principle for control problems in age-dependent population dynamics. J. Math. Biol., 23(1), 75-101., (1985).

[23] Morozov, A., Kuzenkov, O. A., Arashkevich, E. G. Modelling optimal behavioural strategies in structured populations using a novel theoretical framework. Scientiﬁc Reports, 9, (2019).

[24] Ekeland, I. On the variational principle. J. Math. Anal. Appl., 47(2), 324-353., (1974).