Designing profitable supply chains for lithium-ion battery recycling in the United States

1.   Introduction

Optimal control (OC) emerged in the 1950s as an extension of the calculus of variations. A significant breakthrough in this field was the development and proof of Pontryagin's maximum principle (PMP) by Pontryagin and his collaborators ^[34]. This principle established the necessary conditions for an optimal solution and has become one of the most powerful tools for addressing OC problems. This formalization of OC theory raised several new questions, particularly in the context of differential equations. It led to the introduction of new generalized solution concepts and produced significant results concerning the existence of trajectories. The applications of OC are extensive, encompassing fields such as mathematics ^[41], physics ^[9], engineering ^[17], robotics ^[11], biology ^[24], and economics ^[39], among others. The primary objective of an OC problem is to identify, from a set of possible solutions (referred to as the admissible set) the one that minimizes or maximizes a given functional. That is, the goal is to determine an OC trajectory and the corresponding state trajectory. The system's dynamics, modeled by state variables, are influenced by control variables that enter the system of equations, thereby affecting its behavior.

Fractional calculus (FC) is a well-known theory that extends the ideas of integration and differentiation to non-integer orders, providing a more flexible approach to modeling complex systems with non-local memory. Numerous publications have been devoted to exploring FC and fractional differential equations (FDEs) from various perspectives, covering both its fundamental definitions and principles (see, e.g., ^[22,38] and the references therein), and its diverse applications in areas such as biology ^[6], medicine ^[23], epidemiology ^[28], electrochemistry ^[31], physics ^[33], mechanics ^[35], and economy ^[37]. Since the begining of fractional calculus in 1695, numerous definitions of fractional derivatives have been introduced, including the Riemann-Liouville, Caputo, Hadamard, Gr$\ddot{\mbox {u}}$nwald-Letnikov, Riesz, Erdélyi-Kober, and Miller-Ross derivatives, among others. Every form of fractional derivative comes with its own benefits and is selected based on the specific needs of the problem being addressed.

In this work, we will explore recent advancements in fractional derivatives as introduced by ^[10]. Specifically, we will focus on the novel concepts of distributed-order fractional derivatives with respect to an arbitrary kernel in the Riemann-Liouville and Caputo senses. In their study, the authors established necessary and sufficient conditions within the context of the calculus of variations and derived an associated Euler–Lagrange equation. Building upon their results, we extend the analysis to the setting of optimal control theory and derive a corresponding generalization of the Pontryagin maximum principle.

A key advantage of distributed-order fractional derivatives ^[8] lies in their ability to describe systems where the memory effect varies over time. Unlike classical or constant-order fractional derivatives, distributed-order derivatives are defined through an integral over a range of orders, weighted by a given distribution function. This provides a more flexible and realistic modeling tool for systems whose dynamic behavior cannot be adequately captured by a single, fixed fractional order. Distributed-order fractional derivatives have proven especially effective in modeling systems characterized by a broad spectrum of dynamic behaviors. This includes, for instance, viscoelastic materials exhibiting a wide range of relaxation times, transport phenomena influenced by multiple temporal and spatial scales, and control systems affected by diverse time delays. Notable applications span across fields such as viscoelasticity, anomalous diffusion, wave propagation, and the design of fractional-order Proportional-Integral-Derivative (PID) controllers ^[13,18].

In parallel, the concept of fractional derivatives with respect to arbitrary kernels ^[3,38] further improves the model. These operators generalize classical definitions (such as the Riemann–Liouville or Caputo derivatives) by introducing a kernel function that dictates the memory behavior of the system. By selecting an appropriate kernel, one can recover various classical forms or construct entirely new operators. This kernel-based approach allows for incorporating diverse memory effects, singularities, or fading influence in time, which are often observed in real-world processes but difficult to model with standard tools.

These two features—variable order and general kernels—collectively enable the modeling of complex phenomena with nonlocal and history-dependent behavior, particularly in contexts where the memory characteristics are not uniform or stationary. This is especially relevant in systems governed by internal friction, relaxation, hereditary stress-strain relations, or anomalous transport, where the effects of the past on the present are not constant and may depend on the nature of past events in intricate ways.

Building on these two ideas, we consider OC problems that involve a generalized fractional derivative constraint. OC provides a rigorous mathematical framework for determining the best possible strategy to steer a dynamical system from an initial state to a desired final state while minimizing (or maximizing) a given cost functional. The inclusion of fractional operators—particularly of distributed-order or arbitrary-kernel types—within this framework enhances the ability to model systems with complex dynamics and long-range temporal dependencies. It enables decision-makers to formulate strategies that are optimal with respect to cost, energy, time, or other performance metrics, all while satisfying dynamic constraints such as differential equations, control bounds, and state limitations. This is crucial for achieving efficiency, such as reducing fuel consumption in aerospace trajectories, minimizing energy use in industrial processes, or shortening recovery time in medical treatments.

Pioneering works, such as those in ^[1,2], developed necessary conditions for optimality with respect to the classical Caputo fractional derivative. These studies were further expanded in ^[14,15] with the formulation of the fractional Noether's theorem. Since then, numerous studies have emerged on fractional OC problems, addressing various fractional operators such as the Caputo ^[5,20], Riemann–Liouville ^[19], and distributed-order derivatives ^[30], as well as aspects like delays in systems ^[16] and fuzzy theory ^[32]. Numerical methods for solving fractional optimal control problems are available in the literature. For example, ^[21] uses a generalized differential transform method, ^[25,42] employ second-order and third-order numerical integration methods, and ^[27] applies a Legendre orthonormal polynomial basis.

The primary objective of this paper is to establish a generalized fractional PMP applicable to OC problems involving the left-sided distributed-order Caputo fractional derivative with respect to an arbitrary smooth kernel. Additionally, we will present sufficient conditions for optimality based on PMP.

The rest of the paper is structured as follows. Section 2 introduces some fundamental concepts and results from fractional calculus required for our work. In Section 3, we present the OC problem $(\mathcal{P}_{\mathcal{OC}})$ under study, some fundamental lemmas, PMP, and sufficient conditions for the global optimality of the problem $(\mathcal{P}_{\mathcal{OC}})$. Section 4 provides two examples that illustrate the practical relevance of our research. Finally, in Section 5, we summarize our findings and present ideas for future work.

2.   Preliminaries of fractional calculus

This section provides essential definitions and results with respect to the distributed-order Riemann-Liouville and Caputo derivatives depending on a given kernel ^[10]. We assume that readers are already acquainted with the definitions and properties of the classical fractional operators ^[22,38].

Throughout the paper, we suppose that $a < b$ are two real numbers, $J = [a, b]$, $\gamma\in \mathbb{R}^+$, and $[\gamma]$ denotes the integer part of $\gamma$. First, we introduce essential concepts relevant to our work.

Definition 2.1. ^[38] (Fractional integrals in the Riemann–Liouville sense) Let $z\in L^1(J, \mathbb{R})$, and let $\sigma \in C^1(J, \mathbb{R})$ be a continuously differentiable function with a positive derivative.

For $t > a$, the left-sided $\sigma$-R-L fractional integral of the function $z$ of order $\gamma$ is given by

and for $t < b$, the right-sided $\sigma$-R-L fractional integral is

Definition 2.2. ^[38] (Fractional derivatives in the Riemann–Liouville sense) Let $z\in L^1(J, \mathbb{R})$, and let $\sigma \in C^n(J, \mathbb{R})$ be such that $\sigma '(t) > 0$, given $t \in J$.

The left-sided $\sigma$-R-L fractional derivative of the function $z$ of order $\gamma$ is given by

and the right-sided $\sigma$-R-L fractional derivative is

where $n = [ \gamma ] + 1$.

Definition 2.3. ^[3] (Fractional derivatives in the Caputo sense) For a fixed $\gamma \in \mathbb{R}^+$, define $n \in \mathbb{N}$ as: $n = [\gamma ] + 1$ if $\gamma \notin \mathbb{N}$, and $n = \gamma$ if $\gamma \in \mathbb{N}$. Moreover, consider $z, \sigma \in C^n(J, \mathbb{R})$ where $\sigma '(t) > 0$, $t \in J$. The left- and right-sided $\sigma$-C fractional derivatives of the function $z$ of order $\gamma$ are given by

and

respectively.

We remark that, for $0 < \gamma < 1$:

and

If $\gamma = 1$, we get

and

Next, we recall some properties that will be useful in our proofs. For $\gamma\in(0, 1]$ and $z \in C^1(J, \mathbb{R})$, the following relations hold:

and

For further details, we refer the reader to ^[3].

