Solving functional integrodifferential equations with Liouville-Caputo fractional derivatives by fixed point techniques

Abbreviations

FDE, fractional differential equation

LC, Liouville-Caputo

RL, Riemann-Liouville

BS, Banach space

FFIDE, fractional functional integrodifferential equation

BVP, boundary value problems

1.   Introduction

The theory of fractional differential equations (FDEs) generalizes classical differential equations by introducing fractional derivatives, enabling the modeling of complex phenomena exhibiting non-locality, memory, and power-law behavior ^[1,2]. FDEs have been extensively developed using various fractional derivatives, such as Riemann-Liouville (RL), Caputo, and Grü nwald-Letnikov ^[3]. These equations describe anomalous diffusion, relaxation, and oscillation processes, making them suitable for modeling real-world problems in physics, engineering, and biology. FDEs have diverse applications across disciplines, including viscoelasticity ^[4], chaotic dynamics ^[5], image processing ^[6], model financial systems ^[7], population dynamics ^[8], and electrical circuits ^[9]. Recent studies have explored fractional-order controllers in robotics ^[10] and biomedical signal processing^[11]. These applications demonstrate the versatility and potential of FDEs in describing complex systems.

Fractional calculus and fixed point theory have emerged as powerful tools in addressing optimization and inverse problems across various scientific and engineering disciplines. Fractional derivatives and integrals, with their inherent nonlocal properties, are particularly well-suited for modeling systems exhibiting memory effects and long-range dependencies, which are often encountered in optimization problems involving complex systems. Furthermore, fixed point (FP) methods, especially those tailored for non-smooth or set-valued mappings, provide robust frameworks for solving inverse problems, including those arising in image processing and signal reconstruction. The combination of fractional calculus and FP theory offers a synergistic approach, enabling the analysis and solution of challenging optimization and inverse problems that are often intractable using classical techniques. See ^{[12,13,14,15,16]} for more information.

Recent years have witnessed significant interest in boundary value problems (BVPs) of FDEs, encompassing various boundary conditions such as the existence and uniqueness of solutions to fractional boundary value problems ^{[17,18,19,20]}, the existence and uniqueness of solutions to hybrid fractional systems under multi-point, periodic, and anti-periodic boundary conditions ^[21,22], and the stability of mixed integral fractional delay dynamic system equations and pantograph differential equations under impulsive effects and nonlocal conditions ^[23,24]. Integral boundary conditions, in particular, have far-reaching implications in applied fields like heat conduction, electric power networks, elastic stability, telecommunications and electric railway systems. Multi-point BVPs, arising from practical applications, also warrant attention. For example, the existence results for FDEs are established in ^{[25,26,27,28,29]}. Alos, the existence of solutions to fractional functional differential equations ^[30], semilinear fractional differential inclusions ^[31], Hadamard fractional integro-differential equations ^[32], systems of multi-point boundary value problems ^[33], fractional hybrid delay differential equations ^[34], and nonlinear Atangana-Baleanu-type fractional differential equations ^{[35,36,37,38]} have been established. The theory of fractional functional BVPs remains underdeveloped, necessitating further research in mathematical modeling, numerical methods, and computational simulations to address the unresolved aspects.

Benchohra et al. ^[26] proved the existence of a solution via Leray-Schauder nonlinear alternative and uniqueness via Banach's FP theorem for fractional functional differential equations with infinite delay. Chauhan et al. ^[20] explored existence solutions for fractional integro-differential equations with impulses, infinite delay, and integral boundary conditions. Dabas and Gautam ^[39] examined existence results for impulsive neutral fractional integrodifferential equations featuring state-dependent delays and integral boundary conditions.

Inspired by the contributions of ^[20,26,39], first, our study focuses on establishing existence and uniqueness results for a fractional functional integrodifferential equation (FFIDE) featuring infinite delay. It takes the form

where $^{LC}D^{p}$ is the Liouville-Caputo (LC) fractional derivative with order $p.$ Assume that $\mho = \left\{ (\varsigma, \vartheta):0\leq \vartheta \leq \varsigma \leq \sigma \right\}$, $\Xi$ is a BS, and $\Theta$ is a phase space. Then $g:U\times \Theta \times \Xi \rightarrow \Xi,$ $\Omega, \Upsilon :\mho \times \Theta \rightarrow \Xi$ are continuous functions and $\psi \in \Theta.$ Furthermore, $I_{0^{+}}^{q_{j}}$ refers to the RL fractional integral of order $q_{j} > 0,$ and $b_{j}$ represents suitable real constants for $j = 1, 2, \dots, u.$

Supposing that $\varpi :(-\infty, \sigma ]\rightarrow \Theta$ and $\varsigma \in U$, we denote $\varpi _{\varsigma }\in \Theta$ as an element defined by

Throughout this manuscript, we suppose that $\varpi _{\varsigma }(.)$ is the historical state trajectory from time $-\infty$ to $\varsigma$, and $\varpi _{\varsigma }\in \Theta$, where $\Theta$ is an abstract phase space.

The second main result here is to investigate the existence and uniqueness of solutions to the neutral FFIDE with BVPs. It takes the form

where $h, g:U\times \Theta \times \Xi \rightarrow \Xi,$ $\Omega _{1}, \Upsilon _{1}, \Omega _{2}$, and $\Upsilon _{2}:\mho \times \Theta \rightarrow \Xi$ are continuous functions.

● This paper provides a systematic exploration of fractional functional integrodifferential equations.

● Section 2 lays the groundwork by establishing the foundational definitions, notation, and preliminary results.

● Existence and uniqueness criteria for FFIDEs are developed in Section 3, employing both Krasnoselskii's FP and Banach's FP theorem.

● Building upon these findings, Section 4 extends the existence and uniqueness results to neutral FFIDEs with BVPs.

● The applicability and practicality of the theoretical framework are demonstrated through illustrative examples provided in Section 5.

2.   Preliminary work

This section presents fundamental definitions, notation, and lemmas essential for the subsequent analysis. Let $\Xi$ denote a BS equipped with the norm $\left\Vert.\right\Vert$. Furthermore, $C(U, \Xi)$ represents the BS of continuous functions from the interval $U$ to $\Xi$, endowed with a uniform convergence topology and the norm $\left\Vert.\right\Vert _{C}$.

Definition 2.1. ^[2] For the function $g\in L^{1}(\mathbb{R} _{+})$

(ⅰ) The RL fractional integral of order $p > 0$ is given by

whenever the integral exists.

(ⅱ) The LC fractional derivative of order $p\in (v-1, v]$ is described as

where $g$ has absolutely continuous derivatives up to order $(v-1)$.

Remark 2.2. It should be noted that, if we take $v = 1$ in Definition 2.1 (ⅱ), we have $0 < p\leq 1$ and

where $g^{^{\prime }}(\vartheta) = \frac{dg(\vartheta)}{d\vartheta }.$

Now, for simplicity, we denote $^{LC}D_{\varsigma }^{p}$ and $I_{0^{+}}^{p}$ by $^{LC}D^{p}$ and $I^{p},$ respectively.

Lemma 2.3. ^[2] Assume that $p, q\geq 0$, and $g\in L^{1}[b, c].$ Then, for all $\varsigma \in \lbrack b, c],$ we have

(i) $I^{q}I^{p}g(\varsigma) = I^{q+p}g(\varsigma) = I^{p}I^{q}g(\varsigma);$

(ii) $^{LC}D_{\varsigma }^{p}I^{q}g(\varsigma) = g(\varsigma).$

Theorem 2.4. ^[40] (Krasnoselskii's theorem) Assume that $\Lambda \neq \emptyset$ is a closed and convex subset of a BS $\Xi$ and that $\Im, \Re$ are two operators satisfying

(i) for $\varpi, \varrho \in \Lambda,$ $\Im \varpi +\Re \varrho \in \Lambda,$

(ii) $\Im$ is continuous and compact,

(iii) $\Re$ is a contraction.

Then, $w\in \Lambda$ exists such that $w = \Im w+\Re w.$

This paper considers a seminormed linear state space $(\Theta, \left\Vert.\right\Vert _{\Theta })$ of functions from $(-\infty, 0]$ to $\Xi$ satisfying the following hypotheses of Hale and Kato ^[41]:

(H$_{1}$) On the interval $(-\infty, \sigma ],$ if $\varpi :(-\infty, \sigma ]\rightarrow \Xi$ is continuous and $\varpi _{0}\in \Theta$, then for $\varsigma \in U,$ we have the following stipulations:

(1) $\varpi _{\varsigma }\in \Theta,$

(2) $\left\Vert \varpi (\varsigma)\right\Vert _{\Theta }\leq \kappa \left\Vert \varpi _{\varsigma }\right\Vert _{\Theta },$ where $\kappa$ is a non-negative constant and is independent of $\varpi (.),$

(3) There is a continuous function $N_{1}:[0, \infty)\rightarrow \lbrack 0, \infty)$ and a locally bounded function $N_{2}:[0, \infty)\rightarrow \lbrack 0, \infty)$ in order that

where $N_{1}$ and $N_{2}$ are independent of $\varpi (.).$

(H$_{2}$) The space $\Theta$ is complete.

(H$_{3}$) On the interval $U,$ $\varpi _{\varsigma }$ is a B-valued continuous function, where $\varpi (.)$ is described in (H$_{1}$).

Here, we consider $N_{1}^{\ast } = \sup_{\varsigma \in U}N_{1}(\varsigma)$ and $N_{2}^{\ast } = \sup_{\varsigma \in U}N_{2}(\varsigma).$

3.   Solving the FFIDE

This section is devoted to investigating the existence and uniqueness of solution to the considered problem (1.1) by applying Krasnoselskii's and Banach's FP theorems.

Assume the space

and select $P\varpi (\varsigma) = \int_{0}^{\varsigma }\Omega \left(\varsigma, \vartheta, \varpi _{\vartheta }\right) d\vartheta,$ and $Q\varpi (\varsigma) = \int_{0}^{\sigma }\Upsilon \left(\varsigma, \vartheta, \varpi _{\vartheta }\right) d\vartheta.$

Definition 3.1. We say that the function $\varpi \in \widetilde{\mho }$ is a solution to the FFIDE (1.1) if it fulfills the problem

We initiate our analysis of the nonlinear problem (1.1) by examining its linear counterpart, thereby obtaining a foundational solution.

Lemma 3.2. Assume that $\varpi (\varsigma)\in C\left(U, \Xi \right) \ $ satisfies the following problem:

Then the unique solution of the fractional BVP (3.1) can be written as

where $B = \sigma -\sum\limits_{j = 1}^{u}\frac{b_{j}\;\lambda _{j}^{q_{i}+1}}{\Gamma \left(q_{i}+2\right) }\neq 0,$ provided that $\sum\limits_{j = 1}^{u}\frac{b_{j}\;\lambda _{j}^{q_{i}}}{\Gamma \left(q_{i}+1\right) } > 1.$

Proof. Suppose that $\rho _{0}, \rho _{1}\in \Xi$ are vector constants. Based on ^[2], the solution of (3.1) takes the form

Applying the condition $\varpi (\varsigma) = \psi \left(\varsigma \right),$ we get

Using the condition $\varpi (\sigma) = \sum\limits_{j = 1}^{{u} }b_{j}\left(I^{q_{j}}\varpi \right) \left(\lambda _{j}\right),$ we have

From (3.3) and (3.4) in (3.2), we can write

□

After that, we need the following assertions:

(A$_{1}$) For all $\varsigma, \vartheta \in U,$ $\psi _{1}, \psi _{2}\in \Theta$ and $\varpi _{1}, \varpi _{2}, \widetilde{ \varpi }_{1}, \widetilde{\varpi }_{2}\in \Xi,$ $\ell _{g}, \ell _{P}, \ell _{Q}$ exist in order that

(A$_{2}$) For all $\left(\varsigma, \psi, \varpi _{1}, \varpi _{2}\right) \in U\times \Theta \times \Xi \times \Xi$ and $\left(\varsigma, \vartheta, \psi \right) \in \mho \times \Theta,$ $V_{j}\in L^{1}(U, \mathbb{R} _{+})$ $\left(j = 1, 2, 3, 4, 5\right)$ exists such that

(A$_{3}$) We consider $S = \ell _{g}N_{1}^{\ast }\left\{ \xi _{1}+\xi _{2}\left(\ell _{P}+\ell _{Q}\right) \right\} < 1,$ where

Now, the first main result in this part is as follows:

Theorem 3.3. Under Assertions (A$_{1}$) and (A$_{2}$), the BVP (1.1) has at least one solution on $(-\infty, \sigma ]$, provided that

Proof. The FP technique involves equating a given operator to the problem at hand and seeking a unique FP, which corresponds to the problem's unique solution. Therefore, we convert the BVP (1.1) to an FP problem. Define the operator $M: \widetilde{\mho }\rightarrow \widetilde{\mho }$ as

Assume that $\varrho (.):(-\infty, \sigma ]\rightarrow \Xi$ is a function described as

It is clear that $\varrho _{0} = 0.$ For every $\omega \in C(U, \Xi)$ with $\omega (0) = 0,$ we select

If $\varpi (.)$ fulfills (3.5), then we decompose $\varpi (.)$ as $\varpi (\varsigma) = \varrho (\varsigma)+\widetilde{\omega }(\varsigma),$ which leads to $\varpi _{\varsigma } = \varrho _{\varsigma }+\widetilde{\omega }_{\varsigma }$ for all $\varsigma \in U$, and $\omega (.)$ satisfies

Put $G_{0} = \left\{ \omega \in C(U, \Xi):\omega _{0} = 0\right\}$ and consider $\left\Vert.\right\Vert _{G_{0}}$ to be the seminorm in $G_{0}$ given by

Hence, $\left(G_{0}, \left\Vert.\right\Vert _{G_{0}}\right)$ is a BS. Describe the operator $\Phi :G_{0}\rightarrow G_{0}$ as

The existence of an FP for an operator $M$ is equivalent to the operator $\Phi$ having an FP. Hence, we focus on establishing the existence of an FP for $\Phi$.

Consider the set $H_{s} = \left\{ \omega \in G_{0}:\left\Vert \omega \right\Vert _{G_{0}}\leq s\right\}.$ Hence, $H_{s}$ is a bounded, closed, and convex subset of $G_{0}.$ Assume that there is a positive constant $\varepsilon$ such that $\varepsilon < s,$ where

Now, we decompose $\Phi$ as $\Phi _{1}+\Phi _{2}$ on $H_{s}$, where

and

Now, if we let $\omega, \omega ^{\ast }\in H_{s}$ and $\varsigma \in U$, we get

Hence,

where $V(\varsigma) = \max \left\{ V_{1}(\varsigma), V_{2}(\varsigma), V_{3}(\varsigma), V_{4}(\varsigma)\right\}$ and

It follows from (3.6) that $\Phi _{1}\omega +\Phi _{2}\omega ^{\ast }\in H_{s}.$ Now, we show that $\Phi _{2}$ is a contraction. For this, assume that $\omega, \omega ^{\ast }\in H_{s}$, and $\varsigma \in U$. We then

It follows that

Since $\ell < 1,$ then, $\Phi _{2}$ is contraction. Because $g$, $P$ and $Q$ are continuous, and thus $\Phi _{1}$ is continuous. Consider

This proves that $\Phi _{1}$ is uniformly bounded on $H_{s}.$ Finally, we prove that $\Phi _{1}$ is compact. Indeed, we claim that $\Phi _{1}$ is equicontinuous. For $\varsigma _{1}, \varsigma _{2}\in U,$ with $\varsigma _{1} < \varsigma _{2}$ and $\omega \in H_{s},$ one can write

Clearly, $\left\Vert \left(\Phi _{1}\omega \right) \left(\varsigma _{2}\right) -\left(\Phi _{1}\omega \right) \left(\varsigma _{2}\right) \right\Vert _{\Xi }\rightarrow 0$ as $\varsigma _{1}\rightarrow \varsigma _{2}.$ Consequently, $\Phi _{1}$ is equicontinuous. Applying the Arzelà -Ascoli theorem, we establish that $\Phi _{1}$ is compact on $H_{s}$. Consequently, invoking Krasnoselskii's FP theorem, we prove the existence of an FP $\omega \in G_{0}$, satisfying $\Phi \omega = \omega$, thereby yielding a solution to the fractional BVP (1.1). □

Now, for the uniqueness, we apply Banach's FP theorem as follows:

Theorem 3.4. Via Assertions (A$_{1}$) and (A$_{3}$), the BVP (1.1) owns a unique solution on $(-\infty, \sigma ]$.

Proof. Recall the set $H_{s} = \left\{ \omega \in G_{0}:\left\Vert \omega \right\Vert _{G_{0}}\leq s\right\}$ and assume that $\omega, \omega ^{\ast }\in G_{0}.$ For $\varsigma \in U,$ one has

Hence,

By (A$_{3}$), $S < 1,$ $\Phi$ is a contraction. By Banach's FP theorem, $\Phi$ possesses a unique FP, which is a unique solution to the problem (1.1) on the interval $(-\infty, \sigma ]$. □

4.   Solving the neutral FFIDE

In this section, we discuss the existence and uniqueness of solution to the considered problem (1.2) by applying Krasnoselskii's and Banach's FP theorems.

Assume that the space $\widetilde{\mho }$ is defined as in the section above and choose

Definition 4.1. We say that the function $\varpi \in \widetilde{\mho }$ is a solution to the FFIDE (1.1) if it fulfills the problem

With the aid of Lemma 3.2, the solution of the neutral FFIDE (1.2) takes the form

where $B = \sigma -\sum\limits_{j = 1}^{u}\frac{b_{j}\text{ }\lambda _{j}^{q_{i}+1}}{\Gamma \left(q_{i}+2\right) }\neq 0.$

To accomplish our main task here, we needs the following assertions:

(A$_{4}$) For all $\varsigma, \vartheta \in U,$ $\psi _{1}, \psi _{2}\in \Theta$ and $\varpi _{1}, \varpi _{2}, \widetilde{ \varpi }_{1}, \widetilde{\varpi }_{2}\in \Xi,$ $\ell _{g}, \ell _{h}, \ell _{P_{1}}, \ell _{Q_{1}}, \ell _{P_{2}}, \ell _{Q_{2}}$ exist such that

(A$_{5}$) For all $\left(\varsigma, \psi, \varpi _{1}, \varpi _{2}\right) \in U\times \Theta \times \Xi \times \Xi$ and $\left(\varsigma, \vartheta, \psi \right) \in \mho \times \Theta,$ $V_{j}\in L^{1}(U, \mathbb{R} _{+})$ $\left(j = 1, 2, 3, 4, 5\right)$ exists in order that

(A$_{6}$) Assume that $S^{\ast } = \ell _{g}N_{1}^{\ast }\left\{ \xi _{1}+\xi _{2}\left(\ell _{P_{2}}+\ell _{Q_{2}}\right) \right\} +\ell _{h}N_{1}^{\ast }\left\{ \xi _{1}+\xi _{2}\left(\ell _{P_{1}}+\ell _{Q_{1}}\right) \right\} < 1,$ where

Theorem 4.2. Under Assertions (A$_{4}$) and (A$_{5}$), the neutral BVP (1.2) has at least one solution on $(-\infty, \sigma ]$, provided that

Proof. Define the operator $\Game :\widetilde{\mho }\rightarrow \widetilde{\mho }$ as

Analogous to Theorem 3.3, define the operator $\Re :G_{0}\rightarrow G_{0}$ as

Describe the set $H_{\widehat{s}}$ as $H_{\widehat{s}} = \left\{ \omega \in G_{0}:\left\Vert \omega \right\Vert _{G_{0}}\leq \widehat{s}\right\}.$ Let there be a positive constant $\varepsilon ^{\ast }$ such that $\varepsilon ^{\ast } < \widehat{s},$ where

where

Now, we decompose $\Re$ as $\Re _{1}+\Re _{2}$ on $H_{\widehat{s}}$, where

and

Now, for $\omega, \omega ^{\ast }\in H_{\widehat{s}}$ and $\varsigma \in U$, we have

where $V^{\ast }(\varsigma) = \max \left\{ V_{1}(\varsigma), V_{2}(\varsigma), V_{3}(\varsigma), V_{4}(\varsigma), V_{5}(\varsigma), V_{6}(\varsigma)\right\}$ and

Hence,

Thus, $\Re _{1}\omega +\Re _{2}\omega ^{\ast }\in H_{\widehat{s}}.$ Now, we prove that $\Re _{2}$ is a contraction. Let $\omega, \omega ^{\ast }\in H_{ \widehat{s}}$ and $\varsigma \in U$. We then

It follows that

Since $\ell ^{\ast } < 1,$ then, $\Re _{2}$ is contraction. Since $g$, $P_{1}, P_{2}, Q_{1}$ and $Q_{2}$ are continuous, then $\Re _{1}$ is continuous. Furthermore,

Hence, $\Re _{1}$ is uniformly bounded on $H_{s}.$ Ultimately, we claim that $\Re _{1}$ is compact. Indeed, we prove that $\Re _{1}$ is equicontinuous. For $\varsigma _{1}, \varsigma _{2}\in U,$ with $\varsigma _{1} < \varsigma _{2}$ and $\omega \in H_{\widehat{s}},$ one has

Therefore, $\left\Vert \left(\Re _{1}\omega \right) \left(\varsigma _{2}\right) -\left(\Re _{1}\omega \right) \left(\varsigma _{2}\right) \right\Vert _{\Xi }\rightarrow 0$ as $\varsigma _{1}\rightarrow \varsigma _{2}.$ Hence, $\Phi _{1}$ is equicontinuous. By the Arzelà-Ascoli theorem, we establish that $\Re _{1}$ is compact on $H_{\widehat{s}}$. Consequently, invoking Krasnoselskii's FP theorem, $\omega \in G_{0}$ exists such that $\Re \omega = \omega$, which is a solution the neutral BVP (1.2). □

For the uniqueness, we have the following theorem:

Theorem 4.3. Via Assertions (A$_{4}$) and (A$_{6}$), the neutral BVP (1.2) has a unique solution on $(-\infty, \sigma ]$.

Proof. Define the set $H_{\widehat{s}} = \left\{ \omega \in G_{0}:\left\Vert \omega \right\Vert _{G_{0}}\leq s\right\}$ and assume that $\omega, \omega ^{\ast }\in G_{0}.$ For $\varsigma \in U,$ we get

Hence,

By (A$_{6}$), $S^{\ast } < 1.$ Thus, $\Re$ is a contraction. By Banach's FP theorem, $\Re$ has a unique FP, which is a unique solution to the problem (1.2) on $(-\infty, \sigma ]$. □

5.   Supportive examples

This section is devoted to testing the conditions of the proposed systems and their effectiveness, which leads to supporting and enhancing the theoretical results we obtained.

Example 5.1. Assume the following FFIDE:

where $\nu > 0$ is a real constant and define the set $H_{\nu }$ as

under the norm

Assume that $\varpi :(-\infty, \sigma ]\rightarrow \Xi$ in order that $\varpi _{0} = \psi \in H_{\nu }.$ Then

Therefore, $\omega _{\varsigma }\in H_{\nu }.$ Select $N_{1} = N_{2} = \kappa = 1.$ Hence, we show the condition

Clearly, $\left\vert \omega _{\varsigma }(\phi)\right\vert = \left\vert \omega (\varsigma +\phi)\right\vert.$ If $\varsigma +\phi \leq 0,$ we get

In the case of $\varsigma +\phi \geq 0,$ we have

Hence, if $\varsigma +\phi \in \lbrack 0, 1],$ we can write

which implies that

Furthermore, the pair $(H_{\nu }, \left\Vert \omega \right\Vert)$ is a BS and $H_{\nu }$ is a phase space. Here, $p = \frac{5}{2},$ $u = 3,$ and we choose

By simple calculation, we have

From (5.1), we have

where

Now, for $\varpi _{\varsigma }, \varrho _{\varsigma }\in H_{\nu },$ we have

and

It follows from (5.2)–(5.7) that $\ell _{g} = \frac{1}{64},$ $\ell _{P} = \frac{1}{16},$ $\ell _{Q} = \frac{1}{20},$ $V_{1}(\varsigma) = \frac{1 }{64},$ $V_{2}(\varsigma) = \frac{1}{32},$ $V_{3}(\varsigma) = \frac{1}{64},$ $V_{4}(\varsigma) = \frac{1}{32},$ $V_{5}(\varsigma) = \frac{1}{40},$ and $N_{1}^{\ast } = 1.$ Hence

and

Therefore, all the requirements of Theorems 3.3 and 3.4 are satisfied. Then, the considered problem (5.1) has a unique solution on $(-\infty, \sigma ].$

Example 5.2. Assume the following neutral FFIDE:

Assume that $H_{\nu }$ is the phase space, which is defined in Example 5.1, where $p = \frac{5}{2},$ $u = 3$, and

By simple calculation, we have

From (5.8), one can write

where

Now, for $\varpi _{\varsigma }, \varrho _{\varsigma }\in H_{\nu },$ we have

and

From (5.9)–(5.20), we have $\ell _{g} = \ell _{P_{1}} = V_{7}(\varsigma) = \frac{1}{5},$ $\ell _{h} = \ell _{Q_{2}} = V_{10}(\varsigma) = \frac{1}{4},$ $\ell _{P_{2}} = \ell _{Q_{1}} = V_{8}(\varsigma) = \frac{1}{6},$ $V_{1}(\varsigma) = V_{3}(\varsigma) = \frac{1}{35},$ $V_{2}(\varsigma) = \frac{1}{15},$ $V_{4}(\varsigma) = \frac{1 }{20},$ $V_{5}(\varsigma) = \frac{1}{160},$ $V_{6}(\varsigma) = \frac{1}{100}, ,$ $V_{9}(\varsigma) = \frac{1}{60},$ $N_{1}^{\ast } = 1.$ Thus, we can write

and

Hence, all the assertions of Theorems 3.4 and 4.2 are fulfilled. Therefore, the supposed problem (5.8) has a unique solution on $(-\infty, \sigma ].$

6.   Conclusions and future work

The study of FFIDEs presents a formidable challenge due to the inherent complexities arising from the interplay of fractional-order derivatives, functional arguments, and integral operators. Traditional methods often fall short in addressing these equations due to the non-local nature of fractional derivatives and the intricate dependence on past states introduced by functional arguments. Overcoming these difficulties requires the development and application of sophisticated mathematical tools, including specialized FP theorems tailored for fractional settings, careful treatment of infinite delay, and the construction of appropriate function spaces that accommodate the combined effects of these operators. Furthermore, the presence of multi-term fractional integral boundary conditions adds another layer of complexity, demanding innovative techniques for handling the non-local and distributed nature of the boundary constraints. Successfully navigating these hurdles necessitates a deep understanding of fractional calculus, functional analysis, and operator theory, ultimately paving the way for a more comprehensive understanding of the dynamics governed by FFIDEs. This paper investigates the existence and uniqueness of solutions for a class of hybrid fractional-order functional and neutral functional integrodifferential equations, featuring infinite delay and multi-term fractional integral boundary conditions. A rigorous mathematical framework is developed, leveraging FP theorems, to analyze these complex equations. The LC definition of fractional derivatives is employed, facilitating a comprehensive study of nonlocal dynamics. Illustrative examples are provided to demonstrate the applicability and practical relevance of the theoretical results. Future work includes exploring more complex equations (e.g., variable-order, generalized functional arguments), investigating stability, controllability, and numerical methods, and applying these equations to real-world problems. Developing new fixed point theorems tailored for fractional functional integrodifferential equations and studying associated inverse problems are also promising research avenues. Finally, we also look forward to extending the study period outside the proposed period ^[2,3].

Author contributions

Manal Elzain Mohamed Abdalla: Writing–review-editing, formal analysis, funding acquisition; Hasanen A. Hammad: Writing–original draft, conceptualization, investigation, methodology. All authors have read and approved the final version of the manuscript for publication.

Use of Generative-AI tools declaration

The authors declare they have not used artificial intelligence (AI) tools in the creation of this article.

Data availability statement

Data sharing is not applicable to the article as no data sets were generated or analyzed during the current study.

Acknowledgments

The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through the Large Research Project under grant number R.G.P. 2/217/45.

Conflict of interest

All authors confirm that they have no conflict of interest.