New classes of reverse super edge magic graphs

Kotte Amaranadha Reddy; S Sharief Basha; Kotte Amaranadha Reddy; S Sharief Basha

doi:10.3934/math.2022198

AIMS Mathematics

2022, Volume 7, Issue 3: 3590-3602. doi: 10.3934/math.2022198

Previous Article Next Article

Research article Special Issues

New classes of reverse super edge magic graphs

Kotte Amaranadha Reddy ,
S Sharief Basha ^,

Department of Mathematics, Vellore Institute of Technology, Vellore, India

Received: 23 June 2021 Accepted: 30 November 2021 Published: 03 December 2021
MSC : 37E25, 05C38

A reverse edge magic (REM) labeling of a graph $G(V, E)$ with $p$ vertices and $q$ edges is a bijection $f:V\left(G \right)\cup E\left(G \right)\to \{1, 2, \cdot \cdot \cdot, p+q\}$ such that $k = f\left(uv \right)-\{f\left(u \right)+f\left(v \right)\}$ is a constant $k$ for any edge $uv\in E\left(G \right).$ A REM labeling $f$ is called reverse super edge magic (RSEM) labeling if $f(V(G)) = \; \{1, 2, 3, 4, 5, \ldots, v\}$ and $f(E(G)) = \{v+1, v+2, v+3, v+4, v+5, \ldots, v+e\}.$ In this paper, we find some new classes of RSEM labeling and the investigation of the connection between the RSEM labeling and different classes of labeling.

Keywords:

Citation: Kotte Amaranadha Reddy, S Sharief Basha. New classes of reverse super edge magic graphs[J]. AIMS Mathematics, 2022, 7(3): 3590-3602. doi: 10.3934/math.2022198

Related Papers:

[1]	Yayun Fu, Mengyue Shi . A conservative exponential integrators method for fractional conservative differential equations. AIMS Mathematics, 2023, 8(8): 19067-19082. doi: 10.3934/math.2023973
[2]	Yong-Chao Zhang . Least energy solutions to a class of nonlocal Schrödinger equations. AIMS Mathematics, 2024, 9(8): 20763-20772. doi: 10.3934/math.20241009
[3]	Tingting Ma, Yuehua He . An efficient linearly-implicit energy-preserving scheme with fast solver for the fractional nonlinear wave equation. AIMS Mathematics, 2023, 8(11): 26574-26589. doi: 10.3934/math.20231358
[4]	Karmina K. Ali, Resat Yilmazer . Discrete fractional solutions to the effective mass Schrödinger equation by mean of nabla operator. AIMS Mathematics, 2020, 5(2): 894-903. doi: 10.3934/math.2020061
[5]	Erdal Bas, Ramazan Ozarslan . Theory of discrete fractional Sturm–Liouville equations and visual results. AIMS Mathematics, 2019, 4(3): 593-612. doi: 10.3934/math.2019.3.593
[6]	Dengfeng Lu, Shuwei Dai . On a class of three coupled fractional Schrödinger systems with general nonlinearities. AIMS Mathematics, 2023, 8(7): 17142-17153. doi: 10.3934/math.2023875
[7]	Mubashir Qayyum, Efaza Ahmad, Hijaz Ahmad, Bandar Almohsen . New solutions of time-space fractional coupled Schrödinger systems. AIMS Mathematics, 2023, 8(11): 27033-27051. doi: 10.3934/math.20231383
[8]	Xiaojun Zhou, Yue Dai . A spectral collocation method for the coupled system of nonlinear fractional differential equations. AIMS Mathematics, 2022, 7(4): 5670-5689. doi: 10.3934/math.2022314
[9]	Zunyuan Hu, Can Li, Shimin Guo . Fast finite difference/Legendre spectral collocation approximations for a tempered time-fractional diffusion equation. AIMS Mathematics, 2024, 9(12): 34647-34673. doi: 10.3934/math.20241650
[10]	Xiao-Yu Li, Yu-Lan Wang, Zhi-Yuan Li . Numerical simulation for the fractional-in-space Ginzburg-Landau equation using Fourier spectral method. AIMS Mathematics, 2023, 8(1): 2407-2418. doi: 10.3934/math.2023124

Abstract

1. Introduction

Fractional calculus is a popular subject because of having a lot of application areas of theoretical and applied sciences, like engineering, physics, biology, etc. Discrete fractional calculus is more recent area than fractional calculus and it was first defined by Diaz–Osler ^[1], Miller–Ross ^[2] and Gray–Zhang ^[3]. More recently, the theory of discrete fractional calculus have begun to develop rapidly with Goodrich–Peterson ^[4], Baleanu et al. ^[5,6], Ahrendt et al. ^[7], Atici–Eloe ^[8,9], Anastassiou ^[10], Abdeljawad et al. ^{[11,12,13,14,15,16]}, Hein et al. ^[17] and Cheng et al. ^[18], Mozyrska ^[19] and so forth ^{[20,21,22,23,24,25]}.

Fractional Sturm–Liouville differential operators have been studied by Bas et al. ^[26,27], Klimek et al.^[28], Dehghan et al. ^[29]. Besides that, Sturm–Liouville differential and difference operators were studied by ^{[30,31,32,33]}. In this study, we define DFHA operators and prove the self–adjointness of DFHA operator, some spectral properties of the operator.

More recently, Almeida et al. ^[34] have studied discrete and continuous fractional Sturm–Liouville operators, Bas–Ozarslan ^[35] have shown the self–adjointness of discrete fractional Sturm–Liouville operators and proved some spectral properties of the problem.

Sturm–Liouville equation having hydrogen atom potential is defined as follows

$\frac{d^{2}R}{dr^{2}}+\frac{a}{r}\frac{dR}{dr}-\frac{\ell \left( \ell +1\right) }{r^{2}}R+\left( E+\frac{a}{r}\right) R = 0\qquad \left( 0 \lt r \lt \infty \right) .$

In quantum mechanics, the study of the energy levels of the hydrogen atom leads to this equation. Where $R$ is the distance from the mass center to the origin, $\ell$ is a positive integer, $a$ is real number $E$ is energy constant and $r$ is the distance between the nucleus and the electron.

The hydrogen atom is a two–particle system and it composes of an electron and a proton. Interior motion of two particles around the center of mass corresponds to the movement of a single particle by a reduced mass. The distance between the proton and the electron is identified $r$ and $r$ is given by the orientation of the vector pointing from the proton to the electron. Hydrogen atom equation is defined as Schrödinger equation in spherical coordinates and in consequence of some transformations, this equation is defined as

$y^{\prime \prime }+\left( \lambda -\dfrac{l\left( l+1\right) }{x^{2}}+\dfrac{ 2}{x}-q\left( x\right) \right) y = 0.$

Spectral theory of hydrogen atom equation is studied by ^[39,40,41]. Besides that, we can observe that hydrogen atom differential equation has series solution as follows (^[39], p.268)

$\begin{eqnarray} y(x) = &a_0x^{l+1}\left\{1-\dfrac{k-l-1}{1!(2l+2)}.\dfrac{2x}{k}+\dfrac{(k-l-1)(k-l-2)}{2!(2l+2)(2l+3)}\left(\dfrac{2x}{k} \right)^2+\ldots \right.\\ &\left.+(-1)^n\dfrac{(k-l-1)(k-l-2)\ldots3.2.1}{(k-1)!(2l+2)(2l+3)\ldots(2l+n)}\left(\dfrac{2x}{k} \right)^n \right\}, k = 1, 2, \ldots \end{eqnarray}$

(1.1)

Recently, Bohner and Cuchta ^[36,37] studied some special integer order discrete functions, like Laguerre, Hermite, Bessel and especially Cuchta mentioned the difficulty in obtaining series solution of discrete special functions in his dissertation (^[38], p.100). In this regard, finding series solution of DFHA equations is an open problem and has some difficulties in the current situation. For this reason, we study to obtain solutions of DFHA eq.s in a different way with representation of solutions.

In this study, we investigate DFHA equation in Riemann–Liouville and Grü nwald–Letnikov sense. The aim of this study is to contribute to the spectral theory of DFHA operator and behaviors of eigenfunctions and also to obtain the solution of DFHA equation.

We investigate DFHA equation in three different ways;

$i)$ (nabla left and right) Riemann–Liouville (R–L)sense,

$L_{1}x\left( t\right) = \nabla _{a}^{\mu }\left( _{b}\nabla ^{\mu }x\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2}{t} +q\left( t\right) \right) x\left( t\right) = \lambda x\left( t\right) , \text{ }0 \lt \mu \lt 1,$

$ii)$ (delta left and right) Grünwald–Letnikov (G–L) sense,

$L_{2}x\left( t\right) = \Delta _{-}^{\mu }\left( \Delta _{+}^{\mu }x\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2}{t} +q\left( t\right) \right) x\left( t\right) = \lambda x\left( t\right) , \text{ }0 \lt \mu \lt 1,$

$iii)$ (nabla left) Riemann–Liouville (R–L)sense,

$L_{3}x\left( t\right) = \nabla _{a}^{\mu }\left( \nabla _{a}^{\mu }x\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2}{t} +q\left( t\right) \right) x\left( t\right) = \lambda x\left( t\right) , \text{ }0 \lt \mu \lt 1.$

2. Preliminaries

Definition 2.1. ^[42] Falling and rising factorial functions are defined as follows respectively

$\begin{equation} t^{\underline{\alpha }} = \frac{\Gamma \left( t +1\right) }{\Gamma \left( t-\alpha +1\right) }, \end{equation}$

(2.1)

$\begin{equation} t^{\overline{\alpha}} = \frac{\Gamma \left( t+\alpha \right) }{\Gamma \left( t\right) }, \end{equation}$

(2.2)

where $\Gamma$ is the gamma function, $\alpha \in \mathbb{R}$ .

Remark 2.1. Delta and nabla operators hold the following properties

$\begin{eqnarray} \Delta t^{\underline{\alpha }} & = &\alpha t^{\underline{\alpha -1}}, \\ \nabla t^{\overline{\alpha}} & = &\alpha t^{\overline{\alpha -1}}. \end{eqnarray}$

(2.3)

Definition 2.2. ^[2,8,11] Nabla fractional sum operators are given as below,

$\left(i\right)$ The left fractional sum of order $\mu > 0$ is defined by

$\begin{equation} \nabla _{a}^{-\mu }x\left( t\right) = \frac{1}{\Gamma \left( \mu \right) } \sum\limits_{s = a+1}^{t}\left( t-\rho \left( s\right) \right) ^{\overline{\mu -1} }x\left( s\right) , \text{ }t\in \mathbb{N} _{a+1}, \end{equation}$

(2.4)

$\left(ii\right)$ The right fractional sum of order $\mu > 0$ is defined by

$\begin{equation} _{b}\nabla ^{-\mu }x\left( t\right) = \frac{1}{\Gamma \left( \mu \right) } \sum\limits_{s = t}^{b-1}\left( s-\rho \left( t\right) \right) ^{\overline{\mu -1} }x\left( s\right) , \text{ }t\in \text{ }_{b-1} \mathbb{N} , \end{equation}$

(2.5)

where $\rho \left(t\right) = t-1$ is called backward jump operators, $\mathbb{N} _{a} = \left \{ a, a+1, ...\right \},$ $_{b} \mathbb{N} = \left \{ b, b-1, ...\right \}$ .

Definition 2.3. ^[12,14] Nabla fractional difference operators are as follows,

$\left(i\right)$ The left fractional difference of order $\mu > 0$ is defined by

$\begin{equation} \nabla _{a}^{\mu }x\left( t\right) = \nabla ^{n}\nabla _{a}^{-\left( n-\mu \right) }x\left( t\right) = \frac{\nabla ^{n}}{\Gamma \left( n-\mu \right) } \sum\limits_{s = a+1}^{t}\left( t-\rho \left( s\right) \right) ^{\overline{n-\mu -1} }x\left( s\right) , \text{ }t\in \mathbb{N} _{a+1}, \end{equation}$

(2.6)

$\left(ii\right)$ The right fractional difference of order $\mu > 0$ is defined by

$\begin{equation} _{b}\nabla ^{\mu }x\left( t\right) = \left( -1\right) ^{n}\nabla ^{n}\nabla _{a}^{-\left( n-\mu \right) }x\left( t\right) = \frac{\left( -1\right) ^{n}\Delta ^{n}}{\Gamma \left( n-\mu \right) }\sum\limits_{s = t}^{b-1}\left( s-\rho \left( t\right) \right) ^{\overline{n-\mu -1}}x\left( s\right) , \text{ }t\in \text{ }_{b-1} \mathbb{N} . \end{equation}$

(2.7)

Fractional differences in $\left(2.6-2.7\right)$ are called the Riemann–Liouville (R–L) definition of the $\mu$ -th order nabla fractional difference.

Definition 2.4. ^[1,18] Fractional difference operators are given as follows

$\left(i\right)$ The delta left fractional difference of order $\mu,$ $0 < \mu \leq 1,$ is defined by

$\begin{equation} \Delta _{-}^{\mu }x\left( t\right) = \frac{1}{h^{\mu }}\sum\limits_{s = 0}^{t}\left( -1\right) ^{s}\frac{\mu \left( \mu -1\right) ...\left( \mu -s+1\right) }{s!} x\left( t-s\right) , \text{ }t = 1, ..., N. \end{equation}$

(2.8)

$\left(ii\right)$ The delta right fractional difference of order $\mu,$ $0 < \mu \leq 1,$ is defined by

$\begin{equation} \Delta _{+}^{\mu }x\left( t\right) = \frac{1}{h^{\mu }}\sum\limits_{s = 0}^{N-t}\left( -1\right) ^{s}\frac{\mu \left( \mu -1\right) ...\left( \mu -s+1\right) }{s!} x\left( t+s\right) , \text{ }t = 0, .., N-1, \end{equation}$

(2.9)

fractional differences in $\left(2.8-2.9\right)$ are called the Grünwald–Letnikov (G–L) definition of the $\mu$ -th order delta fractional difference.

Definition 2.5 ^[14] Integration by parts formula for R–L nabla fractional difference operator is defined by, $u$ is defined on $_{b} \mathbb{N}$ and $v$ is defined on $\mathbb{N} _{a}$ ,

$\begin{equation} \sum\limits_{s = a+1}^{b-1}u\left( s\right) \nabla _{a}^{\mu }v\left( s\right) = \sum\limits_{s = a+1}^{b-1}v\left( s\right) _{b}\nabla ^{\mu }u\left( s\right) . \end{equation}$

(2.10)

Definition 2.6. ^[34] Integration by parts formula for G–L delta fractional difference operator is defined by, $u$ , $v$ is defined on $\left \{ 0, 1, ..., n\right \}$ , then

$\begin{equation} \sum\limits_{s = 0}^{n}u\left( s\right) \Delta _{-}^{\mu }v\left( s\right) = \sum\limits_{s = 0}^{n}v\left( s\right) \Delta _{+}^{\mu }u\left( s\right) . \end{equation}$

(2.11)

Definition 2.7. ^[17] $f: \mathbb{N} _{a}\rightarrow \mathbb{R}$ , $s\in \Re,$ Laplace transform is defined as follows,

$\mathcal{L} _{a}\left \{ f\right \} \left( s\right) = \sum\limits_{k = 1}^{\infty }\left( 1-s\right) ^{k-1}f\left( a+k\right) ,$

where $\Re = \mathbb{C} \backslash \left \{ 1\right \}$ and $\Re$ is called the set of regressive (complex) functions.

Definition 2.8. ^[17] Let $f, g: \mathbb{N} _{a}\rightarrow \mathbb{R}$ , all $t\in \mathbb{N} _{a+1},$ convolution of $f$ and $g$ is defined as follows

$\left( f\ast g\right) \left( t\right) = \sum\limits_{s = a+1}^{t}f\left( t-\rho \left( s\right) +a\right) g\left( s\right) ,$

where $\rho \left(s\right)$ is the backward jump function defined in ^[42] as

$\rho \left( s\right) = s-1.$

Theorem 2.1. ^[17] $f, g: \mathbb{N} _{a}\rightarrow \mathbb{R},$ convolution theorem is expressed as follows,

$\mathcal{L} _{a}\left \{ f\ast g\right \} \left( s\right) = \mathcal{L} _{a}\left \{ f\right \} \mathcal{L} _{a}\left \{ g\right \} \left( s\right) .$

Lemma 2.1. ^[17] $f: \mathbb{N} _{a}\rightarrow \mathbb{R},$ the following property is valid,

$\mathcal{L} _{a+1}\left \{ f\right \} \left( s\right) = \frac{1}{1-s} \mathcal{L} _{a}\left \{ f\right \} \left( s\right) -\frac{1}{1-s}f\left( a+1\right) .$

Theorem 2.2. ^[17] $f: \mathbb{N} _{a}\rightarrow \mathbb{R},$ $0 < \mu < 1,$ Laplace transform of nabla fractional difference

$\mathcal{L} _{a+1}\left \{ \nabla _{a}^{\mu }f\right \} \left( s\right) = s^{\mu }\mathcal{L} _{a+1}\left \{ f\right \} \left( s\right) -\frac{ 1-s^{\mu }}{1-s}f\left( a+1\right) , \mathit{\text{}}t\in \mathbb{N} _{a+1}.$

Definition 2.9. ^[17] For $\left \vert p\right \vert < 1$ , $\alpha > 0,$ $\beta \in \mathbb{R}$ and $t\in \mathbb{N} _{a}$ , Mittag–Leffler function is defined by

$E_{p, \alpha , \beta }\left( t, a\right) = \sum\limits_{k = 0}^{\infty }p^{k}\frac{\left( t-a\right) ^{\overline{\alpha k+\beta }}}{\Gamma \left( \alpha k+\beta +1\right) }.$

Theorem 2.3. ^[17] For $\left \vert p\right \vert < 1$ , $\alpha > 0,$ $\beta \in \mathbb{R}$ , $\left \vert 1-s\right \vert < 1$ and $\left \vert s\right \vert ^{\alpha } > p$ , Laplace transform of Mittag–Leffler function is as follows,

$\mathcal{L} _{a+1}\left \{ E_{p, \alpha , \beta }\left( ., a\right) \right \} \left( s\right) = \frac{s^{\alpha -\beta -1}}{s^{\alpha }-p}.$

3. Main results

3.1. Discrete fractional hydrogen–atom equations

Let us consider equations in three different forms;

$i)$ $L_{1}$ DFHA operator $L_{1}$ is defined in (nabla left and right) R–L sense,

$\begin{equation} L_{1}x\left( t\right) = \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }x\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}- \dfrac{2}{t}+q\left( t\right) \right) x\left( t\right) = \lambda x\left( t\right) , \text{ }0 \lt \mu \lt 1, \end{equation}$

(3.1)

where $l$ is a positive integer or zero, $q\left(t\right) +\dfrac{2}{t}- \dfrac{l\left(l+1\right) }{t^{2}}$ are named potential function., $\lambda$ is the spectral parameter, $t\in \left[a+1, b-1\right],$ $x\left(t\right) \in l^{2}\left[a+1, b-1\right],$ $a > 0.$

$ii)$ $L_{2}$ DFHA operator $L_{2}$ is defined in (delta left and right) G–L sense,

$\begin{equation} L_{2}x\left( t\right) = \Delta _{-}^{\mu }\left( p\left( t\right) \Delta _{+}^{\mu }x\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2} }-\dfrac{2}{t}+q\left( t\right) \right) x\left( t\right) = \lambda x\left( t\right) , \text{ }0 \lt \mu \lt 1, \end{equation}$

(3.2)

where $p, q, l, \lambda$ is as defined above, $t\in \left[1, n\right],$ $x\left(t\right) \in l^{2}\left[0, n\right].$

$iii)$ $L_{3}$ DFHA operator $L_{3}$ is defined in (nabla left) R–L sense,

$\begin{equation} L_{3}x\left( t\right) = \nabla _{a}^{\mu }\left( \nabla _{a}^{\mu }x\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2}{t} +q\left( t\right) \right) x\left( t\right) = \lambda x\left( t\right) , \text{ }0 \lt \mu \lt 1, \end{equation}$

(3.3)

$p, q, l, \lambda$ is as defined above, $t\in \left[a+1, b-1\right],$ $a > 0.$

Theorem 3.1. DFHA operator $L_{1}$ is self–adjoint.

Proof.

$\begin{eqnarray} u\left( t\right) L_{1}v\left( t\right) &\bf{ = }&u\left( t\right) \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }v\left( t\right) \right) +u\left( t\right) \left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2 }{t}+q\left( t\right) \right) v\left( t\right) , \end{eqnarray}$

(3.4)

$\begin{eqnarray} v\left( t\right) L_{1}u\left( t\right) &\bf{ = }&v\left( t\right) \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }u\left( t\right) \right) +v\left( t\right) \left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2 }{t}+q\left( t\right) \right) u\left( t\right) . \end{eqnarray}$

(3.5)

Subtracting $\left(16-17\right)$ from each other

$u\left( t\right) L_{1}v\left( t\right) -v\left( t\right) L_{1}u\left( t\right) = u\left( t\right) \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }v\left( t\right) \right) -v\left( t\right) \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }u\left( t\right) \right)$

and applying definite sum operator to both side of the last equality, we have

$\begin{equation} \sum\limits_{s = a+1}^{b-1}\left( u\left( s\right) L_{1}v\left( s\right) -v\left( s\right) L_{1}u\left( s\right) \right) = \sum\limits_{s = a+1}^{b-1}u\left( s\right) \nabla _{a}^{\mu }\left( p\left( s\right) _{b}\nabla ^{\mu }v\left( s\right) \right) -\sum\limits_{s = a+1}^{b-1}v\left( s\right) \nabla _{a}^{\mu }\left( p\left( s\right) _{b}\nabla ^{\mu }u\left( s\right) \right) . \end{equation}$

(3.6)

Applying the integration by parts formula $\left(2.10\right)$ to right hand side of $\left(18\right),$ we have

$\begin{eqnarray*} \sum\limits_{s = a+1}^{b-1}\left( u\left( s\right) L_{1}v\left( s\right) -v\left( s\right) L_{1}u\left( s\right) \right) & = &\sum\limits_{s = a+1}^{b-1}p\left( s\right) _{b}\nabla ^{\mu }v\left( s\right) _{b}\nabla ^{\mu }u\left( s\right) \\ &&-\sum\limits_{s = a+1}^{b-1}p\left( s\right) _{b}\nabla ^{\mu }u\left( s\right) _{b}\nabla ^{\mu }v\left( s\right) \\ & = &0, \end{eqnarray*}$

$\left \langle L_{1}u, v\right \rangle = \left \langle u, L_{1}v\right \rangle .$

The proof completes.

Theorem 3.2. Eigenfunctions, corresponding to distinct eigenvalues, of the equation $\left(3.2\right)$ are orthogonal.

Proof. Assume that $\lambda _{\alpha }$ and $\lambda _{\beta }$ are two different eigenvalues corresponds to eigenfunctions $u\left(n\right)$ and $v\left(n\right)$ respectively for the equation $\left(3.1 \right)$ ,

$\begin{eqnarray*} \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }u\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2}{t}+q\left( t\right) \right) u\left( t\right) -\lambda _{\alpha }u\left( t\right) & = &0, \\ \nabla _{a}^{\mu }\left( p\left( t\right) _{b}\nabla ^{\mu }v\left( t\right) \right) +\left( \dfrac{l\left( l+1\right) }{t^{2}}-\dfrac{2}{t}+q\left( t\right) \right) v\left( t\right) -\lambda _{\beta }v\left( t\right) & = &0, \end{eqnarray*}$

Multiplying last two equations to $v\left(n\right)$ and $u\left(n\right)$ respectively, subtracting from each other and applying sum operator, since the self–adjointness of the operator $L_{1},$ we get

$\left( \lambda _{\alpha }-\lambda _{\beta }\right) \sum\limits_{s = a+1}^{b-1}r\left( s\right) u\left( s\right) v\left( s\right) = 0,$

since $\lambda _{\alpha }\neq \lambda _{\beta, }$

$\begin{eqnarray*} \sum\limits_{s = a+1}^{b-1}r\left( s\right) u\left( s\right) v\left( s\right) & = &0, \\ \left \langle u\left( t\right) , v\left( t\right) \right \rangle & = &0, \end{eqnarray*}$

and the proof completes.

Theorem 3.3. All eigenvalues of the equation $\left(3.1\right)$ are real.

Proof. Assume $\lambda = \alpha +i\beta,$ since the self–adjointness of the operator $L_{1},$ we have

$\begin{eqnarray*} \left \langle L_{1}u, u\right \rangle & = &\left \langle u, L_{1}u\right \rangle , \\ \left \langle \lambda u, u\right \rangle & = &\left \langle u, \lambda u\right \rangle , \end{eqnarray*}$

$\left( \lambda -\overline{\lambda }\right) \left \langle u, u\right \rangle = 0$

Since $\left \langle u, u\right \rangle _{r}\neq 0,$

$\lambda = \overline{\lambda }$

and hence $\beta = 0.$ So, the proof is completed.

Self–adjointness of $L_{2}$ DFHA operator G–L sense, reality of eigenvalues and orthogonality of eigenfunctions of the equation 3.2 can be proven in a similar way to the Theorem 3.1–3.2–3.3 by means of Definition 2.5.

3.2. Representation of solution for discrete fractional hydrogen atom equation

Theorem 3.4.

(3.7)

$\begin{equation} x\left( a+1\right) = c_{1}, \mathit{\text{}}\nabla _{a}^{\mu }x\left( a+1\right) = c_{2}, \end{equation}$

(3.8)

where $p\left(t\right) > 0,$ $r\left(t\right) > 0,$ $q\left(t\right)$ is defined and real valued, $\lambda$ is the spectral parameter $.$ The sum representation of solution of the problem $\left(3.7\right) -\left(3.8\right)$ is given as follows,

$\begin{eqnarray} x\left( t\right) & = &c_{1}\left( \left( 1+\dfrac{l\left( l+1\right) }{\left( a+1\right) ^{2}}-\dfrac{2}{a+1}+q\left( a+1\right) \right) E_{\lambda , 2\mu , \mu -1}\left( t, a\right) -\lambda E_{\lambda , 2\mu , 2\mu -1}\left( t, a\right) \right) \\ &&+c_{2}\left( E_{\lambda , 2\mu , 2\mu -1}\left( t, a\right) -E_{\lambda , 2\mu , \mu -1}\left( t, a\right) \right) \\ &&-\sum\limits_{s = a+1}^{t}E_{\lambda , 2\mu , 2\mu -1}\left( t-\rho \left( s\right) +a\right) \left( \dfrac{l\left( l+1\right) }{s^{2}}-\dfrac{2}{s}+q\left( s\right) \right) x\left( s\right) . \end{eqnarray}$

(3.9)

Proof. Taking Laplace transform of the equation $\left(3.7 \right)$ by Theorem 2.2 and take $\left(\dfrac{l\left(l+1\right) }{t^{2}}- \dfrac{2}{t}+q\left(t\right) \right) x\left(t\right) = g\left(t\right),$

$\begin{eqnarray*} \mathcal{L} _{a+1}\left \{ \nabla _{a}^{\mu }\left( \nabla _{a}^{\mu }x\right) \right \} \left( s\right) +\mathcal{L} _{a+1}\left \{ g\right \} \left( s\right) & = &\lambda \mathcal{L} _{a+1}\left \{ x\right \} \left( s\right) , \\ \begin{array}{c} = \end{array} s^{\mu }\mathcal{L} _{a+1}\left \{ \nabla _{a}^{\mu }x\right \} \left( s\right) -\frac{1-s^{\mu }}{1-s}c_{2} & = &\lambda \mathcal{L} _{a+1}\left \{ x\right \} \left( s\right) -\mathcal{L} _{a+1}\left \{ g\right \} \left( s\right) , \\ \begin{array}{c} = \end{array} s^{\mu }\left( s^{\mu }\mathcal{L} _{a+1}\left \{ x\right \} \left( s\right) -\frac{1-s^{\mu }}{1-s}c_{1}\right) -\frac{1-s^{\mu }}{1-s}c_{2} & = &\lambda \mathcal{L} _{a+1}\left \{ x\right \} \left( s\right) -\mathcal{L} _{a+1}\left \{ g\right \} \left( s\right) , \end{eqnarray*}$

$\begin{array}{c} = \end{array} \mathcal{L} _{a+1}\left \{ x\right \} \left( s\right) = \frac{1-s^{\mu }}{1-s} \frac{1}{s^{2\mu }-\lambda }\left( s^{\mu }c_{1}+c_{2}\right) -\frac{1}{ s^{2\mu }-\lambda }\mathcal{L} _{a+1}\left \{ g\right \} \left( s\right) .$

Using Lemma 2.1, we have

$\begin{equation} \mathcal{L} _{a}\left \{ x\right \} \left( s\right) = c_{1}\left( \frac{ s^{\mu }-\lambda }{s^{2\mu }-\lambda }\right) -\frac{1-s}{s^{2\mu }-\lambda } \left( \frac{1}{1-s}\mathcal{L} _{a}\left \{ g\right \} \left( s\right) - \frac{1}{1-s}g\left( a+1\right) \right) +c_{2}\left( \frac{1-s^{\mu }}{ s^{2\mu }-\lambda }\right) . \end{equation}$

(3.10)

Now, taking inverse Laplace transform of the equation $\left(3.10 \right)$ and applying convolution theorem, then we have the representation of solution of the problem $\left(3.7\right) -\left(3.8\right),$ $\left \vert \lambda \right \vert < 1$ , $\left \vert 1-s\right \vert < 1$ and $\left \vert s\right \vert ^{\alpha } > \lambda$ from Theorem 2.3., i.e.

$\begin{eqnarray*} \mathcal{L}_{a}^{-1}\left \{ \frac{s^{\mu }}{s^{2\mu }-\lambda }\right \} & = &E_{\lambda , 2\mu , \mu -1}\left( t, a\right), \\ \mathcal{L}_{a}^{-1}\left \{ \frac{1}{s^{2\mu }-\lambda }\right \} & = &E_{\lambda , 2\mu , 2\mu -1}\left( t, a\right), \end{eqnarray*}$

$\begin{equation*} \mathcal{L}_{a}^{-1}\left \{ \frac{1}{s^{2\mu }-\lambda }\mathcal{L} _{a}\left \{ q\left( s\right) x\left( s\right) \right \} \right \} = \sum\limits_{s = a+1}^{t}E_{\lambda , 2\mu , 2\mu -1}\left( t-\rho \left( s\right) +a\right) q\left( s\right) x\left( s\right) . \end{equation*}$

Consequently, we have sum representation of solution for DFHA problem 3.7–3.8

$\begin{eqnarray*} x\left( t\right) & = &c_{1}\left( \left( 1+\dfrac{l\left( l+1\right) }{\left( a+1\right) ^{2}}-\dfrac{2}{a+1}+q\left( a+1\right) \right) E_{\lambda , 2\mu , \mu -1}\left( t, a\right) -\lambda E_{\lambda , 2\mu , 2\mu -1}\left( t, a\right) \right) \\ &&+c_{2}\left( E_{\lambda , 2\mu , 2\mu -1}\left( t, a\right) -E_{\lambda , 2\mu , \mu -1}\left( t, a\right) \right) \\ &&-\sum\limits_{s = a+1}^{t}E_{\lambda , 2\mu , 2\mu -1}\left( t-\rho \left( s\right) +a\right) \left( \dfrac{l\left( l+1\right) }{s^{2}}-\dfrac{2}{s}+q\left( s\right) \right) x\left( s\right) . \end{eqnarray*}$

3.3. Applications and visual results

Presume that $c_{1} = 1, c_{2} = 0, a = 0$ in the representation of solution $\left(3.9\right)$ and hence we may observe the behaviors of solutions in following figures (Figures 1–7) and tables (Tables 1–3);

Figure 1.

$q\left(t\right) = 0, \lambda = 0.01, l = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 2.

$\lambda = 0.01, \mu = 0.45, l = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 3.

$q\left(t\right) = 0, \mu = 0.4, l = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 4.

$q\left(t\right) = 0, l = 1, \lambda = 0.01, t = 4$ .

DownLoad: Full-Size Img PowerPoint

Figure 5.

$q\left(t\right) = 0, \lambda = -0.01, l = 0.01$ .

DownLoad: Full-Size Img PowerPoint

Figure 6.

$q\left(t\right) = 0, \lambda = -0.01, \mu = 0.9$ .

DownLoad: Full-Size Img PowerPoint

Figure 7.

$q\left(t\right) = 0, \mu = 0.9, l = 0.01$ .

DownLoad: Full-Size Img PowerPoint

Table 1.

$q\left(t\right) = 0, \lambda = 0.01, l = 1$ .

$x\left(t\right)$	$\mu =0.3$	$\mu =0.35$	$\mu =0.4$	$\mu =0.45$	$\mu =0.5$
$x\left(1\right)$	$1$	$1$	$1$	$1$	$1$
$x\left(2\right)$	$0.612$	$0.714$	$1.123$	$0.918$	$1.020$
$x\left(3\right)$	$0.700$	$0.900$	$1.515$	$1.370$	$1.641$
$x\left(5\right)$	$0.881$	$1.336$	$2.402$	$2.747$	$3.773$
$x\left(7\right)$	$1.009$	$1.740$	$3.352$	$4.566$	$7.031$
$x\left(9\right)$	$1.099$	$2.100$	$4.332$	$6.749$	$11.461$
$x\left(12\right)$	$1.190$	$2.570$	$5.745$	$10.623$	$20.450$
$x\left(15\right)$	$1.249$	$2.975$	$6.739$	$15.149$	$32.472$
$x\left(16\right)$	$1.264$	$3.098$	$7.235$	$16.793$	$37.198$
$x\left(18\right)$	$1.289$	$3.330$	$8.233$	$20.279$	$47.789$
$x\left(20\right)$	$1.309$	$3.544$	$9.229$	$24.021$	$59.967$

| Show Table

DownLoad: CSV

Table 2.

$\lambda = 0.01, \mu = 0.45, l = 1$ .

$x\left(t\right)$	$q\left(t\right) =1$	$q\left(t\right) =t$	$q\left(t\right) =\sqrt{t}$
$x\left(1\right)$	$1$	$1$	$1$
$x\left(2\right)$	$7.37\ast 10^{-17}$	$4.41\ast 10^{-17}$	$5.77\ast 10^{-17}$
$x\left(3\right)$	$-0.131$	$-0.057$	$-0.088$
$x\left(5\right)$	$-0.123$	$-0.018$	$-0.049$
$x\left(7\right)$	$-0.080$	$-0.006$	$-0.021$
$x\left(9\right)$	$-0.050$	$-0.003$	$-0.011$
$x\left(12\right)$	$-0.028$	$-0.001$	$-0.005$
$x\left(15\right)$	$-0.017$	$-0.0008$	$-0.003$
$x\left(16\right)$	$-0.015$	$-0.0006$	$-0.0006$
$x\left(18\right)$	$-0.012$	$-0.0005$	$-0.002$
$x\left(20\right)$	$-0.010$	$-0.0003$	$-0.001$

| Show Table

DownLoad: CSV

Table 3.

$q\left(t\right) = 0, \mu = 0.4, l = 1$ .

$x\left(t\right)$	$\lambda =0.1$	$\lambda =0.11$	$\lambda =0.12$
$x\left(1\right)$	$1$	$1$	$1$
$x\left(2\right)$	$1$	$1.025$	$1.052$
$x\left(3\right)$	$1.668$	$1.751$	$1.841$
$x\left(5\right)$	$3.876$	$4.216$	$4.595$
$x\left(7\right)$	$7.243$	$8.107$	$9.095$
$x\left(9\right)$	$11.941$	$13.707$	$12.130$
$x\left(12\right)$	$22.045$	$26.197$	$25.237$
$x\left(15\right)$	$36.831$	$45.198$	$46.330$
$x\left(16\right)$	$43.042$	$53.369$	$55.687$
$x\left(18\right)$	$57.766$	$73.092$	$78.795$
$x\left(20\right)$	$76.055$	$98.154$	$127.306$

| Show Table

DownLoad: CSV

4. Conclusion

We have analyzed DFHA equation in Riemann–Liouville and Grü nwald–Letnikov sense. Self–adjointness of the DFHA operator is presented and also, we have proved some significant spectral properties for instance, orthogonality of distinct eigenfunctions, reality of eigenvalues. Moreover, we give sum representation of the solutions for DFHA problem and find the solutions of the problem. We have carried out simulation analysis with graphics and tables. The aim of this paper is to contribute to the theory of hydrogen atom fractional difference operator.

We observe the behaviors of solutions by changing the order of the derivative $\mu$ in and , by changing the potential function $q\left(t\right)$ in , we compare solutions under different $\lambda$ eigenvalues in , and , also we observe the solutions by changing $\mu$ with a specific eigenvalue in and by changing $l$ values in Figure 6.

We have shown the solutions by changing the order of the derivative $\mu$ in , by changing the potential function $q\left(t\right)$ and $\lambda$ eigenvalues in Table 2, Table 3.

Acknowledgment

The authors would like to thank the editor and anonymous reviewers for their careful reading of our manuscript and their many insightful comments and suggestions.

Conflict of interest

The authors declare no conflict of interest.

References

[1]	A. Kotzig, A. Rosa, Magic valuations of finite graphs, Canad. Math. Bull., 13 (1970), 451–461. doi: 10.4153/CMB-1970-084-1. doi: 10.4153/CMB-1970-084-1
[2]	H. Enomoto, A. S. Lladó, T. Nakamigawa, G. Ringel, Super edge-magic graphs, SUT J. Math., 34 (1998), 105–109.
[3]	R. M. Figueroa-Centenoa, R. Ichishimab, F. A. Muntaner-Batle, The place of super edge magic labelings among other classes of labelings, Discrete Math., 231 (2001), 153–168. doi: 10.1016/S0012-365X(00)00314-9. doi: 10.1016/S0012-365X(00)00314-9
[4]	R. M. Figueroa-Centenoa, R. Ichishimab, F. A. Muntaner-Batle, Mgical coronations of graphs, Australas. J. Comb., 26 (2002), 199–208.
[5]	V. Yegnanarayanan, On magic graphs, Utilitas Math., 59 (2001), 181–204.
[6]	R. M. Figueroa-Centenoa, R. Ichishimab, F. A. Muntaner-Batle, On super edge-magic graphs, Ars Comb., 64 (2002), 81–95.
[7]	J. C. Bermond, Graceful graphs, radio antennae and French windmills, In: Proceedings one day combinatorics conference, Research notes in mathematics, Pitman, 34 (1979), 18–37.
[8]	G. Ringel, A. Llado, Another tree conjecture, Bull. Inst. Combin. Appl., 188 (1996), 83–85.
[9]	K. Kathiresan, Subdivisions of ladders are graceful, Indian J. Pure Appl. Math., 23 (1992), 21–23.
[10]	R. Balakrishnan, R. Sampathkumar, Decompositions of regular graphs into $k_{n}^{c} \vee 2k_{2}$ , Discrete Math., 156 (1996), 19–28. doi: 10.1016/0012-365X(94)00026-F. doi: 10.1016/0012-365X(94)00026-F
[11]	R. Balakrishna, A. Selvam, V. Yegnanarayanan, Some results on elegantgraphs, Indian J. Pure Appl. Math., 28 (1997), 905–916.

This article has been cited by:

1.	Erdal Bas, Funda Metin Turk, Ramazan Ozarslan, Ahu Ercan, Spectral data of conformable Sturm–Liouville direct problems, 2021, 11, 1664-2368, 10.1007/s13324-020-00428-6
2.	Tom Cuchta, Dallas Freeman, Discrete Polylogarithm Functions, 2023, 84, 1338-9750, 19, 10.2478/tmmp-2023-0012
3.	B. Shiri, Y. Guang, D. Baleanu, Inverse problems for discrete Hermite nabla difference equation, 2025, 33, 2769-0911, 10.1080/27690911.2024.2431000
4.	Muhammad Sulthan Zacky, Heru Sukamto, Lila Yuwana, Agus Purwanto, Eny Latifah, The performance of space-fractional quantum carnot engine, 2025, 100, 0031-8949, 025306, 10.1088/1402-4896/ada9de

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)