Since the fractional OC problem examined in this paper involves a fractional derivative of order $\gamma$ within the range (0, 1], we will henceforth assume $\gamma \in (0, 1]$.

To define the distributed-order derivatives, we need to fix the order-weighting function, denoted by $\Phi$. Here, $\Phi$ is a continuous function defined on $[0, 1]$ such that $\Phi([0, 1]) \subseteq [0, 1]$ and $\int_0^1 \Phi(\gamma) d\gamma > 0$.

Definition 2.4. ^[10] (Distributional-order fractional derivatives in the Riemann–Liouville sense) Let $z\in L^1(J, \mathbb{R})$. The left- and right-sided $\sigma$-D-R-L fractional derivatives of a function $z$ with respect to the distribution $\Phi$ are defined by

where $\mathrm{D}^{\gamma, \sigma }_{a^+}$ and $\mathrm{D}^{\gamma, \sigma }_{b^-}$ are the left- and right-sided $\sigma$-R-L fractional derivatives of order $\gamma$, respectively.

Definition 2.5. ^[10] (Distributional-order fractional derivatives in the Caputo sense) The left- and right-sided $\sigma$-D-C fractional derivatives of a function $z \in C^1(J, \mathbb{R})$ with respect to the distribution $\Phi$ are defined by

where $^C\mathrm{D}^{\gamma, \sigma }_{a^+}$ and $^C\mathrm{D}^{\gamma, \sigma }_{b^-}$ are the left- and right-sided $\sigma$-C fractional derivatives of order $\gamma$, respectively.

In what follows, we introduce the concepts of $\sigma$-distributional-order fractional integrals:

where $\mathrm{I}^{1-\gamma, \sigma }_{a^+}$ and $\mathrm{I}^{1-\gamma, \sigma }_{b^-}$ are the left- and right-sided $\sigma$-R-L fractional integrals of order $1-\gamma$, respectively.

It is evident from the definitions that distributed-order operators are linear.

In the following, we present a generalized fractional integration by parts formula which is useful to demonstrate some of our results.

Lemma 2.6. (Generalized fractional integration by parts) ^[10] Given $w\in C(J, \mathbb R)$ and $z\in C^1(J, \mathbb R)$, the following holds:

We now present a general form of Gronwall's inequality, a crucial tool for comparing solutions of FDEs that involve $\sigma$-fractional derivatives.

Lemma 2.7. (Generalized fractional Gronwall inequality) ^[40] Let $u, v\in L^1(J, \mathbb R)$, $\sigma \in C^1(J, \mathbb{R})$ with $\sigma '(t) > 0$ for all $t\in J$, and $h\in C(J, \mathbb R)$. Assume also that $u, v, h$ are nonnegative and $h$ is nondecreasing. If

then

Remark 2.8. In Lemma 2.7, the assumption that $h$ is nondecreasing is often satisfied in practical contexts where $h$ represents quantities such as cumulative costs, energy consumption, or memory effects that naturally increase or remain constant over time.

Next, we recall the definition of the Mittag-Leffler function (with one parameter), which generalizes the standard exponential function and plays a fundamental role in the study of differential equations of fractional order.

Definition 2.9. The Mittag-Leffler function with parameter $\gamma > 0$ is defined by

To conclude this section, we present the following result, which is a corollary of Lemma 2.7 and will be useful in the proof of Lemma 3.3.

Corollary 2.10. (cf. ^[40]) Under the hypotheses of Lemma 2.7, if $v$ is a nondecreasing function, then, for all $t\in J$,

3.   Main results

The aim of this section is to establish necessary and sufficient optimality conditions for a fractional OC problem involving the left-sided Caputo distributed-order fractional derivative with respect to a smooth kernel $\sigma$.

In what follows, we use the standard notations: $PC(J, \mathbb{R})$ denotes the set of all real-valued piecewise continuous functions defined on $J$, and $PC^1(J, \mathbb{R})$ denotes the set of piecewise smooth functions defined on $J$. The fractional OC problem is defined by the following formulation:

Problem $(\mathcal{P_{OC}})$: Determine $x\in PC^1(J, \mathbb{R})$ and $u\in PC(J, \mathbb{R})$ that extremizes the functional

under the following restrictions: the FDE

and the initial condition

where $L \in C^1(J \times\mathbb{R}^2, \mathbb{R})$ and $f\in C^2(J \times\mathbb{R}^2, \mathbb{R}).$

Remark 3.1. There are several equivalent ways to formulate an optimal control problem, notably in the Mayer, Lagrange, and Bolza forms. Through appropriate changes of variables and auxiliary state transformations, one can show that these formulations are theoretically equivalent (see, e.g., Chapter 3 in ^[26]). In this paper, we adopt the Lagrange form for consistency and simplicity.

We will now prove a result that establishes a relationship between an optimal state trajectory of the OC problem $(\mathcal{P_{OC}})$ and the state solution of a perturbed version of $(\mathcal{P_{OC}})$.

Lemma 3.2. Let the pair $(\overline{u}, \overline{x})$ be an optimal solution to $(\mathcal{P_{OC}})$. For $t\in J$, consider a variation of $\overline{u}$ of the form $\overline{u}+\xi\eta$, where $\eta\in PC(J, \mathbb{R})$ and $\xi\in\mathbb{R}$. Denote by $u^\xi(t)$ such variation and $x^\xi$ the solution of

Then, as $\xi$ tends to zero, $x^\xi\to \overline{x}$.

Proof. By hypothesis, we know that

where $x^\xi(a) = \overline{x}(a) = x_a$. From the linearity and the definition of the operator $^C\mathrm{D}^{ \Phi(\gamma), \sigma }_{a^+}$, we may conclude that

Using the mean value theorem in the integral form, we conclude that there exists $\beta\in[0, 1]$ such that

Denoting $\omega: = \int_0^1 \Phi(\gamma)d\gamma$, we conclude that

Thus,

Let $k_1$ and $k_2$ be two positive real numbers such that

and denote by $k$ their maximum. Therefore,

Since $u^\xi(t)-\overline{u}(t) = \xi \eta(t)$, we have

for all $t\in J$.

Now, applying Lemma 2.7, we obtain that

Thus, by applying the mean value theorem once more, we deduce that there exists some $\overline\tau\in J$ such that

So, the proof is finished by taking the limit when $\xi\to0$. □

Before presenting our next result, we extend the one presented in [12, Theorem 3.4] for FDEs with the $\sigma$-Caputo fractional derivative.

Lemma 3.3. Consider the following two FDEs of order $\gamma\in(0, 1]$:

with the same initial condition. Suppose the functions $f$ and $\tilde f$ are Lipschitz continuous with respect to $x$, with Lipschitz constant $L > 0$. If $y$ and $z$ are the unique solutions to the first and second equations, respectively, then there exists a constant $C > 0$, independent of $y$ and $z$, such that

Proof. Given $t\in J$, we have:

Applying Corollary 2.10, with

we obtain

Thus, we conclude that

which proves the desired result. □

Throughout the following, we use the standard Big-$O$ notation: $a(x) = O(b(x))$ means that there exists a constant $C > 0$ such that $|a(x)|\leq C\, |b(x)|$ for all $x$ sufficiently close to 0.

Lemma 3.4. Suppose we are in the conditions of Lemma 3.2. Then, there exists a mapping $\nu:J \to\mathbb R$ such that

Proof. Given that $f\in C^2$, the function $f$ can be expanded as

For all $t\in J$, we have $u^{\xi}(t)-\overline{u}(t) = \xi \eta(t)$. Additionally, as proven in Lemma 3.2, $x^{\xi}(t)-\overline{x}(t) = \xi \chi(t)$, for some finite function $\chi$. Thus, the term $O(|x^{\xi}(t)-\overline{x}(t)|^2, |u^{\xi}(t)-\overline{u}(t)|^2)$ simplifies to $O(\xi^2)$. Hence, we can express the fractional derivative of $x^\xi$ as:

Therefore, for $\xi\neq0$ and $t\in J$,

Consider the following two FDEs:

and

subject to $y(a) = 0$. The existence and uniqueness of solutions for distributed-order FDEs can be ensured using standard results for FDEs involving the $\sigma$-Caputo fractional derivative (see, e.g., ^[4]). This is achieved by using the relation

for some $\overline\gamma\in [0, 1]$ and with $\omega = \int_0^1 \Phi(\gamma)\, d\gamma$. Equation (3.1) has the solution $y_1: = \frac{x^\xi-\overline{x}}{\xi}$, and let $y_2$ be the solution of Eq (3.2). By Lemma 3.3, we obtain that $y_2(t) = \lim_{\xi\to0} y_1(t)$, proving the desired result with $\nu: = y_2$. □

We are now ready to present the main result of our work:

Theorem 3.5. (PMP for $(\mathcal{P_{OC}})$) Let $(\overline{x}, \overline{u})$ be an optimal pair to problem $(\mathcal{P_{OC}})$. Then, there exists a mapping $\lambda\in PC^1(J, \mathbb{R})$ such that:

● The two following FDEs hold:

and

● The transversality condition

Proof. Suppose that $(\overline{x}, \overline{u})$ is a solution to problem $(\mathcal{P_{OC}})$. Consider a variation of $\overline{u}$ of the form $\overline{u}+\xi\eta$, denoted by $u^{\xi}$, where $\eta\in PC(J, \mathbb{R})$ and $\xi\in\mathbb R$. Also, let $x^{\xi}$ be the state variable satisfying

Observe that, as $\xi$ goes to zero, $u^{\xi}$ goes to $\overline{u}$ on $J$, and that

By Lemma 3.2 we can conclude that $x^{\xi}\to \overline{x}$ on $J$ when $\xi\to 0$. Furthermore, by Lemma 3.4, the partial derivative $\frac{\partial x^{\xi}(t)}{\partial\xi}\Big|_{\xi = 0} = \nu(t)$ exists for each $t$. We are considering the following functional

Let $\lambda\in PC^1(J, \mathbb{R})$, to be explained later. From Lemma 2.6, we conclude that

Then,

Since $(\overline{x}, \overline{u})$ is a solution to problem $(\mathcal{P_{OC}})$, we conclude that

where the partial derivatives of $L$ and $f$ are evaluated at the point $(t, \overline{x}(t), \overline{u}(t))$. Now replacing (3.7), we get

Rearranging the terms we have

Introducing the Hamiltonian function

the expression simplifies to:

where $\partial H/\partial x$ and $\partial H/\partial u$ are evaluated at $(t, \overline{x}(t), \overline{u}(t), \lambda(t))$. Choosing the function $\lambda$ as the solution of the system

we get

Since $\eta$ is arbitrary, from the Du Bois-Reymond lemma (see, e.g., ^[7]), we obtain

for all $t\in J$, proving the desired result. □

Definition 3.6. The Hamiltonian is defined as $H(t, x(t), u(t), \lambda(t)) = L(t, x(t), u(t)) + \lambda(t) f(t, x(t), u(t)),$ where $L$ is the Lagrange function, $f$ is the system dynamics, and $\lambda$ is the adjoint variable.

Remark 3.7. We observe that:

$(1)$ Taking $\sigma$ as the identity function, we recover the results presented in ^[29], while also correcting minor inaccuracies in their proofs. This allows us to refine and extend the results previously established for the classical distributed-order fractional derivative. Furthermore, by introducing a general kernel, we not only broaden the scope of the original framework but also provide a more general formulation that encompasses and extends the findings of the aforementioned study.

$(2)$ If $f(t, x(t), u(t)) = u(t), \, t\in J$, our problem $(\mathcal{P_{OC}})$ reduces to the following fractional problem of the calculus of variations:

subject to the initial condition $x(a) = x_a$, where $L\in C^1$. From Theorem 3.5 we deduce that if $\overline{x}(t)$ is an extremizer, then there exists $\lambda\in PC^1(J, \mathbb{R})$ such that:

● $\lambda(t) = - \partial_3 L (t, \overline{x}(t), ^C\mathrm{D}^{ \Phi(\gamma), \sigma }_{a^+} \overline{x}(t)), \quad t\in J,$

● $\left(\mathrm{D}^{ \Phi(\gamma), \sigma }_{b^-}\frac{\lambda(t)}{\sigma '(t)}\right)\sigma '(t) = \partial_2 L (t, \overline{x}(t), ^C\mathrm{D}^{ \Phi(\gamma), \sigma }_{a^+} \overline{x}(t)), \quad t\in J,$

● $I^{1- \Phi(\gamma), \sigma }_{b^-}\frac{\lambda(b)}{\sigma '(b)} = 0,$

where $\partial_i L$ denotes the partial derivative of $L$ with respect to its ith coordinate. Hence, we obtain the Euler-Lagrange equation:

and the tranversality condition:

proved in ^[10].

We conclude this section by establishing sufficient optimality conditions for our OC problem.

Definition 3.8. We say that a triple $(\overline{x}, \overline{u}, \lambda)$ is a Pontryagin extremal of Problem $(\mathcal{P_{OC}})$ if it satisfies conditions (3.3)–(3.5).

Theorem 3.9. (Sufficient conditions for global optimality I) Let $(\overline{x}, \overline{u}, \lambda)$ be a Pontryagin extremal of Problem $(\mathcal{P_{OC}})$ with $\lambda(t)\geq0$, for all $t\in J$.

● If $L$ and $f$ are convex functions, then $(\overline{x}, \overline{u})$ is a global minimizer of functional $\mathcal{J}$;

● If $L$ and $f$ are concave functions, then $(\overline{x}, \overline{u})$ is a global maximizer of functional $\mathcal{J}$.

Proof. We will only prove the case where $L$ and $f$ are convex; the other case is analogous. If $L$ is a convex function, then (see, e.g., ^[36])

for any control $u$ and associate state $x$. Hence,

Using the adjoint equation (3.4) and the optimality condition (3.3), we obtain

Rearranging the terms, we have

Using the generalized integration by parts formula (Lemma 2.6) in the first integral, we obtain

Since

and

we get

Since $f$ is convex, we have

for any admissible control $u$ and its associated state $x$. Moreover, since $\lambda(t)\geq0$, for all $t\in J$, it follows that $\mathcal{J}(\overline{x}, \overline{u})-\mathcal{J}(x, u) \leq 0$, proving the desired result. □ Similarly, we can prove the following result.

Theorem 3.10. (Sufficient conditions for global optimality II) Let $(\overline{x}, \overline{u}, \lambda)$ be a Pontryagin extremal of Problem $(\mathcal{P_{OC}})$ with $\lambda(t) < 0$, for all $t\in J$.

● If $L$ and $f$ are convex functions, then $(\overline{x}, \overline{u})$ is a global maximizer of functional $\mathcal{J}$;

● If $L$ and $f$ are concave functions, then $(\overline{x}, \overline{u})$ is a global minimizer of functional $\mathcal{J}$.

4.   Illustrative examples

We now present two examples in order to illustrate our main result.

Example 4.1. Consider the following problem:

subject to

The order-weighting function is $\phi:[0, 1]\rightarrow [0, 1]$ defined by

The Hamiltonian in this case is given by

From $(3.3)-(3.5)$, we get

and

Note that the triple $(\overline{x}, \overline{u}, \lambda)$ given by:

satisfies the necessary optimality conditions given by the Pontryagin maximum principle. Since the Lagrangian $L$ and $f(t, x(t), u(t)) = u(t)$ are convex, then $(\overline{x}, \overline{u})$ is a global minimizer of functional $\mathcal{J}$.

Example 4.2. Consider the optimal control problem:

subject to the restriction

The corresponding Hamiltonian is given by

From the stationarity condition (3.3), we obtain

The remaining necessary conditions, as given by (3.4) and (3.5), are

and

5.   Conclusions and future work

In this work, we explored necessary and sufficient conditions for optimality in fractional OC problems involving a novel fractional operator. Specifically, we addressed the challenge of determining the optimal processes for a fractional OC problem, where the regularity of the optimal solution is assured. The solution to this question is provided by an extension of the PMP, for which we have established the basic version applicable to fractional OC problems, involving a new generalized fractional derivative with respect to a smooth kernel. Additionally, we presented sufficient conditions of optimality for this class of problems using the newly defined fractional derivative. Future research will aim to extend these theorems to OC problems with bounded control constraints. Additionally, it is crucial to develop numerical methods for determining the optimal pair. One potential approach is to approximate the fractional operators using a finite sum, thereby transforming the fractional problem into a finite-dimensional problem, which can then be more effectively solved numerically.

Author contributions

Fátima Cruz, Ricardo Almeida and Natália Martins: Formal analysis, investigation, methodology, and writing – review $&$ editing, of this manuscript. All authors have read and approved the final version of the manuscript for publication.

Use of Generative-AI tools declaration

The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work is supported by Portuguese funds through the CIDMA - Center for Research and Development in Mathematics and Applications, and the Portuguese Foundation for Science and Technology (FCT-Fundação para a Ciência e a Tecnologia), within project UIDB/04106/2025.

Conflict of interest

Professor Ricardo Almeida is an editorial board member for AIMS Mathematics and was not involved in the editorial review and/or the decision to publish this article.

All authors declare no conflicts of interest.