Research article Special Issues

Numerical analysis of multi-dimensional time-fractional diffusion problems under the Atangana-Baleanu Caputo derivative


  • Received: 12 December 2022 Revised: 01 February 2023 Accepted: 07 February 2023 Published: 27 February 2023
  • This paper presents the Elzaki homotopy perturbation transform scheme (EHPTS) to analyze the approximate solution of the multi-dimensional fractional diffusion equation. The Atangana-Baleanu derivative is considered in the Caputo sense. First, we apply Elzaki transform (ET) to obtain a recurrence relation without any assumption or restrictive variable. Then, this relation becomes very easy to handle for the implementation of the homotopy perturbation scheme (HPS). We observe that HPS produces the iterations in the form of convergence series that approaches the precise solution. We provide the graphical representation in 2D plot distribution and 3D surface solution. The error analysis shows that the solution derived by EHPTS is very close to the exact solution. The obtained series shows that EHPTS is a very simple, straightforward, and efficient tool for other problems of fractional derivatives.

    Citation: Muhammad Nadeem, Ji-Huan He, Hamid. M. Sedighi. Numerical analysis of multi-dimensional time-fractional diffusion problems under the Atangana-Baleanu Caputo derivative[J]. Mathematical Biosciences and Engineering, 2023, 20(5): 8190-8207. doi: 10.3934/mbe.2023356

    Related Papers:

    [1] Mehmet F. Dicle . US Implied Volatility as A predictor of International Returns. Quantitative Finance and Economics, 2017, 1(4): 388-402. doi: 10.3934/QFE.2017.4.388
    [2] Jie Xiong, Geng Deng, Xindong Wang . Extension of SABR Libor Market Model to handle negative interest rates. Quantitative Finance and Economics, 2020, 4(1): 148-171. doi: 10.3934/QFE.2020007
    [3] Hammad Siddiqi . Financial market disruption and investor awareness: the case of implied volatility skew. Quantitative Finance and Economics, 2022, 6(3): 505-517. doi: 10.3934/QFE.2022021
    [4] Elyas Elyasiani, Luca Gambarelli, Silvia Muzzioli . The Information Content of Corridor Volatility Measures During Calm and Turmoil Periods. Quantitative Finance and Economics, 2017, 1(4): 454-473. doi: 10.3934/QFE.2017.4.454
    [5] Mara Madaleno, Elisabete Vieira . Volatility analysis of returns and risk: Family versus nonfamily firms. Quantitative Finance and Economics, 2018, 2(2): 348-372. doi: 10.3934/QFE.2018.2.348
    [6] Didier Sornette, Peter Cauwels, Georgi Smilyanov . Can we use volatility to diagnose financial bubbles? lessons from 40 historical bubbles. Quantitative Finance and Economics, 2018, 2(1): 486-590. doi: 10.3934/QFE.2018.1.1
    [7] Sebastian Ferrando, Andrew Fleck, Alfredo Gonzalez, Alexey Rubtsov . Trajectorial asset models with operational assumptions. Quantitative Finance and Economics, 2019, 3(4): 661-708. doi: 10.3934/QFE.2019.4.661
    [8] David Melkuev, Danqiao Guo, Tony S. Wirjanto . Applications of random-matrix theory and nonparametric change-point analysis to three notable systemic crises. Quantitative Finance and Economics, 2018, 2(2): 413-467. doi: 10.3934/QFE.2018.2.413
    [9] Elvira Caloiero, Massimo Guidolin . Volatility as an Alternative Asset Class: Does It Improve Portfolio Performance?. Quantitative Finance and Economics, 2017, 1(4): 334-362. doi: 10.3934/QFE.2017.4.334
    [10] Chikashi Tsuji . The historical transition of return transmission, volatility spillovers, and dynamic conditional correlations: A fresh perspective and new evidence from the US, UK, and Japanese stock markets. Quantitative Finance and Economics, 2024, 8(2): 410-436. doi: 10.3934/QFE.2024016
  • This paper presents the Elzaki homotopy perturbation transform scheme (EHPTS) to analyze the approximate solution of the multi-dimensional fractional diffusion equation. The Atangana-Baleanu derivative is considered in the Caputo sense. First, we apply Elzaki transform (ET) to obtain a recurrence relation without any assumption or restrictive variable. Then, this relation becomes very easy to handle for the implementation of the homotopy perturbation scheme (HPS). We observe that HPS produces the iterations in the form of convergence series that approaches the precise solution. We provide the graphical representation in 2D plot distribution and 3D surface solution. The error analysis shows that the solution derived by EHPTS is very close to the exact solution. The obtained series shows that EHPTS is a very simple, straightforward, and efficient tool for other problems of fractional derivatives.



    The study of Diophantine equations plays a very important role in number theory, and the integer solutions of Diophantine equations are widely used in cryptography and coding theory. Silverman [1] studied the parametric solution of equation

    X3+Y3=A,

    Li and Yuan [2] proved that the simultaneous Pell equations possess at most one positive integer solution under certain conditions.

    A Diophantine equation of the form

    a1x1+a2x2++atxt=s (1.1)

    is called a multivariate linear Diophantine equation, where s,a1,a2,,at are nonzero integers and t2. It is well known that for any given nonzero integers a and b, there are two integers u and v such that

    ua+vb=(a,b),

    where (a,b) represents the greatest common divisor of a and b. Now we introduce some symbols associated with Eq (1.1) as follows:

    (a1,a2)=d2,(d2,a3)=d3,,(dt1,at)=dt.

    That is, there exist integers u1,u2,,ut and v2,v3,,vt1 such that

    {a1u1+a2u2=d2,d2v2+a3u3=d3,dt1vt1+atut=dt. (1.2)

    Li [3] gave the structure of the general solution of the multivariate linear Diophantine equation.

    Theorem 1.1. ([3]) The multivariate linear Diophantine Eq (1.1) has solutions if and only if dt|s. Furthermore, if dt|s and t4, then the general solutions of Diophantine Eq (1.1) are

    {x1=u1(δ2it1vi+4jtˉaj2ij2visj1+ˉa3s2)+ˉa2s1,x2=u2(δ2it1vi+4jtˉaj2ij2visj1+ˉa3s2)a1d12s1,x3=u3(δ3it1vi+5jtˉaj3ij2visj1+ˉa4s3)d2d13s2,xt1=ut1(δvt1+ˉatst1)dt2d1t1st2,xt=utδdt1d1tst1, (1.3)

    where si(1it1) are arbitrary integers and δ=sd1t,ˉaj=ajd1j for 2jt.

    Let Fq be a finite field of q elements with characteristic p, and Fqn be its extension of degree n, where p is a prime number and n 2 is an integer. Zhu et al. [4] obtained an explicit formula for the number of solutions to the equation

    f(x1)++f(xn)=a

    over Fq. Zhao et al. [5] found an explicit formula for the number of solutions of the two-variable diagonal quartic equation

    x41+x42=c

    over Fq.

    A basis of Fqn over Fq of the form {α,αq,,αqn1} is called a normal basis of Fqn over Fq, and α is called a normal element of Fqn over Fq. An irreducible polynomial f(x) Fq[x] is called a normal polynomial if all the roots of f(x) are normal elements of Fqn over Fq. The trace of α is defined as

    Tr(α)=α+αq++αqn1

    and the trace of f(x) is defined to be the coefficient of xn1. Theorems 1.1 and 1.2 below give a simple criterion to check when an irreducible polynomial is a normal polynomial.

    Theorem 1.2. ([6]) Let n=pe with e1. Then an irreducible polynomial

    f(x)=xn+a1xn1++anFq[x]

    is a normal polynomial if and only if a10.

    Theorem 1.3. ([7]) Let n be a prime different from p, and let q be a primitive root modulo n. Then an irreducible polynomial

    f(x)=xn+a1xn1++anFq[x]

    is a normal polynomial if and only if a10.

    In 2001, Chang([8]) et al. furthermore proved that the conditions in Theorems 1.1 and 1.2 are also necessary.

    Theorem 1.4. ([8]) If every irreducible polynomial

    f(x)=xn+a1xn1++anFq[x]

    with a10 is a normal polynomial, then n is either a power of p or a prime different from p, and q is a primitive root modulo n.

    In 2018, Huang et al. [9] presented a unified proof of Theorems 1.1–1.3 by comparing the number of normal polynomials and that of irreducible polynomials over Fq.

    The factorization of xn1 and its irreducible factors are closely related to the normal elements in Fqn over Fq (see [10, Section 2]). Denote by F(xn1) the set of all distinct monic irreducible factors of xn1 in a given finite field, Φr(x) a r-th cyclotomic polynomial, and φ() the Euler function. Write

    n=mpe,

    where e0 is an integer, p is the characteristic of Fq, and pm. Below are the known results for

    |F(xn1)|=1,2.

    Theorem 1.5. ([11]) The following statements are equivalent:

    (a) |F(xn1)|=1.

    (b) F(xn1) = {x1}.

    (c) n=pe.

    Theorem 1.6. ([11]) The following statements are equivalent:

    (a) |F(xn1)|=2.

    (b) F(xn1) = {x1,1+x++xm1}.

    (c) m is a prime different from p, and q is a primitive root modulo m.

    We summarize the five theorems above into the theorem below:

    Theorem 1.7. The following statements are equivalent:

    (a) Every irreducible polynomial of degree n over Fq with a nonzero trace is a normal polynomial.

    (b) F(xn1){x1,1+x++xn1}.

    (c) (c1) n=pe, or

    (c2) n is a prime different from p, with q being a primitive root modulo n.

    Cao [11] presented a new and unified proof of Theorem 1.6 and also extended Theorem 1.6. In this paper, we give the necessary and sufficient condition for the polynomial xn1 to have s different irreducible factors for a given positive integer s.

    Lemma 2.1 indicates that factorization of xn1 in finite fields is closely related to the cyclotomic polynomials.

    Lemma 2.1. ([12]) Let Fq be a finite field of characteristic p, and let n be a positive integer not divisible by p. Then

    xn1=d|nΦd(x).

    Lemma 2.2. ([13]) Let l be the order of a modulo m, an1 (mod m), then ln.

    Lemma 2.3. ([12]) Let Fq be a finite field and n a positive integer with (q,n)=1. Then the cyclotomic polynomial Φn(x) factors into φ(n)d distinct monic irreducible polynomials in Fq[x] of the same degree d, where d is the order of q modulo n.

    Lemma 2.4 is the well-known theorem about the existence of primitive roots.

    Lemma 2.4. ([14]) Let n be a positive integer. Then n possesses primitive roots if and only if n is of the form 2, 4, pα, or 2pα, where p is an odd prime and α is a positive integer.

    By Lemmas 2.3 and 2.4, we have the following lemma:

    Lemma 2.5. The cyclotomic polynomial Φn(x) is irreducible over Fq if and only if n=2,4,pα,2pα, and q is a primitive root modulo n.

    We define vp(x) to be the greatest power in which a prime p divides x, that is, if vp(x)=α, then pα|x but pα+1x. The following lemma is called the lifting the exponent lemma (LTE):

    Lemma 2.6. ([15]) Let x and y be (not necessarily positive) integers, n be a positive integer, and p be an odd prime such that p|xy and none of x and y is divisible by p. We have

    vp(xnyn)=vp(xy)+vp(n).

    Lemma 2.7. Let p be an odd prime and g be a primitive root modulo p2, with (g,p)=1. Then g is a primitive root modulo pl(l1).

    Proof. We first prove that g is a primitive root modulo pl(l2) by induction on l. Let g be a primitive root modulo pl. The order of g modulo pl+1 is d. We have

    φ(pl)|d,d|φ(pl+1),

    which shows that

    d=pl1(p1)

    or

    d=pl(p1).

    We next prove that

    dpl1(p1).

    According to Euler's theorem, we have

    gpl2(p1)1(mod pl1),

    there exists an integer k such that

    gpl2(p1)=1+kpl1,

    since g is a primitive root modulo pl, we have

    gpl2(p1)1(mod pl),plkpl1,(k,p)=1.

    Obviously, for l2, we have

    2l1l+1

    and

    3(l1)l+1,

    which shows that

    gpl1(p1)=(1+kpl1)p=1+kpl+k2p(p1)2p2(l1)+tp3(l1)+1+kpl(modpl+1),

    where t is an integer and (k,p)=1. It follows that

    gpl1(p1)1+kpl1(mod pl+1).

    Therefore

    dpl1(p1),d=pl(p1),

    and g is a primitive root modulo pl+1.

    We next use the LTE to prove that g is a primitive root modulo p. The Euler's theorem shows that

    gφ(p)1(modp),

    since g is a primitive root modulo p2, we have

    vp(gφ(p)1)=1.

    Let h be the order of g modulo p, then h|φ(p). Similarly, we can obtain

    vp(gh1)=1.

    Let x=gh, y=1, then (p,gh)=1, p|gh1. By the LTE, we have

    vp(gph1)=vp(gh1)+vp(p)=1+1=2,

    which shows that

    p2|gph1,gph1(mod p2),

    if h<φ(p), then ph<φ(p2), a contradiction. Thus h=φ(p), g is a primitive root modulo p.

    Lemma 2.8. Let g be a primitive root modulo pl. Then g is a primitive root modulo 2pl, where g is odd and p and l are the same as mentioned above.

    Proof. Let s be the order of g modulo 2pl. Then

    gs1(mod 2pl),s|φ(2pl).

    So we have

    gs1(mod pl).

    Since g is a primitive root modulo pl, we have φ(pl)|s and hence φ(2pl)|s. So

    s=φ(2pl)

    and g is a primitive root modulo 2pl.

    Combining Lemmas 2.1 and 2.3, we can calculate the number of different irreducible factors for xn1 over Fq. Let m be a positive integer,

    m=pα11pα22pαll

    be its prime decomposition. If

    n=mpe,

    where e0 is an integer, and p is the characteristic of Fq with pm. Then we can calculate that

    |F(xn1)|=mmφ(m)dm,

    where dm denotes the order of q modulo m. Note that there are lk=1(αk+1) items in the summation, where lk=1(αk+1) is the number of factors of m.

    Now assume that s and m are given positive integers, m1,,mt are the t factors of m. Thus, we have

    xn1=ti=1Φmi(x)pe.

    If xn1 factors into s distinct irreducible polynomials in Fq[x], then

    |F(xn1)|=ti=1φ(mi)dmi=s, (3.1)

    there are t items in the summation, the necessary and sufficient condition for

    |F(xn1)|=s

    is determined as follows:

    Observe φ(mi)dmi (i=1,2,,t), and φ(mi) are known; with the difference of q, the values of dmi will also change, that is, the values of φ(mi)dmi will change in the different Fq. Therefore, the t items in (3.1) can be regarded as t variables, and (3.1) can be regarded as the Diophantine equation with t variables

    x1+x2++xt=s. (3.2)

    Remark 3.1. Combining Lemma 2.2 and Euler's theorem, we know that φ(mi)dmi (i=1,2,,t) are positive integers. So we only need to consider the positive integer solutions of (3.2).

    Remark 3.2. The positive integers s and t satisfy ts. Otherwise, if t>s, then it follows from Lemmas 2.1 and 2.3 that

    |F(xn1)|>s.

    Remark 3.3. As we all know

    Φ1(x)=x1

    is a factor of xn1, and x1 is irreducible over Fq; the order of q modulo 1 is d1=1. Thus, at least one positive integer solution of (3.2) whose value is 1.

    We can find the positive integer solutions of (3.2). Without loss of generality, we suppose that

    xi=φ(mi)dmi=ki (i=1,2,,t) (3.3)

    is the positive integer solution of (3.2), thus we have

    dmi=φ(mi)ki (i=1,2,,t). (3.4)

    If there exists q such that the order of q modulo mi is

    dmi=φ(mi)ki  (i=1,2,,t),

    then it follows from Lemma 2.3 that Φmi(x) factors into

    φ(mi)φ(mi)ki=ki

    distinct monic irreducible polynomials over Fq of the same degree

    dmi=φ(mi)ki.

    Therefore, we have

    |F(xn1)|=k1+k2++kt=s,

    that is, xn1 factors into s distinct irreducible polynomials over Fq.

    In conclusion, we have the following result:

    Theorem 3.4. (Main result) Let Fq denote the finite field of q elements with characteristic p. Let p be a prime. Let

    n=mpe

    with e0, pm. Let s be a positive integer. The t factors of m are m1,m2,,mt and dmi denotes the order of q modulo mi. Then

    |F(xn1)|=s,

    if and only if

    xi=φ(mi)dmi  (i=1,2,,t)

    is a solution to the Diophantine equation

    x1+x2++xt=s.

    Proof. We first assume that

    xi=φ(mi)dmi=ki  (i=1,2,,t)

    is the solution of the Diophantine equation

    x1+x2++xt=s.

    By Lemmas 2.1 and 2.3, we have

    xn1=ti=1Φmi(x)pe

    and Φmi(x) factors into

    φ(mi)φ(mi)ki=ki

    distinct monic irreducible polynomials in Fq[x] of the same degree

    dmi=φ(mi)ki  (i=1,2,,t).

    Hence, we obtain

    |F(xn1)|=ti=1ki=s.

    Suppose

    |F(xn1)|=s.

    According to Lemma 2.3, we have

    |F(xn1)|=ti=1φ(mi)dmi=s.

    So

    xi=φ(mi)dmi  (i=1,2,,t)

    is the solution of the Diophantine equation

    x1+x2++xt=s.

    We apply Theorem 3.4 to deduce Theorem 1.6. Note that Theorem 1.5 is trivial. Recall that

    n=mpe,

    where e0 is an integer and p is the characteristic of Fq with pm.

    For Theorem 1.6

    |F(xn1)|=2,

    we know that m has two factors; thus, m is a prime different from p,

    xn1=(Φ1(x)Φm(x))pe.

    The unique positive integer solution of the Diophantine equation

    x1+x2=2

    is

    x1=x2=1.

    Since

    φ(1)d1=1,

    and by Theorem 3.4, we have

    φ(m)dm=1,

    that is, dm=φ(m), q is a primitive root modulo m. It follows from Lemma 2.5 that the cyclotomic polynomial Φm(x) is irreducible over Fq. Thus

    |F(xn1)|=2.

    The necessary and sufficient conditions for

    |F(xn1)|=3,4,5

    are given, respectively, below: for

    |F(xn1)|=3,

    we know that t3, where t is the number of factors of m.

    Case 1. If m has three factors, then m=r2, where r is a prime different from p.

    xn1=(Φ1(x)Φr(x)Φr2(x))pe.

    The unique positive integer solution of the Diophantine equation

    x1+x2+x3=3

    is

    x1=x2=x3=1.

    By Theorem 3.4, we have

    φ(r)dr=1,  φ(r2)dr2=1,

    that is,

    dr=φ(r),  dr2=φ(r2),

    q is a primitive root modulo r and r2. Recall Lemma 2.7: If q is a primitive root modulo r2, then q is a primitive root modulo r. Cyclotomic polynomials Φr(x) and Φr2(x) are irreducible over Fq. Thus

    |F(xn1)|=3.

    Case 2. If m has two factors, then m=r,

    xn1=(Φ1(x)Φr(x))pe.

    The positive integer solution of the Diophantine equation

    x1+x2=3

    is x1=1, x2=2. Hence, we obtain

    φ(r)dr=2,

    that is,

    dr=φ(r)2,

    the order of q modulo r is φ(r)2. According to Lemma 2.3, the cyclotomic polynomial Φr(x) factors into

    φ(r)φ(r)2=2

    distinct monic irreducible polynomials in Fq[x] of the same degree

    dr=φ(r)2.

    Thus

    |F(xn1)|=3.

    In conclusion, we have the following result:

    Theorem 4.1. The following statements are equivalent:

    (a) |F(xn1)|=3.

    (b) (b1) F(xn1) = {x1,f1(x),f2(x)}, where m=r,

    f1(x)f2(x)=Φr(x),degf1=degf2=φ(r)2.

    (b2) F(xn1) = {x1,Φr(x),Φr2(x)}, where m=r2.

    (c) (c1) m=r, and the order of q modulo r is φ(r)2.

    (c2) m=r2, and q is a primitive root modulo r2.

    For

    |F(xn1)|=4,

    we know that t4, where t is the number of factors of m. In the remaining part of this paper, we always assume that r is an odd prime different from p.

    If m has four factors, then the possible values of m are r3, 2r, p1p2 or 8, where p1 and p2 are odd primes different from p. The unique positive integer solution of the Diophantine equation

    x1+x2+x3+x4=4

    is

    x1=x2=x3=x4=1.

    Case 1. If m=r3, then

    xn1=(Φ1(x)Φr(x)Φr2(x)Φr3(x))pe,

    and we have

    φ(r)dr=φ(r2)dr2=φ(r3)dr3=1,

    that is, q is a primitive root modulo rl,l=1,2,3, which requires that q is a primitive root modulo r2. The cyclotomic polynomials Φr(x), Φr2(x), and Φr3(x) are irreducible over Fq. Thus

    |F(xn1)|=4.

    Case 2. If m=2r, then

    xn1=(Φ1(x)Φ2(x)Φr(x)Φ2r(x))pe,

    and we have

    φ(2)d2=φ(r)dr=φ(2r)d2r=1,

    that is, q is a primitive root modulo r and 2r. Recall Lemma 2.8: If q is a primitive root modulo r, then q is a primitive root modulo 2r. The cyclotomic polynomials Φr(x) and Φ2r(x) are irreducible over Fq. Obviously, the order of q modulo 2 is d2=1, and

    Φ2(x)=x+1

    is irreducible over Fq. Thus

    |F(xn1)|=4.

    Case 3. If m=p1p2, then

    xn1=(Φ1(x)Φp1(x)Φp2(x)Φp1p2(x))pe,

    and we have

    φ(p1)dp1=φ(p2)dp2=φ(p1p2)dp1p2=1.

    It follows from Lemma 2.4 that p1p2 has no primitive root, which contradicts

    dp1p2=φ(p1p2).

    The cyclotomic polynomial Φp1p2(x) is reducible over Fq. Thus

    |F(xn1)|>4.

    Case 4. If m=8, then

    xn1=(Φ1(x)Φ2(x)Φ4(x)Φ8(x))pe.

    Since 8 has no primitive root, Φ8(x) is reducible over Fq. Thus

    |F(xn1)|>4.

    If m has three factors, then the possible values of m are r2 or 4. The positive integer solution of the equation

    x1+x2+x3=4

    is

    x1=x2=1,x3=2.

    Case 1. If m=r2, then

    xn1=(Φ1(x)Φr(x)Φr2(x))pe,

    and we have

    φ(r)dr=2,  φ(r2)dr2=1

    or

    φ(r)dr=1,  φ(r2)dr2=2.

    For the former, the order of q modulo r is

    dr=φ(r)2,

    and q is a primitive root modulo r2, which is impossible by Lemma 2.7. For the latter, q is a primitive root modulo r, and the order of q modulo r2 is φ(r2)2. It follows that

    qφ(r)1(mod r),  qφ(r2)2=qrφ(r)21(mod r2),  qrφ(r)21(mod r).

    Thus φ(r)|rφ(r)2. Since r is an odd prime, r2 is not an integer, and the divisibility is not valid.

    Case 2. If m=4, then

    xn1=(Φ1(x)Φ2(x)Φ4(x))pe,

    and we have

    φ(4)d4=2,

    that is, the order of q modulo 4 is

    d4=φ(4)2=1.

    The cyclotomic polynomial

    Φ4(x)=1+x2

    factors into

    φ(4)φ(4)2=2

    distinct monic irreducible polynomials in Fq[x] of the same degree, d4=1. Thus

    |F(xn1)|=4.

    If m has two factors, then m=r,

    xn1=(Φ1(x)Φr(x))pe.

    The positive integer solution of the Diophantine equation

    x1+x2=4

    is x1=1, x2=3. Hence, we have

    φ(r)dr=3,

    that is, the order of q modulo r is φ(r)3. The cyclotomic polynomial Φr(x) factors into three distinct monic irreducible polynomials in Fq[x] of the same degree φ(r)3. Thus

    |F(xn1)|=4.

    In conclusion, we have the following result:

    Theorem 4.2. The following statements are equivalent:

    (a) |F(xn1)|=4.

    (b) (b1) F(xn1) = {x1,f1(x),f2(x),f3(x)}, where m=r,

    f1(x)f2(x)f3(x)=Φr(x),degf1=degf2=degf3=φ(r)3.

    (b2) F(xn1) = {x1,Φr(x),Φr2(x),Φr3(x)}, where m=r3.

    (b3) F(xn1) = {x1,x+1,Φr(x),Φ2r(x)}, where m=2r.

    (b4) F(xn1) = {x1,x+1,x+e1,x+e2}, where m=4, e1 and e2 are integers.

    (c) (c1) m=r, and the order of q modulo r is φ(r)3.

    (c2) m=r3, and q is a primitive root modulo r2.

    (c3) m=2r, and q is a primitive root modulo r.

    (c4) m=4, and the order of q modulo 4 is 1.

    For

    |F(xn1)|=5,

    we know that t5, where t is the number of factors of m.

    If m has five factors, then the possible values of m are r4 or 16. The unique positive integer solution of the equation

    x1+x2+x3+x4+x5=5

    is

    x1=x2=x3=x4=x5=1.

    Case 1. If m=r4, then

    xn1=(Φ1(x)Φr(x)Φr2(x)Φr3(x)Φr4(x))pe,

    and we have

    φ(r)dr=φ(r2)dr2=φ(r3)dr3=φ(r4)dr4=1,

    that is, q is a primitive root modulo rl,l=1,2,3,4, which requires that q is a primitive root modulo r2. The cyclotomic polynomials Φr(x),Φr2(x),Φr3(x), and Φr4(x) are irreducible over Fq. Thus

    |F(xn1)|=5.

    Case 2. If m=16, then

    xn1=(Φ1(x)Φ2(x)Φ4(x)Φ8(x)Φ16(x))pe.

    Since 8 and 16 have no primitive root, we know that Φ8(x) and Φ16(x) are reducible over Fq. Thus

    |F(xn1)|>5.

    If m has four factors, then the possible values of m are r3,8,2r or p1p2. The positive integer solution of the equation

    x1+x2+x3+x4=5

    is

    x1=x2=x3=1,x4=2.

    Case 1. If m=r3, then

    xn1=(Φ1(x)Φr(x)Φr2(x)Φr3(x))pe.

    It follows from Lemma 2.7 that

    φ(r2)dr2=2

    and

    φ(r)dr=φ(r3)dr3=1,

    that is, q is a primitive root modulo r3, the order of q modulo r2 is φ(r2)2. It follows that

    qφ(r3)1(mod r3),  qφ(r2)21(mod r2),  qφ(r2)21(mod r).

    Since q is a primitive root modulo r3, we have

    vr(qφ(r2)21)=2.

    By Lemma 2.6, we have

    vr(qrφ(r2)21)=vr(qφ(r2)21)+vr(r)=2+1=3,
    qφ(r3)21(modr3),

    which is a contradiction. Thus, if m=r3, then

    |F(xn1)|5.

    Case 2. If m=8, then

    xn1=(Φ1(x)Φ2(x)Φ4(x)Φ8(x))pe,

    and we have

    φ(4)d4=1,  φ(8)d8=2,

    that is, q is a primitive root modulo 4, which means q is congruent to 3 modulo 4, therefore the order of q modulo 8 is

    φ(8)2=2.

    Φ4(x) is irreducible over Fq, and Φ8(x) factors into

    φ(8)φ(8)2=2

    distinct monic irreducible polynomials in Fq[x] of the same degree 2. Thus

    |F(xn1)|=5.

    Case 3. If m=2r, then

    xn1=(Φ1(x)Φ2(x)Φr(x)Φ2r(x))pe.

    It follows from Lemma 2.8 that

    φ(r)dr=2,  φ(2r)d2r=1,

    that is, q is a primitive root modulo 2r, and the order of q modulo r is φ(r)2, Φ2r(x) is irreducible over Fq and Φr(x) factors into

    φ(r)φ(r)2=2

    distinct monic irreducible polynomials in Fq[x] of the same degree φ(r)2. Thus

    |F(xn1)|=5.

    Case 4. If m=p1p2, then

    xn1=(Φ1(x)Φp1(x)Φp2(x)Φp1p2(x))pe.

    Since p1p2 has no primitive root, we have

    φ(p1)dp1=φ(p2)dp2=1,  φ(p1p2)dp1p2=2,

    that is, q is a primitive root modulo p1 and p2, and the order of q modulo p1p2 is φ(p1p2)2, Φp1(x) and Φp2(x) are irreducible over Fq, and Φp1p2(x) factors into

    φ(p1p2)φ(p1p2)2=2

    distinct monic irreducible polynomials in Fq[x] of the same degree φ(p1p2)2. Thus

    |F(xn1)|=5.

    If m has three factors, then the possible values of m are 4 or r2. If m=4, there are at most four distinct irreducible factors for xn1. Thus

    |F(xn1)|<5.

    If m=r2, then

    xn1=(Φ1(x)Φr(x)Φr2(x))pe.

    The positive integer solutions of the equation

    x1+x2+x3=5

    are

    x1=1,  x2=x3=2,orx1=x2=1,  x3=3.

    For the former, we have

    φ(r)dr=φ(r2)dr2=2,

    that is, the order of q modulo r is φ(r)2 and the order of q modulo r2 is φ(r2)2, Φr(x), and Φr2(x) factor into 2 distinct monic irreducible polynomials in Fq[x]. Thus

    |F(xn1)|=5.

    For the latter, it follows from Lemma 2.7 that

    φ(r)dr=1andφ(r2)dr2=3,

    that is, q is a primitive root modulo r, and the order of q modulo r2 is φ(r2)3, Φr(x) is irreducible over Fq and Φr2(x) factors into 3 distinct monic irreducible polynomials in Fq[x]. Thus

    |F(xn1)|=5.

    If m has two factors, then

    m=r,  xn1=(Φ1(x)Φr(x))pe.

    The positive integer solution of the equation x1+x2=5 is x1=1,x2=4. Hence, we have

    φ(r)dr=4,

    that is, the order of q modulo r is φ(r)4. The cyclotomic polynomial Φr(x) factors into four distinct monic irreducible polynomials in Fq[x]. Thus

    |F(xn1)|=5.

    In conclusion, we obtain the following result:

    Theorem 4.3. The following statements are equivalent:

    (a) |F(xn1)|=5.

    (b) (b1) F(xn1) = {x1,Φr(x),Φr2(x),Φr3(x),Φr4(x)}, where m=r4.

    (b2) F(xn1) = {x1,x+1,Φ4(x),f1(x),f2(x)}, where m=8,

    f1(x)f2(x)=Φ8(x),degf1=degf2=2.

    (b3) F(xn1) = {x1,x+1,g1(x),g2(x),Φ2r(x)}, where m=2r,

    g1(x)g2(x)=Φr(x),degg1=degg2=φ(r)2.

    (b4) F(xn1) = {x1,Φp1(x),Φp2(x),h1(x),h2(x)}, where m=p1p2,

    h1(x)h2(x)=Φp1p2(x),degh1=degh2=φ(p1p2)2.

    (b5) F(xn1) = {x1,k1(x),k2(x),r1(x),r2(x)}, where m=r2,

    k1(x)k2(x)=Φr(x),degk1=degk2=φ(r)2andr1(x)r2(x)=Φr2(x),degr1=degr2=φ(r2)2.

    Or F(xn1) = {x1,Φr(x),u1(x),u2(x),u3(x)}, where m=r2,

    u1(x)u2(x)u3(x)=Φr2(x),degu1=degu2=degu3=φ(r2)3.

    (b6) F(xn1) = {x1,v1(x),v2(x),v3(x),v4(x)}, where m=r,

    v1(x)v2(x)v3(x)v4(x)=Φr(x),degv1=degv2=degv3=degv4=φ(r)4.

    (c) (c1) m=r4, and q is a primitive root modulo r2.

    (c2) m=8, q is a primitive root modulo 4, and the order of q modulo 8 is 2.

    (c3) m=2r, the order of q modulo r is φ(r)2 and q is a primitive root modulo 2r.

    (c4) m=p1p2, the order of q modulo p1p2 is φ(p1p2)2, and q is a primitive root modulo p1 and p2.

    (c5) m=r2, the order of q modulo r is φ(r)2 and the order of q modulo r2 is φ(r2)2; or q is a primitive root modulo r and the order of q modulo r2 is φ(r2)3.

    (c6) m=r, and the order of q modulo r is φ(r)4.

    In this final section, we provide two examples.

    Example 5.1. Let Fq denote the finite field of q elements with characteristic p. Let p be a prime, let r be a prime different from p, and let n=rlpe, with l1, e0. Denote by F(xn1) the set of all distinct monic irreducible factors of xn1 in a given finite field. Given a positive integer s, we consider the special case for

    |F(xn1)|=s=l+1.

    Since n=rlpe,

    xn1=(Φ1(x)Φr(x)Φr2(x)Φrl(x))pe.

    The unique positive integer solution of the equation

    x0+x1+x2++xl=l+1

    is

    xi=φ(ri)dri=1  (i=0,1,,l),

    where dri denotes the order of q modulo ri (i=0,1,,l). That is, q is a primitive root modulo rj (j=1,2,,l). Recall Lemma 2.7 that if q is a primitive root modulo r2, then q is a primitive root modulo rj (j=1,2,,l). Cyclotomic polynomials Φrj(x) (j=1,2,,l) are irreducible over Fq.

    In conclusion, if q is a primitive root modulo r2, then

    |F(xrlpe1)|=l+1.

    Example 5.2. We factor polynomial x251 into distinct monic irreducible polynomials over F7. Since 25=52,

    x251=x521=Φ1(x)Φ5(x)Φ25(x).

    We first calculate

    72(mod5),  7241(mod5),  741(mod5),

    where 7 is a primitive root modulo 5, thus

    Φ5(x)=1+x+x2+x3+x4

    is irreducible over F7.

    We next calculate

    77(mod25),  72241(mod25),  741(mod25),

    the order of  7 modulo 25 is 4, thus

    Φ25(x)=1+x5+x10+x15+x20

    factors into

    φ(25)4=5

    distinct monic irreducible polynomials over F7 of the same degree 4,

    f1(x)=1+2x+4x2+2x3+x4,f2(x)=1+4x+4x3+x4,f3(x)=1+4x+3x2+4x3+x4,f4(x)=1+5x+5x2+5x3+x4,f5(x)=1+6x+5x2+6x3+x4,

    respectively.

    Thus

    x251=((x1)(1+x+x2+x3+x4)(1+2x+4x2+2x3+x4)(1+4x+4x3+x4)(1+4x+3x2+4x3+x4)(1+5x+5x2+5x3+x4)(1+6x+5x2+6x3+x4))

    and

    |F(x251)|=7.

    {Let Fq be the finite field of q elements, and Fqn be its extension of degree n. Denote by F(xn1) the set of all distinct monic irreducible factors of the polynomial xn1 in the finite field Fq. Given a positive integer s, we use the properties of cyclotomic polynomials in finite fields and results from the Diophantine equations to provide the sufficient and necessary condition for

    |F(xn1)|=s.

    As an application, we also obtain the sufficient and necessary conditions for

    |F(xn1)|=3,4,5.

    Weitao Xie: the first draft of the manuscript; Jiayu Zhang: preliminaries collection and analysis; Wei Cao: originally raised the problem and commented on previous versions of the manuscript. All authors contributed to the study conception and design. All authors read and approved the final manuscript.

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

    The authors thank the anonymous referees for their helpful comments that improved the quality of the manuscript. This work was jointly supported by the Fujian Provincial Natural Science Foundation of China (Grant No. 2022J02046), and Fujian Key Laboratory of Granular Computing and Applications (Minnan Normal University).

    All authors declare no conflicts of interest in this paper.



    [1] A. Arikoglu, I. Ozkol, Solution of fractional differential equations by using differential transform method, Chaos Solitons Fractals, 34 (2007), 1473–1481.
    [2] Y. Gao, J. G. Liu, T. Luo, Y. Xiang, Revisit of the Peierls-Nabarro model for edge dislocations in Hilbert space, Discrete Contin. Dyn. Syst. B, 26 (2021), 3177–3207.
    [3] N. Sene, A novel fractional-order system described by the Caputo derivative, its numerical discretization, and qualitative properties, in Handbook of Fractional Calculus for Engineering and Science, Chapman and Hall/CRC, (2022), 205–240.
    [4] H. Dong, Y. Gao, Existence and uniqueness of bounded stable solutions to the Peierls-Nabarro model for curved dislocations, Calculus Var. Partial Differ. Equations, 60 (2021), 1–26.
    [5] S. Sharma, R. Bairwa, A reliable treatment of iterative laplace transform method for fractional telegraph equations, Ann. Pure Appl. Math., 9 (2015), 81–89.
    [6] M. S. Mohamed, Y. S. Hamed, Solving the convection–diffusion equation by means of the optimal q-homotopy analysis method, Results Phys., 6 (2016), 20–25.
    [7] M. Yavuz, N. Sene, Approximate solutions of the model describing fluid flow using generalized ρ-Laplace transform method and heat balance integral method, Axioms 9 (2020), 123.
    [8] N. Sene, A. N. Fall, Homotopy perturbation ρ-Laplace transform method and its application to the fractional diffusion equation and the fractional diffusion-reaction equation, Fractal Fractional, 3 (2019), 14.
    [9] L. Wang, Y. Ma, Z. Meng, Haar wavelet method for solving fractional partial differential equations numerically, Appl. Math. Comput., 227 (2014), 66–76.
    [10] M. Khader, N. Sweilam, On the approximate solutions for system of fractional integro-differential equations using chebyshev pseudo-spectral method, Appl. Math. Modell., 37 (2013), 9819–9828.
    [11] M. Osman, H. Almusawa, K. U. Tariq, S. Anwar, S. Kumar, M. Younis, et al., On global behavior for complex soliton solutions of the perturbed nonlinear schrödinger equation in nonlinear optical fibers, J. Ocean Eng. Sci., 7 (2022), 431–443.
    [12] N. Heymans, I. Podlubny, Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives, Rheol. Acta, 45 (2006), 765–771.
    [13] M. Al-Refai, A. M. Jarrah, Fundamental results on weighted Caputo–Fabrizio fractional derivative, Chaos Solitons Fractals, 126 (2019), 7–11.
    [14] F. Jarad, T. Abdeljawad, Z. Hammouch, On a class of ordinary differential equations in the frame of Atangana–Baleanu fractional derivative, Chaos Solitons Fractals, 117 (2018), 16–20.
    [15] J. J. Nieto, A. Ouahab, V. Venktesh, Implicit fractional differential equations via the Liouville–Caputo derivative, Mathematics, 3 (2015), 398–411.
    [16] K. A. Abro, A. Atangana, A comparative analysis of electromechanical model of piezoelectric actuator through Caputo–Fabrizio and Atangana–Baleanu fractional derivatives, Math. Methods Appl. Sci., 43 (2020), 9681–9691.
    [17] M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Prog. Fractional Differ. Appl., 1 (2015), 73–85.
    [18] M. Toufik, A. Atangana, New numerical approximation of fractional derivative with non-local and non-singular kernel: Application to chaotic models, Eur. Phys. J. Plus, 132 (2017), 1–16.
    [19] Y. Gao, L. Jin, T. J. Xiao, A new method to obtain uniform decay rates for multidimensional wave equations with nonlinear acoustic boundary conditions, SIAM J. Control Optim., 56 (2018), 1303-1320.
    [20] Y. Gao, J. G. Liu, Z. Liu, Existence and rigidity of the vectorial Peierls-Nabarro model for dislocations in high dimensions, Nonlinearity, 34 (2021), 7778.
    [21] I. Jaradat, M. Alquran, R. A. Muhsen, An analytical framework of 2D diffusion, wave-like, telegraph, and Burgers' models with twofold Caputo derivatives ordering, Nonlinear Dyn., 93 (2018), 1911–1922.
    [22] M. Dehghan, F. Shakeri, Application of He's variational iteration method for solving the Cauchy reaction–diffusion problem, J. Comput. Appl. Math., 214 (2008), 435–446.
    [23] B. K. Singh, V. K. Srivastava, Approximate series solution of multi-dimensional, time fractional-order (heat-like) diffusion equations using FRDTM, R. Soc. Open Sci., 2 (2015), 140511.
    [24] K. Shah, H. Khalil, R. A. Khan, Analytical solutions of fractional order diffusion equations by natural transform method, Iran. J. Sci. Technol. Trans. A Sci., 42 (2018), 1479–1490.
    [25] D. Kumar, J. Singh, S. Kumar, Analytical modeling for fractional multi-dimensional diffusion equations by using Laplace transform, Commun. Numer. Anal., 1 (2015), 16–29.
    [26] J. H. He, Homotopy perturbation method: a new nonlinear analytical technique, Appl. Math. Comput., 135 (2003), 73–79.
    [27] J. H. He, M. L. Jiao, K. A. Gepreel, Y. Khan, Homotopy perturbation method for strongly nonlinear oscillators, Math. Comput. Simul., 204 (2023), 243–258.
    [28] Z. Odibat, S. Momani, Modified homotopy perturbation method: application to quadratic riccati differential equation of fractional order, Chaos Solitons Fractals, 36 (2008), 167–174.
    [29] T. M. Elzaki, The new integral transform elzaki transform, Global J. Pure Appl. Math., 7 (2011), 57–64.
    [30] S. Rashid, S. Sultana, B. Kanwal, F. Jarad, A. Khalid, Fuzzy fractional estimates of swift-hohenberg model obtained using the Atangana-Baleanu fractional derivative operator, AIMS Math., 7 (2022), 16067–16101.
    [31] H. Yasmin, N. Iqbal, A comparative study of the fractional coupled burgers and Hirota–Satsuma KDV equations via analytical techniques, Symmetry, 14 (2022), 1364.
    [32] M. Nadeem, J. H. He, A. Islam, The homotopy perturbation method for fractional differential equations: Part 1 mohand transform, Int. J. Numer. Methods Heat Fluid Flow, 31 (2021), 3490–3504.
    [33] V. Morales-Delgado, J. Gómez-Aguilar, S. Kumar, M. Taneco-Hernández, Analytical solutions of the keller-segel chemotaxis model involving fractional operators without singular kernel, Eur. Phys. J. Plus, 133 (2018), 200.
    [34] N. Iqbal, S. Alshammari, T. Botmart, Evaluation of regularized long-wave equation via caputo and caputo-fabrizio fractional derivatives, AIMS Math., 7 (2022), 20401–20419.
    [35] N. A. Shah, I. Dassios, J. D. Chung, Numerical investigation of time-fractional equivalent width equations that describe hydromagnetic waves, Symmetry, 13 (2021), 418.
  • This article has been cited by:

    1. Lingxi Gao, Bochi Sun, Ziqing Du, Guangming Lv, How Wealth Inequality Affects Happiness: The Perspective of Social Comparison, 2022, 13, 1664-1078, 10.3389/fpsyg.2022.829707
    2. Wenhui Li, Qi Zhu, Fenghua Wen, Normaziah Mohd Nor, The evolution of day-of-the-week and the implications in crude oil market, 2022, 106, 01409883, 105817, 10.1016/j.eneco.2022.105817
  • Reader Comments
  • © 2023 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)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1668) PDF downloads(66) Cited by(15)

Figures and Tables

Figures(6)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog