#### Abstract

Excessive use of social media is a developing concern in the twenty-first century. This issue needs to be addressed before it has any more significant consequences than what we are currently experiencing. As a preventive technique, advertisements and awareness-raising campaigns about the detrimental impact of digital technologies are used. The application of novel mathematical techniques and terminologies in this field of study will have significant potential to enhance healthy living by preventing certain ailments. This is the most compelling justification for conducting a new study with the most up-to-date techniques at our disposal. This study investigates clear and concise transmission in order to generate a deterministic mathematical model of social media addiction using the fractal-fractional (FF) derivative operator. Also, the analysis of the model in terms of the invariant domain, the existence of a positive invariant solution, and equilibria assumptions are stated in a detailed manner. Besides that, the basic reproduction number is computed, demonstrating that the proposed methodology is more efficacious. The Atangana-Baleanu FF differential operators are recently defined in FF differential operators that are applied to characterize the model’s mathematical algorithm. We investigated the numerical behaviour of the in three ways: (i) changing the fractional-order as well as the fractal-dimension ; (ii) changing while keeping constant; and (iii) changing while keeping constant. Our examined visualizations and simulation studies using MATLAB for the numerical modelling of the aforementioned system showed that the novel developed Atangana-Baleanu FF differential operators produce remarkable results when compared to the classical frame.

#### 2. Preliminaries

In what follows, it is vital to investigate some fundamental FF operator notions prior to continuing on to the formal model. Consider there be a function , which is continuous and fractal differentiable on and has fractal-dimension and fractional-order , in addition to the specifications in [21, 24].

Definition 2.1. (see [21, 24]). Suppose there be a FF operator of having power law kernel in terms of Riemann–Liouville (RL) can be expressed in the form:where and .

Definition 2.2. (see [21, 24]). Suppose there be a FF operator of having exponential kernel in terms of RL can be expressed in the form:such that with .

Definition 2.3. (see [21, 24]). Suppose there be a FF operator of with Mittag-Leffler kernel in terms of RL can be expressed in the form:such that with .

Definition 2.4. (see [21, 24]). The corresponding FF integral form of (3) is described as:

Definition 2.5. (see [21, 24]). The corresponding FF integral version of (4) is described as:

Definition 2.6. (see [21, 24]). The corresponding FF integral form of (5) is described as:

Definition 2.7. (see [14]). Let and the Atangana-Baleanu-Caputo derivative operator is defined as:where represents the normalization function.

#### 3. Qualitative analysis of

In this section, we demonstrate that presented in (1) is epidemiologically viable by ensuring that the system’s corresponding model parameters are non-negative for every time-step . This is based on the more straightforward argument that the model with non-negative ICs becomes non-negative for all . The preceding is a lemma.

Lemma 3.1. Suppose there be the initial data , where . Thus the system (1) are positive for all. Also,having.

Proof. Assume that . Therefore, , the first equation of the framework (1) consists of followingwith , then (10) diminishes toIt follows thatConsequently, we haveNote thatThus, we can obtain for any by repeating the previous methods for the leftover equations of model (1).
Now adding the cohorts lead to the subsequentIf there is no death to the , thenTherefore, we havewhich is the desired proof.
Further, in order to show the invariant region for the proposed system (1), suppose

Lemma 3.2. The domain represented by is positively invariant for the system (1) along with non-negative ICsfor all.

Proof. In view of (12), then we haveTherefore, we have , if . Thus, we havewhich shows that, the domain presented by is positively invariant. Moreover, if or approaches to asymptotically. Hence, the domains presented by capture all of the possibilities in .

##### 3.1. Existence and nonnegativity of the solution

Further, we investigate the existence and nonnegativity of the system (2).

Theorem 3.3. If there be a unique solution of the system (2) and capture the solution in.

Proof. In order to prove the solution of the system (2) is positive, we havewhich indicates that the system (2) solution exist in . Summing up all cohorts in (2), we haveMoreover, we haveand hence, the biologically viable domain for the system (2) can be represented by

The framework for mentioned-above (2) in the FF operator in the Atangana-Baleanu sense is implemented to provide the outcomes in the next subsection.

##### 3.2. Stability result for disease free case

The stability outcomes for the framework introduced by at disease free equilibrium (DFE) are explored in this section. We can get the respective formulas by changing the right side terms of the (2) to zero, as

we have the following as follows

Furthermore, the fundamental reproduction number , which may be determined by applying the next generation methodology for the scheme, can be used to examine the robustness of DFE at . The infectious cohorts in the system (2) are , and the matrices and are obtained as follows:

Hence, the fundamental reproductive number can be calculated as

##### 3.3. Strength number

Following the work [28], we will present the strength number.In recent years, the concept of reproduction in a specific infectious problem has been extensively used in epidemiology modelling. As predicted by the concept, two components, and , will be identified in (27)

will be employed to generate the reproductive number [29]. The nonlinear portion of the classes that are infected is how the component , which is quite intriguing, gets derived

Again, we have

At DFE we have

In this case, we have the followingThen leads toTherefore, the dispersion will have a single magnitude and fade out if there is no regeneration mechanism, which is indicated by a value of . Furthermore, denotes a strength that will trigger a renewing mechanism, indicating that the expansion will contain multiple waves. However, biologists will give a precise explanation of the aforementioned number.

Theorem 3.4. The DFE at for model (2) is locally asymptotically stable whensatisfying the assumption

Proof. To illustrate the provided hypothesis, we must first acquire the Jacobian matrix by evaluating system (2) at the DFE , we haveThe negative eigenvalues are and remaining eigenvalues can be achieved from the following expression .the coefficients and are positive, for the DFE case, the value of should be less than 1. So that the Rough-Hurtwiz condition is satisfied for the assumptions presented if and only if and . Thus, the Rough-Hurtwiz condition promise the local asymptotic stability of the system presented (2) at DFE . The aforesaid results were achieved utilizing the FF framework as described and utilized in [30].

##### 3.4. Endemic equilibrium and their stability

The endemic equilibria of the (2) designated by and the outcome indicated in (2) are presented in this subsection as

Theorem 3.5. The presented by (2) has the following assertions:(a)if , then system (2) exhibits unique endemic equilibrium.(b)if , then system (2) has forward bifurcation.(c)if , then system (2) does not contain the endemic equilibrium point or has backward bifurcation.

#### 4. The fractal-fractional model

Here, we used the novel FF methodology in this section to reassemble the classical integer-order system with a non-singular and nonlocal kernel (2). The framework that ensues when the FF operator is taken into account is (2).

##### 4.1. Existence-uniqueness outcomes of FF-

Now, the existence-uniqueness of the obtained in the FF operator are succinctly discussed in (2). To do so, we shall use a FF derivative to generate the generic Cauchy problem:

In view of Definition (12), the right hand side of (38) yields:

Considering the implementation of the appropriate integral, the following conclusions are drawn as:

Employing the Picard-Lindelöf method, we havewhere .

Accordingly, surmise that

Furthermore, the norm is written as follows:and consider the operationsdescribed as

The essential objective is to illustrate that the aforementioned operator can convert a completely empty metric space onto itself. We also aim to illustrate that it has the potential to map contractions. First and foremost, we show that

Inserting , then produces the foregoing

Therefore,

Then, surmising . To obtain at the following result, apply the Banach fixed point theorem:where .

Owing to the contraction mapping , we have

Consequently, we have

If the supposition made is correct, thenthen the contraction criterion is achieved, i.e.,

In a nutshell, the proof is completed by demonstrating that there is only one solution.

In the next, we describes the numerical solutions for the proposed system.

##### 4.2. Newton polynomial approach

Here, we configure a comprehensive analysis of the numerical approach, which relies on an efficient Newton polynomial method. This methodology, which was also originally envisioned in [34], is more effective than some of the previous methods available in the analysis. To continue further with the approach, we apply the equation:

Integrating (55) with respect to , produces

Taking , then (56) reduces to

At , we have

Therefore, we have

To estimate the mapping, employ the Newton polynomial , we have

Substituting (60) into (57), yields

Simple computations yield

Note that

Using the fact that

Furthermore, we have

Therefore, a general approximate solution of the is as follows:

##### 4.3. New numerical technique for FF-AB derivative model

The objective of this task is to provide a structured approach technique for interacting with the (1) social media framework, using the FF operator in the Atangana-Baleanu context. converting the (2) system to the FF-Atangana-Baleanu derivative configuration as follows:

Employing the AB fractional integral operator, the preceding conclusions were made as