Research article Special Issues

Self-mastery among Chinese Older Adults in the Greater Chicago Area

  • Background: Self-mastery is an important psychological resource to cope with stressful situations. However, we have limited understanding of self-mastery among minority aging populations. Objective: This study aims to examine the presence and levels of self-mastery among U.S. Chinese older adults. Methods: Data were drawn from the PINE study, a population-based survey of U.S. Chinese older adults in the Greater Chicago area. Guided by a community-based participatory research approach, a total of 3,159 Chinese older adults aged 60 and above were surveyed. A Chinese version of the Self-Mastery Scale was used to assess self-mastery. Results: Out of the 7-item Chinese Self-Mastery Scale, approximately 42.8% to 87.5% of Chinese older adults experienced some degree of self-mastery in their lives. Older adults with no formal education and the oldest-old aged 85 and over had the lowest level of self-mastery in our study. A higher mastery level was associated with being married, having fewer children, better self-reported health status, better quality of life, and positive health changes. Conclusion: Although self-mastery is commonly experienced among the Chinese aging population in the Greater Chicago area, specific subgroups are still vulnerable. Future longitudinal studies are needed to improve the understanding of risk factors and outcomes associated with self-mastery among Chinese older adults.

    Citation: Xinqi Dong, Manrui Zhang, Melissa A. Simon. Self-mastery among Chinese Older Adults in the Greater Chicago Area[J]. AIMS Medical Science, 2014, 1(1): 57-72. doi: 10.3934/medsci.2014.1.57

    Related Papers:

    [1] Gaixiang Cai, Fengru Xiao, Guidong Yu . The identification numbers of lollipop graphs. AIMS Mathematics, 2025, 10(4): 7813-7827. doi: 10.3934/math.2025358
    [2] Ali N. A. Koam, Ali Ahmad, Azeem Haider, Moin A. Ansari . Computation of eccentric topological indices of zero-divisor graphs based on their edges. AIMS Mathematics, 2022, 7(7): 11509-11518. doi: 10.3934/math.2022641
    [3] Bilal A. Rather, M. Aijaz, Fawad Ali, Nabil Mlaiki, Asad Ullah . On distance signless Laplacian eigenvalues of zero divisor graph of commutative rings. AIMS Mathematics, 2022, 7(7): 12635-12649. doi: 10.3934/math.2022699
    [4] Mohammad Hamidi, Irina Cristea . Hyperideal-based zero-divisor graph of the general hyperring $ \mathbb{Z}_{n} $. AIMS Mathematics, 2024, 9(6): 15891-15910. doi: 10.3934/math.2024768
    [5] Meijin Luo, Qiutao Qin . Exponents of a class of special three-colored primitive digraphs with $ n $ vertices in graph theory. AIMS Mathematics, 2025, 10(4): 9415-9434. doi: 10.3934/math.2025435
    [6] Yunfeng Tang, Huixin Yin, Miaomiao Han . Star edge coloring of $ K_{2, t} $-free planar graphs. AIMS Mathematics, 2023, 8(6): 13154-13161. doi: 10.3934/math.2023664
    [7] Hongyu Chen, Li Zhang . A smaller upper bound for the list injective chromatic number of planar graphs. AIMS Mathematics, 2025, 10(1): 289-310. doi: 10.3934/math.2025014
    [8] Jin Cai, Shuangliang Tian, Lizhen Peng . On star and acyclic coloring of generalized lexicographic product of graphs. AIMS Mathematics, 2022, 7(8): 14270-14281. doi: 10.3934/math.2022786
    [9] Xiaoxue Hu, Jiangxu Kong . An improved upper bound for the dynamic list coloring of 1-planar graphs. AIMS Mathematics, 2022, 7(5): 7337-7348. doi: 10.3934/math.2022409
    [10] Jianxin Luo, Jiangxu Kong . An upper bound for the semistrong chromatic index of Halin graphs. AIMS Mathematics, 2025, 10(7): 15811-15820. doi: 10.3934/math.2025708
  • Background: Self-mastery is an important psychological resource to cope with stressful situations. However, we have limited understanding of self-mastery among minority aging populations. Objective: This study aims to examine the presence and levels of self-mastery among U.S. Chinese older adults. Methods: Data were drawn from the PINE study, a population-based survey of U.S. Chinese older adults in the Greater Chicago area. Guided by a community-based participatory research approach, a total of 3,159 Chinese older adults aged 60 and above were surveyed. A Chinese version of the Self-Mastery Scale was used to assess self-mastery. Results: Out of the 7-item Chinese Self-Mastery Scale, approximately 42.8% to 87.5% of Chinese older adults experienced some degree of self-mastery in their lives. Older adults with no formal education and the oldest-old aged 85 and over had the lowest level of self-mastery in our study. A higher mastery level was associated with being married, having fewer children, better self-reported health status, better quality of life, and positive health changes. Conclusion: Although self-mastery is commonly experienced among the Chinese aging population in the Greater Chicago area, specific subgroups are still vulnerable. Future longitudinal studies are needed to improve the understanding of risk factors and outcomes associated with self-mastery among Chinese older adults.


    Some human pathogens, such as hepatitis B virus (HBV), hepatitis C virus (HCV) and human immunodeficiency virus (HIV), possess the ability to suppress the immune response and damage the immune system, ultimately leading to the development of certain chronic diseases and even tumors [1]. In recent years, researchers have extensively studied numerous models pertaining to these infections [2,3,4]. In recent years, there has been notable progress in treating these infections. For example, the latest direct-acting antiviral drugs (DAA) for HCV can cure the majority of cases quickly, with cure rates reaching up to 95%[5,6]. For HBV, the optimal approach involves long-term medication to manage the condition, reduce viral transmission and stabilize host immune system[7,8]. Similarly, antiretroviral therapy (ART) is commonly used to suppress HIV replication[9]. Notably, a small group of elite controllers can naturally manage HIV without treatment[10]. However, challenges persist, including drug resistance and the risk of viral replication or recurrence[5,6,7,8,9]. In the past few years, immune checkpoint therapy, a method that shows significant efficacy in completely eliminating tumor cells for certain types of tumors, has been discovered, which activates the host's own immune system to inhibit tumor growth and spread [11,12,13]. Enhancing antiviral treatment by strengthening specific immune responses has become a hot topic among researchers.

    In 2003, Komarova et al. [1] first established an immunosuppressive infection model to investigate the relationship between drug efficacy and treatment duration, the model is given as follows:

    $ {dydt=ry(1yK)aypyz,dzdt=cyz1+ηybzqyz:=zf(y),
    $
    (1.1)

    where $ y $, $ z $ denote the population size of viruses and immune cells at time $ t $, respectively. Viral reproduction follows logistic growth, $ r $ represents the viral replication rate at a viral load $ K $. It is hypothesized that the virus population undergoes decay at a rate $ a $, and immune cells eliminate the viruses at a rate $ pyz $. Furthermore, immune cells are suppressed by the viruses at a rate $ qyz $ and die at a rate $ b $. In the context of antigenic stimulation, the proliferation of immune cells relies on both the viruses and immune cells, and it can be quantified by the expansion function

    $ \Phi(y, z) = \frac{cyz}{1+\eta y}, $

    where $ \eta $ represents the inhibitory effects of the viruses on the proliferation of immune cells.

    In 2014, Shu et al. [14] introduced a time delay in the immune expansion function of the model described in [1] and studied the existence of Hopf bifurcation and periodic solutions. In 2017, Tian et al. [15] incorporated diffusion into the model presented in [14] and studied the direction and stability of Hopf bifurcations. However, viruses require a certain amount of time to effectively stimulate immune cells. We refer to this time period as the virus stimulation delay, denoted as $ \tau $. The immune cells at time $ t $ are essentially activated by the viruses at time $ t-\tau $. Therefore, in 2021, Li et al. [16] studied the function as follows:

    $ Φ(yτ,z)=cy(tτ)z1+ηy(tτ).
    $

    However, it has been noted by Wodarz (see pp 31, Chapter 2 in [17]) that the growth rate of immune cells, denoted as $ zf(y) $, in all of these models was assumed to be directly proportional to $ z $, but nonlinearity can better capture the realistic dynamics. Thus, a more realistic model, proposed by De Boer and Perelson [18,19] in 1995–1997, has been suggested to provide a better representation of the actual dynamics. In this model, the immune expansion function is described by the well-known Beddington-DeAngelis functional response [20,21], the form of which is as follows:

    $ Φ(y,z)=cyz1+ηy+sz.
    $

    It was initially used in predator-prey models [22,23,24]. The parameter $ s $ in the function represents the strength of competition between immune responses. This immune expansion function provides a more biologically realistic description of the interactions and competitive relationships between virus populations and immune cell populations within a host [25]. Thus the Beddington-DeAngelis functional response has also been widely discussed in immunosuppressive infection models [26,27,28].

    In 2023, Chen et al. [29] introduced the Beddington-DeAngelis functional response into the immunosuppressive infection model established by Komarova et al. [1], the model is as follows:

    $ {dydt=ry(1yK)aypyz,dzdt=cy(tτ)z1+ηy(tτ)+szbzqyz.
    $
    (1.2)

    In [29], Chen et al. discussed the local and global stability of equilibrium points in the absence of delay, as well as the stability and Hopf bifurcation of model (1.2) when delay is present. They designed antiviral drug treatment strategies from the perspectives of adaptivity and feasibility of immune competition strength and viral suppression intensity.

    In practical problems, the spread of viruses within a host is possible. Thus we introduce the diffusion into the model proposed by Chen et al. in [29]. As far as we know, this task has not been attempted yet. Therefore, we study the reaction-diffusion system under the Neumann boundary condition as follows:

    $ {y(x,t)t=d1Δy(x,t)+ry(x,t)(1y(x,t)K)ay(x,t)py(x,t)z(x,t),z(x,t)t=d2Δz(x,t)+cy(x,tτ)z(x,t)1+ηy(x,tτ)+sz(x,t)bz(x,t)qy(x,t)z(x,t),νy(x,t)=νz(x,t)=0,xΩ,t>0,y(x,t)=ψ1(x,t),z(x,t)=ψ2(x,t),xΩ,t[τ,0],
    $
    (1.3)

    where $ \Omega = (0, l \pi) $ with $ l > 0 $, $ y = y(x, t) $ and $ z = z(x, t) $ denoting the population densities of viruses and immune cells at location $ x $ and time $ t $, respectively. The parameters $ s $ and $ \tau $ stand for the competition intensity between immune responses and the viral stimulation delay, respectively. The diffusion rates of viruses and immune cells are $ d_1 $ and $ d_2 $, respectively. The same parameters as in model (1.1) have the same biological meanings. Note that all parameters in (1.3) remain nonnegative.

    The remaining sections of the paper are organized as follows: Section 2 discusses the local asymptotic stability of the positive constant steady state solution and the existence of Hopf bifurcation. In Section 3, we derive the normal form of the Hopf bifurcation by using the multiple time scales (MTS) method. Section 4 presents numerical simulations to illustrate the theoretical analysis, and the conclusion is given in Section 5.

    In this section, we concentrate on studying the local asymptotic stability of the positive constant steady state solution and the existence of Hopf bifurcation for system (1.3). From a biological perspective, we only focus on the positive constant steady state solution $ (y_*, z_*) $. It is evident that $ (y_*, z_*) $ satisfies the following equations:

    $ ry(1yK)aypyz=0,cyz1+ηy+szbzqyz=0.
    $
    (2.1)

    Then we can write $ y_* $ in terms of $ z_* $ as

    $ y=K(1a+pzr),
    $
    (2.2)

    substituting Eq (2.2) into Eq (2.1), we obtain

    $ Az2+Bz+C=0,
    $

    where

    $ A=qηK2p2+Kpqrs,B=2qη(ra)pK2+[r(bηc+q)p(ra)rqs]Kbr2s,C=qη(ra)2K2+r(bηc+q)(ar)Kbr2.
    $
    (2.3)

    Denote

    $ Δ1=B24AC.
    $
    (2.4)

    Lemma 1. For the positive constant steady state solution of system (1.3), we have the following results.

    (i) If $ \Delta_1 > 0 $, $ AC > 0, \; AB < 0, \mathit{\text{and}} a+pz_* < r $ holds, then system (1.3) has two positive constant steady state solutions, $ E_{*_{1}} = (y_{*_{1}}, z_{*_{1}}) $ and $ E_{*_{2}} = (y_{*_{2}}, z_{*_{2}}) $. Here $ z_{*_{1}} = \frac{-B + \sqrt{\Delta_1}}{2A} $, $ z_{*_{2}} = \frac{-B - \sqrt{\Delta_1}}{2A} $, with $ y_{*_{1}} = K\left(1 - \frac{a + pz_{*_{1}}}{r}\right) $ and $ y_{*_{2}} = K\left(1 - \frac{a + pz_{*_{2}}}{r}\right) $;

    (ii) If $ \Delta_1 > 0 $, $ AC < 0, \mathit{\text{and}} a+pz_* < r $ holds, then system (1.3) has one positive constant steady state solution $ E_{*_{3}} = (y_{*_{3}}, z_{*_{3}}) $. Here $ z_{*_{3}} = \frac{-B + \sqrt{\Delta_1}}{2A} $, with $ y_{*_{3}} = K\left(1 - \frac{a + pz_{*_{3}}}{r}\right) $;

    (iii) If $ \Delta_1 = 0 $, $ AB < 0, \mathit{\text{and}} a+pz_* < r $ holds, then system (1.3) has one positive constant steady state solution $ E_{*_{4}} = (y_{*_{4}}, z_{*_{4}}) $. Here $ z_{*_{4}} = \frac{-B}{2A} $, with $ y_{*_{4}} = K\left(1 - \frac{a + pz_{*_{4}}}{r}\right) $.

    Without loss of generality, we denote $ E_{*_{1}} $, $ E_{*_{2}} $, $ E_{*_{3}} $, and $ E_{*_{4}} $ as $ E_* $. Therefore, when the parameters of system (1.3) satisfy any one of the conditions (ⅰ)–(ⅱ) of Lemma 1, the positive constant steady state solution $ E_* = (y_*, z_*) $ exists.

    Now, we analyze the stability of $ E_* = (y_*, z_*) $. Set $ U(x, t) = (y(x, t), z(x, t))^{\mathrm{T}} $, then the linearization of system (1.3) at $ E_* $ can be written as

    $ U(x,t)t=DΔU(x,t)+AU(x,t)+BU(x,tτ),
    $
    (2.5)

    where

    $ D=(d100d2),A=(ryKpyqzzs(b+qy)2cy),B=(00z(b+qy)[cη(b+qy)]cy0).
    $

    Therefore, the characteristic equation of the linearization of system (1.3) at $ E_* $ can be obtained as

    $ λ2+Tnλ+Nn+C0(eλτb0q)=0,nN0={0,1,2,3,...},
    $
    (2.6)

    with

    $ Tn=n2l2d1+n2l2d2+ryK+zs(b+qy)2cy,Nn=n4l4d1d2+n2l2d1zs(b+qy)2cy+n2l2d2ryK+rzs(b+qy)2Kc,C0=pyz,b0=(b+qy)[cη(b+qy)]cy.
    $

    Lemma 2. When $ \tau = 0 $, the positive constant steady state solution $ E_* $ of system (1.3) is locally asymptotically stable if $ b_0-q > 0 $.

    Proof. When $ \tau = 0 $, characteristic Eq (2.6) becomes

    $ λ2+Tnλ+Nn+C0(b0q)=0,nN0.
    $
    (2.7)

    It is clear that $ T_n > 0 $, $ N_n > 0 $, thus if $ b_0-q > 0 $ holds, the real parts of the roots of Eq (2.6) are negative. Thus, $ E_* $ is locally asymptotically stable.

    This completes the proof.

    Suppose that $ \lambda = $i$ \omega\; (\omega > 0) $ is a root of Eq (2.6), then we have

    $ {C0b0cos(ωτ)=ω2+C0qNn,C0b0sin(ωτ)=Tnω.
    $
    (2.8)

    Square and add above equations, then denote $ g\mathrel{\overset{\triangle}{ = }}\omega^2. $ We obtain

    $ φ(g)=g2+(2C0q2Nn+T2n)g+(C0qNn)2C20b20=0.
    $
    (2.9)

    Next, we will discuss the existence of $ \omega $. For the sake of convenience in our discussion, we present the following assumed conditions:

    $ (H1)Δ2=0,2C0q2Nn+T2n<0,(H2)Δ2>0,(C0qNn)2C20b20<0,(H3)Δ2>0,2C0q2Nn+T2n<0,(C0qNn)2C20b20>0,(H4)Δ2<0,(H5)Δ2>0,2C0q2Nn+T2n>0,(C0qNn)2C20b20>0,
    $

    where $ \Delta_2 = (2C_0 q - 2 N_n+T_n^2)^2-4\left[(C_0 q-N_n)^2-C_0^2 b_0^2\right] $, and $ N_n, T_n, b_0, C_0 $ are given in Eq (2.6).

    Lemma 3. When $ \tau > 0 $, concerning the cases of positive roots of Eq (2.9), we have the following results:

    (i) Equation (2.9) has a unique positive root if $ \mathrm{(H1)} $ or $ \mathrm{(H2)} $ is satisfied, then the positive constant steady state solution of system (1.3) will not exhibit stability switches,

    (ii) Equation (2.9) has two positive roots if $ \mathrm{(H3)} $ is satisfied, then the positive constant steady state solution of system (1.3) will exhibit stability switches,

    (iii) Equation (2.9) has no positive root if $ \mathrm{(H4)} $ or $ \mathrm{(H5)} $ is satisfied.

    Proof. If the condition $ \mathit{\text{(H1)}} $ is satisfied, we have $ g_1 = g_2 $, $ g_1+g_2 = -(2C_0 q - 2 N_n+T_n^2) > 0 $. If the condition $ \mathit{\text{(H2)}} $ is satisfied, we have $ g_1\cdot g_2 = (C_0 q-N_n)^2-C_0^2 b_0^2 < 0 $. Both of them imply that there is a unique positive root $ g $ satisfying Eq (2.9). This proves (ⅰ).

    Similarly, we conclude (ⅱ)–(ⅲ).

    Considering the scenario with only a pair of purely imaginary roots, the unique positive $ \omega $ that satisfies Eq (2.9) can be described as follows:

    $ ω+n=(2C0q2Nn+T2n)+Δ22,nN0.
    $
    (2.10)

    Denote

    $ Rn=sin(ω+nτ)=Tnω+nC0b0,Qn=cos(ω+nτ)=(ω+n)2+C0qNnC0b0.
    $
    (2.11)

    Now, define a set

    $ I1={nN0|such that (H1) or (H2) holds}.
    $
    (2.12)

    From Eq (2.11), we obtain

    $ τ+n,j={1ω+n(arccos(Qn)+2jπ),forRn0,1ω+n(2πarccos(Qn)+2jπ),forRn<0.nI1,jN0.
    $
    (2.13)

    From Eq (2.13), we find that $ \{\tau_{n, j}^+\}_{j = 0}^{\infty} $ is increasing about $ j $ for some fixed $ n \in I_1 $. Thus, for fixed $ n $, one has $ \tau_{n, 0}^+ = \min _{j \in \mathbb{N}_0}\left\{\tau_{n, j}^ + \right\} $. To investigate the stability of $ E_* $, we define the smallest critical value of time delay

    $ τc=min{τ+n,0|nI1}.
    $
    (2.14)

    Theorem 1. If $ b_0-q > 0 $ and one of the two conditions $ {{\text{(H1)}}} $ or $ {{\text{(H2)}}} $ is satisfied, then the positive constant steady state solution $ E_* $ of system (1.3) is locally asymptotically stable when $ \tau \in [0, \tau_c) $, and unstable when $ \tau \in \left(\tau_c, \infty\right) $. Moreover, system (1.3) undergoes the Hopf bifurcation at $ E_* $ when $ \tau = \tau_c $, where $ \tau_c $ is given in Eq (2.14).

    Proof. Differentiating the both sides of the characteristic Eq (2.6) with respect to time delay $ \tau $ gives rise to

    $ (dλdτ)1=2λ+Tnλ[C0q(λ2+Tnλ+Nn)]τλ.
    $

    As a result, one has

    $ (Redλdτ)1|τ=τ+n,j=Δ2T2n(ω+n,j)2+[NnC0q(ω+n,j)2]2>0,nI1,jN0.
    $
    (2.15)

    This means that the transversality condition is satisfied.

    This completes the proof.

    We have obtained that under certain conditions, the characteristic Eq (2.6) of system (1.3) has a unique pair of purely imaginary roots $ \lambda = \pm $i$ \omega_n^+ $, then system (1.3) undergoes the Hopf bifurcation at $ E_* $ when $ \tau = \tau_c $, with $ \tau_c $ as given at $ (2.14) $. In this section, we will utilize the method of MTS to derive the normal form of the Hopf bifurcation.

    Denote that

    $ f(y(x,tτ),z(x,t))=cy(x,tτ)z(x,t)1+ηy(x,tτ)+sz(x,t)
    $
    (3.1)

    and

    $ L=1+ηy+sz,α=1+sz,β=1+ηy.
    $
    (3.2)

    Perform a Taylor expansion of $ f(y(x, t-\tau), z(x, t)) $ at $ E_* $, then we obtain

    $ f(y(x,tτ),z(x,t))=cyzL+czαL2(y(x,tτ)y)+cyβL2(z(x,t)z)cηzαL3(y(x,tτ)y)2+c(2ηyzs+ηy+α)L3(y(x,tτ)y)(z(x,t)z)csyβL3(z(x,t)z)2+cη2zαL4(y(x,tτ)y)3+cη[3szαL41+2szL3](y(x,tτ)y)2(z(x,t)z)+cs[1zs(1+2ηy)+η2y2]L4(y(x,tτ)y)(z(x,t)z)2+cs2yβL4(z(x,t)z)3+.
    $
    (3.3)

    Now, we set $ \hat{y}(x, t) = y(x, \tau t)-y_* $, $ \hat{z}(x, t) = z(x, \tau t)-z_* $, and still denote $ \hat{y}(x, t) $ and $ \hat{z}(x, t) $ by $ y(x, t) $ and $ z(x, t) $, respectively. Then system (1.3) takes the form

    $ {y(x,t)t=τ[d1Δy(x,t)+r(y(x,t)+y)(1y(x,t)+yK)a(y(x,t)+y)p(y(x,t)+y)(z(x,t)+z)],z(x,t)t=τ[d2Δz(x,t)+czαL2y(x,t1)+cyβL2z(x,t)cηzαL3y2(x,t1)+c(2ηyzs+ηy+α)L3y(x,t1)z(x,t)csyβL3z2(x,t)+cη2zαL4y3(x,t1)+cη[3szαL41+2szL3]y2(x,t1)z(x,t)+cs[1zs(1+2ηy)+η2y2]L4y(x,t1)z2(x,t)+cs2yβL4z3(x,t)bz(x,t)q[y(x,t)z(x,t)+y(x,t)z+z(x,t)y]],νy(x,t)=νz(x,t)=0,xΩ,t>0,y(x,t)=ψ1(x,t),z(x,t)=ψ2(x,t),xΩ=(0,lπ),t[τ,0].
    $
    (3.4)

    Let $ h = \left(h_{11}, h_{12}\right)^{\mathrm{T}} $ be the eigenvector of the characteristic matrix of the linearization of system Eq (3.4) at $ E_* $ corresponding to the eigenvalue $ \mathrm{i} \omega \tau $ and $ h^* = d\left(h_{21}, h_{22}\right)^{\mathrm{T}} $ be the eigenvector of the adjoint matrix corresponding to the eigenvalue $ -\mathrm{i} \omega \tau $. They satisfy the inner product

    $ \langle h^*, h\rangle = \bar{h}^*\mathrm{^T}\cdot h = 1. \quad $

    By calculating, we have

    $ h11=1,h12=1py(iω+n2l2d1+rKy),h21=1py[iω+n2l2d2+zs(b+qy)2cy],h22=1,d=(h21ˉh11+h22ˉh12)1.
    $
    (3.5)

    Consider the delay $ \tau $ as a bifurcation parameter, let $ \tau = \tau_c+\varepsilon \mu $, where $ \tau_c $ is the critical point of the Hopf bifurcation, given as Eq (2.14). $ \mu $ is a small perturbation parameter, and $ \varepsilon $ is the characteristic scale used for non-dimensionalization. The solution of Eq (3.4) is assumed to take the form:

    $ U(x,t)=U(x,T0,T1,T2,)=+k=1εkUk(x,T0,T1,T2,),
    $
    (3.6)

    where

    $ U(x,T0,T1,T2,)=(y(x,T0,T1,T2,),z(x,T0,T1,T2,))T,Uk(x,T0,T1,T2,)=(yk(x,T0,T1,T2,),zk(x,T0,T1,T2,))T.
    $
    (3.7)

    The derivative with respect to $ t $ is now transformed into

    $ t=T0+εT1+ε2T2+=D0+εD1+ε2D2+,
    $
    (3.8)

    with the differential operator $ D_i = \frac{\partial}{\partial T_i}, \; i \in \mathbb{N}_0. $ Let

    $ yj=yj(x,T0,T1,T2,),yj,1=yj(x,T01,T1,T2,),zj=zj(x,T0,T1,T2,),zj,1=zj(x,T01,T1,T2,),
    $
    (3.9)

    where $ j\in \mathbb{N} = \{1, 2, 3, ...\} $. According to Eqs (3.7) and (3.8), we get

    $ Ut=εD0U1+ε2D1U1+ε3D2U1+ε2D0U2+ε3D1U2+ε3D0U3+.
    $
    (3.10)

    To deal with the delayed terms, we expand $ y(x, t-1) $ at $ y(x, T_0-1, T_1, T_2, \cdots) $, then we have

    $ y(x,t1)=εy11+ε2y21+ε3y31ε2D1y11ε3D1y21ε3D2y11+,
    $
    (3.11)

    where $ y_{j1} = y_j\left(x, T_0-1, T_1, T_2, \cdots\right), \; j\in \mathbb{N} $.

    Substituting Eqs (3.6)–(3.11) into Eq (3.4) will lead to some perturbation equations. By collecting the coefficients with respect to different powers of $ \varepsilon $, we can obtain a set of ordered linear differential equations. First, for the $ \varepsilon $-order terms, we have

    $ {D0y1τc[ry1ay1+d1Δy12ryy1Kpy1zpyz1]=0,D0z1τc[d2Δz1+czαy11L2+cyβz1L2bz1q(yz1+y1z)]=0,
    $
    (3.12)

    where $ L, \alpha, \beta $ are given by Eq (3.2). Thus the solution of Eq (3.12) can be expressed in the form as follows:

    $ y1=G(T1,T2,)eiωτcT0h11cos(nx)+c.c.,z1=G(T1,T2,)eiωτcT0h12cos(nx)+c.c.,
    $
    (3.13)

    where $ h_{11}, h_{12} $ are given by Eq (3.5), and $ c. c. $ represents the complex conjugate of the previous term. Next, for the $ \varepsilon^2 $-order terms, we obtain

    $ {D0y2τc[ry2ay2+d1Δy22ryy2Kpy2zpyz2]=D1y1τcry21Kμay1+μry1+μd1Δy1μpy1zμpyz12μryy1Kτcpy1z1,D0z2τc[d2Δz2+czαy21L2+cyβz2L2bz2q(yz2+y2z)]=τc[czαD1y11L2cηzαy211L3+c(2ηyzs+ηy+α)y11z1L3cysβz21L3qy1z1]+μ[d2Δz1+czαy11L2+cyβz1L2bz1q(yz1+y1z)]D1z1,
    $
    (3.14)

    where $ L, \alpha, \beta $ are given by Eq (3.2). We substitute solution Eq (3.13) into the right side of Eq (3.14), and denote the coefficient vector of $ \mathrm{e}^{\mathrm{i} \omega \tau_c T_0} $ as $ m_1 $, which satisfies the solvability condition

    $ \langle h^*, (m_1,\cos(nx))\rangle = 0. $

    Then one has

    $ GT1=MμG,
    $
    (3.15)

    where

    $ M=ˉh21M1+ˉh22M2h11ˉh21+h12ˉh22+τczcαL2eiωτch11ˉh22,
    $
    (3.16)

    with

    $ M1=ah11+rh11n2d1h11pzh11pyh122ryKh11,M2=n2d2h12+czαL2eiωτch11+cyβL2h12bh12qyh12qzh11.
    $

    Suppose that the solution of Eq (3.14) is

    $ {y2=+k=0(η0,kGˉG+η1,ke2iωτcT0G2+ˉη1,ke2iωτcT0ˉG2)cos(kx),z2=+k=0(ζ0,kGˉG+ζ1,ke2iωτcT0G2+ˉζ1,ke2iωτcT0ˉG2)cos(kx).
    $
    (3.17)

    Let $ C_k = \langle\cos (n x) \cos (n x), \cos (k x)\rangle = \int_0^{l\pi} \cos (n x) \cos (n x) \cos (k x) \mathrm{d} x $. Besides, we have

    $ lπ0cos(nx)cos(kx)dx={0kn,lπ2k=n0,lπk=n=0.
    $

    We substitute solution Eqs (3.13) and (3.17) into Eq (3.14), then make the inner product with $ \cos (kx) $ on both sides of the obtained equation, where $ k\in \mathbb{N}_0 $. By adding the obtained results, we can compare the coefficients in front of terms $ G^2 $ and $ G \bar{G} $, then we get

    $ η1,k=Ck(W(J1+2iω)J2py)lπ0cos(kx)cos(kx)dx[(2iω+J1)(l1+2iω)pyczαL2e2iωτcqz],η0,k=Ck(SJ1l2py)lπ0cos(kx)cos(kx)dx[J1l1pyczαL2qz],ζ1,k=Cklπ0cos(kx)cos(kx)dxJ2pyl1+2iωpyη1,k,ζ0,k=Cklπ0cos(kx)cos(kx)dxl2pyl1pyη0,k,
    $
    (3.18)

    where

    $ W=cηzαL3e2iωτch211+c(2ηyzs+ηy+α)L3eiωτccysβL3h212qh11h12,S=cηzαL32h11ˉh11+c(2ηyzs+ηy+α)L3(eiωτch11ˉh12+eiωτcˉh11h12)cysβL32h12ˉh12q(h11ˉh12+h12ˉh11),J1=d2k2cyβL2+b+qy,J2=rKh211ph11h12,l1=ar+2ryK+pz+d1k2,l2=rK2h11ˉh11p(h11ˉh12+h12ˉh11),L=1+ηy+sz,α=1+sz,β=1+ηy.
    $

    Specifically, if $ n = 0 $, the condition for $ C_k \neq0 $ is $ k = 0 $. Thus we consider $ n = k = 0 $, then we have $ C_k = l \pi $, $ \int_0^{l\pi} \cos (n x) \cos (k x) \mathrm{d} x = l \pi $. Thus we can simply get the value of $ \eta_{1, 0} $, $ \eta_{0, 0} $, $ \zeta_{1, 0} $, $ \zeta_{0, 0} $. For the $ \varepsilon^3 $-order terms, we get

    $ {D0y3+τc(ar+2ryK+pz)y3d1τcΔy3+τcpyz3=μd1Δy2μ(ar+2ryK+pz)y2D2y1D1y2μpy1z1μpyz2μry21K2τcry2y1Kτcpy1z2τcpy2z1,D0z3τcd2Δz3τcczαL2y31(τccyβL2bqy)z3=D2z1D1z2+τc[czαD1y21czαD2y11L22cηzαD1y211+2cηzαy11y21L3+c(2ηyzs+ηy+α)y11z2L3+c(2ηyzs+ηy+α)(D1y11z1+y21z1)L32cysβz1z2L3+c(1+z(2ηy1)s)sy11z21L4+cys2βz31L4+czη2αy311L4+cη[3szαL42sz+1L3]y211z1qy1z2qy2z1]+μ[d2Δz2+czαD1y11+czαy21L2+(cyβL2bqy)z2cηzαy211L3+c(2ηyzs+ηy+α)y11z1L3cyβsz21L3qy1z1qy2z],
    $
    (3.19)

    where $ L, \alpha, \beta $ are given by Eq (3.2). We substitute solution Eqs (3.13) and (3.17) into Eq (3.19), then denote the coefficient vector of $ \mathrm{e}^{\mathrm{i} \omega \tau_c T_0} $ as $ m_2 $, which satisfies the solvability condition

    $ \langle h^*, (m_2,\cos(nx))\rangle = 0. $

    Since $ \mu $ is the disturbance parameter and $ \mu^2 $ has very little influence on the result, the items containing $ \mu^2 G $ can be ignored. Then we obtain

    $ \frac{\partial G}{\partial T_2} = \chi G^2 \bar{G}, $

    where

    $ χ=RCk+QV,
    $
    (3.20)

    with

    $ R=ˉh21[2τcrK(h11+k=0η0,k+ˉh11+k=0η1,k)τcp(h11+k=0ζ0,k+ˉh11+k=0ζ1,k+h12+k=0η0,k+ˉh12+k=0η1,k)]+ˉh22τc[2cηzαL3(eiωτch11+k=0η0,k+eiωτcˉh11+k=0η1,k)+c(2ηyzs+ηy+α)L3(eiωτch11+k=0ζ0,k+eiωτcˉh11+k=0ζ1,k+h12+k=0η0,k+ˉh12+k=0η1,k)2cysβL3(h12+k=0ζ0,k+ˉh12+k=0ζ1,k)],Q=ˉh22τc[c(1+z(2ηy1)s)sL4(2eiωτch11h12ˉh12+eiωτcˉh11h212)+cyβs2L4(2h212ˉh12)+czη2αL4(2eiωτch211ˉh11)+cη[3szαL42sz+1L3](e2iωτch211ˉh12+2h11ˉh11h12)q(h11+k=0ζ0,k+ˉh11+k=0ζ1,k+h12+k=0η0,k+ˉh12+k=0η1,k)]lπ0cos4(nx)dx,V=[h11ˉh21+h12ˉh22+τcczαL2eiωτch11ˉh22]lπ0cos(nx)cos(nx)dx,L=1+ηy+sz,α=1+sz,β=1+ηy.
    $

    Therefore, the normal form of the Hopf bifurcation for system Eq (1.3) reduced on the center manifold is

    $ ˙G=MμG+χG2ˉG,
    $
    (3.21)

    where $ M $ and $ \chi $ are given by Eqs (3.16) and (3.20), respectively. Now introduce $ G = r \mathrm{e}^{\mathrm{i} \theta } $ and put it into the normal form Eq (3.21), then one has

    $ {˙r=Re(M)μr+Re(χ)r3,˙θ=Im(M)μ+Im(χ)r2.
    $
    (3.22)

    For the stability of the periodic solutions of system (1.3), we have the following result.

    Theorem 2. When $ \frac{\mathrm{Re}(M)\mu}{\mathrm{Re}(\chi)} < 0 $, system Eq (1.3) exists as periodic solutions near the positive constant steady state solution $ E_* $. Moreover,

    (i) If $ \mathrm{Re}(M)\mu < 0 $, the bifurcating periodic solutions reduced on the center manifold are unstable, and when $ \mu > 0 $ ($ \mu < 0 $), the direction of bifurcation is forward (backward);

    (ii) If $ \mathrm{Re}(M)\mu > 0 $, the bifurcating periodic solutions reduced on the center manifold are stable, and when $ \mu > 0 $ ($ \mu < 0 $), the direction of bifurcation is forward (backward).

    In this section, we will perform numerical simulations for different values of $ \tau $. Based on the work of Komarova et al. [1], we fix

    $ a=3.0day1,p=1.0mm3cells1day1,q=1.0mm3virus1day1,r=6.0day1.
    $

    For the remaining parameters, we consider the actual meaning of the parameters and choose

    $ d1=0.4mm3day1,K=10virusmm3,b=2day1,d2=0.1mm3day1,c=6.5mm3virus1day1,η=0.35mm1virus,s=0.02mm3immunecells1,l=2.
    $

    It is easy to check that the condition (ⅱ) of Lemma 1 is satisfied, which implies that the system (1.3) has only one positive constant steady state solution,

    $ E_* = (y_*,z_*) = (0.4600,2.7239). $

    In addition, we can also find that this set of parameters satisfies the condition $ \text{(H2)} $ of Lemma 3. Furthermore, by Eqs (2.10), (2.12) and (2.14), we obtain

    $ ω+0=2.1289,I1={n=0,1,2,3},τc=τ+0,0=0.0667.
    $

    Therefore, by Theorem 2.1, $ E_* $ is locally asymptotically stable when $ \tau \in [0, \tau_c) = [0, 0.0667) $, and unstable when $ \tau \in \left(0.0667, +\infty\right) $. Take $ \tau = 0.055 \in [0, 0.0667) $ as an example, then $ E_* $ is locally asymptotically stable as is illustrated by the computer simulations (see Figure 1).

    Figure 1.  Choose the initial function as ($ 0.4600+0.0001\cos x, 2.7239+0.002\cos x $), then the constant steady state solution $ E_* = (0.4600, 2.7239) $ of system (1.3) is locally asymptotically stable when $ \tau = 0.055 < \tau_c $.

    Biologically, when the virus stimulation delay does not exceed the threshold, immune cells can control the virus in the long-term at low levels. Although the virus is not completely eliminated, it does not cause serious disease deterioration, this is the ideal scenario that we hope for.

    According to Eqs (3.16) and (3.20), we can get

    $ M = 0.3874 -10.6020\mathrm{i},\; \; \; \; \chi = -0.4221 - 0.2652\mathrm{i}, $

    which implies that $ \mathrm{Re}(M) > 0 $, and $ \mathrm{Re}(\chi) < 0 $. From Theorem 2, when $ \mu > 0 $, system (1.3) has stable bifurcating periodic solutions. In addition, the direction of the Hopf bifurcation is forward. Take $ \tau = 0.067\in \left(0.0667, +\infty\right) $ as an example, then system (1.3) generates a spatially homogeneous stable periodic solution (see Figure 2, where Figure 2(a), (c) depict simulations of the virus and immune cell dynamics, respectively, while Figure 2(b), (d) represent their respective local zoom-in plots).

    Figure 2.  Choose the initial function as ($ 0.4600+0.0001\cos x $, $ 2.7239+0.002\cos x $), then system (1.3) has a spatially homogeneous stable periodic solution when $ \tau = 0.067 > \tau_c $.

    Biologically, when the virus stimulation delay slightly exceeds the threshold $ \tau_c $, the densities of immune cells and virus exhibits periodic oscillations. As the number of immune cells increases, they effectively control the virus and reduce the viral load. After a while, an immunoevasion mechanism for the virus may emerge and cause the number of immune cells to drop, allowing the virus to remultiply. This scenario is not ideal, but it is still manageable.

    Once $ \tau $ significantly exceeds the threshold, take $ \tau = 0.076\in \left(0.0667, +\infty\right) $ as an example, the virus proliferates rapidly (see Figure 3). This is the worst-case scenario, as the virus could spread rapidly and pose significant risks to the host.

    Figure 3.  Choose the initial function as ($ 0.4600+0.0001\cos x $, $ 2.7239+0.002\cos x $), then the constant steady state solution $ E_* $ of system (1.3) is unstable when $ \tau = 0.076 > \tau_c $.

    In this paper, we incorporated diffusion into an immunosuppressive infection model with virus stimulation delay and Beddington-DeAngelis functional response. We studied the stability of the constant steady state solution of the model and provided conditions for the occurrence of Hopf bifurcation. Moreover, we analyzed the stability and direction of the Hopf bifurcation, and further derived the normal form. By analyzing the model, we obtained a threshold $ \tau_c $ for the virus stimulation delay $ \tau $ and discovered that when $ \tau $ was less than this threshold, the constant steady state solution was locally asymptotically stable. However, once $ \tau $ exceeded the threshold, the constant steady state solution lost stability, and the model exhibited stable periodic solutions. This implied that when the virus could quickly stimulate immune cells, the host possessed strong immunity, which could inhibit virus growth and spread. Conversely, rapid virus proliferation posed significant risks to the host. The discovery of the threshold could help biologists predict, judge, and control the growth of virus cells and select effective treatment strategies.

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

    This study was funded by Fundamental Research Funds for the Central Universities of China (Grant No.2572022DJ07).

    The authors declare there is no conflicts of interest.

    [1] Pearlin LI, Schooler C. (1978) The structure of coping. J Health Soc Behav 19: 2-21. doi: 10.2307/2136319
    [2] Pearlin LI, Pioli MF. (2003) Personal control: Some conceptual turf and future directions. New York: Springer.
    [3] Skinner EA. (1996) A guide to constructs of control. J Person Soc Psychol 71: 549. doi: 10.1037/0022-3514.71.3.549
    [4] Roepke SK, Grant I. (2011) Toward a more complete understanding of the effects of personal mastery on cardiometabolic health. Health Psychol 30: 615. doi: 10.1037/a0023480
    [5] Skaff MM. (2007) Sense of control and health. Handbook of health psychology and aging 186-209.
    [6] Lightsey OR. (1997) Stress buffers and dysphoria: A prospective study. J Cogn Psychother 11:263-277.
    [7] Miller LM, West RL. (2010) The effects of age, control beliefs, and feedback on self-regulation of reading and problem solving. Exp Aging Res 36: 40-63.
    [8] Bisschop MI, Kriegsman DM, Beekman AD, et al. (2004) Chronic diseases and depression: the modifying role of psychosocial resources. Social Sci Med 59: 721-733. doi: 10.1016/j.socscimed.2003.11.038
    [9] Ward MM. (2013) Sense of control and self-reported health in a population-based sample of older Americans: assessment of potential confounding by affect, personality, and social support. Int J Behav Med 20: 140-147. doi: 10.1007/s12529-011-9218-x
    [10] Bowen R, South M, Fischer D, et al. (1994) Depression, mastery and number of group sessions attended predict outcome of patients with panic and agoraphobia in a behavioural/medication program. Can J Psychiatry 39: 283-288.
    [11] Lachman ME. (1986) Locus of control in aging research: a case for multidimensional and domain-specific assessment. Psychol Aging 1: 34. doi: 10.1037/0882-7974.1.1.34
    [12] Reker GT, Peacock EJ, Wong PT. (1987) Meaning and purpose in life and well-being: A life-span perspective. J Gerontol 42: 44-49. doi: 10.1093/geronj/42.1.44
    [13] Lachman ME, Firth KM. (2004) The adaptive value of feeling in control during midlife. In: How healthy are we 320-349.
    [14] Wolinsky FD, Wyrwich KW, Babu AN, et al. (2003) Age, aging, and the sense of control among older adults: A longitudinal reconsideration. J Gerontol Series B: Psychol Sci Social Sci 58: S212-S220. doi: 10.1093/geronb/58.4.S212
    [15] Pearlin LI, Schieman S, Fazio EM, et al. (2005) Stress, health, and the life course: Some conceptual perspectives. J Health Soc Behav 46: 205-219. doi: 10.1177/002214650504600206
    [16] Dong X, Beck T, Simon MA. (2009) Loneliness and mistreatment of older Chinese women: does social support matter? J Women Aging 21: 293-302. doi: 10.1080/08952840903285252
    [17] Baltes PB, Baltes MM, Baltes PB, et al. (1990) Psychological perspectives on successful aging: The model of selective optimization with compensation. Success Aging Persp Behave Sci 1: 1-34.
    [18] US Census Bureau. (2010) American Community Survey. Washington: US Census Bureau.
    [19] Liu EY, Mencken FC. (2010) Fatalistic voluntarism and life happiness in post-socialist China. Sociol Spect 30: 270-288. doi: 10.1080/02732171003635455
    [20] Terrill R. (1979) China enters the 1980s. Foreign Aff 58: 920.
    [21] Ho DY. (1994) Filial piety, authoritarian moralism, and cognitive conservatism in Chinese societies. Gene Soc Gen Psychol Monogr 120: 349-365.
    [22] Dong X, Wong E, Simon MA. (2014) Study design and implementation of the PINE Study. J Health Aging 0898264314526620.
    [23] Dong X, Chang ES, Wong E, et al. (2011) Working with culture: lessons learned from a community-engaged project in a Chinese aging population. Aging Health 7: 529-537. doi: 10.2217/ahe.11.43
    [24] Simon MA, Chang E, Rajan K, et al. (2014) Demographic characteristics of U. S. Chinese older adults in the greater Chicago area: Assessing the representativeness of the PINE study. J Aging Health. In Press.
    [25] Forbes DA. (2001) Enhancing mastery and sense of coherence: Important determinants of health in older adults. Geriatr Nurs 22: 29-32. doi: 10.1067/mgn.2001.113532
    [26] Penninx B, Tilburg T, Kriegsman M, et al. (1997) Effects of social support and personal coping resources on mortality in older age: the longitudinal aging study amsterdam. Am J Epidemiol 146:510-519. doi: 10.1093/oxfordjournals.aje.a009305
    [27] Skaff MM, Pearlin LI, Mullan JT. (1996) Transitions in the caregiving career: effects on sense of mastery. Psychol Aging 11: 247-257. doi: 10.1037/0882-7974.11.2.247
    [28] Dong X, Simon MA, Wilson RS, et al. (2010) Decline in cognitive function and risk of elder self-neglect: finding from the Chicago health aging project. J Am Geriatr Society 58: 2292-2299. doi: 10.1111/j.1532-5415.2010.03156.x
    [29] Dong X, Simon M, de Leon CM, et al. (2009) Elder self-neglect and abuse and mortality risk in a community-dwelling population. JAMA: J Am Med Assoc 302: 517-526. doi: 10.1001/jama.2009.1109
    [30] Brenner S. (1995) Why women rule the roost: rethinking Javanese ideologies of gender and self-control. In: Bewitching women pious men: Gender and body politics in Southeast Asia 19-50.
    [31] Dong X, Beck T, Simon MA. (2010) The associations of gender, depression and elder mistreatment in a community-dwelling Chinese population: the modifying effect of social support. Arch Gerontol Geriatr 50: 202-208. doi: 10.1016/j.archger.2009.03.011
    [32] Pearlin L, Nguyen K, Schieman S, et al. (2007) The life-course origins of mastery among older people. J Health Soc Behav 48: 164-179. doi: 10.1177/002214650704800205
    [33] Marshall GN, Lang EL. (1990) Optimism, self-mastery, and symptoms of depression in women professionals. J Per Soc Psychol 59: 132-139. doi: 10.1037/0022-3514.59.1.132
    [34] Jonker AGC. (2010) Health decline and well-being in old age: the need of coping. Tijdschrift voor Gerontologie en Geriatrie 42: 102-104.
    [35] Marcus P. (2008) Victory through vegetables: self-mastery through a vegetarian way of life. Psychoanal. Rev 95: 61-77. doi: 10.1521/prev.2008.95.1.61
    [36] Paquet C, Dub L, Gauvin L, et al. (2010) Sense of mastery and metabolic risk: moderating role of the local fast-food environment. Psychosom Med 72: 324-331. doi: 10.1097/PSY.0b013e3181cdf439
    [37] Pudrovska T, Schieman S, Pearlin L, et al. (2005) The sense of mastery as a mediator and moderator in the association between economic hardship and health in late life. J Aging Health17: 634-660.
    [38] Jonker A, Comijs C, Knipscheer K, et al. (2009) The role of coping resources on change in well-being during persistent health decline. J Aging Health 21: 1063-1082. doi: 10.1177/0898264309344682
    [39] Jang Y, Haley WE, Small B, et al. (2002) The role of mastery and social resources in the associations between disability and depression in later life. The Gerontologist 42: 807-813. doi: 10.1093/geront/42.6.807
    [40] Dong X, Simon MA. "Is greater social support a protective factor against elder mistreatment?. " Gerontology 54: 381-388.
    [41] Dong X, Simon MA, Gorbien M. (2007) Elder abuse and neglect in an urban Chinese population. J Elder Abuse Neglect 19(3-4): 79-96.
    [42] Dong X, Simon MA, Odwazny R, et al. (2008) Depression and elder abuse and neglect among a community-dwelling Chinese elderly population. J Elder Abuse Negl 20(1): 25-41.
    [43] Dong X, de Leon CFM, Evans DA. (2009) Is greater self-neglect severity associated with lower levels of physical function? J Aging Health 21(4):596-610.
    [44] Dong X, Simon MA, Gorbien M, et al. (2007) Loneliness in older Chinese adults: a risk factor for elder mistreatment. J Am Geriatr Society 55: 1831-1835. doi: 10.1111/j.1532-5415.2007.01429.x
    [45] Dong X, Chang ES, Simon M, et al. (2011) Sustaining Community-University Partnerships: Lessons learned from a participatory research project with elderly Chinese. Gateways: Int J Comm Res Engag 4: 31-47. doi: 10.5130/ijcre.v4i0.1767
    [46] Dong X, Li Y, Chen R, et al. (2013) Evaluation of community health education workshops among Chinese older adults in Chicago: a community-based participatory research approach. J Educ Train Stud 1(1): 170-181.
  • This article has been cited by:

    1. K. Ishwarya, K. G. Rani, K. Appathurai, R. Surendran, R. Selvanarayanan, C. Y. Lau, 2024, 3161, 0094-243X, 020301, 10.1063/5.0229271
    2. Ishwarya K, Saranya G, Appathurai K, Surendran R, 2023, Exploring the Significance and Paradoxical Nature of the Zero Value in Mathematics Using Artificial Intelligence, 979-8-3503-0088-8, 861, 10.1109/ICOSEC58147.2023.10276215
    3. Shiying Wang, Hongmei Li, Lina Zhao, The m-Component Connectivity of Leaf-Sort Graphs, 2024, 12, 2227-7390, 404, 10.3390/math12030404
    4. Krittawit Limkul, Sayan Panma, On the Independence Number of Cayley Digraphs of Clifford Semigroups, 2023, 11, 2227-7390, 3445, 10.3390/math11163445
    5. S. T. Vikram, S. Balaji, An Analysis of the Factors Influencing the Strong Chromatic Index of Graphs Derived by Inflating a Few Common Classes of Graphs, 2023, 15, 2073-8994, 1301, 10.3390/sym15071301
    6. J. Anitha, Indra Rajasingh, Hossein Rashmanlou, Propagating Sets in Sierpiński Fractal Graphs, 2025, 59, 0988-3754, 1, 10.1051/ita/2024016
  • Reader Comments
  • © 2014 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(7183) PDF downloads(1105) Cited by(4)

Figures and Tables

Tables(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog