Research article

Comparative analysis of copper and zinc based agrichemical biocide products: materials characteristics, phytotoxicity and in vitro antimicrobial efficacy

  • Received: 18 May 2016 Accepted: 20 July 2016 Published: 26 July 2016
  • In the past few decades, copper based biocides have been extensively used in food crop protection including citrus, small fruits and in all garden vegetable production facilities. Continuous and rampant use of copper based biocides over decades has led to accumulation of this metal in the soil and the surrounding ecosystem. Toxic levels of copper and its derivatives in both the soil and in the run off pose serious environmental and public health concerns. Alternatives to copper are in great need for the agriculture industry to produce food crops with minimal environmental risks. A combination of copper and zinc metal containing biocide such as Nordox 30/30 or an improved version of zinc-only containing biocide would be a good alternative to copper-only products if the efficacy can be maintained. As of yet there is no published literature on the comparative study of the materials characteristics and phyto-compatibility properties of copper and zinc-based commercial products that would allow us to evaluate the advantages and disadvantages of both versions of pesticides. In this report, we compared copper hydroxide and zinc oxide based commercially available biocides along with suitable control materials to assess their efficacy as biocides. We present a detailed material characterization of the biocides including morphological studies involving electron microscopy, molecular structure studies involving X-ray diffraction, phytotoxicity studies in model plant (tomato) and antimicrobial studies involving surrogate plant pathogens (Xanthomonas alfalfae subsp. citrumelonis, Pseudomonas syringae pv. syringae and Clavibacter michiganensis subsp. michiganensis). Zinc based compounds were found to possess comparable to superior antimicrobial properties while exhibiting significantly lower phytotoxicity when compared to copper based products thus suggesting their potential as an alternative.

    Citation: Parthiban Rajasekaran, Harikishan Kannan, Smruti Das, Mikaeel Young, Swadeshmukul Santra. Comparative analysis of copper and zinc based agrichemical biocide products: materials characteristics, phytotoxicity and in vitro antimicrobial efficacy[J]. AIMS Environmental Science, 2016, 3(3): 439-455. doi: 10.3934/environsci.2016.3.439

    Related Papers:

    [1] Nalin Fonseka, Jerome Goddard Ⅱ, Alketa Henderson, Dustin Nichols, Ratnasingham Shivaji . Modeling effects of matrix heterogeneity on population persistence at the patch-level. Mathematical Biosciences and Engineering, 2022, 19(12): 13675-13709. doi: 10.3934/mbe.2022638
    [2] Robert Stephen Cantrell, Chris Cosner, William F. Fagan . The implications of model formulation when transitioning from spatial to landscape ecology. Mathematical Biosciences and Engineering, 2012, 9(1): 27-60. doi: 10.3934/mbe.2012.9.27
    [3] Nazanin Zaker, Christina A. Cobbold, Frithjof Lutscher . The effect of landscape fragmentation on Turing-pattern formation. Mathematical Biosciences and Engineering, 2022, 19(3): 2506-2537. doi: 10.3934/mbe.2022116
    [4] Zhilan Feng, Robert Swihart, Yingfei Yi, Huaiping Zhu . Coexistence in a metapopulation model with explicit local dynamics. Mathematical Biosciences and Engineering, 2004, 1(1): 131-145. doi: 10.3934/mbe.2004.1.131
    [5] James T. Cronin, Jerome Goddard II, Amila Muthunayake, Ratnasingham Shivaji . Modeling the effects of trait-mediated dispersal on coexistence of mutualists. Mathematical Biosciences and Engineering, 2020, 17(6): 7838-7861. doi: 10.3934/mbe.2020399
    [6] Nariyuki Nakagiri, Hiroki Yokoi, Yukio Sakisaka, Kei-ichi Tainaka . Population persistence under two conservation measures: Paradox of habitat protection in a patchy environment. Mathematical Biosciences and Engineering, 2022, 19(9): 9244-9257. doi: 10.3934/mbe.2022429
    [7] Kunquan Lan, Wei Lin . Population models with quasi-constant-yield harvest rates. Mathematical Biosciences and Engineering, 2017, 14(2): 467-490. doi: 10.3934/mbe.2017029
    [8] Minjuan Gao, Lijuan Chen, Fengde Chen . Dynamical analysis of a discrete two-patch model with the Allee effect and nonlinear dispersal. Mathematical Biosciences and Engineering, 2024, 21(4): 5499-5520. doi: 10.3934/mbe.2024242
    [9] Robert Stephen Cantrell, Chris Cosner, William F. Fagan . Edge-linked dynamics and the scale-dependence of competitive. Mathematical Biosciences and Engineering, 2005, 2(4): 833-868. doi: 10.3934/mbe.2005.2.833
    [10] Kuang-Hui Lin, Yuan Lou, Chih-Wen Shih, Tze-Hung Tsai . Global dynamics for two-species competition in patchy environment. Mathematical Biosciences and Engineering, 2014, 11(4): 947-970. doi: 10.3934/mbe.2014.11.947
  • In the past few decades, copper based biocides have been extensively used in food crop protection including citrus, small fruits and in all garden vegetable production facilities. Continuous and rampant use of copper based biocides over decades has led to accumulation of this metal in the soil and the surrounding ecosystem. Toxic levels of copper and its derivatives in both the soil and in the run off pose serious environmental and public health concerns. Alternatives to copper are in great need for the agriculture industry to produce food crops with minimal environmental risks. A combination of copper and zinc metal containing biocide such as Nordox 30/30 or an improved version of zinc-only containing biocide would be a good alternative to copper-only products if the efficacy can be maintained. As of yet there is no published literature on the comparative study of the materials characteristics and phyto-compatibility properties of copper and zinc-based commercial products that would allow us to evaluate the advantages and disadvantages of both versions of pesticides. In this report, we compared copper hydroxide and zinc oxide based commercially available biocides along with suitable control materials to assess their efficacy as biocides. We present a detailed material characterization of the biocides including morphological studies involving electron microscopy, molecular structure studies involving X-ray diffraction, phytotoxicity studies in model plant (tomato) and antimicrobial studies involving surrogate plant pathogens (Xanthomonas alfalfae subsp. citrumelonis, Pseudomonas syringae pv. syringae and Clavibacter michiganensis subsp. michiganensis). Zinc based compounds were found to possess comparable to superior antimicrobial properties while exhibiting significantly lower phytotoxicity when compared to copper based products thus suggesting their potential as an alternative.


    A molecular graph is a graph such that its vertices correspond to the atoms and the edges to the bonds. Chemical Graph Theory is a branch of Mathematical Chemistry, which has an important effect on the development of the Chemical Sciences. A topological index is a numerical parameter mathematically derived from the graph structure. Several such topological indices have been considered in Theoretical Chemistry and have found some applications, especially in QSPR/QSAR study see [1,2].

    Let G be a finite, simple, connected graph with vertex set V(G) and edge set E(G). The degree dG(u) of a vertex u is the number of vertices adjacent to u. Let Δ(G) and δ(G) denote the maximum and minimum degree, respectively, among the vertices of G. The Revan vertex degree of a vertex uV(G) is defined as rG(u)=Δ(G)+δ(G)dG(u). The edge connecting the vertices u,vV(G) will be denoted by uv. We refer to [3] for undefined notations and terminologies.

    The first and second Revan indices of a graph G were introduced by Kulli in [4], and they are defined as

    R1(G)=uvE(G)(rG(u)+rG(v)),R2(G)=uvE(G)rG(u)rG(v).

    The F-Revan index of a graph is defined in [5] as

    FR(G)=uvE(G)(rG(u)2+rG(v)2).

    The Sombor index was introduced by Gutman in [6], defined as

    SO(G)=uvE(G)dG(u)2+dG(v)2.

    Recently, a large number of studies on Sombor indices have been published, see e.g., [7,8,9,10,11,12,13,14,15].

    The Revan Sombor index was proposed by Kulli and Gutman in [16] as

    RSO(G)=uvE(G)rG(u)2+rG(v)2.

    The first and second Revan (a,b)-KA indices of a graph G are defined (generalizing the Revan Sombor index), respectively, as

    RKA1a,b(G)=uvE(G)(rG(u)a+rG(v)a)b,RKA21,c(G)=uvE(G)(rG(u)rG(v))c.

    Recently, some (a,b)-KA indices have been used to analyze the properties of chemical compounds and nanostructures. In particular, they have been applied to the study of benzenoid systems and phenylenes, see for example [17,18,19,20].

    Therefore, motivated both by their practical applications as well as by their theoretical properties, mainly because of the close relationship with the well-known Sombor index [15], in this paper we perform analytical and statistical studies on the recently introduced Revan Sombor indices.

    Theorem 1. If G is a graph with maximum degree Δ and minimum degree δ, and cR{0}, then

    2c/2ΔcRKA21,c(G)RKA12,c/2(G)2c/2δcRKA21,c(G) if c>0,2c/2δcRKA21,c(G)RKA12,c/2(G)2c/2ΔcRKA21,c(G) if c<0,

    and the equality in each bound is attained for each c if and only if G is a regular graph.

    Proof. If c>0, then

    RKA12,c/2(G)=uvE(G)(rG(u)2+rG(v)2)c/2=uvE(G)(rG(u)rG(v))c(1rG(u)2+1rG(v)2)c/2uvE(G)(rG(u)rG(v))c(1Δ2+1Δ2)c/2=2c/2ΔcRKA21,c(G),

    and we obtain the converse inequality if c<0. The equality in each bound is attained if and only if rG(u)=rG(v)=Δ, i.e., dG(u)=dG(v)=δ for every uvE(G), that is, G is a regular graph.

    Also, if c>0, then

    RKA12,c/2(G)=uvE(G)(rG(u)rG(v))c(1rG(u)2+1rG(v)2)c/2uvE(G)(rG(u)rG(v))c(1δ2+1δ2)c/2=2c/2δcRKA21,c(G),

    and we obtain the converse inequality if c<0. The equality in each bound is attained if and only if rG(u)=rG(v)=δ, i.e., dG(u)=dG(v)=Δ for every uvE(G), that is, G is a regular graph.

    Notice that the limit as c approaches zero, the bounds, both upper and lower, are all equal to m (number of edges of G).

    If we choose c=1 in Theorem 1, then we obtain the following result.

    Corollary 2. If G is a graph, then

    2ΔR2(G)RSO(G)2δR2(G)

    and the equality in each bound is attained if and only if G is a regular graph.

    Theorem 3. If G is a graph with m edges and maximum degree Δ, and c>0, then

    RKA12,c/2(G)mΔcRKA11,c(G),

    and the equality in the bound is attained for each c if and only if G is a regular graph.

    Proof. Cauchy-Schwarz inequality gives

    RKA12,c/2(G)=uvE(G)(rG(u)2+rG(v)2)c/2uvE(G)Δc/2(rG(u)+rG(v))c/2(uvE(G)Δc)(uvE(G)(rG(u)+rG(v))c)=mΔcRKA11,c(G).

    If the equality in the bound is attained for some c>0, then rG(u)=rG(v)=Δ for every uvE(G) and so, G is a regular graph.

    If G is regular, then

    mΔcRKA11,c(G)=mΔcm2cΔc=m2c/2Δc=RKA12,c/2(G).

    If we choose c=1 in Theorem 3, then we obtain the following result.

    Corollary 4. If G is a graph with m edges and maximum degree Δ, then

    RSO(G)mΔR1(G)

    and the equality in the bound is attained if and only if G is a regular graph.

    We need the following well known inequalities for x1,x2>0 [21]:

    xα1+xα2<(x1+x2)α2α1(xα1+xα2)if α>1,2α1(xα1+xα2)(x1+x2)α<xα1+xα2if 0<α<1,(x1+x2)α2α1(xα1+xα2)if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if x1=x2.

    These inequalities allow to obtain the following result.

    Theorem 5. If G is a graph and a,b,cR{0}, then

    RKA1ab/c,c(G)<RKA1a,b(G)2bcRKA1ab/c,c(G) if b>c,bc>0,2bcRKA1ab/c,c(G)RKA1a,b(G)<RKA1ab/c,c(G) if b<c,bc>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b<0,c>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b>0,c<0,

    and the equality in each non-strict inequality is attained for each a,b,c if and only if each connected component of G is a regular graph.

    Proof. If α=b/c, x1=rG(u)a and x2=rG(v)a, then the previous inequalities give

    rG(u)ab/c+rG(v)ab/c<(rG(u)a+rG(v)a)b/c2b/c1(rG(u)ab/c+rG(v)ab/c)if b/c>1,2b/c1(rG(u)ab/c+rG(v)ab/c)(rG(u)a+rG(v)a)b/c<rG(u)ab/c+rG(v)ab/cif 0<b/c<1,(rG(u)a+rG(v)a)b/c2b/c1(rG(u)ab/c+rG(v)ab/c)if b/c<0,

    and the equality in each non-strict inequality is attained if and only if rG(u)=rG(v), i.e., dG(u)=dG(v).

    Hence, we obtain

    (rG(u)ab/c+rG(v)ab/c)c<(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif b/c>1,c>0,2bc(rG(u)ab/c+rG(v)ab/c)c(rG(u)a+rG(v)a)b<(rG(u)ab/c+rG(v)ab/c)cif b/c>1,c<0,2bc(rG(u)ab/c+rG(v)ab/c)c(rG(u)a+rG(v)a)b<(rG(u)ab/c+rG(v)ab/c)cif 0<b/c<1,c>0,(rG(u)ab/c+rG(v)ab/c)c<(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif 0<b/c<1,c<0,(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif b<0,c>0,(rG(u)a+rG(v)a)b2bc(rG(u)ab/c+rG(v)ab/c)cif b>0,c<0,

    and the equality in each non-strict inequality is attained if and only if dG(u)=dG(v).

    If we sum on uvE(G) these inequalities, then we obtain the desired inequalities. Furthermore, the equality in each non-strict inequality is attained for each a,b,c if and only if dG(u)=dG(v) for every uvE(G); and this happens if and only if each connected component of G is a regular graph.

    Remark 6. The excluded case b=c in Theorem 5 is not interesting, since RKA1ab/c,c(G)=RKA1a,b(G) if b=c.

    If we take a=2, b=1/2 and c=1/α in Theorem 5, then we obtain the following result for the Revan Sombor index.

    Corollary 7. If G is a graph and αR{0}, then

    RKA1α,1/α(G)<RSO(G)21/21/αRKA1α,1/α(G) if α>2,21/21/αRKA1α,1/α(G)RSO(G)<RKA1α,1/α(G) if 0<α<2,21/21/αRKA1α,1/α(G)RSO(G) if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if each connected component of G is a regular graph.

    The above bounds (in Theorem 5 and Corollary 7) are achieved for several families of graphs, e.g., the complete graph Kn, the complete bipartite graph Km,n, the cycle graph Cn, etc.

    Recall that one of the most studied topological indices is the first Zagreb index, defined by

    M1(G)=uvE(G)(dG(u)+dG(v))=uV(G)dG(u)2.

    Our next result relates M1 and RKA1α,1/α.

    Theorem 8. If G is a graph with m edges, maximum degree Δ and minimum degree δ, and αR{0}, then

    2(Δ+δ)m211/αRKA1α,1/α(G)M1(G)<2(Δ+δ)mRKA1α,1/α(G) if α>1,2(Δ+δ)mRKA1α,1/α(G)<M1(G)2(Δ+δ)m211/αRKA1α,1/α(G) if 0<α<1,M1(G)2(Δ+δ)m211/αRKA1α,1/α(G) if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if each connected component of G is a regular graph.

    Proof. If we take a=b=1 and c=1/α in Theorem 5, then we obtain the following result for RKA11,1(G):

    RKA1α,1/α(G)<RKA11,1(G)211/αRKA1α,1/α(G)if α>1,211/αRKA1α,1/α(G)RKA11,1(G)<RKA1α,1/α(G)if 0<α<1,211/αRKA1α,1/α(G)RKA11,1(G)if α<0,

    and the equality in each non-strict inequality is attained for each α if and only if each connected component of G is a regular graph.

    These inequalities and the identity

    RKA11,1(G)=uvE(G)(rG(u)+rG(v))=2(Δ+δ)mM1(G)

    finish the proof.

    The first and second (a,b)-KA indices of a graph G are defined in [17], respectively, as

    KA1a,b(G)=uvE(G)(dG(u)a+dG(v)a)b,KA2a,b(G)=uvE(G)(dG(u)adG(v)a)b.

    The following result relates the KA1a,b and RKA1a,b indices.

    Theorem 9. If G is a graph with maximum degree Δ and minimum degree δ, and a,bR{0}, then

    δabΔabKA1a,b(G)RKA1a,b(G)ΔabδabKA1a,b(G) if ab>0,ΔabδabKA1a,b(G)RKA1a,b(G)δabΔabKA1a,b(G) if ab<0,

    and the equality in each inequality is attained for each a,b if and only if G is a regular graph.

    Proof. Let us define the function f:[δ,Δ]×[δ,Δ](0,) by

    f(x,y)=(Δ+δx)a+(Δ+δy)axa+ya.

    Note that if a>0, the numerator of f is strictly decreasing on [δ,Δ] on each variable and its denominator is strictly increasing on [δ,Δ] on each variable. Hence, f is a strictly decreasing function on [δ,Δ] on each variable, and so,

    δaΔa=f(Δ,Δ)f(x,y)f(δ,δ)=Δaδa.

    If x=dG(u) and y=dG(v), then

    δaΔa(dG(u)a+dG(v)a)rG(u)a+rG(v)aΔaδa(dG(u)a+dG(v)a).

    Notice that if b>0, then the above inequality is preserved; thus

    δabΔab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bΔabδab(dG(u)a+dG(v)a)b.

    In contrast, if b<0, then the inequality above is not preserved; thus

    Δabδab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bδabΔab(dG(u)a+dG(v)a)b.

    If a<0, then f is a strictly increasing function on [δ,Δ] on each variable, and so,

    Δaδa(dG(u)a+dG(v)a)rG(u)a+rG(v)aδaΔa(dG(u)a+dG(v)a).

    Thus,

    δabΔab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bΔabδab(dG(u)a+dG(v)a)b

    if b<0, and

    Δabδab(dG(u)a+dG(v)a)b(rG(u)a+rG(v)a)bδabΔab(dG(u)a+dG(v)a)b

    if b>0.

    If we sum on uvE(G) these inequalities, then we obtain the desired inequalities. Furthermore, since f is a strictly monotone function, the equality in each inequality is attained for each a,b if and only if dG(u)=dG(v)=δ for every uvE(G) or dG(u)=dG(v)=Δ for every uvE(G); and this happens if and only if G is a regular graph.

    The Revan Sombor and the Sombor indices have the same value for every regular graph. The following result, which is a consequence of Theorem 9 with a=2 and b=1/2, relates these two indices for every graph.

    Corollary 10. If G is a graph with maximum degree Δ and minimum degree δ, then

    δΔSO(G)RSO(G)ΔδSO(G)

    and the equality in each inequality is attained if and only if G is a regular graph.

    Another remarkable topological index is the harmonic index, defined in [22] as

    H(G)=uvE(G)2du+dv.

    This index has attracted a great interest in the lasts years (see, e.g., [23,24,25,26,27]).

    The next result relates the RKA1a,1/a and the harmonic indices.

    Theorem 11. If G is a graph with maximum degree Δ and minimum degree δ, and aR, then

    21/aΔδH(G)RKA1a,1/a(G)<12(Δ+δ)2H(G) if a>1,2ΔδH(G)<RKA1a,1/a(G)22+1/a(Δ+δ)2H(G) if 0<a<1,RKA1a,1/a(G)22+1/a(Δ+δ)2H(G) if a<0.

    The equality in the first inequality is attained for each a if and only if each connected component of G is Δ-regular or δ-regular. The equality in the two last inequalities is attained for each a if G is a regular graph.

    Proof. Recall that

    (xa1+xa2)1/a<x1+x2211/a(xa1+xa2)1/aif a>1,211/a(xa1+xa2)1/ax1+x2<(xa1+xa2)1/aif 0<a<1,211/a(xa1+xa2)1/ax1+x2if a<0,

    for x1,x2>0, and the equality in each non-strict inequality is attained for each a if and only if x1=x2.

    Let us define the function g:[δ,Δ]×[δ,Δ](0,) by

    g(x,y)=((Δ+δx)a+(Δ+δy)a)1/a(x+y).

    If a>1, then

    g(x,y)<((Δ+δx)+(Δ+δy))(x+y)=(2(Δ+δ)(x+y))(x+y)(2(Δ+δ)(Δ+δ))(Δ+δ)=(Δ+δ)2,

    since x+y[2δ,2Δ]. Thus,

    (rG(u)a+rG(v)a)1/a<(Δ+δ)21dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)<12(Δ+δ)2H(G).

    Also,

    g(x,y)21+1/a((Δ+δx)+(Δ+δy))(x+y)=21+1/a(2(Δ+δ)(x+y))(x+y)21+1/a4Δδ

    since x+y[2δ,2Δ]. The equality in the previous bound is attained if and only if we have either x=y=δ or x=y=Δ. Therefore,

    (rG(u)a+rG(v)a)1/a21/aΔδ2dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)21/aΔδH(G).

    The equality in this bound is attained for each a>1 if and only if for each edge uvE(G) we have either dG(u)=dG(v)=Δ or dG(u)=dG(v)=δ, that is, each connected component of G is Δ-regular or δ-regular.

    Assume now a<1, then

    g(x,y)21+1/a((Δ+δx)+(Δ+δy))(x+y)=21+1/a(2(Δ+δ)(x+y))(x+y)21+1/a(2(Δ+δ)(Δ+δ))(Δ+δ)=21+1/a(Δ+δ)2,

    since x+y[2δ,2Δ]. Thus,

    (rG(u)a+rG(v)a)1/a22+1/a(Δ+δ)22dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)22+1/a(Δ+δ)2H(G).

    If G is a regular graph with m edges, then

    RKA1a,1/a(G)=21/aΔm=22+1/a(2Δ)2mΔ=22+1/a(Δ+δ)2H(G).

    If 0<a<1, then we have

    g(x,y)>((Δ+δx)+(Δ+δy))(x+y)=(2(Δ+δ)(x+y))(x+y)4Δδ

    since x+y[2δ,2Δ]. Therefore,

    (rG(u)a+rG(v)a)1/a>2Δδ2dG(u)+dG(v)

    for every uvE(G), and so,

    RKA1a,1/a(G)>2ΔδH(G).

    If we choose a=2 in Theorem 11, then we obtain the following result for the Revan Sombor index.

    Corollary 12. If G is a graph with maximum degree Δ and minimum degree δ, then

    2ΔδH(G)RSO(G)<12(Δ+δ)2H(G)

    and the equality in the first inequality is attained if and only if each connected component of G is Δ-regular or δ-regular.

    Recently, within a statistical approach to degree-based topological indices (TIs) on random graphs (see e.g., [28,29,30]), it has been shown that the average values of TIs like the average Randic index, R(G), the average harmonic index, H(G), and the average of some Sombor indices are highly correlated with the average Shannon entropy of the eigenvectors of the corresponding adjacency matrix [31,32]. Here, denotes the average over an ensemble of random graphs. This is a notable result because it puts forward the application of degree-based TIs beyond mathematical chemistry (it is relevant to add that random graphs have also been studied by means of eigenvalue-based TIs such as the Estrada index, the Laplacian Estrada index, and Rodriguez-Velazquez indices, see e.g., [33,34].) Specifically, on the one hand, certain average TIs can provide equivalent information than traditional, spectrum-based, random matrix theory measures. On the other hand, TIs may be used to predict spectral properties of random graphs.

    Therefore, motivated by potential applications, in what follows we apply for the first time (to our knowledge) Revan degree-based indices on random graphs. Indeed, since the inequalities obtained in Section 2 are not restricted to any particular type of graph, we first quantify them for random graphs and later we extend some of them to index average values, as needed in statistical studies of random graphs.

    Below we consider two prominent models of random graphs: Erdös-Rényi (ER) graphs and random geometric (RG) graphs. ER graphs [35,36] GER(n,p) are formed by n vertices connected independently with probability p[0,1]. While RG graphs [37,38] GRG(n,r) consist of n vertices uniformly and independently distributed on the unit square, where two vertices are connected by an edge if their Euclidean distance is less or equal than the connection radius r[0,2].

    In the Theorems presented in Section 2 we know the graph properties needed for the equalities involved have to be attained, however we do not really know how strong the inequalities could be for arbitrary graphs. Thus, in this subsection we quantify the inequalities of those Theorems on random graphs.

    To ease the quantification of the six Theorems presented in Section 2 we:

    (i) write the left inequality of Corollary 2 as

    0RSO(G)2ΔR2(G), (3.1)

    (ii) write Corollary 4 as

    0mΔR1(G)RSO(G), (3.2)

    (iii) write the right inequality in the second line of Corollary 7, for α=1/2, as

    0<RKA11/2,2(G)RSO(G), (3.3)

    (iv) write the left inequality in the second line of Theorem 8, for α=1/2, as

    0<M1(G)2(Δ+δ)m+RKA11/2,2(G), (3.4)

    (v) write the left inequality of Corollary 10 as

    0RSO(G)δΔSO(G), (3.5)

    and

    (vi) write the left inequality of Corollary 12 as

    0RSO(G)2δΔH(G). (3.6)

    Then, in Figure 1 [in Figure 2] we plot the right hand side of Eqs (3.1)–(3.6) as a function of the probability p [the connection radius r] of ER [RG] graphs of size n. For each value of p [r] we present results for 10 graph realizations.

    Figure 1.  Right hand side of (a) Eq (3.1), (b) Eq (3.2), (c) Eq (3.3), (d) Eq (3.4), (e) Eq (3.5) and (f) Eq (3.6) as a function of the probability p of Erdős-Rényi graphs GER(n,p) of size n. Each symbol corresponds to a single graph realization. Dashed lines in (c) [(d)] correspond to Eq (3.7) [(3.8)] with dG=(n1)p.
    Figure 2.  Right hand side of (a) Eq (3.1), (b) Eq (3.2), (c) Eq (3.3), (d) Eq (3.4), (e) Eq (3.5) and (f) Eq (3.6) as a function of the connection radius r of random geometric graphs GRG(n,r) of size n. Each symbol corresponds to a single graph realization. Dashed vertical lines in (a, b) [(e, f)] mark the maximum data values which occur at r0.8 [r0.64]. Dashed lines in (c) [(d)] correspond to Eq (3.7) [Eq (3.8)] with dG=(n1)r2(π8r/3+r2/2).

    For both random graph models, ER and RG graphs, we observe two different behaviors for the right hand side of Eqs (3.1)–(3.6): On the one hand, the right hand side of Eqs (3.1), (3.2), (3.5) and (3.6) are close to zero for n0 and p0 or r0, then grow with p or r, approach a maximum value at p=pmax or r=rmax, and finally decrease abruptly and become zero at p=1 or r=2. For ER graphs pmax is very close to 1, so the decrease of the right hand side of Eqs (3.1), (3.2), (3.5) and (3.6) is not displayed in Figure 1; while for RG graphs rmax0.8 for Eqs (3.1) and (3.2) and rmax0.64 for Eqs (3.5) and (3.6), as can be seen in Figure 2. On the other hand, the right hand side of Eqs (3.3) and (3.4) is close to zero for n0 and p0 or r0, then grow with p or r approaching their maximum values at p=1 or r=2. That is, the right hand side of Eqs (3.3) and (3.4) are increasing functions of n and p or r.

    Moreover, we can estimate the right hand side of Eqs (3.3) and (3.4) as follows. As it is shown in Ref. [39], for ER and RG graphs in the dilute limit (i.e., when dG1) we have that δΔdGrG; where dG and rG are the average degree and average Revan degree, respectively. Therefore, for ER and RG graphs in the dilute limit, we can write the right hand side of Eq (3.3) as

    RKA11/2,2(G)RSO(G)=uvE(G)(rG(u)1/2+rG(v)1/2)2uvE(G)rG(u)2+rG(v)2uvE(G)(2rG1/2)2uvE(G)2rG2=4rGm2rGm4dGm2dGm=422ndG2; (3.7)

    while the right hand side of Eq (3.4) can be written as

    M1(G)2(Δ+δ)m+RKA11/2,2(G)=uvE(G)(dG(u)+dG(v))2(Δ+δ)m+uvE(G)(rG(u)1/2+rG(v)1/2)2uvE(G)2dG4dGm+uvE(G)(2rG1/2)2=2dGm+4rGm2dGm+4dGm=2dGm=ndG2. (3.8)

    Indeed, in panels (c) and (d) of Figures 1 and 2 we plot (as dashed lines) Eqs (3.7) and (3.8), respectively, and observe a very good correspondence with the computational data, even for relatively small values of p or r. In Figure 1(c), (d) we used dG=(n1)p for ER graphs; while in Figure 2(c), (d) we used dG=(n1)f(r), with f(r)=r2(π8r/3+r2/2), for RG graphs.

    Since there is a growing interest in the average values of TIs, in particular when applied to random graphs (see e.g., [28,29,30,31]), it should be desirable to state inequalities for them. Fortunately, due to linearity, Theorem 5 can be straightforwardly extended to index average values; it reads as:

    Corollary 13. If {G} is an ensemble of graphs and a,b,cR{0}, then

    RKA1ab/c,c(G)<RKA1a,b(G)2bcRKA1ab/c,c(G) if b>c,bc>0,2bcRKA1ab/c,c(G)RKA1a,b(G)<RKA1ab/c,c(G) if b<c,bc>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b<0,c>0,RKA1a,b(G)2bcRKA1ab/c,c(G) if b>0,c<0.

    Moreover, from Theorem 8 we make the following conjecture:

    Conjecture 14. If {G} is an ensemble of graphs, with average maximum degree Δ, average minimum degree δ, average number of edges m, and αR{0}, then

    2(Δ+δ)m211/αRKA1α,1/α(G)M1(G) <2(Δ+δ)mRKA1α,1/α(G) if α>1,2(Δ+δ)mRKA1α,1/α(G)<M1(G) 2(Δ+δ)m211/αRKA1α,1/α(G) if 0<α<1,M1(G)2(Δ+δ)m211/αRKA1α,1/α(G) if α<0.

    Now, to validate Conjecture 14 we split and write the second inequality, with α=1/2, as

    0<M1(G)2(Δ+δ)m+RKA11/2,2(G) (3.9)

    and

    02(Δ+δ)m21RKA11/2,2(G)M1(G). (3.10)

    Then, in Figure 3 [in Figure 4] we plot the right hand side of Eqs (3.9) and (3.10) as a function of the probability p [the connection radius r] of ER [RG] graphs of size n. The averages are computed over ensembles of 106 random graphs. As can be clearly observed in these figures, the right hand side of Eqs (3.9) and (3.10) is larger than zero for all the combinations of n and p [r] used in this work; thus we validate Conjecture 14 on both ER and RG graphs. In Figure 3(c), (d) we compute m as m=n(n1)p/2, while in Figure 4(c), (d) we use m=n(n1)f(r)/2 with f(r)=r2(π8r/3+r2/2).

    Figure 3.  Right hand side of (a) Eq (3.9) and (b) Eq (3.10) as a function of the probability p of Erdős-Rényi graphs GER(n,p) of size n. Each data value was computed by averaging over 106 random graphs GER(n,p). Dashed lines in (a) are Eq (3.11) with dG=(n1)p. The dashed vertical line in (b) marks the maximum of the curves which occur at p0.5.
    Figure 4.  Right hand side of (a) Eq (3.9) and (b) Eq (3.10) as a function of the connection radius r of random geometric graphs GRG(n,r) of size n. Each data value was computed by averaging over 106 random graphs GRG(n,r). Dashed lines in (a) are Eq (3.11) with dG=(n1)r2(π8r/3+r2/2). The dashed vertical line in (b) marks the maximum of the curves which occur at r0.64.

    Also, from Figure 3(a) [from Figure 4(a)] we observe that the right hand side of Eq (3.9) is close to zero for n0 and p0 [r0], then grows with p [r] approaching their maximum value at p=1 [r=2]. While from Figure 3(b) [from Figure 4(b)] we note that the right hand side of Eq (3.10) is close to zero for n0 and p0 [r0], then grows with p [r], approaches a maximum value at p0.5 [r0.64], and finally decrease abruptly and become zero at p=1 [r=2].

    By the use of the same approximations that allowed us to get Eqs (3.7) and (3.8) in the previous subsection, we can also estimate the right hand side of Eq (3.9) for both ER and RG graphs in the dilute limit. This results in

    M1(G)2(Δ+δ)m+RKA11/2,2(G)ndG2. (3.11)

    Therefore, in Figure 3(a) [from Figure 4(a)] we plot Eq (3.11), as dashed lines, and observe good correspondence with the computational data.

    In this work, we study the recently introduced Revan Sombor index and the first and second Revan (a,b)-KA indices. First, we present new relations to provide bounds on these Revan indices of the Sombor family which also serve to relate them with other Revan indices (such as R1,2, which are Revan versions of the first and second Zagreb indices) and with standard degree-based indices (such as the Sombor index, the first and second (a,b)-KA indices, the first Zagreb index and the Harmonic index). We also quantify our relations on two models of random graphs: Erdös-Rényi graphs and random geometric graphs. Then, motivated by the growing interest in the application of topological indices in the analysis of random graphs, we extend some of our relations to index average values, so they can be effectively used for the statistical study of ensembles of random graphs. It is important to stress that our results on index average values are essentially obtained by averaging over the graph ensemble in the statistical sense rather than in the rigorous probabilistic sense, which is a topic we plan to explore in the near future.

    We hope our results may motivate further analytical as well as computational studies of Revan Sombor indices. Indeed, since the average Sombor indices have shown to work as complexity measures of random graphs and networks (equivalent to random matrix theory measures) [32], we plan to explore wether Revan Sombor indices could also be used in this direction.

    J. A. M. -B. thanks support from CONACyT (Grant No. 286633), CONACyT-Fronteras (Grant No. 425854), and VIEP-BUAP (Grant No. 100405811-VIEP2022), Mexico. The research of J. M. R. and J. M. S. was supported by a grant from Agencia Estatal de Investigación (PID2019-106433GBI00/AEI/10.13039/501100011033), Spain. J. M. R. was supported by the Madrid Government (Comunidad de Madrid-Spain) under the Multiannual Agreement with UC3M in the line of Excellence of University Professors (EPUC3M23), and in the context of the V PRICIT (Regional Programme of Research and Technological Innovation).

    The authors declare there is no conflict of interest.

    [1] Wen A, Balogh B, Momol MT, et al. (2009) Management of bacterial spot of tomato with phosphorous acid salts. Crop prot 28: 859-863.
    [2] Graham J, Gruber B, Bock C (2014) Research progress for integrated canker management. Citrus Industry 95: 20-24.
    [3] Arias-Estévez M, López-Periago E, Martínez-Carballo E, et al. (2008) The mobility and degradation of pesticides in soils and the pollution of groundwater resources. Agr Ecosyst Environ 123: 247-260.
    [4] Hoang TC, Rogevich EC, Rand GM, et al. (2008) Copper desorption in flooded agricultural soils and toxicity to the Florida apple snail (Pomacea paludosa): Implications in Everglades restoration. Environ Poll 154: 338-347.
    [5] Behlau F, Belasque J, Graham J, et al. (2010) Effect of frequency of copper applications on control of citrus canker and the yield of young bearing sweet orange trees. Crop Prot 29: 300-305. doi: 10.1016/j.cropro.2009.12.010
    [6] Schuler LJ, Hoang TC, Rand GM (2008) Aquatic risk assessment of copper in freshwater and saltwater ecosystems of South Florida. Ecotoxicol 17: 642-659. doi: 10.1007/s10646-008-0236-7
    [7] Behlau F, Canteros BI, Jones JB, et al. (2012) Copper resistance genes from different xanthomonads and citrus epiphytic bacteria confer resistance to Xanthomonas citri subsp. citri. Eur j plant pathol 133: 949-963. doi: 10.1007/s10658-012-9966-8
    [8] Behlau F, Canteros BI, Minsavage GV, et al. (2011) Molecular characterization of copper resistance genes from Xanthomonas citri subsp. citri and Xanthomonas alfalfae subsp. citrumelonis. Appl environ microb 77: 4089-4096.
    [9] Ed Etxeberria PG, Priyanka Bhattacharya, Parvesh Sharma, et al. (2016) Determining the size exclusion for nanoparticles in citrus leaves. Hortic Sci 51: 1-6.
    [10] Hendricks KE, Donahoo RS, Roberts PD, et al. (2013) Effect of copper on growth characteristics and disease control of the recently introduced Guignardia citricarpa on Citrus in Florida. Am J Plant Sci 4: 282.
    [11] Jorgensen JH, Turnidge JD (2015) Susceptibility Test Methods: Dilution and Disk Diffusion Methods*.
    [12] Young M, Santra S (2014) Copper (Cu)–Silica Nanocomposite Containing Valence-Engineered Cu: A New Strategy for Improving the Antimicrobial Efficacy of Cu Biocides. J agr food chem 62: 6043-6052. doi: 10.1021/jf502350w
    [13] Hans M, Erbe A, Mathews S, et al. (2013) Role of copper oxides in contact killing of bacteria. Langmuir 29: 16160-16166. doi: 10.1021/la404091z
    [14] Balogh B, Jones J, Momol M, et al. (2003) Improved efficacy of newly formulated bacteriophages for management of bacterial spot on tomato. Plant Dis 87: 949-954. doi: 10.1094/PDIS.2003.87.8.949
    [15] Baron M, Arellano JB, Gorge JL (1995) Copper and photosysten II-A controversial relationship. Physiol Plant 94: 174-180. doi: 10.1111/j.1399-3054.1995.tb00799.x
    [16] Kuepper H, Goetz B, Mijovilovich A, et al. (2009) Complexation and Toxicity of Copper in Higher Plants. I. Characterization of Copper Accumulation, Speciation, and Toxicity in Crassula helmsii as a New Copper Accumulator. Plant Physiol 151: 702-714.
    [17] Gupta N, Ram H, Kumar B (2016) Mechanism of Zinc absorption in plants: uptake, transport, translocation and accumulation. Rev Environ Sci Bio-Technol 15: 89-109. doi: 10.1007/s11157-016-9390-1
    [18] Haslett BS, Reid RJ, Rengel Z (2001) Zinc mobility in wheat: Uptake and distribution of zinc applied to leaves or roots. Ann Bot 87: 379-386. doi: 10.1006/anbo.2000.1349
    [19] Aslani F, Bagheri S, Julkapli NM, et al. (2014) Effects of Engineered Nanomaterials on Plants Growth: An Overview. Scientific World J 2014: 641759.
    [20] Chichiricco G, Poma A (2015) Penetration and Toxicity of Nanomaterials in Higher Plants. Nanomater 5: 851-873. doi: 10.3390/nano5020851
  • Reader Comments
  • © 2016 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(10841) PDF downloads(1472) Cited by(19)

Figures and Tables

Figures(9)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog