#### 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.

#### 1. Introduction

With the popularity of networking communication comes an increase in the psychoactive properties of this technological innovation. Numerous investigations have linked compulsive media platform use to undesirable outcomes, including decreased efficiency, adverse social interactions, and decreased life contentment. The misuse of social networking sites has increased at an alarming rate. For example, Facebook has 68 percent and 73 percent of the adolescent community in the United States, respectively [1]. Over prescription of social media addiction is connected with impaired effectiveness appraisal, problematic interpersonal interactions, insomnia, reduced levels of happiness, and sentiments of hostility, stress, and melancholy, see [2–7].To be consistent with the preponderance of the research, we included the terms “social media addiction” or “addicted social networking usage” (in a non-clinical meaning) in the entirety of that kind of analysis, while acknowledging the concerns surrounding the terminology [8]. Whenever it becomes more appropriate to employ certain concepts, exclusions are created. For the purposes of this evaluation, we characterize excessive social media adoption as becoming extremely apprehensive regarding social media, intrinsically incentivized, and dedicating a significant portion of time and vitality to incorporating social media to the juncture where an individual’s group opportunities, intimate communication, surveys, and/or general well-being and well-being are jeopardized [9].There is now a discrepancy in theorizing , particularly how it progresses. There have been multiple evaluations of conceptualizations to elucidate , but limited research has addressed the mathematical frameworks implemented in actual research on . One evaluation, in particular, barely discussed three extensively acknowledged internet dependency scenarios. One questionnaire range of interventions includes broad scientific viewpoints (for example, neurological aspects and behavioural insights) without addressing relevant paradigms. Although some alternative analyses reviewed various relevant conceptual models, analyzed their flaws, and proposed hypothetical modifications, they did not describe how scientific findings about implement these conceptualizations and the various attributes in these structures. Furthermore, new scientific findings on are appearing at an astonishing rate; there is a need for a comprehensive overview of essential mathematical approaches. Alemneh and Alemu [10] constructed and examined a mathematical formulation for the dissemination and prevention of in global society in 2021, as continues to follow:The current community is categorized into five compartments depending on addiction level in the scheme (1). Compartment 1: individuals that are likely to be exposed are those that are not obsessed but are accessible to is represented by . Compartment 2: individuals who utilize SM relatively regularly yet do not become obsessed are referred to as exposed individuals . Compartment 3: compulsive populace are persons who are obsessed to SM and spend a significant amount of time doing it . Compartment 4: retrieved individuals are persons who have progressed from , is the number of individuals who have benefited from . Compartment 5: represents anyone that doesn’t use or stop using SM perpetually. The number of individuals in the community is: . The model’s suppositions are as follows: the dissemination of the issue occurs in a confined area and is not reliant on intimate relations, racial group, or sentient socialist system. Representatives combining uniformly, and social networking dependent individuals will convey to non-addictive people when they are in contact with the compression of habit forming. Furthermore, the mathematical expressions of this scheme are amalgamated by employing the social networking addiction process, which begins with the introduction of vulnerable people into the community at a rate of . They are prompted to transfer to the revealed condition by addicted humans having a stress interaction incidence of and a probabilistic transference ratio of . At a proportion of , certain vulnerable persons migrate to a set of individuals who never access digital platforms. The challenged people are divided into two classes: one that develops obsession and advances to the addictive class at a pace , and another that recovers using therapy at a pace . Many dependent individuals transfer to the recovery class at a pace of or die as a result of the overuse of dependence on social networks at a ratio of . People who have overcome become vulnerable again at a ratio of or eventually discontinue access to digital platforms at a ratio of . Eventually, all of the inhabitants in each group have a spontaneous mortality rate of . Alemneh and Alemu [10] further examined the robustness of the equilibria and applied Pontryagin’s maximal principle to design the best monitoring mechanism.In past times, fractional calculus has been consistently proven to be an outstanding approach for illustrating the heredity properties of diverse formations. Furthermore, fractional differential operators have been used in a plethora of distinctive manifestations, such as hydroelectricity, fluid dynamics, chromatography, commerce, and viruses, utilising numerous fractional strategies such as Caputo, Hadamard, Riemann-Liouville, Katugumpola, Caputo–Fabrizio, Atangana-Baleanu and fractal–fractional, see [11–13]. This combo has suddenly accumulated a wealth of traction, partly because fractional differential equations have been shown to be superb tools for displaying a few incredibly challenging phenomena in a wide range of various and endless academic disciplines; assessors are encouraged to [14–20]. In 2019, the African mathematician Abdon-Atangana [21] proposed novel differential operators with fractal derivatives based on the combination of index law, modified Mittag-Leffler rule, and exponential decay law. Novel formulations, according to [20], offer additional non-local physical challenges with fractal characteristics simultaneously. Chen et al. [21] investigated the classification of anomalous propagation by deriving the underlying component of the numerical scheme using fractal derivative. For the purpose of evaluating processing performance and dispersion velocity, they compared fractional and fractal derivatives. Wei et al. [22] employed a scalability transformation technique to offer fractal modelling of elastic deformation. Sania et al. [24] introduced differential operators with fractional order and fractal dimension to describe additional chaotic complexity, wherein various kinds of orientations were mentioned: King Cobra, Shilnikov, Thomas cyclically symmetric, Langford, and Rössler. Rashid et al. [25] discussed the complex oscillatory behaviour of a human liver model via FF derivative operator technique. According to current evidence analysis of the data collected for this study, none of the works have addressed the mathematical framework of involving fractal-fractional derivatives. The novelty of this article is the addition of the FF fractional derivative to the model. The improved transmitting system incorporating the fractal-fractional derivative in the Atangana-Baleanu sense is proposed in the following framework:having initial conditions , where denotes the FF operator in terms of Atangana-Baleanu derivative operator, whilst presents the fractional-order and denotes the fractal-dimension, respectively. Furthermore, Table 1 contains explanations of all attributes. Recently, Alemneh and Alemu [10] demonstrated the mathematical modeling with optimal control analysis of . Kongson et al. [27] expounded the estimates for model pertaining to the Atnagana-Baleanu fractional derivative operator in the Caputo context. Leveraging the aforementioned proclivity, we employ a newly formed arbitrary order derivative in the model. The main purpose of this research is to investigate the model via a pioneering FF derivative operator in the Atangana-Baleanu sense and to characterize the complexities of the uniqueness and existence of the aforementioned framework response by utilizing the Picard-Lindlöf and contraction mapping techniques. Furthermore, to the best of our knowledge, there is no previous paper related to qualitative aspects of fractal–fractional derivative based on the Atangana-Baleanu context. As a result, the main objective of this paper is to bridge that gap. Graphical results demonstrate that while fixing and varying fractional-order and fractal-dimension how they affect the model’s characterizations.This article’s entire tasks are organized into five portions, which are provided as they continue to be implemented: Section 2 summarizes and presents the core notions and relevant formulations of FF derivatives in the Atangana-Baleanu interpretations. Section 3 describes the analysis of the model that corresponds to the positivity, bounded domain, invariant region, and equilibria points via the FF derivative operator technique. Section 4 illustrates and performs a convergence and uniqueness investigation for the above-mentioned with the aid of the FF derivative operator and the numerical formulation of the Newton polynomial approach. Section 5 deals with the numerical outcomes and description of the suggested model’s solution. The prime focus of the planned study is on which we will exhibit simulation consequences. Ultimately, Section 6 summarizes the remarks and discusses the promising possibilities of the fractal-fractional model.

#### 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 ICs**for 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 when**satisfying 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