\subsection{Notations}

Consider two sets $A$ and $B$:
\begin{itemize}
\item $|A|$ refers to the number of elements of set $A$.
\item $A-B$ refers to the set of elements present in set $A$ and absent in set $B$.
\item $\triangle$ operator: $A\triangle B$ is defined as $A \cup B - A \cap B$.
For instance if $A=\{1,2,3,4,5\}$ and $B=\{4,5,6,7,8\}$ then $A \triangle B = \{1,2,3,6,7,8\}$
\item Consider $G$ a dynamic graph, $G_t = (V_t,E_t)$ denotes the state of the graph at time $t$, where $V_t$ is the set of vertices and $E_t$ is the set of edges.
\item $|V_t|$ (resp. $|E_t|$) corresponds to the number of vertices (resp. edges) of graph $G_t$.
\item for simplifying notations in the document, $|V_t|$ is often denoted by $n_t$
\item $G=(V,E)$ is said to be a null graph if both $V = \emptyset$ and $E = \emptyset$. In the report such a null graph can also be called an empty graph.
\end{itemize}

\subsection{Position of the work with respect to Temporal Networks}

Current definitions of temporal networks (TN) include \textit{time-varying graphs} \cite{santoro_time-varying_2011}, \textit{temporal networks} \cite{kempe_connectivity_nodate}, \textit{evolving graphs} \cite{ferreira_note_2002} etc.
They all define structures described by sequence of static graphs, ordered by a timestamp (e.g., $G=(G_i=(V_i, E_i))_{i \geqslant 0}$) where $i$ refers to the time step).
It is worth mentioning that TN definitions do not include information about the generative process.
Thus, the way the graph at time $t+1$ is obtained from the graph at time $t$ is not described.
In this report, the emphasis is precisely on the study of generative processes.
This work is therefore positioned upstream of TN.
In the sequel graphs produced by generators are called dynamic graphs or simply graphs.

\subsection{Generalities}

From a general point of view, a dynamic graph generator can be defined as a process with input data, that produces at each time step $t+1$ a new static graph $G_{t+1}$. It is produced from already generated static graphs $\{G_1,\ldots,G_t\}$ and possibly additional information.
Thus, the output of a dynamic graph generator is a flow of static graphs identified by time stamps.
The time stamps may also corresponds to events, and in such a case, the time interval between two time stamps may be different.
However, in this report, for sake of clarity, we consider integer time stamps.
If the flow stops, for whatever reason (e.g. clock has been stopped, evolution process is finished) at step $T$, the set of generated static graphs $\{G_1,G_2,\ldots,G_T\}$ corresponds to a temporal network (TN).

\begin{center}
\centering
\begin{tikzpicture}[-,thick,xscale=1.0,yscale=1.0]
\tikzstyle{process}=[rectangle,draw=black,fill=blue!10,text=black]
\tikzstyle{input}=[rectangle,draw=black,fill=white,text=black]
\tikzstyle{output}=[rectangle,draw=black,fill=blue!50,text=white]

    \node[process, text width=3.2cm] (G) at (5,5) {Evolution Process $\{G_1,G_2,\ldots,G_t\}$};
    \node[input] (I) at (1,5) {Inputs};
    \node[process] (CK) at (4,7) {Clock};
    \node[output] (TN) at (9,5) {$G_{t+1}$};
    \node (TIME) at (3.5,6.2) {$t+1$};
    \node[process] (INFO) at (5,3.5) {Information};

    \draw[->,thick] (I.east) --  (G.west);
    \draw[->,thick] (CK) -- (4,5.6);
    \draw[->,thick] (G) -- (TN);
    \draw[->,thick] (INFO) -- (G);
    \draw[->,>=stealth,thick] (G.east) .. controls (10,6) and (6,8) .. (G.north)
        node[sloped, above,midway] {$G_{t+1}$};

\end{tikzpicture}
\end{center}

\subsection{Sustainability}
\label{sec:sustainability}

The goal of this section is to introduce a novel notion for qualifying the dynamics of a graph.
The only measurement of the order (or of the density) of the graph is not enough for qualifying its dynamics.
For instance, if a dynamic graph becomes static, all vertices remain the same and the graph order does not change.
Conversely, if between two consecutive time steps all vertices are replaced by new ones, the order also remains the same, but the dynamics is different.
Sustainability qualifies a dynamic graph that never becomes null or periodic (which includes static).
A graph owing the sustainability property is said sustainable. \\

\begin{definition}
\underline{Graph sustainability}\\
A dynamic graph $G$ is said sustainable if both Condition 1 and Condition 2 are not verified.
$$
\begin{array}{ll}
   \mbox{Condition 1:} & \exists T \in \mathds{N}, \forall t \geqslant T, G_{t} = (\emptyset,\emptyset)\\
   \mbox{Condition 2:} & \exists T \in \mathds{N}~\mbox{and}~\exists k \in \mathds{N^*}, \forall t \geqslant T, G_t = G_{t+k}
\end{array}
$$
\end{definition}
\vspace{0.3cm}

Some well-known graph generators produce sustainable dynamic graphs.
For instance generators of growing networks.
Indeed, for all $t \in \mathds{N}$, $|V_t| > |V_{t-1}|$, and $G_t \neq (\emptyset,\emptyset)$.
For these generators, graph sustainability is obvious and does not require any analysis.

Unlike these cases, some generators are based on mechanisms making the evolution of the vertices (and edges) more difficult to predict, and the dynamics is worth studying.
For that purpose, we propose to consider a metric enabling a quantification of the dynamics.

\subsection{Nervousness}
\label{sec:nervousness}

This metric provides a way of measuring the dynamics of a graph.
Note that this metric is not new and is known as the Jaccard distance, which can be defined as one minus the coefficient of community as defined in \cite{jaccard_1912}.
However, in the context of dynamic graphs it seems to us more meaningful to call it nervousness.
It is different from the burstiness which is defined at the node/edge level during the lifetime of the graph \cite{goh.barabasi.2008} and aims at representing the frequency of events occurring on each node/edge.
Nervousness metric aim is to capture the dynamics of creation and deletion of nodes and edges between two time steps at graph level. \\

\begin{definition}
\underline{Vertices Nervousness}:\\
Given a dynamic graph $G$, such that at time $t$ $G_t = (V_t,E_t)$.
We call {\bf vertices nervousness} at time $t$ and denoted by \vn{t}, the ratio:
\[
  \mbox{\vn{t}} = \frac{|V_{t+1} \triangle V_t|}{|V_{t+1} \cup V_t|} = \frac{|V_{t+1} \cup V_t - V_t \cap V_{t+1}|}{|V_{t+1} \cup V_t|}
\]
\end{definition}
\vspace{0.2cm}

This metric is complementary with the graph order measure.
Indeed, graph order can remain constant between two consecutive time steps although some vertices change.
If all vertices are replaced, nervousness equals 1.
If all vertices are kept, nervousness is 0.
Similarly we define the edges nervousness as: \\

\begin{definition}
\underline{Edges Nervousness}:\\ Given a dynamic graph $G$, such that at time $t$ $G_t = (V_t,E_t)$.
We call {\bf Edges Nervousness} at time $t$ and denoted by \en{t}, the ratio:
\[
  \mbox{\en{t}}= \frac{|E_{t+1} \triangle E_t|}{|E_{t+1} \cup E_t|}
\]
\end{definition}
\vspace{0.3cm}

At graph level, nervousness is defined as a couple:\\

\begin{definition}
\underline{Graph Nervousness}:\\ Given a dynamic graph $G$, the {\bf Graph Nervousness} at time $t$ is:
\[
 \mbox{\gn{t}}=(\mbox{\vn{t},\en{t}})
\]
\end{definition}
\vspace{0.2cm}

For illustrating these definitions, consider the following cases for a dynamic graph, from $t$ to $t+1$.
We denote $|V_t| = n_t$.  We also assume that between $t$ and $t+1$ the order remains the same, thus $|V_{t+1}| = n_{t+1} = n_t = n$.
\begin{itemize}
    \item if all vertices are replaced: \\
	$$\mbox{\vn{t}}= \frac{|V_t \triangle V_{t+1}|}{|V_t \cup V_{t+1}|} = \frac{2n}{2n} = 1$$
    \item if half of the vertices are replaced: \\
	$$\mbox{\vn{t}} = \frac{|V_t \triangle V_{t+1}|}{|V_t \cup V_{t+1}|} = \frac{3n/2-n/2}{3n/2} = \frac{n}{3n/2} = \frac{2}{3}$$
    \item if the vertices remain the same, the union of the sets is equal to their intersection thus:\\
	$$\mbox{\vn{t}}=0$$
\end{itemize}

When the order changes, for instance if all vertices are duplicated, thus $|V_{t+1}| = 2n_t = 2n$: \\
	$$\mbox{\vn{t}} = \frac{|V_t \triangle V_{t+1}|}{|V_t \cup V_{t+1}|} = \frac{2n-n}{2n} = \frac{1}{2}$$

Same results hold for edge nervousness.
Dynamic graphs keeping the same set of vertices but with varying set of edges have a graph nervousness equal to: \gn{t} = (0,\en{t})\\

If \gn{t} = (\vn{t},0) with \vn{t} $ \neq 0$, then the set $V_{t+1}-V_t$ and the set $V_t - V_{t+1}$ contain only isolated nodes.

\subsection{Sustainability vs Nervousness}
\label{sec:sustainabilityvsnervousness}

Sustainability and nervousness are closely related.
Sustainability describes a dynamic graph property while nervousness enables the measure of the evolution of vertices and edges sets between two consecutive time steps.
When nervousness is null for both sets, the graph is static and thus does not have the sustainability property.
Given a dynamic graph $G$, if for all $t \in \mathds{N}$, \vn{t} $\neq 0~\mbox{or}~\mbox{\en{t}} \neq 0$, then $G$ holds the sustainability property, except if $G$ is periodic.

\subsection{D3G3: definition}

In this section we define a parameterized model generating families of dynamic graphs.
An instance of the generative model is defined by a set of parameters.
For studying the model, we analyze, according to the parameters set, the dynamic graphs families produced and rely on both the sustainability and the nervousness for that purpose.

The generator has two types of inputs: a set of parameters, $S_p$, and an initial graph, called seed graph and denoted $G_0$.
At each time step $t+1$ it produces, from the previous graph $G_t$, a new graph $G_{t+1}$ as illustrated on the figure.

\begin{center}
\centering
\begin{tikzpicture}[-,thick,xscale=1.0,yscale=1.0]
\tikzstyle{process}=[rectangle,draw=black,fill=blue!10,text=black]
\tikzstyle{input}=[rectangle,draw=black,fill=white,text=black]
\tikzstyle{output}=[rectangle,draw=black,fill=blue!50,text=white]

    \node[process,text width=3cm,align=center] (G) at (6,5) {Evolution Process $G_t$};
    \node[input] (G0) at (1,4) {seed graph $G_0$};
    \node[input] (P) at (1,6) {$S_p=\{\mbox{parameters}\}$};
    \node[process] (CK) at (5,7) {Clock};
    \node (TIME) at (4.5,6.2) {$t+1$};
    \node[output] (TN) at (10,5) {$G_{t+1}$};

    \draw[->,thick] (P.east) -- (3,6) -- (3.5,5) -- (G.west);
    \draw[->,thick] (G0.east) -- (3,4) -- (3.5,5) -- (G.west);
    \draw[->,thick] (CK) -- (5,5.5);
    \draw[->,thick] (G) -- (TN);
    \draw[->,>=stealth,thick] (G.east) .. controls (10,6) and (7,8) .. (G.north)
        node[sloped, above,midway] {$G_{t+1}$};
\end{tikzpicture}
\end{center}

Graphs produced by D3G3 are geometric graphs.
A geometric graph is defined by an euclidean space and a threshold $d$.
For this study, without loss of generality we consider a 2D-unit-torus (i.e., a square $[0;1[^2$ where the two opposite sides are connected).
Each vertex is characterized by a set of coordinates, such that given two vertices $u$ and $v$ it is possible to compute their euclidean distance: $dist(u,v)$.
Given $V$ the set of vertices, the set of edges $E$ is defined in the following way:
$E=\{(u, v) \in V^2~|~\mbox{dist}(u, v) \leqslant d\}$\\

Graphs generated by D3G3 are produced thanks to an evolution process.
This mechanism is parameterized by an initial graph (the seed graph) and by two transition rules driving the evolution of the graph between two consecutive time steps.
Apart from a random generator, no external decision or additional information is used by this mechanism.
Rules are based on node degrees only and rely on a random generator for positioning new nodes in the 2D euclidean space.
This leads to the name of the generator: \textit{Degree-Driven Dynamic Geometric Graphs Generator} or D3G3. \\

From now, graphs we are studying are referred to as sequences of static graphs $(G_t=(V_t,E_t))_{t\geqslant0}$, where $t\geqslant0$ is the time step.
The initial graph, $G_0$ ($t=0$) is called the seed graph. \\

\begin{definition}
\underline{Degree Driven Dynamic Geometric Graph Generator}
\\
An instance of D3G3 is defined by an initial graph, a set of parameters and two rules:
\begin{itemize}[topsep=2pt]
\setlength\itemsep{-0.1em}
\item $G_0 \neq (\emptyset,\emptyset)$ the seed graph,
\item parameters:
    \begin{itemize}[topsep=0pt]
    \setlength\itemsep{-0.1em}
    \item $d \in ]0, \frac{\sqrt{2}}{2}[$
    \item $S_s$ a set of non-negative integers
    \item $S_c$ a set of non-negative integers
    \end{itemize}
\item rules applied on $G_t$ leading to $G_{t+1}$:
    \begin{itemize}[topsep=0pt]
    \setlength\itemsep{-0.1em}
    \item if $v \in V_t$, then $v \in V_{t+1}$ iff $\mbox{deg}(v) \in S_S$ (conservation rule)
    \item if $v \in V_t$ and if $\mbox{deg}(v) \in S_C$ then add a new vertex to $V_{t+1}$ with a random position in the unit-torus (creation rule)
    \end{itemize}
\end{itemize}
\end{definition}
\vspace{0.2cm}

Thus, evolution of the graph between two consecutive time steps $t$ and $t+1$, is driven by two rules applied to each vertex $v \in V_t$ simultaneously.
The first rule determines for a vertex $v \in V_t$ whether it is kept at step $t+1$ while the second rule concerns the possibility for a vertex $v \in V_t$ to create a new vertex in $V_{t+1}$ according to its degree. \\

\begin{definition}
\underline{Conserved/Create/Removed/Duplicated nodes}
\\
Let $t \geqslant 0$ and $G=(G_t)$ a D3G. Let $u \in V_t$ and $v \in V_{t+1}$, then
\begin{itemize}
    \item $u$ is said to be a conserved node iff $u \in V_t \cap V_{t+1}$.
    \item $u$ is said to be removed iff $u \in V_t - V_{t+1}$.
    \item $u$ is said to be a creator/creating node iff $\mbox{deg}(u) \in S_C$.
    \item $v$ is said to be a created node iff $v \in V_{t+1} - V_t$.
    \item $u$ is said to duplicate iff it is both a conserved and a creator node.
\end{itemize}
\end{definition}
\vspace{0.2cm}

Once created, a node never change its position.
Positions of created nodes do not depend on the creating nodes positions.
The position of a created node is chosen randomly and uniformly over the unit-torus.

\subsection{Limit cases analyses}

Regarding the generative process, the obtained graphs depend on the threshold $d$, on the seed graph
$G_0$ and on the two sets $S_C$ and $S_S$. Each combination of $d$, $S_C$ and $S_S$ does not necessarily lead to sustainable graphs.
$n_t$ denotes the order of $G_t$, and $n_0$ the order of the seed graph ($G_0$).

Some preliminary remarks and a first theorem help the analysis.

\begin{remark}\label{remark:deathnodenevercomeagain}
As soon as a node is removed from the graph, the same node will never appear again in the graph: \\
$\forall v \in V_{t-1}, v \notin V_t$, $\forall t' > t, v \notin V_{t'}$
\end{remark}

\begin{remark}\label{remark:noperiodic}
Graphs generated by D3G3 are not periodic (period $k > 1$).\\
$\forall G = (V,E), \nexists k>1,~\forall t,  V_{t+k} = V_t$
\end{remark}

Remark \ref{remark:noperiodic} is the direct consequence of Remark \ref{remark:deathnodenevercomeagain}.
From these two remarks, Theorem \ref{theo:conclusionsustainable} can be formulated.

\begin{theorem}\label{theo:conclusionsustainable}
Given $G=(V,E)$ a dynamic graph generated by D3G3, if its order and its nervousness are never equal to 0 then the graph is sustainable.\\
$\forall t>0$, if $n_t > 0$ and \vn{t} $>0$ the graph is sustainable.
\end{theorem}

\subsubsection{Case:  $S_S = \mathds{N}$ and $S_C = \mathds{N}$}

At each time step each vertex of the graph is both a conserved and a creator vertex, it is then duplicated.
    Then, as soon as $G_0 \neq (\emptyset,\emptyset)$, the order of the graph increases exponentially: $n_t = 2^{t}n_0$ and thus
    $\lim_{t \rightarrow \infty} n_t = +\infty$.\\
    Consider an instance of dynamic graph produced by D3G3 with these parameters sets, is this graph sustainable?
    The analysis of $\cal{N}$$^V$~leads to the conclusion.
    Indeed, every vertex is duplicated, thus, $|V_{t+1}| = 2n_t$.
    From that, it comes $\mbox{\vn{t}} = \frac{|V_t \triangle V_{t+1}|}{|V_t \cup V_{t+1}|} = \frac{2n_t-n_t}{2n_t} = \frac{1}{2}$.
    In conclusion, for all $t$, (i) the graph is never null and (ii) \vn{t} $> 0$, thus from Theorem \ref{theo:conclusionsustainable} the graph is sustainable.

\subsubsection{Case:  $S_S = \mathds{N}$ and $S_C = \emptyset$}

 If $S_S = \mathds{N}$ and $S_C = \emptyset$, then for all $t$, whatever the degree of any vertex, it is not a creator, hence no new nodes are added to the graph.
In addition, all vertices have their degree in $S_S$, thus all nodes are conversed between $G_t$ and $G_{t+1}$.
So, $G_t = G_0$.
The nervousness \gn{t} $= (0,0)$, the graph is static: it is not sustainable.

\subsubsection{Case: $S_S = \emptyset$  and $S_C = \mathds{N}$}

This case is the opposite situation of the previous one.
All vertices are creators, but none is conserved.
Hence, for all $t > 0$, if $v \in V_t$, $v \notin V_{t+1}$, but, as $deg(v) \in S_C$, the creation rule generates a new node that replaces $v$.
So, for all $t$, $n_t = n_0$.
The graph is only dynamic, no subgraph is conserved between two consecutive time steps, all nodes are renewed.
From the nervousness point of view, this implies $V_t \cap V_{t+1} = \emptyset$ and $\mbox{\vn{t}} = \frac{|V_t \cup V_{t+1}|}{|V_t \cup V_{t+1}|}=1$, and the graph is sustainable.

\subsubsection{Case: $S_S$ is a non-empty finite set and   $S_C = \emptyset$}

Note first that no vertex will be created in this situation since $S_C = \emptyset$.
Thus for all $t$, $n_t \leq n_0$.
If $G_0$ is not null, then, several cases may occur:
\begin{enumerate}
\item first case: all vertices in $G_t$ have their degree in $S_S$: for all $v \in V_t$, $\mbox{deg(v)} \in S_S$.
Then, all vertices are conserved, the order remains unchanged and the graph is static, $G_{t+1} = G_t$,  hence the graph is not sustainable,
\item second case: opposite to the first one, for all $v \in V_t, \mbox{deg(v)} \notin S_S$, hence $G_{t+1}$ is null and the graph is not sustainable,
\item last case: some vertices have their degree not belonging to $S_S$, they are removed and as a consequence $n_{t+1} < n_t$.
At each occurrence of this case, the order of the graph is strictly decreased by one, then this can happen at most $n_0$ times before the graph becomes empty unless it becomes static before.
\end{enumerate}

As a consequence, for the limit case, the graph is not sustainable.

\subsubsection{Case: $S_S = \emptyset$ and  $S_C$ is a non-empty finite set}

Note first that for all $t$, $V_{t+1} \cap V_t = \emptyset$, no vertex is conserved from $t$ to $t+1$.
In addition, for all $t$, $n_{t+1} < n_t$, unless if for all $v \in V_t, \mbox{deg(v)} \in S_C$, then $n_{t+1} = n_t$.
So the evolution of the graph order is non increasing.
In addition, as $S_S = \emptyset$, \vn{t}$ = 1$, so the graph can never become static.
The only way for the graph to be non sustainable is to become null.
It is then not possible to conclude about the sustainability.
Indeed, consider the case for which at $t$, $0 < n_t \leq k+1$ and $S_C = [0,k]$ then for all $t' > t$, $n_{t'} = n_t$.
Such a dynamic graph is sustainable.
Conversely, if $S_C = \{k\}$, as soon as less than $k+1$ nodes have their degree equal to $k$ the graph becomes empty after two time steps.

\subsubsection{Graph Order Increase}

For some specific parameters sets, the increase of the order of the graph asymptotically tends to zero, which has some effect on graph sustainability.
This is the object of Theorem \ref{theo:low_increase}.

\begin{theorem}\label{theo:low_increase}
Let $G=(G_t)_{t \geqslant 0}$, $S_S$ and $S_C$ such as $S_S$ or $S_C$ equals $\mathds{N}$ while the other one is finite and not empty, then the probability that graph order increases tends toward 0 as the graph order tends to infinity:
$P(n_{t+1} > n_t) \underset{\scriptscriptstyle {n_t \to +\infty}}{\to} 0$.
\end{theorem}

\noindent {\textbf{Proof}}
Let denote $P(n_{t+1} > n_t)$ the probability that graph order increases between step $t$ and step $t+1$ and $D=S_S \cap S_C$.
For a node $u \in V_t$, $u$ duplicates if and only if $\mbox{deg}(u) \in D$.
\begin{align*}
P(n_{t+1} > n_t) &= P(\exists u, \mbox{deg}(u) \in D) \\
&= 1 - P(\forall u, \mbox{deg}(u) \notin D) \\
&= 1 - {(P(\neg D))}^{n_t }
\end{align*}
Where $P(\neg D)$ represents the probability a node does not satisfy the duplicating condition.
As positions of vertices are independent of each others:
\[
P(\neg D)=1-\sum_{k\in D}{n_t - 1 \choose k}p^k{(1-p)}^{n_t-k-1}
\]
Let $M = \max{D}$. For large values of $n_t$ and as $p \in ]0, 1[$:
\[\forall k \leqslant M, 0 < {n_t - 1 \choose k}p^k{(1-p)}^{n_t-1-k} \leqslant {n_t - 1 \choose M} {(1-p)}^{n_t-1-M}\]
As $0 < |D| \leqslant M+1$:
\begin{align} \label{eq:inequality_growth}
P(\neg D) \geqslant 1-\left({(M+1){n_t-1 \choose M}{(1-p)}^{n_t-1-M}}\right)^{n_t}
\end{align}
Knowing that
\[
\forall N, \forall k \leqslant N, {N \choose k} = \frac{N!}{k!(N-k)!} = \frac{1}{k!}{\displaystyle \prod_{i=N-k}^{N}(N-i)}
\]
and that $M$ is a constant not depending on $n_t$:
\[
{n_t - 1 \choose M}  \underset{\scriptscriptstyle {n_t \to +\infty}}{\sim} \frac{{n_t}^M}{M!}
\]
Rewriting the right side of \ref{eq:inequality_growth} leads to:
\begin{equation*}
A{n_t-1 \choose M}(1-p)^{n_t} \underset{\scriptscriptstyle {n_t \to +\infty}}{\sim} \frac{A}{M!}n_t^M(1-p)^{n_t}\textnormal{~~~where $A=\frac{M+1}{{(1-p)}^{M+1}}$}
\end{equation*}
As $1-p \in ]0;1[$ is not depending on $n_t$,
$n_t^M(1-p)^{n_t}$ tends toward 0 as $n_t$ tends to infinity. Thus, $A{n_t - 1 \choose M}(1-p)^{n_t} \to 0$, thereby:
\begin{align*}
(P(\neg D))^{n_t} &\underset{\scriptscriptstyle {n_t \to +\infty}}{\sim}
\exp{(n_t\ln{(1-\frac{A}{M!}{n_t}^M(1-p)^{n_t}}))}\\
& \underset{\scriptscriptstyle {n_t \to +\infty}}{\sim}
\exp{(\frac{A}{M!}{n_t}^{M+1}(1-p)^{n_t})}
\end{align*}
From this, $(P(\neg D))^{n_t} \to 1$,
and hence we deduce the wanted theorem. $\square$

\subsubsection{Case: $S_S = \mathds{N}$ and $S_C$ is a non-empty finite set}

As $S_S = \mathds{N}$, between two time steps, all nodes are conserved, thus, for all $t$, $V_t \subseteq V_{t+1}$.
The first consequence is that if $G_0$ is not a null graph, for all $t$, $G_t$ is not a null graph.
Let us analyze the evolution of the order of the graph.
Given $t$ such that, there exists a vertex $v$ such that $\mbox{deg(v)} \in S_C$, then $n_{t+1} > n_t$.
But, if for some $t'$, all $v \in V_{t'} , \mbox{deg(v)} \notin S_C$, then for all $t > t'$, $n_t = n_{t'}$ and the graph will be static.
From Theorem \ref{theo:low_increase}
$P(n_{t+1} > n_t) \underset{\scriptscriptstyle {n_t \to +\infty}}{\to} 0$ which implies $P(n_{t+1} - n_t = 0) \underset{\scriptscriptstyle {n_t \to +\infty}}{\to} 1$.
As all vertices are conserved between two time steps:
\[
P(|V_{t+1}-V_t| = 0) \underset{\scriptscriptstyle {n_t \to +\infty}}{\to} 1
\]
This leads to this limit on the probability the graph becomes static:
\[
 P(G_t~\mbox{is static}) \underset{\scriptscriptstyle {n_t \to +\infty}}{\to} 1\\
\]

Hence, for this limit case, we can conclude that the graph is asymptotically non sustainable.

\subsubsection{Case: $S_S$ is a non-empty finite set and   $S_C = \mathds{N}$}

For all $t$, for all $v \in V_t$, either $v$ is duplicated (if $\mbox{deg(v)} \in S_S$) or $v$ is not conserved but is a creator.
In both cases, between two consecutive time steps each vertex generates the creation of a new one.
Thus for all $t$,  \vn{t}~$>~0$ and $n_t~>~0$, hence, from remark 3, the graph is sustainable.

\subsubsection{Summary of results}

\begin{table}[!h]
    \centering
    \begin{tabular}{|c|c|c|c|} \hline
    \diagbox{$S_S$}{$S_C$} & $\mathds{N}$ & finite set & $\emptyset$ \\ \hline
     \multirow{3}{*}{$\mathds{N}$}
                & $\forall t, n_t = 2^tn_0$        & $\forall t, n_{t+1} \geqslant n_t$         & $\forall t, G_t = G_0$ \\
                & $\forall t, \mbox{\vn{t}} = 0.5$ & $\forall t, 0 \leq \mbox{\vn{t}} \leq 0.5$ & $\forall t, \mbox{\gn{t}} = (0, 0)$ \\
                &                                  & $\lim_{n_t \rightarrow \infty} P(n_{t+1} > n_t)$ = 0 & \\ \hline

    \multirow{2}{*}{finite}
        & $\forall t, n_{t+1} \geqslant n_t$ & General
        & $\forall t, n_{t+1} \leqslant n_t$ \\
       & $\forall t, 0.5 \leq \mbox{\vn{t}} \leq 1$ & cases & $\lim_{t \rightarrow \infty} n = \mbox{constant}$\\
       set & $\lim_{n_t \rightarrow \infty} P(n_{t+1} > n_t) = 0$ & (see Section \ref{sec:generalcases})  & $\lim_{t \rightarrow \infty} \mbox{\gn{t}} = (0, 0)$ \\ \hline

    \multirow{3}{*}{$\emptyset$} & $\forall t, n_{t+1} = n_t$  & $\forall t, n_{t+1} \leqslant n_t$ & \multirow{3}{*}{$\forall G_0, G_1=(\emptyset, \emptyset)$} \\
    & $\forall t, \mbox{\gn{t}} = (1, 1)$ &  $\forall t, V_t \neq \emptyset \implies \mbox{\gn{t}} = (1,1)$ & \\  \hline
    \end{tabular}
    \caption{Order and Nervousness evolution for the different cases. $n_t$ denotes the order of graph $G_t$, \vn{t} its vertices nervousness and \gn{t} the graph nervousness.}
    \label{tab:cases}
\end{table}
Stemmed from these results, sustainability property of dynamic graphs can be established.
The results are reported on Table \ref{tab:casessustainability}.

\begin{center}
    \begin{table}[!h]
    \centering
    \begin{tabular}{|c|c|c|c|} \hline
    \diagbox{$S_S$}{$S_C$} & $\mathds{N}$ & finite set & $\emptyset$ \\ \hline
     \multirow{2}{*}{$\mathds{N}$} & Sustainable & Asymptotically  & Non sustainable \\
	                                            &                     & non sustainable & \\ \hline
    \multirow{1}{*}{finite}       & Sustainable & General cases  & Non sustainable \\
                             set          &                    & (see Section \ref{sec:generalcases}) & \\ \hline
    \multirow{2}{*}{$\emptyset$} & Sustainable & Depends on  & Non sustainable\\
    &  & the parameters  & \\ \hline
    \end{tabular}
    \caption{Sustainability of dynamic graphs according to parameters sets of D3G3.}
    \label{tab:casessustainability}
\end{table}
\end{center}

\subsubsection{\texorpdfstring{$S_s \cap S_c = \emptyset$}{Sc cap Ss = 0}}

\begin{theorem}
{Disjoint sets}
\label{theo:disjoint}
\\
Let $t\geqslant0$ and $G_t = (V_t,E_t)$ a graph and $S_s$ and $S_c$ two sets of positive integers. If $S_s\cap S_c = \emptyset$, then the series $(|V_t|)_{t \geqslant 0}$ is decreasing.
\end{theorem}

\noindent {\bf Proof}

Let $(G_t=(V_t,E_t))_{t \geqslant 0}$ a graph and $S_s$ and $S_c$ two disjoint sets.
To prove the result, it is enough to prove that for all $t \geqslant 0$, $|V_{t+1}| \leqslant |V_t|$.

Let $t \geqslant 0$ and $u \in V_t$.
As $S_s\cap S_c = \emptyset$, nodes in $V_t$ can be divided in three groups:
\begin{itemize}
\item $V_s = \{u, \mbox{deg}(u) \in S_s~\mbox{and deg(u)} \notin S_c\}$,
\item $V_c = \{u, \mbox{deg}(u) \in S_c~\mbox{and deg(u)} \notin S_s\}$
\item $V_{others} = \{u, \mbox{deg}(u) \notin S_c~\mbox{and deg(u)} \notin S_s\}$
\end{itemize}

From the rules, each node in $V_s$ is conserved, each node in $V_c$ is removed but generates a new node and each node in $V_{others}$ is removed, thus:

$|V_{t+1}| = |V_s|+|V_c| \leq |V_s|+|V_c|+|V_{others}| = |V_t|$
$\square$

\subsubsection{\texorpdfstring{$S_s \cup S_c = \mathds{N}$}{Ss cup Sc = N}}

\begin{theorem}{Union set}\label{theo:union}
\\
Let $S_s$ and $S_c$ subsets of $\mathds{N}$. If $S_s \cup S_c = \mathds{N}$, then the series $(|V_t|)_{t \geqslant 0}$ is increasing.
\end{theorem}

\noindent {\bf Proof}

Let $(G_t=(V_t,E_t))_{t \geqslant 0}$ a graph and $S_s$ and $S_c$ two non disjoint sets.
To prove the result, it is sufficient to prove that for all $t \geqslant 0$, $|V_{t+1}| \geqslant |V_t|$.

As for previous case, let $t \geqslant 0$ and $u \in V_t$.
As $S_s\cup S_c = \mathds{N}$, nodes in $V_t$ can be divided in three groups:

\begin{itemize}
\item $V_s = \{u, \mbox{deg}(u) \in S_s~\mbox{and deg(u)} \notin S_c\}$,
\item $V_c = \{u, \mbox{deg}(u) \in S_c~\mbox{and deg(u)} \notin S_s\}$
\item $V_{others} = \{u, \mbox{deg}(u) \in S_c~\mbox{and deg(u)} \in S_s\}$
\end{itemize}

From the rules it comes: each node in $V_s$ is conserved, each node in $V_c$ is removed but generates a new node and each node in $V_{others}$ is conserved and generates a new node, thus:

$|V_{t+1}| = |V_s|+|V_c|+2|V_{others}| \geq |V_s|+|V_c|+|V_{others}| = |V_t|$
$\square$

\subsubsection{Partition sets}

In this section, $S_S$ and $S_C$ are considered to be a partition of $\mathds{N}$. This means $S_S \cap S_C = \emptyset$ and $S_S \cup S_C = \mathds{N}$. From theorems \ref{theo:disjoint} and \ref{theo:union}, one can say that for every graph $G_t=(V_t,E_t)$, the series $(|V_t|)_{t_ \geqslant 0}$ remains steady.
Two cases rises from that situations:

\begin{itemize}
    \item $S_S=\mathds{N}$ \textbf{and} $S_C=\emptyset$: in that case the graph is constant ($\forall t, G_t = G_0$).
    \item $S_S=\emptyset$ \textbf{and} $S_C=\mathds{N}$: the series of static graphs $(G_t)_{t \geqslant 0}$ is a series of independent random geometric graph with a constant number of nodes ($\forall t, n_t=n_0$).
\end{itemize}

\subsection{Model and conjecture}

In this section parameters $S_S$ and $S_C$ are limited to equal sets of consecutive integers.
Both sets are such that $S_S=S_C=[m,M]$ (called segments), where $m, M \in {\mathds{N}}^2$ and referred to as $S$ in the following.
The evolution of graph order for different values of parameters $m$ and $M$ is investigated.
Some statements and properties are theoretically and experimentally proved for the special case $S=\{0\}$.
A relationship between graph order at a step $t+1$ and at step $t$ and an upper bound for $n_t$ ($t > 0$) are given.
Then, a theoretical analysis of the general case is provided, and a new concept named sustainable interval is introduced.
In the last part of this section, vertices nervousness of graphs is studied through experimentation.
It is shown to be equal in average to $\frac{2}{3}$.
The reason behind this particular value will be explained in this last part.

\subsection{\texorpdfstring{$S=\{0\}$}{S=0}}

The case $S_S=S_C=S=\{0\}$ is considered in this section.
The seed graph, $G_0$, is supposed to be a random geometric graph whose order is arbitrarily chosen.
We first derive a bound on the graph order.
Then the transition from step 0 to step 1 is studied.
Finally, the mean value of graph order is estimated.
An approximation for small values of the distance threshold $d$ is provided.

\begin{theorem}{(Bounded graph order)}
\\
Let $S=\{0\}$, $d \in ]0; \frac{1}{2}[$ and $G_0=(V_0, E_0)$ such that there exists at least one node $u \in V_0$ being isolated (i.e., $\mbox{deg}(u)=0$), then, for all $t > 0$, $n_t \leqslant \frac{8}{\pi d^2}$
\end{theorem}

\noindent {\bf Proof}

The unit torus and $0 < d < \frac{1}{2}$ are considered.
Let $S=\{0\}$ and $G_0=(V_0,E_0)$.
Let $t > 0$ such that $G_t$ is not empty.
The size of $V_t$ is maximized as soon as there is not enough free space on the torus for adding a new isolated vertex.
If we consider an empty torus, a vertex can be put anywhere.
The area covered by this node is equal to $\pi d^2$.
For the rest of the proof, let consider, for each vertex $u$ on the torus, the disk of radius $d/2$ and center $u$ and referred it as $D(u)$.
From this, the condition for two nodes $u$ and $v$ to be non connected is:
$D(u)\cap D(v) = \emptyset$.
After the addition of the first vertex $u_1$, a second vertex $u_2$ can be added to the torus if it satisfies $D(u_1) \cap D(u_2) = \emptyset$.
The area occupied by the two disks is then $2\times \pi (d/2)^2$, and the remaining free area is thus $1-2\times \pi (d/2)^2$.
Assuming $N$ vertices with non intersecting disks are already present in the torus, the free area is then $1-N\pi \frac{d^2}{4}$.

This quantity has to be positive so $N$ must be lower than or equal to $\frac{4}{\pi d^2}$.
This last quantity is an upper bound for the number of isolated vertices.
However, the rules says that isolated vertices duplicate.
Thus, the number of vertices at one step $t$ can not exceed twice the upper bound of isolated vertices.
Hence, for all $t>0$, $n_t \leqslant \frac{8}{\pi d^2}$.$\square$ \\

This theorem provides an upper bound for graph order for $S=\{0\}$, ensuring graph order cannot exceed a certain value.
However, this does not provide any information about graph order evolution, which is the purpose of the two following theorems.

\begin{theorem}\label{theo:RGGOrderS1}
{(Expected graph order at step 1)}
\\
Let $S=\{0\}$, $d>0$ and $G_0=(V_0, E_0)$ be a random geometric graph of order $n_0$, then $\frac{n_1}{2} \sim B(n_0, (1-p(d))^{n_0 - 1})$, where $p(d)$ is the area of a circle of radius $d$ on the torus. An expectation value for $n_1$ is therefore $2n_0\cdot(1-p(d))^{n_0-1}$.
\end{theorem}

\noindent {\bf Proof}

For sake of clarity, in the remaining part of this proof, $p(d)$ will be referred to as $p$.
At the very first step $t=0$, $G_t=G_0$ is a random geometric graph and its nodes are uniformly distributed over the unit torus.
Let $(u, v)\in {V_0}^2$.
For a fixed threshold $d$, let consider the probability that $u$ and $v$ are connected.
$v$ is connected to $u$ if and only if $\mbox{dist}(u,v) \leqslant d$.
It means $v$ is in the disk of center $u$ and radius $d$.
If we denote by $X(u,v)$ the event ``$u$ and $v$ are connected", the wanted probability is the ratio of the area of the surface of a disk of radius $d$ over the area of the unit torus.
For all $(u,v)\in {V_0}^2$, $X(u,v) \sim B(p)$.

Let's study the degree distribution of a node $u\in V_0$.
As the position of every point is independent one from the others, variables $X(u,x)$ are independent for all $x\neq u$.
More over the degree of $u$ is the number of connections $u$ has to other nodes :
\[
\mbox{deg}(u)=\sum_{u\neq v}X(u,v)
\]
All variables $X(u,v)$ being independent for all $v\neq u$,
$\mbox{deg}(u) \sim B(|V_0|-1, p)$ as a sum of $|V_0|-1=n_0-1$ independent Bernoulli variables of same parameter p.
Let's consider $Y_0(u)$ the event ``$u$ is conserved at step 1" ($Y_0(u)=1$ if and only if $u$ is conserved and 0 otherwise), then, as $S=\{0\}$:
\[
P(Y_0(u)=1)=P(\mbox{deg}(u)=0)=(1-p)^{n_0-1}
\]
This means ${Y_0(u)} \sim B((1-p)^{n_0-1}) $.
Thus the number of conserved vertices at step 1 is: \[
Y_0=\sum_{u\in V_0}{Y_0(u)}
\]
As all points have an independent position, $\frac{n_1}{2} = Y_0 \sim B(n_0, (1-p)^{n_0-1})$.$\square$

\begin{theorem}\label{theo:RGGOrder}
{(Expected value of graph order)}
\\
Let $S=\{0\}$, $d>0$ and $G_0=(V_0, E_0)$ such that there exists at least one node $u \in V_0$ being isolated (i.e., $\mbox{deg}(u)=0$), then
either the graph becomes empty, or the average number of conserved nodes is $l(d) = 1-\frac{\log{(\frac{\sqrt{1+4\alpha}-1}{2}})}{\log{\alpha}}$ with $\alpha=\frac{1}{1-p}$ and $p=p(d)$.

\end{theorem}

\noindent {\bf Proof}

Let $t \geqslant 1$. Two cases are to be discussed:  the case of conserved vertices from step $t-1$ to step $t$ ($V_t \cap V_{t-1}$) and the case of created nodes at step $t$ ($V_t - V_{t-1}$).
As the number of created nodes is the same as the number of conserved nodes from $t-1$ to $t$, we set $c_t=|V_t \cap V_{t-1}|=|V_t - V_{t-1}|$.

First let's study the number of conserved vertices from step $t$ to step $t+1$ among those conserved from step $t-1$ to step $t$.
$c_{t+1}^{\mbox{\footnotesize conserved}}$ denotes this number.
Let $u\in V_t \cap V_{t-1}$. The probability for $u$ to be conserved is the probability that its degree to created nodes remains equal to 0.
\[
\mbox{deg}(u)=\sum_{v\in V_t-V_{t-1}}X(u,v)
\]
Let $v \in V_t-V_{t-1}$.
As in the previous section, $X(u,v) \sim B(p)$ and $\mbox{deg}(u) \sim B(c_t, p)$ as a sum of independent Bernoulli variables of same parameter $p$.
$Y_t(u)$ denotes the event ``$u$ is conserved at step $t+1$".
The probability that $u$ survives is $P(Y_t(u)=1)=P(\mbox{deg}(u)=0)=(1-p)^{c_t}$, thus: $Y_t(u) \sim B((1-p)^{c_t})$.
Therefore, the number of conserved vertices at step $t+1$ among those conserved at step $t$ is:
\[
c_{t+1}^{\mbox{\footnotesize conserved}}=\sum_{u\in V_t \cap V_{t-1}}Y_t(u)
\]
As the position of created nodes are independent from themselves and from conserved vertices, $Y_t(u)$ are independent for all $u\in V_t \cap V_{t-1}$, $c_{t+1}^{\mbox{\footnotesize conserved}} \sim B(c_t, (1-p)^{c_t})$.

Let's study the number of conserved vertices among created nodes.
$c_{t+1}^{\mbox{\footnotesize created}}$ denotes this number.
Let $u \in V_t-V_{t-1}$. To study the degree of $u$, two cases must be studied.
The first one is the number of connections between $u$ and all other created nodes (denoted as $\mbox{deg}^C(u)$). The second one is the number of connections to already present nodes (denoted as $\mbox{deg}^S(u)$).
$\mbox{deg}^C(u)$ and $\mbox{deg}^S(u)$ can be obtained using the following formulas:

\begin{align*}
    \mbox{deg}^C(u) &= \sum_{v\in V_t - V_{t-1}, u \neq v }X(u,v) \\
    \mbox{deg}^S(u) &= \sum_{v \in V_t \cap V_{t-1}}X(u,v)
\end{align*}

As the position of created points on the torus are independent one from the others, $\mbox{deg}^C(u)$ is a sum of independent Bernoulli variables and therefore,
$\mbox{deg}^C(u) \sim {B(c_t -1, p)}$. For $\mbox{deg}^S(u)$, connections between a created node and an already present node are not independent from each other: knowing $u$ is connected to an already present node means it is close to it and as other conserved nodes are farther than $d$, it implies that $\mbox{deg}^S(u)$ is not a sum of independent Bernoulli variables.
However, as a first approximation, this quantity will be considered as a sum of independent Bernoulli variables.

Thus, the computation of the expectation of $c_{t+1} = c_{t+1}^{\mbox{\footnotesize conserved}}+c_{t+1}^{\mbox{\footnotesize created}}$ gives:
\[
c_{t+1}=c_t(1-p)^{c_t}+c_t(1-p)^{2c_t-1}
\]
By looking for a limit to this series gives $l \geqslant 0$ satisfying:
\[
l=l(1-p)^l+l(1-p)^{2l-1}
\]
Solving this equation gives $l=0$ or :
\begin{eqnarray*}
    l = 1 -  \frac{\log{\left({\frac{\sqrt{1+4\alpha}-1}{2}}\right)}}{\log{\alpha}} & & \textnormal{with $\alpha=\frac{1}{1-p}$}
\end{eqnarray*}
$\square$

Experiments have been run to see if this relationship holds.
The results are summarized on figure \ref{fig:expected_relationship_isolated}.
This figure shows the accuracy of the expectation of graph order given in theorem \ref{theo:RGGOrder}.
Indeed, the blue curve, close to the dashed line, highlights that theoretical expectation and experimental results are close to be equal.
This result being proved leads to an approximation for small values of threshold $d$.

\begin{corollary}
Let $d > 0$ and $l(d)$ as defined in the previous theorem (\ref{theo:RGGOrder}). Then for small values of $d$:
\[l(d) \sim -\frac{\log{\left(\frac{\sqrt{5}-1}{2}\right)}}{\pi d^2} = \frac{\log{\phi}}{\pi d^2}\]
where $\phi$ is the golden ratio $\left(\frac{1+\sqrt{5}}{2}\right)$.
\end{corollary}

\noindent{\textbf{Proof}}
Let $d > 0$ be small. Thus, applying Taylor expansion gives $\frac{1}{1-\pi d^2} \sim 1 + \pi d^2$ and $\log{\left(\frac{1}{1-\pi d^2}\right)} \sim \pi d^2$.
The numerator comes from $4 \cdot \frac{1}{1-\pi d^2} \simeq 4$.
The golden ratio is obtained using operations on $\log$ and by noticing that $\frac{2}{\sqrt{5}-1} = \frac{2(\sqrt{5}+1)}{4} = \phi$, the golden ratio.
Combining these results leads to the statement of the corollary.
$\square$\\

It is therefore possible to state that, in the case where $S=\{0\}$, it is possible to theoretically get an expectation of graph order as well as to get an upper bound for graph order depending on parameter $d$.

\begin{figure}[h!]
  \begin{center}
    \begin{tikzpicture}
      \begin{axis}[
          width=\linewidth,
          grid=major,
          grid style={dashed,gray!30},
          xlabel=$c_t$,
          ylabel=Expected value,

          title style={yshift=0.5ex, align=left},
          title={Relationship between the average value of $c_t$ and\\the expected value mentioned in theorem \ref{theo:RGGOrder}.},
          xmin = 0, xmax = 155000,
          ymin = 0, ymax = 155000
        ]

        \addplot[color=blue!50!cyan,only marks]
        table[x=Expected,y=Survivors,col sep=comma] {NodeCountExpected.csv};

        \addplot[color=black,smooth,dashed,tension=0.7,very thick] coordinates{(0,0)(155000,155000)};
        \legend{$c_t=f($Expected average value$)$, Identity ($y=x$)}
      \end{axis}
    \end{tikzpicture}
    \caption{Relationship between the average value of $c_t$ and the expected value. Each point correspond to a single threshold $d$. $d$ is ranging from 0.001 to 0.01 with a step of 0.0005 and from 0.01 to 0.2 with a step of 0.005}
    \label{fig:expected_relationship_isolated}
  \end{center}
\end{figure}

\subsection{The general case}

Now the focus is on $S=[m, M]$ for every $m$ and $M$ integers.
The goal is to provide a tool aiming at stating, for given parameters $m, M$ and $d$, whether the graph is likely to be sustainable or not.
This part mainly focuses on a simpler model.
This model is studied as it helps understanding the evolution of graph order.

\input{R4-RedistributedModel.tex}

\subsection{Vertex nervousness}

The goal is to highlight a characterization aspect of the segment family using the vertex nervousness metric.
As edge nervousness will not be studied for that case, vertex nervousness will be referred to as nervousness in this section.
As in this particular configuration, survivors are the same as created nodes, it is possible to state particular results about the value of nervousness:
\begin{theorem}\label{theo:nervousness_segment}
Let $S$ be a segment set of non-negative integer and $d \in ]0, \frac{1}{2}[$. Let $G$ be a generated graph of order $n_t$ at step $t$ and number of survivor from step $t$ to step $t+1$ referred to as $s_t$. Then:
\[
    \mbox{VN}_t = \frac{n_t}{n_t+s_t}
\]
\end{theorem}

\noindent{\textbf{Proof}}
To prove this result, it is sufficient to notice that $n_{t+1} = 2s_t$, as $S_S=S_C$, which means the number of survivors is the same as the number of created nodes.
Thus, applying some basic result about set sizes and noticing that $s_t = |V_t \cap V_{t+1}|$, leads to:
\begin{align*}
    |V_t \cup V_{t+1}| &= n_t + n_{t+1} - s_t = n_t + s_t \\
    |V_t \triangle V_{t+1}| &= n_t + n_{t+1} - 2|V_t \cap V_{t+1}| = n_t
\end{align*}
It follows that vertex nervousness is well equal to $\frac{n_t}{n_t+s_t}$.
$\square$\\
Result about the nervousness observed in generated graphs parameterized with a segment set $S$ is stated in the following conjecture:
\begin{conjecture}
Let $m, M \in \mathds{N}$. Let $S = [m, M]$ and $d > 0$ be parameters of the graph $G=(G_t)_{t \leqslant 0}$. Then the nervous of the graph is in average equal to $\frac{2}{3}$.
\end{conjecture}
Although this conjecture has not been proved theoretically, experimentation have been performed.
They all highlight this conjecture telling that the average nervousness of generated graphs is roughly equal to $\frac{2}{3}$.
Results of this experimentation are gathered on picture \ref{fig:nervousess_mean}.
A possible interpretation of this conjecture and performed experimentation relies on the result stated in theorem \ref{theo:nervousness_segment} and on results from last part.
Indeed, if vertex nervousness is close to $\frac{2}{3}$, it means $s_t\simeq\frac{n_t}{2}$.
Then, as $n_{t+1}=2s_t$, it comes $n_{t+1}\simeq n_t$, which means that graph order is close to a fixed point of the relationship $f_{S,d}$ mentioned in the previous section.

\begin{figure}
    \centering
\begin{tikzpicture}
\begin{axis}[
          view = {45}{10},
          width=0.7\linewidth,
          xlabel=$m$,
          ylabel=$M$,
          zlabel=$\bar{\mathcal{N}^v_t}$,
          legend style={at={(0.05.,.98)}},

          title style={yshift=0.5ex},
          title={Average nervousness},
          zmin=0.57, zmax=0.7,
          only marks
        ]

    \addplot3+
    [
        samples= 9,
        on layer=background,
        z filter/.expression={z<0.66 ? z : nan}
    ]
    table [x=mS, y=MS, z=Nervousness,col sep=semicolon]{nervousness_mean.csv};

    \addplot3+
    [
        samples= 9,
        fill=red,
        on layer=foreground,
        z filter/.expression={z>0.66 ? z : nan}
    ]
    table [x=mS, y=MS, z=Nervousness,col sep=semicolon]{nervousness_mean.csv};

    \addplot3 [
    fill opacity=0.75,
    shader=faceted,
    domain=20:65,
    domain y = 25:145,
    samples = 8,
    samples y = 8,
    surf,
    fill=yellow,
    faceted color = teal,
    on layer=main] {0.6666666666};

\end{axis}
\end{tikzpicture}
    \caption{Mean value of nervousness got from experimentation.
    Points represent the average over 20 run and 30000 time steps for a single $m$ and $M$.
    The yellow surface is the plan of equation $z=\frac{2}{3}$.
    For all these parameters, $d$ is set to $0.05$.
    Red points represent nervousness of value greater than $\frac{2}{3}$.
    Blue points represent nervousness of value lower than $\frac{2}{3}$.
    }
    \label{fig:nervousess_mean}
\end{figure}

\subsubsection{Study of graph evolution}

In this Section we aim at estimating the evolution of the graph order during graph dynamics.
However, in the D3G3 model, between two time steps, non-conserved nodes are removed from the graph and conserved nodes are located at the same position, which entails a remanent graph.
This remanent graph induces a structure influencing the computation of graph order.
More precisely, nodes that are about to be removed connected to conserved ones interfere in the probability that conserved nodes at time $t$ are still conserved at time $t+1$.
This is linked to computing the degree of the neighbors of a node $u$ knowing the degree of node $u$.
To our knowledge, this is a difficult question.
For that purpose, a relaxed version of the D3G3  model is considered enabling analytical study of this evolution.
In this model, conserved nodes are moved (i.e., their position are changed) such that obtained graph is a new random geometric graph at each step.
We call this model "the redistributed model".
This will help us proving the following theorem:

\begin{theorem}\label{theo:FRGGP_law}
Let $G=(G_t)$ be a dynamic graph obtained with the redistributed model, then at every step $t$, $\frac{n_{t+1}}{2} \sim B(n_t, p(S, d, n_t))$, where $p(S, d, n_t)$ is the probability that a node is conserved between step $t$ and $t+1$:
\[
p(S, d, n_t) = \sum_{k = m}^M{n_t - 1 \choose k}p^k(1-p)^{n_t - 1 - k}
\]
Here, $p(d)$ refers to the probability for two different nodes to be connected (i.e., the probability that the distance between them is lower than or equal to $d$), which is, for $d \leqslant \frac{1}{2}$, $\pi d^2$.
\end{theorem}

\noindent{\bf Proof}

{In the redistributed model, at time step $t$ a RGG ($G_t$) is built.
If the graph order at time $t$ is equal to $n_t$, the graph order at $t+1$ is equal to twice the number of surviving nodes at time $t$.}
As every node has an independent position {in the torus}, this probability is the same for all nodes.
Let's denote it $p(S, d, n_t)$.
Let $u \in V_t$.
Then:
\begin{align}
p(S, d, n_t) &= P(\mbox{deg}(u) \in S) = \sum_{k = m}^MP(\mbox{deg}(u) = k) \\
p(S, d, n_t) &= \sum_{k = m}^M{n_t - 1 \choose k}p^k(1-p)^{n_t - 1 - k} \label{eq:probaSurviving}
\end{align}
Assuming one node is a conserved node does not affect the probability of conservation for other nodes.
The number of conserved nodes can be computed summing independent Bernoulli's events of parameter $p(S, d, n_t)$.
This gives $\frac{n_{t+1}}{2}$ follows a binomial distribution of parameter $n_t$ and $p(S, d, n_t)$.
$\square$ \\

Computing expectation for a binomial distribution leads to an expectation for $n_{t+1}$ knowing $n_t$.
Indeed, this expectation is $2n_tp(S, d, n_t)$.
For a fixed set $S$, this provides a relationship between $n_t$ and $n_{t+1}$:
\begin{definition}
\underline{Expectation of graph order}:\\
Let $m$, $M$ and $d$ be parameters for the redistributed model.
Let $G=(G_t)$ be an obtained graph with such parameters.
Then, the expectation of graph order at step $t+1$ ($n_{t+1}$) knowing graph order at step $t$ ($n_t$) is $f_{S,d}(n_t)$ and satisfies $n_{t+1} = f_{S,d}(n_t) = 2n_t p(S,d,n_t)$, and then:
\begin{equation}
\forall n \in \mathds{N}, f_{S,d}(n) = 2n p(S,d,n)
\end{equation}
\end{definition}

This quantity is referred to as the relationship in the sequel.
Studying the relation for every value of $m$, $M$ and $d$ turns out to be a difficult problem.
However some results may be conjectured.
A first conjecture concerns the variations of the relationship:
\begin{conjecture}\label{cj:increase_decrease}
Let $m, M$ and $d$ be parameters of the model. Let $S=[m, M]$ and $f_{S,d}$ the relationship as defined above. Then there exists $n_* \in \mathds{N}$ such that $f_{S,d}$ is increasing on $[0, n_*]$ and decreasing on $[n_*+1, +\infty[$.
\end{conjecture}
This conjecture is difficult to prove due to the sum involved in the computation of $f_{S,d}$.
However, it is not necessary to study the relationship for all integers.
It is possible to perform the study on a limited interval.
This is the purpose of theorem \ref{theo:increase_decrease} (below).
But before proving this theorem, it is necessary to provide another formulae computing variations of $f_{S,d}$:
\begin{lemma}
\label{lm:variation}
Let $m$, $M$ and $d$ be parameters of the model.
Let $\Delta f_{S,d}$ defined as the variation of $f_{S,d}$: for $n \in \mathds{N}, \Delta f_{S,d}(n)= f_{S,d}(n+1) - f_{S,d}(n)$.
Then:
\[
\forall n \in \mathds{N}, \Delta f_{S,d}(n) = 2\sum_{k=m}^M {(k+1){n \choose k}p^k(1-p)^{n-1-k}}\left(1-\frac{n+1}{k+1}p\right)
\]
\end{lemma}

\noindent \textbf{Proof}:
Let $m$, $M$ and $d$ be parameters of the model.
Let $n$ a be non-negative integer.
This proof only focuses on the terms of the sum of $\Delta f_{S,d}$:
\begin{align*}
\Delta f_{S,d}(n) &= 2\left(\sum_{k=m}^M(n+1){n \choose k}p^k(1-p)^{n-k}-n{n-1 \choose k}p^k(1-p)^{n-1-k}\right) \\
&= 2\sum_{k=m}^M{p^k(1-p)^{n-1-k}\left((n+1){n \choose k}(1-p) - {n-1 \choose k}\right)}
\end{align*}
Let $k \in \mathds{N}$ such that $m \leqslant k \leqslant M$.
Every term of the sum of $\Delta f_{S,d}$ can be expressed as follow only using results on binomial coefficients:
\begin{align*}
(n+1){n \choose k}(1-p) - n{n-1 \choose k} &= (k+1){n+1 \choose k+1}(1-p) - (k+1){n \choose k+1} \\
&= (k+1)\left({n+1 \choose k+1}-p{n+1 \choose k+1} - {n \choose k+1}\right) \\
&= (k+1)\left({n \choose k} - p{n+1 \choose k+1}\right) \\
&= (k+1)\left({n \choose k} - p\frac{n+1}{k+1}{n \choose k}\right) \\
&= (k+1){n \choose k}\left(1-p\frac{n+1}{k+1}\right) \\
\end{align*}
\vspace{-0.1in}
This last equality leads to the following form of $\Delta f_{S,d}$:
\begin{equation*}
\Delta f_{S,d}(n) = 2\sum_{k=m}^M {(k+1){n \choose k}p^k(1-p)^{n-1-k}}\left(1-\frac{n+1}{k+1}p\right)
\end{equation*}
$\square$

It is now possible to state the following theorem about variations of $f_{S,d}$:
\begin{theorem}
\label{theo:increase_decrease}
Let $m, M$ and $d$ be the parameters of the model. Let $S=[m, M]$ and $f_{S,d}$ the relationship as defined above. Let $p=p(d)$ be the probability for two different nodes to be connected. Then, $f_{S,d}$ is increasing between 0 and $\frac{m+1}{p}-1$ and decreasing from $\frac{M+1}{p}-1$ to infinity.
\end{theorem}

\noindent \textbf{Proof}:
The goal is to prove that $\Delta f_{S,d}(n)$ is positive for $n < \frac{m+1}{p}-1$ and negative for $n > \frac{M+1}{p}-1$.
To understand this, $\Delta f_{S,d}(n)$ can be rewritten as follow (lemma \ref{lm:variation}):
\[
\forall n \in \mathds{N}, \Delta f_{S,d}(n) = 2\sum_{k=m}^M {(k+1){n \choose k}p^k(1-p)^{n-1-k}}\left(1-\frac{n+1}{k+1}p\right)
\]
It is sufficient to notice that, for all $k \in S$, the sign of every single term of the sum is the sign of $\left(1-\frac{n+1}{k+1}p\right)$.
For fixed $k$, the term is positive if and only if $n$ is lower than $\frac{k+1}{p}-1$.
As this last term is an increasing function of $k$,
all terms of the sum are therefore positive if $n$ is lower than $\frac{m+1}{p}-1$ and negative if $n$ is greater than $\frac{M+1}{p}-1$.
Hence, the relationship is increasing from 0 to $\frac{m+1}{p}-1$ and decreasing from $\frac{M+1}{p}-1$ to infinity.
$\square$\\

\noindent Thanks to theorem \ref{theo:increase_decrease}, conjecture \ref{cj:increase_decrease} is proved for intervals $[0,x_m]$ and $[x_M, \infty[$ $x_m=\frac{m+1}{p}-1$ and $x_M=\frac{M+1}{p}-1$.
At this stage, quantifying more precisely the evolution of the graph order is not achievable.
However, a study of the fixed points of $f_{S,d}$ enables to draw some conclusion about generated graphs sustainability.

\subsubsection{Graph evolution and sustainability}

First note that knowing the variations of $f_{S,d}$ is not enough to deal with graphs sustainability.
Indeed, as claimed by the following theorem, big graphs are not sustainable.
\begin{theorem}
\label{theo:big_graphs_collapse}
\underline{Non-sustainability of big graphs}:\\
Let $m$, $M$ and $d$ be parameters of the model.
Let $f_{S,d}$ be the relationship.
Then, there exists $N>0$ such that for all $n > N, f_{S,d}(n) < 1$.
\end{theorem}

\noindent \textbf{Proof}:
For this proof, it is sufficient to prove that $f_{S,d}(n) \to 0$ when $n \to +\infty$.
Let $n$ such that $n \geqslant 2M+1$.
In this situation, for all $k \leqslant M$, binomial coefficient ${n-1 \choose k} \leqslant {n-1 \choose M}$.
Moreover, as $(1-p) < 1$, $x \longmapsto (1-p)^x$ is decreasing.
Therefore, for all $k \leqslant M$, $(1-p)^{n-1-k} \leqslant (1-p)^{n-1-M}$.
It is thus possible to get the following inequality for all $k \leqslant M$:
\begin{equation*}
{n-1 \choose k}(1-p)^{n-1-k}p^k \leqslant {n-1 \choose M}(1-p)^{n-1-M}p^k
\end{equation*}
Noticing $p < 1$ and $f_{S,d}$ is a sum of $M-m+1$ elements, $f_{S,d}(n)$ can be bounded as follow
\begin{equation*}
f_{S,d}(n) \leqslant 2n(M-m+1)\left({n-1 \choose M}(1-p)^{n-1-M}\right)
\end{equation*}
As $M$ is fixed, the binomial coefficient ${n-1 \choose M}$ is equivalent to a polynomial of degree $M$ as $n$ grows to infinity:
\begin{equation*}
{n-1 \choose M} \sim \frac{n^M}{M!}
\end{equation*}
Therefore, $f_{S,d}(n)$ is equivalent to the product of a polynomial and an exponential function converging towards 0.
This implies $f_{S,d}(n)$ converges towards 0 as $n$ tends to infinity.
$\square$ \\
This theorem says that there always exists a graph order limit such that graphs whose order are greater than this limit are likely to become empty.
Therefore, it is not possible to obtain sustainable graphs with a large amount of nodes.

A new mathematical concept is now introduced aiming at classifying parameters into three classes.
This concept is referred to fixed point and is defined as follow:
\begin{definition}
\underline{Fixed Point}:\\
Let $m$, $M$ and $d$ be parameters of the model.
A fixed point for the relationship $f_{S,d}$ is an non-negative integer $n$ such that:
$$\left\{
  \begin{array}{ c l }
    & f_{S,d}(n) \leqslant n ~\mbox{and} ~f_{S,d}(n+1) > n+1 \\
    \text{or} & f_{S,d}(n) \geqslant n ~\mbox{and} ~f_{S,d}(n+1) < n+1
  \end{array}
\right.$$
\end{definition}
Such fixed points characterize variation of graph order.
Indeed, graph of order $n$ for $n$ taken between two consecutive fixed points is either always decreasing or increasing.
From experiment performed on the redistributed model as well as on D3G3, three different cases appear and are conjectured as follow:
\begin{conjecture}
\label{cj:three_fixed_point}
For all $m$, $M$ and $d$ being parameters of the model, the relationship $f_{S,d}$ has either one, two or three fixed points.
\end{conjecture}
This conjecture is the main tool aiming at studying sustainability in the segment case.
Indeed, in the three different cases, it is possible to answer whether a given set of parameters is sustainable or not.
However, their is no characterisation about parameters value that may help founding which case parameters lead to.
The only one claim that can be made is that $d$ does have an influence on this case.

The conjecture \ref{cj:three_fixed_point} is assumed in this subsection.
This section aims at stating about sustainability in the three different cases.
This is illustrated by a description of the behavior of the relationship $f_{S,d}$ in every case.

\subsubsubsection{One fixed point}

First let's consider the case where the relationship has only one fixed point.
When it has only one fixed point, this point is $0$.
This comes from $f_{S,d}(0) = 0$.
Moreover, for all $n$, $f_{S,d}(n) < n$.
As for a snapshot graph of order $n_t$ at step $t$, $f_{S,d}(n_t)$ gives the expectation value of $n_{t+1}$ at step $t+1$.
Graph orders of generated graphs are decreasing in average.
Graphs obtained in this case are therefore not sustainable.
This is illustrated by Figure \ref{sub:One_fixed_point}.
\begin{figure}[!h]
    \centering
    \includegraphics[width=0.9\textwidth]{NodeCount_ntp1_S_3___4_d_0.04.png}
     \caption{One fixed point}
     \label{sub:One_fixed_point}
\end{figure}

\subsubsubsection{Two fixed points}

For the two fixed points case, 0 is also a fixed point.
The argument is also because $f_{S,d}(0)=0$.
The other one is greater than zero.
The case where $f_{S,d}$ has two fixed points is assumed to happen if and only if $0 \in S$ and is stated in the following conjecture.
\begin{conjecture}{(Characterisation of the two fixed points)}
\label{cj:caracterisation_of_two_fp}
The relationship $f_{S,d}$ has two fixed points if and only if $m=0$.
\end{conjecture}
An argument is that a snapshot graph with one vertex becomes empty if and only if $0 \notin S$, that is $m=0$.
Moreover, for $n=0$, the first term of the sum defining $f_{S,d}$ is equal to 1 so $f_{S,d}(0) = 2$.
Graphs generated in such configurations are therefore sustainable as long as their graph order does not exceed a limit.
This limit is a consequence of theorem \ref{theo:big_graphs_collapse}.
In this case, graphs whose order exceeds the limit are likely to become empty.
This is illustrated by figure \ref{sub:Two_fixed_points}.
\begin{figure}[!h]
    \centering
    \includegraphics[width=0.9\textwidth]{NodeCount_ntp1_S_0___5_d_0_04.png}
     \caption{Two fixed points}
     \label{sub:Two_fixed_points}
\end{figure}

\subsubsubsection{Three fixed points}

For the last case, the goal is to show that graph order is likely to remain bounded.
Deeply looking at this case raises the question of values of graph order for which the size is not too large and not too small so that it does not collapse.
{For that purpose we define an interval, called \textit{sustainable interval}, such that, if the graph order remains within that interval, this ensures the persistence of the graph.}
This sustainable interval is considered as a tool to study graph sustainability.
It concerns expectation of graph order evolution through time.
It says that if the image of the function $f_{S,d}$ for all integers within the interval does not exceed the upper bound, then the graph is likely not to collapse.
Let's define more precisely this concept:
\begin{definition}
\underline{sustainable interval}:\\
Let $m$, $M$ and $d$ be parameters of the model.
Let consider $f_{S,d}$ set such that it has three fixed points.
Let $N_m$ be the first positive fixed point and  $N'_m$ the smallest integer greater than $N_m$ such that $f_{S,d}(N'_m) \geqslant N_m$ and $f_{S,d}(N'_m + 1) < N_m$.
The \textit{sustainable interval} associated to $m$, $M$ and $d$ is defined as the interval $[N_m, N'_m]$.
\end{definition}

This definition is illustrated through figure \ref{sub:Three_fixed_points}.
\begin{figure}[!h]
    \centering
    \includegraphics[width=0.9\textwidth]{NodeCount_ntp1_S_3___8_d_0_02.png}
     \caption{Three fixed points}
     \label{sub:Three_fixed_points}
\end{figure}
Such an interval satisfies a property about the values $f_{S,d}$ takes when it is restricted to it:
\begin{theorem}
\underline{Sustainability in the sustainable interval}:\\
Let $m$, $M$ and $d$ be parameters of the model.
Let assume the relationship $f_{S,d}$ has three fixed points and that $[N_m, N'_m]$ is its associated sustainable interval.
Then, the relationship satisfies:
\begin{equation*}
\forall n \in [N_m, N'_m], f_{S,d}(n) \geqslant N_m
\end{equation*}
Moreover, if the relationship does not exceed $N'_m$, then the relationship satisfies:
\begin{equation*}
\forall n \in [N_m, N'_m], f_{S,d}(n) \in [N_m, N'_m]
\end{equation*}
\end{theorem}

\noindent \textbf{Proof:}
The first part of the theorem comes from definition of the sustainable interval.
$\square$\\

\noindent Main interpretation of that theorem is graphs are sustainable in probability in the sustainable interval if and only if there are no values of $f_{S,d}$ that exceed the upper bound of the sustainable interval.\\

The following paragraphs provides arguments aiming at obtaining the sustainable interval.
They also provide arguments to check whether the relationship exceeds the upper bound of the interval.
The theorem \ref{theo:increase_decrease} clearly gives bounds to find out the maximum of the relationship $f_{S,d}$.
Three algorithms are sufficient to answer both questions:
an algorithm to compute the argument of the maximum of the relationship $f_{S,d}$, an algorithm to find its fixed point between 0 and the argument of the maximum and an algorithm to solve $f_{S,d}(n)=y$ for $n$ greater than the argument of the maximum and $y>0$ lower than or equal to the maximum.
In the following, these algorithms are first implemented.
It is then explained how to use them to answer questions about the sustainable interval.\\

\noindent \underline{\textbf{The argument maximum}}: To compute the argument maximum of the relationship, it is sufficient to study $f_{S,d}$ on the interval $[x_m, x_M]$ for $x_m$ and $x_M$ as defined above.
This is a consequence of theorem \ref{theo:increase_decrease}. Let's denote it $N_*$.\\

\noindent \underline{\textbf{The first positive fixed point}}: To find the fixed point of $f_{S,d}$ mentioned in the definition of the sustainable interval, it is sufficient to compute the argument maximum of it. The previous algorithm answers this question.
Then, as the relationship is increasing from 0 to $N_*$, it is sufficient to iterate and find an integer $n$ such that $f_{S,d}(n) \leqslant n$ and $f_{S,d}(n + 1) > n + 1$.\\

\noindent \underline{\textbf{The solution of the equation}}: For the last algorithm, the goal is to find an integer $n$ such that $n$ is greater than $N_*$ of $f_{S,d}$, $f_{S,d}(n) \geqslant y$ and $f_{S,d}(n) < y$, for a fixed $y$ which is assumed to be positive and lower than the maximum of $f_{S,d}$.\\

\noindent From these algorithms it is possible to implement algorithms stating the existence of the sustainable interval and its bounds.
For the existence or not of the sustainable interval, it is sufficient to check whether the maximum of the relationship is greater than its argument.
This comes from that sustainable interval exists if and only if there are values of the relationship that exceed their argument.
As the relationship is increasing from 0 to $f_{S,d}(N_*)$, then  sustainable interval exists if and only if $f_{S,d}(N_*) > N_*$.
For computing the sustainable interval boundaries, it is sufficient to know the value of the first fixed point $N_m$ (as it provides the lower bound) and to solve the equation $f_{S,d}(x)=N_m$ as finding the corresponding $x$ to this equation provides the upper bound ($N'_m$).
The existence of $N'_m$ is ensured by theorem \ref{theo:big_graphs_collapse}.

\section{Introduction}
\label{sec:intro}

\input{0_introduction}

\section{Definitions and Generative Model and Definitions}
\label{sec:model}

\input{1_model-and-definitions}

\section{Theoretical Analysis}
\label{sec:theoreticalanalysis}

\input{2_analysis}

\subsection{General Cases}
\label{sec:generalcases}

\input{3_general-cases}

\section{Segments}
\label{sec:segments}

\input{4_segments}

\section{Conclusion}
\input{6_conclusion}

\bibliographystyle{plain}


\section{Introduction}
Suppose that $W$ is a non-void set endowed with a metric $d$. The set $CB(W)$ represents the collection of all non-void bounded closed subsets of $W$. A map $H:CB(W)\times CB(W) \to \mathbb{R}$ is described as a Hausdorff metric on $CB(W)$, whenever for $Y$ and $Z$ lies in $CB(W)$,
\begin{center}
	$\displaystyle{	H(Y,Z)= \max\left\{ \sup_{u\in Y}  D(u,Z),\,\,\sup_{v\in Z}D(v,Y)\right\}}$,
\end{center}
whereas the notation $D(u,Z)$ describes the $\inf_{v\in Z}d(u,v)$. Further, $CL(W)$ expresses the collection of all non-void subsets of $W$, which are closed. Given a set-valued map $F:W\to CL(W)$, an element $w^*\in W$ is stated as a fixed point for $F$, if it lie in $\{w\in W: w\in Fw\}$. The occurrence of fixed points for a set-valued contraction first derived by Nadler \cite{nadler} in the year of 1969.
Later, Mizoguchi and Takahashi \cite{mizo} improved the Nadler's \cite{nadler} result for a set-valued non-linear contraction in the year of 1989. The authors \cite{mizo} deduced that if $F:W\to CB(W)$ and for each $v,w$ in $W$ fulfils $H(Fv, Fw)\leq k(d(v,w))d(v,w)$, where $k \in U=\{h:[0,\infty)\to [0,1)| \limsup_{r\to t^+}h(r)<1,\forall \,0\leq t<\infty\}$, then occurrence of a fixed point for the map $F$ is guaranteed whenever $(W,d)$ is complete. There are several generalization of this result in the literature, which can be found in \cite{berinde,mam, mam_feng}. In the year 2006, Feng and Liu \cite{feng} extended the well known Nadler's \cite{nadler} result for a set-valued generalized contraction $F$ from $W$ into $CL(W)$, where the Hausdorff distance function is not involved.
\begin{theorem}\cite{feng}\label{feng}
	For the metric space $(W,d)$, a set-valued map $F:W\to CL(W)$ in order that for each $v\in W$, $w\in Fv$ fulfils
	\begin{equation}
		D(w,Fw)\leq \theta d(v,w)~\textrm{where}~0\leq \theta<1 .
	\end{equation}
	Then there is an element $w^*\in Fw^*$ if the map $h: W\to \mathbb{R}$ given by $h(v)=D(v,Fv)$ is lower semicontinuous and $(W,d)$ is complete.
\end{theorem}

For some more results about the fixed points for set-valued generalized contraction can be found in \cite{azam,klim}.

In the parallel development, many mathematicians scrutinized common fixed point results for a couple of single valued mappings on a metric space, one may see \cite{abbas, Doric}, whereas in \cite{khan,kumam}, the authors studied the common fixed point theory for a couple of single valued and set-valued mappings.
In the year 2007, Kamran \cite{kamran} established a common fixed point result for a single valued and set-valued generalized contractive mappings. This result actually enables the author \cite{kamran} to generalize several fixed point and common fixed point results on metric space available in the literature. Indeed, the author \cite{kamran} established the below stated result.
\begin{theorem}\cite{kamran}\label{thm1}
	For the metric space $(W,d)$, the maps $F:W\to CB(W)$, $f:W\to W$ in order that for each $u\in W$, $Fu\subset fW$ and
	\begin{equation}\label{eqn1}
		H(Fv,Fw)\leq k(d(fv,fw))d(fv,fw)+M\,D(fv,Fw)~\text{for each}~ v,w \in W,
	\end{equation}
	where $k\in U$, $M\geq 0$. Then occurrence of $w^*$ having $fw^*\in Fw^*$ is guaranteed if $fW$ is complete. Moreover, $w^*=fw^*\in Fw^*$, whenever $ffw^*=fw^*$ and $\{ffw^*\}\subset Ffw^*$.
\end{theorem}
\noindent We note that, the above mentioned result deduces Mizoguchi-Takahashi's \cite{mizo} theorem by considering $f$ to be an identity map and $M=0$.

On the other side, Jachymski \cite{Jach} extended the renowned Banach contraction principle via $G$-contraction on a metric space having  graph $G$. With the aid of this result, the author \cite{Jach} enables to extend and unify some fixed point results proved in metric spaces as well as metric spaces having partial order.
Followed by Jachymski \cite{Jach}, an extension of the Theorem \ref{feng} due to Feng and Liu \cite{feng} was derived by Sultana and Vetrivel \cite{mam_feng}  in the year 2017, on a metric space consisting with directed graph $G$.
The authors \cite{mam_feng} demonstrated the below stated result.
\begin{theorem}\cite{mam_feng}\label{mam}
	Let $(W,d)$ along with a directed graph $G$ and a map $F:W\to CL(W)$ in order that for every $v \in W$, $w\in Fv$ with $(v,w)\in E(G)$ fulfils
	\begin{enumerate}
		\item [(i)] $D(w,Fw)\leq k(d(v,w))d(v,w)$ for $k\in U$,
		\item [(ii)] for any $p\in Fw$ having $d(p,w)\leq d(v,w)$ indicates $(p,w)\in E(G)$.
	\end{enumerate}
	Suppose that $w_0\in W$,  $w_1\in Fw_0$ with $(w_0,w_1)\in E(G)$. Moreover, assume that if any sequence $\{w_n\}_n\subset W$ with $(w_{n-1},w_{n})\in E(G)$ for every $n$, and $w_n\to w$ then there is $\{w_{n_r}\}_{r\in \mathbb{N}}$ in order that $D(w,Fw)\leq \liminf_{r\to \infty}D(w_{n_r},F(w_{n_r}))$. Then there is an element $w^*\in Fw^*$ if $(W,d)$ is complete.
\end{theorem}
\noindent As an application of the Theorem \ref{mam}, the authors \cite{mam_feng} improved the Kelisky and Rivlin \cite[Theorem 1]{KR8} result about the convergence of iterates for Bernstein operator on a complete normed linear space. In the literature, there are several fixed point theorems using the notion of graph theory, which can be found in \cite{frigon,khan, mam}.

Inspired by Kamran \cite{kamran} and Sultana-Vetrivel \cite{mam_feng},
in this present manuscript, we derive a common fixed point result on a metric space having graphical structure for a set-valued mapping, which is supposed to be closed valued instead of bounded and closed valued.
The main result is indeed an extension of the Theorem \ref{thm1} due to Kamran \cite{kamran} and the Theorem \ref{mam} due to Sultana-Vetrivel \cite{mam_feng}.
As an implementation of our result,  we establish the convergence of iterates for some nonlinear operators, which enables us to derive the convergence of iterates for a nonlinear $q$-analogue Bernstein operator. As an another application, we derive a sufficient criteria for the presence of a solution for a fractional order differential equation. An invariant best approximation result is also discussed in this article.
\bigskip
\section{Preliminaries}
In this part, we utilize some essential symbols and definitions, which is needful all throughout this article.
An element $w^*$ is described as a coincidence point for a couple of mappings $F:W\to CL(W)$ and $f:W\to W$, if $\{fw^*\}\subset Fw^*$. Furthermore, an element $w^*$ is described as a common fixed point for the map $f$ coupled with $F$ if $w^*=fw^*\in Fw^*$. Additionally, the map $f$ is stated as a $F$-weakly commuting \cite{kamran} at an element $z^*$ in $W$ if $ffz^*\in Ffz^*$. Up to the end, the collection of common fixed points and coincidence points for the maps $f$ and $F$ are expressed by $Fix(f,F)$ and $Coin(f,F)$ respectively.

Now we present some fundamental terminology of graph theory (see \cite{John}), which is mandatory up to the end of this article. For the metric space $(W,d)$, let us assume a directed graph $G$ which is a pair of $(V(G),E(G))$, whereas the collection of all vertices are symbolized by $V(G)$, which is actually whole set $W$ and the collection of edges $E(G)$ includes the set $\Delta=\{(s,s):(s,s)\in W\times W\}$. Further, this reflexive graph $G$ does not possess any parallel edges.

\section{Main results}
Up to the end of this part, we consider a metric space $(W,d)$ along with a directed graph $G(V(G),E(G))$, whereas $V(G)=W$, $\Delta \subseteq E(G)$ and $G$ has no parallel edges. In the below stated theorem, we scrutinize common fixed points on a metric space having a graph $G$ for a set-valued generalized contraction, which is assumed to be only closed-valued.
\begin{theorem}\label{thm3.2}
	For the metric space $(W,d)$, the maps $F:W \to CL(W)$, $f:W\to W$ in order that for each $v\in W$, $fw\in Fv$ with $(fv,fw)\in E(G)$ fulfils,
	\begin{enumerate}\label{main}
		\item[(i)] $D(fw,Fw)\leq k(d(fv,fw))d(fv,fw)$, where $k\in U$,
		\item[(ii)] for any $fp\in Fw$ with $d(fw,fp)\leq d(fv,fw)$, indicates $(fw,fp)\in E(G)$.
	\end{enumerate}
	Suppose that $w_0\in W$, $p_0\in F(w_0)$ with $(fw_0,p_0)\in E(G)$. Also assume that if any $\{fu_n\}_{n}\subset W$ with $(fu_{n-1},fu_{n})\in E(G)$, $fu_n\in Fu_{n-1}$ for every $n\in \mathbb{N}$ and $fu_n\to fu$, then there is $\{fu_{n_r}\}\subseteq\{fu_n\}$ in order that $D(fu,Fu)\leq \liminf_{r\to \infty}D(fu_{n_r},Fu_{n_r})$. Then there is a sequence $\{fw_n\}_{n}$ with $fw_n\in Fw_{n-1}$ converges to $fw^*$, where $fw^*\in Fw^*$
	is guaranteed if for each $u\in W$, $Fu\subset fW$ and $(fW,d)$ is complete.

	Moreover, $f$ and $F$ possess a common fixed point whenever the map $f$ is weakly commuting at $w^*$ along with $fw^*=ffw^*$.

\end{theorem}

	\begin{proof}
		As $p_0\in Fw_0$ and $Fw_0\subset fW$, then there is $w_1\in W$ in order that $p_0=fw_1$. Consequently, we obtain $(fw_0,fw_1)\in E(G)$ and $fw_1\in F(w_0)$.

		Since $D(fw_1,Fw_1)=\inf\{d(fw_1,y): y\in Fw_1\}$, then for $\frac{1}{\sqrt{k(d(fw_0,fw_1))}}>1$ there is $w_2\in W$ with $fw_2\in F(w_1)$ fulfils
		\begin{eqnarray}\label{eqn_1st}
			d(fw_1,fw_2)&\leq& D(fw_1,Fw_1)+ \left(\frac{1}{\sqrt{k(d(fw_0,fw_1))}}-1\right)D(fw_1,Fw_1)\nonumber\\&\leq&\frac{1}{\sqrt{k(d(fw_0,fw_1))}} k(d(fw_0,fw_1))d(fw_0,fw_1)\nonumber\\

			&\leq&\sqrt{k(d(fw_0,fw_1))}d(fw_0,fw_1)< d(fw_0,fw_1).
		\end{eqnarray}

		Due to the fact that $F$ meets the hypothesis $(ii)$ and $(fw_0,fw_1)\in E(G)$ with $fw_1\in Fw_0$, then the previous inequation $(\ref{eqn_1st})$ leads to $(fw_1,fw_2)\in E(G)$.

		Repeating the above methodology, for $\frac{1}{\sqrt{k(d(fw_1,fw_2))}}>1$ we get $fw_3\in F(w_2)$ fulfilling
		\begin{eqnarray*}
			d(fw_2,fw_3)&\leq& D(fw_2,Fw_2)+\left( \frac{1}{\sqrt{k(d(fw_1,fw_2))}}-1\right)D(fw_2,Fw_2)\nonumber\\
			&\leq& \frac{1}{\sqrt{k(d(fw_1,fw_2))}}k(d(fw_1,fw_2))d(fw_1,fw_2)< d(fw_1,fw_2

		\end{eqnarray*}
		Hence by the same argument as before we have $(fw_2,fw_3)\in E(G)$. Following this approach, we are able to construct two sequences $\{f(w_n)\}_{n}$, $\{F(w_n)\}_{n}$ in order that for every $n$, $(fw_{n-1},fw_n)\in E(G)$ and $f(w_n)\in F(w_{n-1})$ fulfils
		\begin{equation}\label{eqn2}
			d(fw_n,fw_{n+1})

			\leq\sqrt{k(d(fw_{n-1},fw_n))}d(fw_{n-1},fw_n)<d(fw_{n-1},fw_n).

		\end{equation}
		Now for $n\geq 1$, we consider $s_n=d(fw_n,fw_{n+1})$. Consequently from $(\ref{eqn2})$, it is clearly observed that for each $n$, $s_n$ are non-negative and $s_{n+1}<s_n$.

		Therefore $\{s_n\}_n$ is convergent and assume $s_n \to q$, where $q\geq 0$. On the account of $\limsup_{t\to q^+}k(t)<1$, there is $0\leq \alpha<1$ and $M\geq 1$ meeting the condition $k(s_n)\leq \alpha$ for each $n\geq M$. Now for every $n> M$, the inequation $(\ref{eqn2})$ leads to
		\begin{eqnarray}\label{bounded}
			d(fw_n,fw_{n+1})&\leq&\alpha^{(n-M)/2}\prod_{j=1}^{M-1}\sqrt{k(d(fw_{j-1},fw_{j}))}d(fw_{0},fw_1)\nonumber\\

			&<&B\alpha^{n/2}d(fw_{0},fw_1),
		\end{eqnarray}
		where $B=\alpha^{-{M/2}}\prod_{j=1}^{M-1}\sqrt{k(d(fw_{j-1},fw_{j}))}$.

		Now using $(\ref{bounded})$ and for any arbitrary positive integer $m$ and $n>M$,
		\begin{equation*}
			\begin{aligned}
				d(fw_n,fw_{n+m})&\leq d(fw_{n},fw_{n+1})+d(fw_{n+1},fw_{n+2})+\cdots+\\
				& \qquad d(fw_{n+m-1},fw_{n+m})\\

				&\leq B\alpha^{n/2}\left[1+\alpha^{1/2}+ \cdots+ \alpha^{(m-1)/2}\right]d(fw_{0},fw_1)\\

				&\leq B\frac{\alpha^{n/2}}{1-\alpha^{1/2}}d(fw_{0},fw_1).
			\end{aligned}
		\end{equation*}
		Thus the sequence $\{fw_n\}_{n}\subset fW$ turns out to be Cauchy, on the account of $\alpha^{n/2}\to 0$ as $n\to \infty$.

		Since $fW$ is complete, then $\{fw_n\}_{n}$ is converges to $fw^*$, for some $w^*\in W$.
		Again for every $n$, $(fw_{n-1},fw_n)\in E(G)$ and $fw_n\in Fw_{n-1}$, it yields from $(i)$ that $D(fw_n,Fw_n)\leq k(d(fw_{n-1},fw_n))d(fw_{n-1},fw_n)<d(fw_{n-1},fw_n)$. Hence $\lim_{n \rightarrow \infty}D(fw_n,Fw_n)=0$.

		On the account of $(fw_{n-1},fw_n)\in E(G)$ with $fw_{n}\in Fw_{n-1}$ for each $n$ and $fw_n\to fw^*$, consequently by the hypothesis there is a subsequence $\{fw_{n_r}\}_{r\in \mathbb{N}}$ having  $D(fw^*,Fw^*)\leq \liminf_{r\to \infty}D(fw_{n_r},F(w_{n_r}))$. Thus we can conclude that $D(fw^*,Fw^*)=0$. Due to the fact that $Fw^*$ is closed, therefore $fw^*\in Fw^*$.

		Now let $a=fw^*$, then the hypothesis $fw^*=ffw^*$ indicates that $fa=a$. Again $f$ and $F$ are weakly commuting at $w^*$, hence $ffw^*\in Ffw^*$ leads to $a=fa\in Fa$.
	\end{proof}

	\begin{remark}
		Take a graph $G$ on the set $W$ in order that $V(G)=W$ and $E(G)= W\times W$. Consequently, the aforementioned Theorem \ref{feng} due to Feng-Liu \cite{feng} is followed from our main Theorem \ref{thm3.2} by considering the map $f$ is an identity map on $W$ and for each $0\leq t<\infty$, $k(t)=\theta$, where $\theta\in [0,1)$.
	\end{remark}

	The succeeding corollary deals with common fixed points for set-valued generalized contraction in a metric space, which is indeed an extension of the Theorem \ref{feng} due to Feng-Liu \cite{feng}. Moreover, this result is followed from the aforementioned Theorem \ref{thm3.2} by choosing a graph $G$ on $W$ in order that $V(G)=W$ and $E(G)=W\times W$.
	\begin{corollary}\label{cor}
		For the metric space $(W,d)$, the maps $F:W \to CL(W)$, $f:W\to W$ in order that for each $v\in W$, $fw\in Fv$ fulfils,
		\begin{equation*}
			D(fw,Fw)\leq k(d(fv,fw))d(fv,fw), ~\textrm{where}~ k\in U.
		\end{equation*}
		Suppose that the function $h: W\to \mathbb{R}$ by $h(v)=D(fv,Fv)$ is lower semi-continuous. Then occurrence of $w^*$ having $fw^*\in Fw^*$ is guaranteed if for each $u\in W$, $Fu\subset fW$ and $(fW,d)$ is complete.

		Moreover, $f$ and $F$ possess a common fixed point whenever the map $f$ is weakly commuting at $w^*$ along with $fw^*=ffw^*$.
	\end{corollary}

	The below-stated example illustrates the preceding Theorem \ref{thm3.2}.

	\begin{example}\label{example}
		Consider a set $W=\{0,1\}\cup\left\{\frac{1}{3^n}:n\in \mathbb{N}\right\}$ together with usual metric. Suppose that a single valued map $f:W\to W$ by
		\begin{center}
			$f(w)=
			\begin{cases}
				~\frac{1}{3^{n+1}} &\mbox{if}~~w=\frac{1}{3^n},~n\in \mathbb{N}\cup\{0\}\\
				~0~&\mbox{if}~~w=0.
			\end{cases}$
		\end{center}
		and set-valued map $F: W\to CL(W)$ is defined as
		\begin{center}
			$F(w)=
			\begin{cases}
				~\left\{0,\frac{1}{3}\right\} &\mbox{if}~~w=0,\\
				~\left\{\frac{1}{3},\frac{1}{3^{n+2}}\right\} &\mbox{if}~~w=\frac{1}{3^n},~\text{where}~ n\in \mathbb{N},\\
				~\left\{0\right\}~&\mbox{if}~~w=1,
			\end{cases}$
		\end{center}
		Take a graph $G$ on $W$ in order that $V(G)=W$ and $E(G)=\{(v,w)\in W\times W: d(v,w)< \frac{1}{9}\}$. Then, $G$ does not have any parallel edges and $E(G)$ contains the set $\Delta$. Now for $v=\frac{1}{3^n}$ and $fw=\frac{1}{3^{n+2}}\in Fv$, we see that for each $n$, $(fv,fw)\in E(G)$. Further,
		\begin{equation*}
			D(fw,Fw)=D\left(\frac{1}{3^{n+2}},\left\{\frac{1}{3},\frac{1}{3^{n+3}}\right\}\right)=\frac{1}{3^{n+2}}-\frac{1}{3^{n+3}}\leq \frac{1}{3}d(fv,fw).
		\end{equation*}
		Again, for $v=0$ and $fw=0\in Fv$, the pair $(fv,fw)\in E(G)$. Also, it is easy to check that $D(fw,Fw)\leq \frac{1}{3}d(fv,fw)$. Thus we conclude that for each $v\in W$, $fw\in Fv$ with $(fv,fw)\in E(G)$ fulfils $D(fw,Fw)\leq k(d(fv,fw))d(fv,fw)$, where $k(t)=\frac{1}{3}$ for each $t\in [0,\infty)$. Furthermore, if $fp\in Fw$ with $d(fw,fp)\leq d(fv,fw)$, then $(fw,fp)\in E(G)$.
		Now for choosing $w_0=\frac{1}{3}$, we observe that $\frac{1}{27}\in Fw_0$ and $(fw_0,\frac{1}{27})\in E(G)$. It is clear to observe that the map $h:W\to \mathbb{R}$ by $h(w)=D(fw,Fw)$ is lower semi-continuous. Further, $(fW,d)$ is complete and for each $u\in W$, $Fu\subset fW$. Thus all the criteria of the Theorem \ref{thm3.2} is fulfilled, which confirm the occurrence of coincidence point $w^*$. We visualize that $0$ is a coincidence point of $f$ and $F$, that is, $f0\in F0$. Again $ff0=f0$ and $ff0\in Ff0$. We note that $0$ is a common fixed point of $f$ and $F$.

		On the other side, if we choose $v=\frac{1}{3^n}$ for each $n\geq 1$ and $fw=\frac{1}{3}\in Fw$, then $d(fv,fw)=\frac{1}{3}-\frac{1}{3^{n+1}}$ and $D(fw,Fw)=\frac{1}{3}$. Therefore for such chosen points $D(fw,Fw)>d(fv,fw)$. This indicates that the Theorem \ref{thm3.2} is indeed an extension of the Corollary \ref{cor}.
	\end{example}

	\begin{remark}
		It is worth to observe that the preceding Theorem \ref{thm3.2} is an extension of the Theorem \ref{thm1} due to Kamran \cite{kamran}. Indeed, consider $W$ incorporates with a graph $G$ having $V(G)=W$ and $E(G)=W\times W$. Let $v\in W$, $fw\in Fv$ with $(fv,fw)\in E(G)$, then it follows
		\begin{eqnarray*}
			D(fw,Fw)\leq H(Fv,Fw)&\leq& k(d(fv,fw))d(fv,fw)+M\,D(fw,Fv)\\
			&=& k(d(fv,fw))d(fv,fw),~\textrm{where}~ k\in U.
		\end{eqnarray*}
		Moreover, for the chosen graph $G$, it always follows that $(fw,fp)\in E(G)$.
		Also we are able to find $w_0,p\in W$ in order that $(fw_0,p)\in E(G)$ and $p\in Fw_0$. Now assume a $\{fw_n\}_n\in W$ in order that $fw_n\in Fw_{n-1}$ for every $n$, and $fw_n\to fw$. Evidently for each $n$, $(fw_n,fw_{n+1})\in E(G)$. Consequently for each $n$ and $fz\in Fw_n$,
		\begin{equation*}
			\begin{aligned}
				D(fw,Fw
				&\leq  d(fw,fw_n)+d(fw_n,fz)+H(Fw_n,Fw)\\
				&\leq d(fw,fw_n)+ D(fw_n,Fw_n)+ d(fw,fw_n)+ M\,D(fw,Fw_n)\\

				&\leq d(fw,fw_n)+ D(fw_n,Fw_n)+ d(fw,fw_n).~ \textrm{[as $fw_n\in F(w_{n-1})]$}
			\end{aligned}
		\end{equation*}

		Now taking $n\to \infty$ it yields, $D(fw,Fw)\leq\liminf_{n\to \infty} D(fw_n,Fw_n)$. Hence all the criteria of the Theorem \ref{thm3.2} is fulfilled.
	\end{remark}

	The below-mentioned example demonstrates that, the Theorem \ref{thm3.2} indeed an extension of Theorem \ref{thm1} due to Kamran \cite{kamran}.

	\begin{example}
		In the Example \ref{example}, choose two points $0$ and $1$ in $W$. Then we see that $H(F0,F1)=\frac{1}{3}$, $d(f0,f1)=\frac{1}{3}$ and $D(f1,F0)=0$. Therefore for all $k\in U$ and $M\geq 0$, $H(F0,F1)=d(f0,f1)> k(d(f0,f1))d(f0,f1)+M\,D(f1,F0)$. This implies that $F$ does not fulfils the equation (\ref{eqn1})  of the Theorem \ref{thm1} due to Kamran.
	\end{example}

	\begin{remark}
		The Theorem \ref{mam} due to Sultana and Vetrivel \cite{mam_feng} follows
		from our Theorem \ref{thm3.2} by choosing $f$ to be an identity map.
	\end{remark}

	Let the set $W$ be endowed with a norm $\|.\|$ and $Q\subseteq W$ with $Q\neq \emptyset$. Now for $z\in W$, a set $B_Q(z)$ contains $y\in Q$ in order that $D(z,Q)=\|y-z\|$. Then the set $B_Q(z)$ is called as a best $Q$-approximates of $z$ over $Q$. Invariant best approximation problems can be solved through the perception of common fixed points for single valued mappings, as one may see in \cite{thagafi,pathak}. Also in \cite{kamran1,kumam}, the authors established invariant best approximation results for a couple of single valued and set-valued mappings. In the succeeding theorem, we investigate invariant best approximation for single valued as well as set-valued maps meeting the condition $(i)$ of the Theorem \ref{thm3.2} on a normed linear space without considering the Hausdorff distance function.

	\begin{theorem}
		Suppose that $(W,\|.\|)$ is a normed linear space and $Q\subseteq W$ such that $Q\neq \emptyset$. Assume that the mappings $f:Q\to Q$ and $F:Q \to CL(Q)$ fulfils the following:
		\begin{enumerate}
			\item[(i)] for every $v\in B_Q(z)$ and $fw\in Fv$ fulfils $D(fw,Fw)\leq k(\|fv-fw\|)\|fv-fw\|$ where $k \in U$,

			\item[(ii)] $f(B_Q(z))=B_Q(z)$ and $B_Q(z)$ is complete,

			\item[(iii)] $\sup_{u\in Fp}\|u-z\|\leq \|fp-z\|$ for each $p\in B_Q(z)$.
		\end{enumerate}
		Then $w^*\in Coin(f,F)\cap B_Q(z)$ if for any $\{fw_n\}_n\in B_Q(z)$ with $fw_n\in Fw_{n-1}$ and $fw_n\to fw$ implies $D(fw,Fw)\leq \liminf_{n\to \infty}D(fw_n,Fw_n)$. Furthermore, $Fix(f,F)\cap B_Q(z)\neq \emptyset$ if $ffw^*=fw^*$ and $ffw^*\in Ffw^*$.
	\end{theorem}

	\begin{proof}
		Consider a graph $G$ on $B_Q(z)$ in order that $V(G)=B_Q(z)$ and $E(G)=B_Q(z)\times B_Q(z)$. Now assume an element $p$ lies in $B_Q(z)$ and $w\in Fp$. Then we have $fp\in B_Q(z)$ for every $p\in B_Q(z)$, on the account of $f(B_Q(z))=B_Q(z)$.
		Therefore we get $\|fp-z\|=D(z,Q)$, by using the definition of $B_Q(z)$. Again the condition $(iii)$ yields that $$\|w-z\|\leq \sup_{u\in Fp}\|u-z\|\leq \|fp-z\|=D(z,Q).$$

		Hence $w\in B_Q(z)$. This indicates that  for each $p \in B_Q(z)$,  $Fp\subseteq B_Q(z)$. As the set $Fp$ is closed for each $p\in Q$, then $Fp$ is also closed for any $p \in B_Q(z)$. Subsequently, $F|_{ B_Q(z)}$ is a set-valued mapping from $B_Q(z)$ into $CL( B_Q(z))$ and $f|_{ B_Q(z)}: B_Q(z)\to  B_Q(z)$. Then eventually
		$$Fix(f|_{ B_Q(z)},F|_{ B_Q(z)})=Fix(f,F)\cap B_Q(z).$$
		Thus this theorem is followed by applying our main Theorem \ref{thm3.2} under the assumption $W=B_Q(z)$.
	\end{proof}

	Jachymski \cite[Theorem 4.1]{Jach}, established the convergence of iterates for some linear operator with the perception of fixed point theory on a normed linear space, which is indeed an extension of Kelisky-Rivlin \cite[Theorem 1]{KR8} theorem about the convergence of iterates for Bernstein operator \cite{KR8}. Later, Sultana and Vetrivel \cite[Theorem 5]{mam_feng} improved this result on a complete normed linear space. The coming theorem deals with the convergence of iterates for certain operators in a normed linear space, where the operators need not to be linear. Essentially, our theorem generalizes the Sultana and Vetrivel's \cite[Theorem 5]{mam_feng} result.

	\begin{theorem}\label{iterates}
		For the complete normed linear space $(W,\|.\|)$, let $W_0\subset W$ be a closed subspace. A map $T:W\to W$ with $(I-T)(W)\subseteq W_0$ and for $w-Tw\in W_0$,
		\begin{equation}\label{convergence}
			\|Tw-T^2w\|\leq k(\|w-Tw\|)\|w-Tw\|\quad\text{where}~k\in U.
		\end{equation}
		Then for every $w\in W$, $\lim_{n \rightarrow \infty}T^nw\in \{w^*\in W: w^*=Tw^*\}\cap (w+W_0)$ if $h:W\to \mathbb{R}$ by $h(u)=D(u,Tu)$ is lower semicontinuous.
	\end{theorem}

	\begin{proof}
		Take a graph $G$ on $W$ in order that $V(G)=W$ and $E(G)$ is the collection of $(v,w)\in W\times W$ having $v-w\in W_0$. Consider $f$ is an identity map on $W$ and a set valued map $F$ is chosen by $Fw=\{Tw\}$, for each $w\in W$. Let $v$ be any arbitrary point in $W$ and $w=fw\in Fv=\{Tv\}$ with $(fv,fw)=(v,w)\in E(G)$. Hence $(v,Tv)\in E(G)$, that is $v-Tv\in W_0$. Consequently by (\ref{convergence}), $\|fw-Fw\|=\|Tv-T(Tv )\|=\|Tv-T^2v\|\leq k(\|v-Tv\|)\|v-Tv\|$ where $k\in U$. On the account of $(I-T)(W)\subseteq W_0$, it yields that for $p=fp \in Fw=\{Tw\}$, we have $w-p=w-Tw\in W_0$, that is $(w,p)=(fw,fp)\in E(G)$. Hence condition (i) and (ii) of Theorem \ref{thm3.2} is fulfilled. Choose an element $w\in W$. Moreover, for $w=w_0\in W$, there is $p_0\in F(w_0)$ with $(fw_0,fp_0)\in E(G)$ due to the fact that $(I-T)(W)\subseteq W_0$.
		Thus for $w\in W$, applying the previous Theorem \ref{thm3.2} we obtain that $\lim_{n\to \infty}T^nw\in \{w^*\in W: w^*=Tw^*\}$.

		As $(I-T)W\subseteq W_0$, then it is easy observe that $w-T^nw\in W_0$ for each $n \geq 1$. On the account of $W_0$ is closed, we obtain $\lim_{n \rightarrow \infty}T^nw\in (w+W_0)$. Hence we conclude that  $\lim_{n \rightarrow \infty}T^nw\in \{w^*\in W: w^*=Tw^*\}\cap (w+W_0)$.
	\end{proof}

	\section{Applications}
	In this part we provide some applications of our main Theorem \ref{thm3.2}. First we discuss about the occurrence of solution of a fractional differential equation and at the end of this article, we deduce the convergence of iterates for a nonlinear $q$-analogue Bernstein operator.

	A fractional differential equation is a generalization of a differential equation, whose order lies in the set of natural numbers. A fractional derivative looks like $D^\beta$, where $\beta>0$ and the symbol ${}^cD^\beta$ denotes the Caputo fractional derivative \cite{baleanu} with order $\beta$. In fact if $\beta \in \mathbb{N}$, then fractional differential equation turns out to be a normal differential equation. The presence of solution for fractional differential equations are deduced through the perception of fixed point theory, (see \cite{baleanu,khan,gopal}).
	In this article, we assume a generalized fractional differential equation
	\begin{equation}\label{eqn6}
		{}^cD^\beta w(b)+g(b,f(w(b)))=0~ \textrm{whenever}~~ w\in C[0,1],~ 0<b<1~\textrm{and}~1<\beta,
	\end{equation}
	having boundary criteria $w(0)=0$ and $w(1)=0$, where the maps $f:C[0,1]\to C[0,1]$ and $g:[0,1]\times C[0,1]\to \mathbb{R}$ are both continuous. The Green function $G(b,a)$ corresponding to the equation $(\ref{eqn6})$ takes the value $(b(1-a))^{\beta-1}-(b-a)^{\beta-1}$ whenever $0\leq a\leq b\leq 1$, and it attains the value $\frac{(b(1-a))^{\beta-1}}{\Gamma(\beta)}$ whenever $0\leq b\leq a\leq1$.

	The collection of every continuous real-valued mappings on $[0,1]$ is represented by $C([0,1])$, which incorporates with $\|.\|=\max_{b\in [0,1]}|w(b)|$. By applying our main Theorem \ref{thm3.2}, in the succeeding result we establish the sufficient criteria for ensuring a solution of the above mentioned fractional differential equation $(\ref{eqn6})$.

	\begin{theorem}\label{thm4.1}
		Take the fractional differential equation as stated in $(\ref{eqn6})$. Assume that the succeeding criteria happen:
		\begin{enumerate}
			\item[(i)] for every $v,w\in C[0,1]$ and an element $b$ lies in $[0,1]$,
			\begin{equation*}
				\left|g(b,f(v(b)))-g(b,f(w(b)))\right|\leq k\left(\left\|fv-fw\right\|\right)\left|fv(b)-fw(b)\right|
			\end{equation*}
		    where $k\in U$,
			\item[(ii)] for each $w\in C[0,1]$, $Fw\subseteq f\left(C[0,1]\right)$ where $F$ is a single valued mapping on $C[0,1]$, which is stated as
			\begin{equation*}
				F(w(b))=\int_0^1G(b,a)g(a,(fw)(a)) da.
			\end{equation*}
		\end{enumerate}
		Then equation $(\ref{eqn6})$ admits a solution if  $ f(C[0,1])$ is closed set in $C[0,1]$.
	\end{theorem}

	\begin{proof}
		We assume that $W=C[0,1]$ incorporate with a graph $G$ having $V(G)=W$ and $E(G)=\{(v,w):(v,w)\in W\times W\}$. Now we visualize that, a function $w$ lies in $W$ meets the equation $(\ref{eqn6})$, if and only if it fulfils the succeeding equation,
		\begin{equation}\label{eqn7}
			w(b)=\int_0^1G(b,a)g(a,(fw)(a))da.
		\end{equation}
		Let $v\in W$ and $fw\in Fv=\int_0^1G(b,a)g(a,(fv)(a))da$. Now
		\begin{eqnarray*}
			|fw(b)-Fw(b)|&= & \left|\int_0^1G(b,a)g(a,(fv)(a))da-\int_0^1G(b,a)g(a,(fw)(a))da\right|\\
			&\leq &\int_0^1\left|G(b,a)\right|\left|g(a,(fv)(a))-g(a,(fw)(a))\right|da\\
			&\leq &k\left(\left\|fv-fw\right\|\right)\int_0^1\left|G(b,a)\right|\left|fv(a)-fw(a)\right|da~\text{[from $(i)$]}\\
			&\leq &k\left(\left\|fv-fw\right\|\right)\|fv-fw\|\int_0^1\left|G(b,a)\right|da\\
			&\leq & k\left(\left\|fv-fw\right\|\right)\|fv-fw\|.
		\end{eqnarray*}
		Therefore $\|fw-Fw\|\leq k(\|fv-fw\|)\|fv-fw\|$. Assume $\{fu_n\}_n \in W$ in order that $fu_n\to fu$. Consequently for each $n$,
		\begin{eqnarray*}
			|fu-Fu|&\leq &|fu-fu_n|+|fu_n-Fu_n|+|Fu_n-Fu|\\
			&\leq & \left|\int_0^1G(b,a)g(a,(fu_n)(a))da-\int_0^1G(b,a)g(a,(fu)(a))da\right|\\  & &~ + |fu-fu_n|+|fu_n-Fu_n|\\
			&\leq & |fu-fu_n|+|fu_n-Fu_n|+k(\|fu_n-fu\|)\|fu_n-fu\|\textrm{[by~ $(i)$]}.
		\end{eqnarray*}
		Therefore for each $n$, $\|fu-Fu\|\leq \|fu-fu_n\|+\|fu_n-Fu_n\|+\|fu_n-fu\|$. Hence $\|fu-Fu\|\leq\liminf_{n\to \infty}\|fu_n-Fu_n\|$. Then using Theorem \ref{thm3.2} we can conclude that $F$ and $f$ has a function $w^*(b)\in W$ in order that $Fw^*(b)=f(w^*(b))$. Hence $f(w^*(b))$ fulfils the equation $(\ref{eqn7})$, therefore the given equation $(\ref{eqn6})$ has a solution.
	\end{proof}

	For $n\in\mathbb{N}$ and $q>0$, Lupa\c{s} \cite{lupas} first presented the $q$-analogue Bernstein operator $L_{n,q}$ on the space $C[0,1]$, in the year $1987$. In fact if $q=1$, then this operator reduces to the classical Bernstein operator \cite{KR8}. For $q>0$ and any $\phi$ lies in $C[0,1]$, we consider a nonlinear $q$-analogue Bernstein operator,
	\begin{equation}\label{eq1_Bern}
		\tilde {L}_{n,q}(\phi(a))=\sum_{i=0}^n \left|\phi\bigg(\frac{[i]_q}{[n]_q}\bigg)\right|b_{n,i}(q,a)\quad\textrm{for}~a\in[0,1],
	\end{equation}
	where for each $i\in \mathbb{Z}^{+}$,
	\begin{equation*}
		b_{n,i}(q,a)=\frac{\left[\begin{array}{l}{n}\\{i}\end{array}\right]_q q^{\frac{i(i-1)}{2}}a^i(1-a)^{n-i}}{\prod_{j=0}^{n-1}\{1-a+q^ja\}},
		\hspace*{0.15cm}
		\left[\begin{array}{l}{n}\\{i}\end{array}\right]_q=\frac{[n]_q!}{[n-i]_q![i]_q!},
		\hspace*{0.15cm}
		[n]_q!=\prod_{i=1}^{n}[i]_q,
	\end{equation*}
	and $[0]_q!=1$, where $[i]_q=\sum_{k=0}^{i-1}q^{k}$ for $i\geq 1$ and $[0]_q=0$.
	Agratini \cite{Agra} first studied the convergence of iterates of $q$-analogue Bernstein operator $L_{n,q}$ on the space $C[0,1]$ in the year 2008. In the upcoming corollary we derive the the convergence of iterates for nonlinear  $q$-analogue Bernstein operator $\tilde{L}_{n,q}$ through Theorem \ref{iterates}.

	\begin{corollary}
		Suppose that $\tilde{L}_{n,q}~(q>0)$ is an operator as mentioned in $(\ref{eq1_Bern})$. Then for each $\varphi\in C[0,1]$ and fixed $n\in \mathbb{N}$, $\lim_{j\to \infty}\tilde{L}_{n,q}^j\varphi \in \{\varphi^{*}\in C[0,1]:\tilde{L}_{n,q}\varphi^{*}=\varphi^{*}\}$. Furthermore,
		for each $\varphi\in W^{'}=\{\varsigma\in C[0,1]$: $\varsigma(1)\geq 0$ and $\varsigma(0)\geq 0$\},
		\begin{equation*}
			\lim_{j\to \infty}(\tilde{L}_{n,q}^j\varphi)(a)=\varphi(0)(1-a)+\varphi(1)a \quad \textrm{whereas}~0\leq a\leq 1.
		\end{equation*}
	\end{corollary}

	\begin{proof}
		Consider $W=C[0,1]$ with supremum norm. Let $W_0\subseteq W$ consist of all $\varphi\in W$ whose zeros are $0$ and $1$. Take a function $\phi$ in $W$ in order that $\phi-\tilde{L}_{n,q}\phi\in W_0$. Then for each $a$ lies in $[0,1]$ we have,
		\begin{eqnarray}

			\big|\tilde{L}_{n,q}\phi(a)-\tilde{L}_{n,q}^2\phi(a)\big|&\leq&\sum_{i=1}^{n-1} b_{n,i}(q,a) \left|(\phi-\tilde{L}_{n,q}\phi)\left(\frac{[i]_q}{[n]_q}\right)\right|\nonumber\\
			&\leq&\left[1-b_{n,n}(q,a)-b_{n,0}(q,a)\right]\|\phi-\tilde{L}_{n,q}\phi\|.\label{eq10_application2}

		\end{eqnarray}
		Consider $h_{n}(a)=q^{n(n-1)/2}a^n+(1-a)^n$, then for $n=1$, the minimum value of $h_{n}(a)$ is $1$. For $n \geq 2$, the minimum value of $h_{n}(a)$  is $\left(q^{n/2}/(1+q^{n/2})\right)^{n-1}$, which occurs at $a=(1+q^{n/2})^{-1}$.
		On the other hand consider $g_{n}(a)=\prod_{s=0}^{n-1}(1-a+q^s a)$, then one can easily verify that maximum of $g_{n}(a)$ over $[0,1]$ is less than or  equal to $\max\{1,q^{(n-1)^2}\}$. Consequently,
		\begin{equation}
			\frac{1}{\max\{q^{(n-1)^2},1\}}\bigg(\frac{q^{n/2}}{1+q^{n/2}}\bigg)^{n-1} \leq \displaystyle{\frac{h_{n}(a)}{g_{n}(a)}}= b_{n,n}(q,a)+b_{n,0}(q,a).

			\label{eq9_application2}
		\end{equation}
		Again we observe that $0<b_{n,q}=\frac{1}{\max\{1,q^{(n-1)^2}\}}\bigg(\frac{q^{n/2}}{1+q^{n/2}}\bigg)^{n-1} \leq 1$. Subsequently from (\ref{eq10_application2}) and (\ref{eq9_application2}), it follows that
		$||\tilde{L}_{n,q} \phi-\tilde{L}_{n,q}^2\phi||\leq k(||\phi-\tilde{L}_{n,q} \phi||)||\phi-\tilde{L}_{n,q} \phi||$, where $\forall\,t \geq 0$, $k(t)=(1-b_{n,q})\in U$.
		Consider a sequence $\{\eta_j\}_j\in W$ in order that $\eta_j\to \eta$. Then for every $a\in[0,1] $ and $j\in \mathbb{Z}^{+}$,
		\begin{equation*}
			\begin{aligned}
			|\eta(a) -\tilde{L}_{n,q}\eta(a)|&\leq|\eta(a)-\eta_j(a)|+|\eta_j(a) -\tilde{L}_{n,q}\eta_j(a)|+|\tilde{L}_{n,q}\eta_j(a)-\tilde{L}_{n,q}\eta(a)|\\
			&\leq \|\eta-\eta_j\|+\|\eta_j -\tilde{L}_{n,q}\eta_j\|+\|\eta_j-\eta\|.
			\end{aligned}
		\end{equation*}
		Therefore $\|\eta -\tilde{L}_{n,q}\eta\|\leq\|\eta-\eta_j\|+\|\eta_j -\tilde{L}_{n,q}\eta_j\|+\|\eta_j-\eta\|$. Now taking $j\to \infty$, we have $\|\eta -\tilde{L}_{n,q}\eta\|\leq\liminf_{j\to \infty}\|\eta_j -\tilde{L}_{n,q}\eta_j\|$. Again it is simple to observe that $(I-\tilde{L}_{n,q})(W)\subseteq W_0$. Then using the Theorem $\ref{iterates}$, for any $\varphi\in W$,
		$\displaystyle{\left\{\lim_{j\to \infty}\tilde{L}_{n,q}^j\varphi\right\}}\subseteq(\varphi+W_0)\cap Fix\,\tilde{L}_{n,q}$. It is simple to verify that the set $(\varphi+W_0)\cap Fix\,\tilde{L}_{n,q}$ is singleton, also $a$ and $1-a$ are fixed points of $\tilde{L}_{n,q}$. Again for each $\varphi\in W^{'}$, $\varphi(0)(1-a)+\varphi(1)a \in Fix\,\tilde{L}_{n,q} $. Furthermore, $\varphi(0)(1-a)+\varphi(1)a\in \varphi+W_0$. Consequently, we conclude that 	$\lim_{j\to \infty}(\tilde{L}_{n,q}^j\varphi)(a)=\varphi(0)(1-a)+\varphi(1)a$ for each $\varphi\in W^{'}$.
	\end{proof}

	\section*{Acknowledgment}
	The first author is thankful to the Ministry of Human Resource Development (MHRD), India for the financial support.