Arxiv:1704.05697V1 [Math.OC]

Home , Variational principle

arXiv:1704.05697v1 [math.OC] 19 Apr 2017 hsi rpito ae hs nladdfiiefr swit is form definite and final whose paper a Fractals of & preprint a is This oet ar,GogoS aen n efi .M Torres M. F. Delfim and Taverna S. Giorgio Garra, Roberto Derivatives Caputo Generalized with Principles Variational Herglotz Fractional 10.1016/ DOI: 2017; Apr 19 publication for accepted Article ine oko atse l 8 ] eew r neetdi h gen the in interested are Herglotz. we by Here introduced 9]. variations [8, ref of al. are calculus et continuous variations Bastos of of calculus work Th fractional pioneer discrete [1]. the problems order in variational in interested discrete proved formulas for d parts equations fractional by such Euler–Lagrange integration of formalism discrete properties field nonsingula and the 1 a studied where + with [1], 1 derivatives see called fractional function of Leffler so versions the discrete Euler–Lagr within For fractional [5]. obtained the There are [5]. equations in clas Hamilton found the be of can generalization formula, fractional with a Bruno equations of Euler–Lagrange use fractional making by finding fo derivatives, for Hamiltonian method the classic A using addition, investigated ma [25]. In a are mechanics. in derivatives classical obtained in fractional are found with one motion the varia of to equations of similar Hamilton calculus ner the the where of [25], refe in problems can given fractional we fract u to formulation, with approach Lagrangian fields subject fractional Hamiltonian classical the a For using 20]. is by 19, tives, applications [3, research their current and strong principles variational Fractional Introduction 1. lt n13 1]gvsavrainlpicpedsrpino non-cons of description principle variational a gives [14] 1930 in glotz Abstract. Keywords. (2010). Classification effects Subject memory Mathematics arbitrary and oscillator elasticity, harmonic and damped mass a depending consider variati fractional we generalized application, Herglotz-type an on for As depending theorem Lagrangians with Noether-like lems a and tions h eeaie aitoa rnil rtypooe yGsa Her Gustav by proposed firstly principle variational generalized The operators. fractional generalized equations, Lagrange SN 9007.Sbitd1 e 06 ril eie 17 revised Article 2016; Dec 16 Submitted 0960-0779. ISSN: , eoti ue–arneeutos rnvraiycondi transversality equations, Euler–Lagrange obtain We rcinlvrainlpicpe,Hrlt rbe,Eul problem, Herglotz principles, variational Fractional 63;4K5 49S05. 49K05; 26A33; j.chaos.2017.04.035. h derivatives. . nlprob- onal ihtime- with ho,Solitons Chaos, ffrne are ifferences oa deriva- ional re othe to ereed o[] An [7]. to r ia Faásical di Mittag- r oobtain to p 2017; Apr readers e neand ange rmalism er– ervative eralized lfields al Caputo in is tions - nder n- - 2 R. Garra, G. S. Taverna and D. F. M. Torres systems even when the Lagrangian is autonomous [28, 30]. It is essentially based on the following problem: find the trajectories x(t), satisfying given boundary conditions, that extremize (minimize or maximize) the terminal value z(b) of the functional z that satisfies the differential equation z˙(t)= L (t, x(t), x˙(t),z(t)) , t ∈ [a,b], subject to the initial condition z(a)= γ. Herglotz proved that the necessary condition for a trajectory to be an extremizer of the generalized variational problem is to satisfy the generalized Euler–Lagrange equation ∂L d ∂L ∂L ∂L − + =0. ∂x dt ∂x˙ ∂z ∂x˙ The main physical motivation for the development of generalized variational methods behind the classical calculus of variations is linked to the inverse variational problem of classical mechanics in the cases where dissipation is not negligible [18]. In [2], Almeida and Malinowska have considered a fractional variational Herglotz principle, where fractionality stands in the dependence of the Lagrangian by the Caputo fractional derivative of the generalized variables. See also the more recent papers [33, 34] on fractional Herglotz variational principles. In [33], new necessary conditions for higher-order generalized variational problems with time delay, which are semi-invariant under a group of transformations that depends on arbitrary functions, is obtained. Fractional variational problems of Herglotz type of variable order are investigated in [34], where necessary optimality conditions, described by fractional differential equations depending on a combined Caputo fractional derivative of variable order, are proved, both for one and several independent variables [34]. Using such results, it is possible to find, by using a variational approach, the equations of motion of a dissipative mechanical system with memory. This is useful, since in many classical cases, memory effects play a relevant role in systems with dissipation. A relevant example is given by the Basset memory force acting on a sphere rotating in a Stokes fluid. It is well-known that this force can be represented by means of Caputo derivatives of order 1/2 (see, e.g., [6] and references therein). However, a limit in this approach stands to the fact that the particular choice of the dependence of the Lagrangian by the Caputo fractional derivative of the generalized variable, implies that it describes equations of motion of systems with power-law memory kernels. On the other hand, in the framework of the fractional calculus of variations (a quite recent topic of research, started from the seminal investigations of Riewe [26] and then developed by many researchers, see for example the recent monographs [3, 4, 16, 19, 20]), Odzijewicz et al. [21, 22] discusses the case of the Lagrangian depending on generalized Caputo-type derivatives with arbitrary completely monotonic kernels [19]. Here we consider a generalized calculus of variations in the sense of Herglotz with Lagrangians depending on generalized Caputo-type operators treated in [12, 19, 21]. Our aim is to find a general and adequate variational approach to describe mechanical systems with arbitrary memory forces. Fractional Herglotz Variational Principles 3

The paper is organized as follows. In Section 2, we recall the necessary deﬁnitions and results from the generalized fractional variational calculus. Our results are then given in Sections 3, 4 and 5: we prove in Section 3 necessary optimality conditions of Euler–Lagrange type (Theorem 3.1) and transversality conditions (Theorem 3.3) to the generalized fractional variational problem of Herglotz; we show in Section 4 how our approach can deal, in an elegant way, with dissipative dynamical systems with memory eﬀects and time-varying mass and elasticity; and we obtain a generalized fractional Herglotz Noether theorem (Theorem 5.2) in Section 5. We end with Section 6 of conclusions and some directions of future work.

2. Preliminaries In this section, we recall the main deﬁnitions of the generalized Riemann– Liouville and Caputo-like operators and their properties, according to the analysis of generalized fractional variational principles developed in [19, 21]. For a general introduction to fractional diﬀerential operators and equations we refer to the classical encyclopedic book [27]. See also [24]. For an introduction to the fractional variational methods, and in particular integration by parts formulas for fractional integrals and derivatives, we refer to the monographs [4, 20]. For computational and numerical aspects see [3, 11, 15].

α Deﬁnition 2.1. The operator KP is given by α α KP [f](x)= KP [t → f(t)](x) x b = p kα(x, t)f(t)dt + q kα(t, x)f(t)dt, Za Zx where P = ha,x,b,p,qi is the parameter set, x ∈ [a,b], p, q ∈ R, and kα(x, t) is a completely monotonic kernel.

For the sake of completeness, we should remark that similar generalizations of the Riemann–Liouville integrals have been considered in the framework of the fractional action-like variational approach (FALVA) [13, 17]. On the other hand, a similar generalization is considered in a probabilistic framework in [35].

Theorem 2.2 (See [21, Theorem 3]). Let kα ∈ L1([a,b]) and

kα(x, t)= kα(x − t).

α Then, the operator KP : L1([a,b]) → L1([a,b]) is a well-deﬁned bounded and linear operator.

Deﬁnition 2.3. Let P be a given parameter set and α ∈ (0, 1). The operator α 1−α A = D ◦ KP is the generalized Riemann–Liouville derivative, where D stands for the conventional integer order derivative operator. 4 R. Garra, G. S. Taverna and D. F. M. Torres

α The corresponding generalized Caputo derivative is deﬁned as BP = 1−α KP ◦ D. A key-role in the following analysis is played by the following theorem that provides the integration by parts formula for the generalized operators deﬁned before.

Theorem 2.4 (See [21, Theorem 11]). Let α ∈ (0, 1) and P = ha,x,b,p,qi. 1−α ∗ If f,KP g ∈ AC([a,b]), where P = ha,x,b,q,pi and kα(·) is a square- integrable function on ∆ = [a,b] × [a,b], then

b a b α 1−α α g(x)B [f](x)dx = f(x)K ∗ [g](x) − f(x)A ∗ [g](x)dx. Z P P Z P a b a

3. Generalized fractional Herglotz variational principles One of the main aims of this work is to prove generalized Euler–Lagrange equations related to the generalized fractional variational principle of Her- glotz. In particular, the generalization is based on the fact that the La- α grangian depends on the generalized Caputo derivative BP . As explained before, by using this approach, we will be able to ﬁnd a variational approach to mechanical systems involving an arbitrary (suitable) memory kernel kα(t). Therefore, let us consider the diﬀerential equation

α z˙(t)= L (t, x(t),BP [x](t),z(t)) , t ∈ [a,b], (3.1) with the initial condition z(a)= za. We moreover assume that n • x(a)= xa, x(b)= xb, xa, xb ∈ R , • α ∈ (0, 1), 1 Rn α 1 Rn • x ∈ C ([a,b], ), BP [x] ∈ C ([a,b], ), • the Lagrangian L :[a,b] × R2n+1 → R is of class C1 and the maps ∂L t → λ(t) α [x, z](t) exist and are continuous on [a,b], where we use ∂BP xj the notations

α [x, z](t) := (t, x(t),BP [x](t),z(t)), t ∂L λ(t) = exp − [x, z](τ)dτ . Za ∂z The generalized fractional Herglotz variational principle is formulated as follows: Let functional z(t) = z[x; t] be given by the diﬀerential equation (3.1) and η ∈ C1([a,b], R) be an arbitrary function such that η(a)= α 1 R η(b)=0 and BP [η] ∈ C ([a,b], ). Then, the value of the functional z has an extremum for the function x if and only if d z[x + ǫη; b] =0. dǫ ǫ=0

Fractional Herglotz Variational Principles 5

We are now able to state and prove a generalized fractional necessary optimality condition of Euler–Lagrange type that, together with the fractional Herglotz Noether theorem (see Theorem 5.2), constitute the central results of the paper. Theorem 3.1 (Generalized fractional Herglotz Euler–Lagrange equations). Let function x be such that z[x; b] in (3.1) attains an extremum. Then x(t), t ∈ [a,b], is a solution to the generalized Euler–Lagrange equations

∂L α ∂L λ(t) [x, z](t)+ AP ∗ λ(t) α [x, z](t) =0, (3.2) ∂xj ∂BP xj j =1,...,n. Proof. The proof uses the generalized integration by parts formula for Caputo- α like operators BP given by Theorem 2.4. Let x be such that the functional z[x; b] attains an extremum. The rate of change of z in the direction η is given by d θ(t)= z[x + ǫη; t] . dǫ ǫ=0

The variation ǫη of the argument in equation (3.1) is given by d z[x + ǫη; t]= L (t, x(t)+ ǫη(t),Bα [x](t)+ ǫBα [η](t),z[x + ǫη; t]) . dt P P Observing that, from equation (3.1), we have d d θ(t)= L t, x(t)+ ǫη(t),Bα [x](t)+ ǫBα [η](t),z[x + ǫη; t] , dt dǫ P P ǫ=0 this gives a diﬀerential equation of the form n dθ(t) ∂L ∂L ∂L α − [x, z](t)θ(t)= [x, z](t)ηj (t)+ α BP [ηj ](t) , dt ∂z ∂xj ∂B xj Xj=1 P whose solution is t n ∂L ∂L α [x, z](t)ηj (t)+ α BP [ηj ](t) λ(t)dt = θ(t)λ(t) − θ(a), Z ∂xj ∂B xj a Xj=1 P

t ∂L where λ(t) = exp − a ∂z [x, z](τ)dτ and θ(a) = 0. For t = b, we get b n R ∂L ∂L α [x, z](t)ηj (t)+ α BP [ηj ](t) λ(t)dt = θ(b)λ(t). Z ∂xj ∂B xj a Xj=1 P Since θ(b) is the variation of z[x; b], if x gives a maximum, also θ(b) = 0, and therefore we get b n ∂L ∂L α [x, z](t)ηj (t)+ α BP [ηj ](t) λ(t)dt =0. Z ∂xj ∂B xj a Xj=1 P Then, by using the integration by parts formula (see Theorem 2.4), and the fact that η(a) = η(b) = 0, we obtain the claimed result by the fundamental lemma of the calculus of variations. 6 R. Garra, G. S. Taverna and D. F. M. Torres

1 −α Remark 3.2. Let kα(x, t) = Γ(1−α) (x − t) , α ∈ (0, 1), and the parameter set be given by P = ha, x, b, 1, 0i. In this particular case, the operator α C α BP coincides with the standard Caputo fractional derivative a Dx , and our Theorem 3.1 gives the Euler–Lagrange equation of [2]. Observe that, in order to determine uniquely the unknown function that satisfies equation (3.1), the following system of differential equations z˙(t)= L[x, z](t),  ∂L α ∂L (3.3) λ(t) [x, z](t)+ AP ∗ λ(t) α [x, z](t) =0, ∂xj ∂BP xj j = 1,...,n, should be studied with the given boundary conditions. In the case xj (b) is not fixed, similar arguments as those used in the proof of The- orem 3.1 allow us to obtain the generalized fractional integral transversality conditions (3.4). Theorem 3.3 (Generalized fractional Herglotz transversality conditions). Let x be such that z(b)= z[x; b] in equation (3.1) attains an extremum. Then x is a solution to the system (3.3). Moreover, if xj (b) is not fixed, j ∈{1,...,n}, then the integral transversality condition

1−α ∂L KP ∗ t → λ(t) α [x, z](t) (b) = 0 (3.4) ∂BP xj holds.

4. An application: generalized damped harmonic oscillator with memory effects As already discussed in the literature, generalized variational methods can be useful to treat the inverse problem for dissipative systems where memory or damping effects are not negligible. For example, in [21], the authors discussed an application of generalized variational problems to the Caldirola– Kanai approach to quantum dissipative systems; while in [6] an application to the inverse problem for the Basset system was discussed. One can argue that the starting point to research on generalized variational problems was given by the Lemma of Bauer [10], stating that the equations of motionof a classical dissipative system with constant coefficients cannot be derived from a classical variational approach. Here we show the peculiarity of the approach considered in this paper, to treat dissipative dynamical systems with memory effects and time-varying mass and elasticity. Let us consider the mechanical system described by the following autonomous Lagrangian: 1 1 L (x(t),Bα [x](t),z(t)) = m (Bα [x](t))2 − kx2(t)+ λ z(t), (4.1) P 2 P 2 0 where we take a = 0, b> 0, p = 1 and q = 0. By using the generalized Euler– Lagrange equations discussed in Theorem 3.1, we have that the equation of Fractional Herglotz Variational Principles 7 motion in this case is given by

α −λ0t α −λ0t mAP ∗ e (BP [x](t)) − ke x(t)=0. (4.2) Equation (4.2) describes a dynamical oscillatory system with exponentially time-decaying mass and elasticity coefficient and arbitrary memory. The gen- α eralized velocity BP [x](t) can be physically interpreted as a time–weighted velocity, where memory effects are induced by viscosity (and clearly depends by the relaxation kernel kα in the definition of the generalized operators α α BP and AP ∗ ). Therefore, the generalized Herglotz approach here considered provides a variational method to treat the inverse problem of a damped harmonic oscillator with time-depending mass and elasticity (as a consequence of the Herglotz approach) and with arbitrary memory effect implying a velocity delay (due to the dependence of the Lagrangian by a generalized fractional integral). Clearly, in the case in which the memory is neglected, that is, α → 1, the Lagrangian (4.1) takes the form 1 1 L (x(t), x˙(t),z(t)) = mx˙ 2(t) − kx2(t)+ λ z(t) (4.3) 2 2 0

d −λ0t −λ0t and we recover from (4.2) the equation m dt e x˙(t) + ke x(t)=0 of an harmonic oscillator with dissipation. Moreover, if λ0 = 0, then we are in the case in which the Lagrangian (4.3) does not depend on z(t) and we have the classical equation of the harmonic oscillator: if α → 1 and λ0 = 0, then the Euler–Lagrange equation (4.2) reduces to mx¨(t)+ kx(t) = 0.

5. Noether theorem for generalized fractional Herglotz variational problems The analysis of possible generalizations of Noether-type theorems to fractional and Herglotz variational principles, has been subject of recent research: we refer, for example, to [30, 31] and the references therein. Noether-like theorems play indeed a key-role in mathematical-physics by giving the relation between the invariance of the action with respect to some parametric transformation and the existence of conserved quantities [36, 37]. Here we consider a one-parameter family of transformations

x¯j = hj (t,x,s), j =1, . . . , n, (5.1)

2 depending on the parameter s ∈ (−ǫ, +ǫ), with hj ∈ C such that

n hj (t, x, 0) = xj for all (t, x) ∈ [a,b] × R .

The Taylor expansion is given by

hj (t,x,δ)= hj (t, x, 0) + sξj (t, x)+ o(s)

= xj + sξj (t, x)+ o(s), 8 R. Garra, G. S. Taverna and D. F. M. Torres where ξ (t, x)= ∂ h (t,x,s) . In this case the linear approximation of the j s j s=0 transformation (5.1) is simply given by

x¯j (t)= hj (t, x(t),s), j =1, . . . , n. (5.2) Let d θ(t)= z¯[x + sξ,t]| ds s=0 be the total variation produced by the transformation (5.1). Definition 5.1. The transformationx ¯ given by (5.1) leaves the functional z (3.1) invariant if θ(t) ≡ 0. By using the generalized fractional Euler–Lagrange equation (3.2) of Theorem 3.1, we are now able to state the following Noether-type result. Theorem 5.2 (Generalized fractional Herglotz Noether theorem). If the functional z in (3.1) is invariant in the sense of Definition 5.1, then n ˆα ∂L O λ(t) α [x, z](t), ξj (t, x(t)) =0 ∂B xj Xj=1 P holds along the solutions of the generalized Euler–Lagrange equation (3.2), where t ∂L λ(t) = exp − [x, z](τ)dτ Z0 ∂z and ˆα α α O [f,g] := fBP g − gAP ∗ f. Proof. By using the transformation (5.2) into equation (3.1), we get d z¯(t)= L (t, x¯(t),Bα x¯(t), z¯(t)) . dt P Differentiating with respect to s and setting s = 0, we obtain n ˙ ∂L ∂L ∂L α θ(t) − θ(t)= ξj + α BP ξj , (5.3) ∂z ∂xj ∂B xj Xj=1 P where we omit, in order to simplify the notation, that the partial derivatives of L are evaluated at [x, z](t) and ξj at (t, x(t)). The solution of (5.3) is given by t n ∂L ∂L α θ(t)λ(t) − θ(a)= ξj + α BP ξj λ(τ)dτ. (5.4) Z ∂xj ∂B xj a Xj=1 P Since, along the solutions of the generalized Euler–Lagrange equation (3.2), we have that ∂L α ∂L λ(t) = −AP ∗ λ(t) α , (5.5) ∂xj ∂BP xj we obtain the claimed result by substituting (5.5) into (5.4). Fractional Herglotz Variational Principles 9

6. Conclusions We introduced the study of fractional variational problems of Herglotz type that depend on generalized fractional operators in the sense of [19, 23]. As a particular case, one gets a generalization of the Herglotz variational principle for non-conservative systems with Caputo derivatives. Main results give a necessary optimality condition of Euler–Lagrange type (Theorem 3.1), integral transversality conditions (Theorem 3.3), and a Noether-type theorem (Theorem 5.2). Our motivation comes from physics, where such variational principles can be used to describe mechanical systems with memory of arbitrary form. As an application, a fractional mechanical system is analyzed with a fractionally generalized velocity that reproduces, for α = 1, the standard Lagrangian of a harmonic oscillator with exponential damping, which also contains the non-damped conservative oscillator. The research here initiated can now be enriched in diﬀerent directions, by trying to bring to the fractional setting the recent results [28, 29, 30, 31, 32] of Santos et al. on Herglotz variational problems.

Acknowledgments Torres was supported by Portuguese funds through CIDMA and FCT, within project UID/MAT/04106/2013. The authors are grateful to two anonymous referees for valuable comments and suggestions, which helped them to im- prove the quality of the paper.

References [1] T. Abdeljawad and D. Baleanu, Discrete fractional differences with nonsingular discrete Mittag-Leffler kernels, Adv. Difference Equ. 2016 (2016), 232, 18 pp. [2] R. Almeida and A. B. Malinowska, Fractional variational principle of Herglotz, Discrete Contin. Dyn. Syst. Ser. B 19 (2014), no. 8, 2367–2381. [3] R. Almeida, S. Pooseh and D. F. M. Torres, Computational methods in the fractional calculus of variations, Imp. Coll. Press, London, 2015. [4] T. M. Atanacković, S. Pilipović, B. Stankovićand D. Zorica, Fractional calculus with applications in mechanics, Mechanical Engineering and Solid Mechanics Series, ISTE, London, 2014. [5] D. Baleanu, About fractional quantization and fractional variational principles, Commun. Nonlinear Sci. Numer. Simul. 14 (2009), no. 6, 2520–2523. [6] D. Baleanu, R. Garra and I. Petras, A fractional variational approach to the fractional Basset-type equation, Rep. Math. Phys. 72 (2013), no. 1, 57–64. [7] D. Baleanu and S. I. Muslih, Lagrangian formulation of classical fields within Riemann-Liouville fractional derivatives, Phys. Scripta 72 (2005), no. 2-3, 119– 121. [8] N. R. O. Bastos, R. A. C. Ferreira and D. F. M. Torres, Necessary optimality conditions for fractional difference problems of the calculus of variations, Discrete Contin. Dyn. Syst. 29 (2011), no. 2, 417–437. arXiv:1007.0594 10 R. Garra, G. S. Taverna and D. F. M. Torres

[9] N. R. O. Bastos, R. A. C. Ferreira and D. F. M. Torres, Discrete-time fractional variational problems, Signal Process. 91 (2011), no. 3, 513–524. arXiv:1005.0252 [10] P. S. Bauer, Dissipative dynamical systems I, Proc. Nat. Acad. Sci. 17 (1931), 311–314. [11] T. Blaszczyk, M. Ciesielski, M. Klimek and J. Leszczynski, Numerical solution of fractional oscillator equation, Appl. Math. Comput. 218 (2011), no. 6, 2480– 2488. [12] L. Bourdin, T. Odzijewicz and D. F. M. Torres, Existence of minimizers for generalized Lagrangian functionals and a necessary optimality condition— application to fractional variational problems, Differential Integral Equations 27 (2014), no. 7-8, 743–766. arXiv:1403.3937 [13] R. A. El-Nabulsi and D. F. M. Torres, Fractional actionlike variational problems, J. Math. Phys. 49 (2008), no. 5, 053521, 7 pp. arXiv:0804.4500 [14] G. Herglotz, Berührungstransformationen, Lectures at the University of Göttingen, Göttingen, 1930. [15] S. Jahanshahi and D. F. M. Torres, A simple accurate method for solving fractional variational and optimal control problems, J. Optim. Theory Appl., in press. DOI:10.1007/s10957-016-0884-3 arXiv:1601.06416 [16] M. Klimek, On solutions of linear fractional differential equations of a variational type, The Publishing Office of Czestochowa University of Technology, Czestochowa, 2009. [17] C. Liu, Action principle for a kind of mechanical system with friction force, Fizika A 17.1 (2008), 29–34. [18] J. Lopuszanski, The inverse variational problem in classical mechanics, World Sci. Publishing, River Edge, NJ, 1999. [19] A. B. Malinowska, T. Odzijewicz and D. F. M. Torres, Advanced methods in the fractional calculus of variations, Springer Briefs in Applied Sciences and Technology, Springer, Cham, 2015. [20] A. B. Malinowska and D. F. M. Torres, Introduction to the fractional calculus of variations, Imp. Coll. Press, London, 2012. [21] T. Odzijewicz, A. B. Malinowska and D. F. M. Torres, Fractional calculus of variations in terms of a generalized fractional integral with applications to physics, Abstr. Appl. Anal. 2012 (2012), Art. ID 871912, 24 pp. arXiv:1203.1961 [22] T. Odzijewicz, A. B. Malinowska and D. F. M. Torres, A generalized fractional calculus of variations, Control Cybernet. 42 (2013), no. 2, 443–458. arXiv:1304.5282 [23] T. Odzijewicz and D. F. M. Torres, The generalized fractional calculus of variations, Southeast Asian Bull. Math. 38 (2014), no. 1, 93–117. arXiv:1401.7291 [24] I. Podlubny, Fractional differential equations, Mathematics in Science and En- gineering, 198, Academic Press, San Diego, CA, 1999. [25] E. M. Rabei, K. I. Nawafleh, R. S. Hijjawi, S. I. Muslih and D. Baleanu, The Hamilton formalism with fractional derivatives, J. Math. Anal. Appl. 327 (2007), no. 2, 891–897. [26] F. Riewe, Nonconservative Lagrangian and Hamiltonian mechanics, Phys. Rev. E (3) 53 (1996), no. 2, 1890–1899. Fractional Herglotz Variational Principles 11

[27] S. G. Samko, A. A. Kilbas and O. I. Marichev, Fractional integrals and derivatives, translated from the 1987 Russian original, Gordon and Breach, Yverdon, 1993. [28] S. P. S. Santos, N. Martins and D. F. M. Torres, Higher-order variational problems of Herglotz type, Vietnam J. Math. 42 (2014), no. 4, 409–419. arXiv:1309.6518 [29] S. P. S. Santos, N. Martins and D. F. M. Torres, An optimal control approach to Herglotz variational problems. In: Optimization in the Natural Sciences, Com- munications in Computer and Information Science, Vol. 499, 2015, pp. 107–117. arXiv:1412.0433 [30] S. P. S. Santos, N. Martins and D. F. M. Torres, Variational problems of Herglotz type with time delay: DuBois-Reymond condition and Noether’s first theorem, Discrete Contin. Dyn. Syst. 35 (2015), no. 9, 4593–4610. arXiv:1501.04873 [31] S. P. S. Santos, N. Martins and D. F. M. Torres, Noether’s theorem for higher- order variational problems of Herglotz type, Proc. 10th AIMS Conference on Dynamical Systems, Differential Equations and Applications, 2015, 990–999. arXiv:1507.05911 [32] S. P. S. Santos, N. Martins and D. F. M. Torres, Higher-order variational problems of Herglotz type with time delay, Pure and Applied Functional Analysis 1 (2016), no. 2, 291–307. arXiv:1603.04034 [33] S. P. S. Santos, N. Martins and D. F. M. Torres, Noether currents for higher- order variational problems of Herglotz type with time delay, Discrete Contin. Dyn. Syst. Ser. S, in press. Preprint: arXiv:1704.00088 [34] D. Tavares, R. Almeida and D. F. M. Torres, Fractional Herglotz variational problems of variable order, Discrete Contin. Dyn. Syst. Ser. S, in press. Preprint: arXiv:1703.09104 [35] B. Toaldo, Convolution-type derivatives, hitting-times of subordinators and time-changed C0-semigroups, Potential Anal. 42 (2015), no. 1, 115–140. [36] D. F. M. Torres, Gauge symmetries and Noether currents in optimal control, Appl. Math. E-Notes 3 (2003), 49–57. arXiv:math/0301116 [37] D. F. M. Torres, Carathéodory equivalence Noether theorems, and Tonelli full- regularity in the calculus of variations and optimal control, J. Math. Sci. (N. Y.) 120 (2004), no. 1, 1032–1050. arXiv:math/0206230

Roberto Garra Dipartimento di Scienze Statistiche, “Sapienza” Universit`adi Roma, Rome, Italy e-mail: [email protected] Giorgio S. Taverna Institute for Climate and Atmospheric Science, University of Leeds, Leeds, UK e-mail: [email protected] Delﬁm F. M. Torres Center for Research and Development in Mathematics and Applications (CIDMA), Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal e-mail: [email protected]