Influence of Topological Properties of Complex Networks on the Effect of Spectral Coarse-Grained Network ()

Lang Zeng^{}, Zhen Jia^{*}, Yingying Wang^{}

College of Science, Guilin University of Technology, Guilin, China.

**DOI: **10.4236/cn.2018.103008
PDF HTML XML
636
Downloads
1,119
Views
Citations

College of Science, Guilin University of Technology, Guilin, China.

Recently, some coarse-graining methods based on network synchronization have been proposed to reduce the network size while preserving the synchronizability of the original network. In this paper, we investigate the effects of the coarse graining process on synchronizability over complex networks under different average path lengths and different degrees of distribution. A large amount of experiments demonstrate a close correlation between the average path length, the heterogeneity of the degree distribution and the ability of spectral coarse-grained scheme in preserving the network synchronizability. We find that synchronizability can be well preserved in spectral coarse-grained networks when the considered networks have a longer average path length or a larger degree of variance.

Keywords

Complex Network, Synchronization, Spectral Coarse-Graining, Average Path Length, Degree Distribution

Share and Cite:

Zeng, L. , Jia, Z. and Wang, Y. (2018) Influence of Topological Properties of Complex Networks on the Effect of Spectral Coarse-Grained Network. *Communications and Network*, **10**, 93-104. doi: 10.4236/cn.2018.103008.

1. Introduction

Synchronization, as an emerging phenomenon of a population of dynamically interacting systems, is ubiquitous in nature and plays an important role within various contexts in biology, chemistry, ecology, sociology, and technology [1] . In the past two decades, synchronization of complex networks has been extensively investigated and discussed, both numerically and theoretically [2] - [8] . Unfortunately, many real-world networks are too huge to be dealt with current tools or algorithms. For example, human brain contains 10^{11} neurons and more than 10^{15} connections, which brings a big challenge to research on such networks. Especially for the dynamics of large-scale coupling nodes, large number of coupled differential equations result in the trouble of the computation and simulation. Many techniques and methods at the level of mesoscale networks are useless in large networks. Therefore, when using model reduction for simulation and analysis, it is crucial to be able to reduce the network size while keeping most of the relevant properties of the initial networks [1] . Essentially, the coarse-graining process is very similar to the problem of cluster finding or community detection in networks [9] [10] [11] [12] [13] .

Therefore, some coarse-grained methods are proposed to try to transform large-scale networks into mesoscale networks [14] [15] [16] . For example, D. Gfeller et al. [14] [15] proposed a spectral coarse-graining method (SCG) aiming at maintaining network synchronizability. Zhou and Jia proposed an improved spectral coarse-grained algorithm (ISCG) [16] based on the SCG method, which reduced the computational complexity and improved the synchronization capability. Chen et al. [1] investigated the effects of the coarse graining process on synchronizability over complex clustered networks. They found that a more prominent cluster structure implies a higher capability of the coarse graining strategy in preserving the synchronizability of the initial networks.

Many articles only studied the maintenance of the topological properties of the initial networks by the coarse-grained method, but few people discussed the influence of the topological properties of the initial network on the effect of the coarse graining method. In this paper, we find that the longer the average path length or the more heterogeneous the degree distribution of the initial network, the spectral coarse-grained algorithm becomes more effective in keeping the synchronization ability.

The rest of the paper is organized as follows. In Section 2, Mathematical basis is introduced. Influence of topological properties of complex networks on the effect of coarse-grained networks are proposed in Section 3. In Section 4, we verify our results through the real world networks. Finally, some conclusions are drawn in Section 5.

2. Mathematical Basis

2.1. Characterization of Network Synchronization Capability

Consider a general complex dynamic network with N nodes. The dynamic equation as follow:

${\stackrel{\dot{}}{x}}_{i}=f\left({x}_{i}\right)-c{\displaystyle \underset{j=1}{\overset{N}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{l}_{ij}H\left({x}_{j}\right),\text{\hspace{0.17em}}i=1,2,\cdots ,N$ (1)

where
${x}_{i}\in {R}^{n}$ is the n-dimensional state variable of the i_{th} node,
$c>0$ is the coupling strength,
$H\mathrm{:}{R}_{n}\to {R}_{n}$ is the inner coupling function, and the Laplacian matrix
$L={\left({l}_{ij}\right)}_{N\times N}$ describes the coupling topology of the network, in which
${l}_{ij}=-1$ if j connects to i (otherwise 0). The matrix L satisfies the dissipation coupling conditions:
${\sum}_{j=1}^{N}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{l}_{ij}=0$ . If the network is undirected and connected, then L is a symmetric and positive semi-definite matrix with nonnegative eigenvalues satisfying
$0={\lambda}_{1}<{\lambda}_{2}\le \cdots \le {\lambda}_{N}$ . If there is
$s\left(t\right)\in {R}^{n}$ and
$t\to \infty $ , then
${x}_{i}\left(t\right)\to s\left(t\right),i=1,2,\cdots ,N$ , the state of all nodes of the Equation (1) is fully (asymptotically) synchronized to
$s\left(t\right)$ and
$s\left(t\right)$ is called a synchronous state, and which
${\xi}_{i}$ is the variation of the i_{th} node, The variational equation is obtained by:

$\stackrel{\dot{}}{\xi}=Df\left(s\right)\xi -cDH\left(s\right)\xi {L}^{\text{T}}$ (2)

where $Df\left(s\right)$ and $DH\left(s\right)$ are the Jacobian matrices of $f\left(s\right)$ and $H\left(s\right)$ with respect to $s\left(t\right)$ , Diagonalizing Equation (2) yields the following form:

$\stackrel{\dot{}}{\eta}=\left[Df\left(s\right)\eta -c{\lambda}_{k}DH\left(s\right)\right]{\eta}_{k},k=2,\cdots ,N$ (3)

where ${\eta}_{k}$ is the eigenmode associated with the eigenvalue ${\lambda}_{k}$ of L. Generalize Equation (3) to get the main stability equation:

$\stackrel{\dot{}}{y}=\left[Df\left(s\right)-\alpha DH\left(s\right)\right]y$ (4)

The largest Lyapunov exponent of this equation is a function of real variables α, It is called the main stable function of network (1) [17] . For node’s dynamics, there exist $c{\lambda}_{2}<{\alpha}_{1}$ , such that the synchronized state is linearly stable. As a consequence, the network can be synchronized if ${\alpha}_{1}<c{\lambda}_{2},c{\lambda}_{N}<{\alpha}_{2}$ , which means that larger ${\lambda}_{2}/{\lambda}_{N}$ or larger ${\lambda}_{2}$ leads to better synchronizability. In this paper, both ${\lambda}_{2}/{\lambda}_{N}$ and ${\lambda}_{2}$ are used to characterize network synchronizability.

2.2. The Spectral Coarse-Graining Method

In 2008, D. Gfeller et al. proposed a spectral coarse graining method based on merging nodes with similar characteristic components in reference [15] . The aim is to obtain a simplified network that maintains the synchronization capability of the initial network. The network synchronization capability is represented by ${\lambda}_{2}/{\lambda}_{N}$ and ${\lambda}_{2}$ . Take the indicator of maintaining synchronization as an example, the algorithm is as follows:

First, determine which nodes merge. Let p^{2} denotes the eigenvector for the smallest nonzero eigenvalue
${\lambda}_{2}$ . Merging nodes that correspond to the same or similar components in p^{2}. Here
${p}_{\mathrm{max}}^{2}$ and
${p}_{\mathrm{min}}^{2}$ are the largest and the smallest components of p^{2}. Divide the elements in p^{2} evenly between
${p}_{\mathrm{max}}^{2}$ and
${p}_{\mathrm{min}}^{2}$ into I intervals. The smaller I is, the smaller the size of the coarse-grained network will be, and the more the size of the network is reduced.

Second, update edges and extract the coarse-grained network. Let the N nodes of the initial network be labeled with $i=1,2,\cdots ,N$ , and the coarse-grained network has $\stackrel{\u02dc}{N}$ nodes corresponding to $\stackrel{\u02dc}{N}$ groups, labeled with $C=1,2,\cdots ,\stackrel{\u02dc}{N}$ . The edges of the coarse-grained network (corresponding to the new Laplacian matrix $\stackrel{\u02dc}{L}$ ) can be updated by the following matrix product

$\stackrel{\u02dc}{L}=KLQ,$ (5)

where, $K\in {R}^{\stackrel{\u02dc}{N}\times N}$ and $Q\in {R}^{N\times \stackrel{\u02dc}{N}}$ are defined by

${K}_{Ci}={\Psi}_{C,{C}_{i}}/\left|C\right|;\text{\hspace{0.17em}}{Q}_{iC}={\Psi}_{C,{C}_{i}}.$ (6)

Here, $\left|C\right|$ is the cardinality of group C; ${C}_{i}$ is the index of the I’th node group, and Ψ is the Kronecker symbol.

2.3. The Random Interchanging Algorithm

To investigate the structural effects on network synchronizability, we use random interchanging algorithm [18] to adjust the average path length while keeping degree distribution unchanged. The process is as follows:

1) Randomly pick two existing edges ${e}_{1}={x}_{1}{x}_{2}$ and ${e}_{2}={x}_{3}{x}_{4}$ , such that ${x}_{1}\ne {x}_{2}\ne {x}_{3}\ne {x}_{4}$ and there is no edge between ${x}_{1}$ and ${x}_{4}$ as well as ${x}_{2}$ and ${x}_{3}$ .

2) Cross reconnect these four nodes., that is, connect ${x}_{1}$ and ${x}_{4}$ as well as ${x}_{2}$ and ${x}_{3}$ , and remove the edges ${e}_{1}$ and ${e}_{2}$ .

3) Ensure that the network is connected and calculate whether this interchange increases/decreases the network average path length. If it does, accept the new configuration, else restore the original network structure.

4) Repeat step 1) unless the desired average path length is achieved.

Because the algorithm is only reconnected, it does not change the degree of any node. So the degree distribution and degree sequence are fixed. Figure 1 provides a sketch map of random interchanging algorithm, which may help us understand the program flow.

3. Influence of Topological Properties of Complex Networks on the Effect of Coarse-Grained Network

3.1. Influence of the Average Path Length

We consider a small-world network with average degree $\langle deg\rangle =4$ , 1000 nodes. By random interchanging, we obtain four small-world networks with the same degree distribution and different average path lengths, respectively 5.5219, 5.5518, 5.9301 and 6.6772. Similarly, we consider the scale-free network with

Figure 1. A schematic diagram of random interchanging algorithm (cited in reference [18] ).

power exponent 2.05, 1000 nodes. Through random interchanging, we obtain four scale-free networks with the same degree distribution and different average path lengths, and the average path lengths are 4.1000, 4.2530, 4.3048, 4.4591.

Based on the spectral coarse-graining method, The initial network of N = 1000 nodes is merged into $\stackrel{\u02dc}{N}$ groups, and ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{N}$ denote the smallest non-zero eigenvalue and the largest non-zero eigenvalue of the new Laplacian matrix, respectively. The relationship between $\stackrel{\u02dc}{N}$ , ${\stackrel{\u02dc}{\lambda}}_{N}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ for different values of average path lengths are shown in Figure 2(a) Figure 2(c) and Figure 3(a) Figure 3(c), respectively. We find that the more the average path length increases, the better the synchronization ability of the network is maintained. This means that the smaller the average path length of the network, the larger the size of the coarse-grained network is needed to keep the values of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ unchanged with the corresponding values of the initial

Figure 2. (a) (c) Evolution of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ with respect to $\stackrel{\u02dc}{N}$ for WS small-world networks with a constant degree of distribution and different average path lengths; (b) (d) Evolution of the corresponding absolute error.

Figure 3. (a) (c) Evolution of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ with respect to $\stackrel{\u02dc}{N}$ for scale free networks with a constant degree of distribution and different average path lengths; (b) (d) Evolution of the corresponding absolute error.

network. On the other hand, Figure 2(b) Figure 2(d) and Figure 3(b) Figure 3(d) show the evolution trend of the absolute error of the eigenvalue ${\stackrel{\u02dc}{\lambda}}_{N}$ and the eigenratio ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ between the initial network and the coarse-grained network as a function of $\stackrel{\u02dc}{N}$ for different average path lengths. It can be seen that when the average path length decreases, the absolute error of ${\stackrel{\u02dc}{\lambda}}_{N}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ between the initial network and coarse-grained networks increases. This result shows that compared with a network with long average path length, a network with short average path length needs more $\stackrel{\u02dc}{N}$ to achieve the same absolute error.

In summary, regardless of the small-world networks or scale-free networks, as the average path length increases, the eigenvalue ${\stackrel{\u02dc}{\lambda}}_{2}$ and the eigenratio ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ of coarse-grained networks become closer to the original values of the initial networks, which is to say that, with the increase of the average path length, the spectral coarse-grained method becomes more effective in keeping the network synchronization.

In Figure 4(a), the initial network is small-world network with 1000 nodes and average degree $\langle deg\rangle =4$ . Different curves indicate the variance of the degree distribution of the generated small-world networks, which are 1.0840, 0.7620 and 0.4180, respectively. In Figure 4(b), the initial network is scale-free network with 1000 nodes and average degree $\langle deg\rangle =3.994$ . The variance of degree distributions expressed by each curve is 26.368, 19.558 and 15.348, respectively. It is obvious that the above two figures show positive correlation between average path length and clustering coefficient. Chen et al. [12] studied that the more obvious the clustering is, the better the coarse-graining effect is, which means that the longer the average path length, the better the ability of spectral coarse-graining method in preserving the network synchronizability.

3.2. Influence of the Degree Distribution

In order to measure the heterogeneity of the degree distribution of the network, we use the degree variance σ to represent the heterogeneity of the degree distribution. The greater the degree variance, the more inhomogeneous the network.

Similarly, we first obtain small-world networks with different four-degree variances by reconnecting probabilities p = 0.1, 0.2, 0.3, and 0.4. We use the random interchanging algorithm to equalize the average path lengths of these four networks. Then we get four networks with the same average path length and different degree distribution. The degree variance is 0.3657, 0.7640, 1.0840, and 1.3200, respectively. Equally, we consider the four scale-free networks with 1000 nodes, uniform average path length and different degree distributions. The degree variances are 14.2060, 15.2419, 18.6693, and 24.4806, respectively.

Figure 4. The relationship between average distance and clustering coefficient in the process of random interchanging.

Figure 5(a) Figure 5(c) and Figure 6(a) Figure 6(c) show the evolution of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ respect to $\stackrel{\u02dc}{N}$ for networks with uniform average path length and different degree variances. We observe that with smaller σ, the more uniformly distributed the network degree, the larger the size of the coarse-grained network needs to be, to remain the values of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ unchanged. In addition, Figure 5(b) Figure 5(d) and Figure 6(b) Figure 6(d) indicate the evolution trends of the absolute error of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ of the initial network and the coarse-grained networks, respectively. It can be seen that as the degree variance decreases, the error of the eigenvalues and eigenvalues ratio corresponding to the coarse-grained network and the initial network become larger. This also shows that the reduction of the variance of degree is not conducive to coarse graining network.

Figure 5. (a) (c) Evolution of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ with respect to $\stackrel{\u02dc}{N}$ for WS small-world networks with the same average path length and different degree of distribution; (b) (d) Evolution of the corresponding absolute error.

Figure 6. (a) (c) Evolution of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ with respect to $\stackrel{\u02dc}{N}$ for scale free networks with the same average path length and different degree of distribution; (b) (d) Evolution of the corresponding absolute error.

The simulation results show, with the increase of degree variance, the distribution of network degree becomes more heterogeneous in both small-world networks and scale-free networks. The ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ of coarse-grained network are more close to the corresponding values of the initial networks, that is to say, the effect of coarse-grained network will be better with the increase of degree variance.

4. Example

To demonstrate our conclusions, we consider two real world networks-scientific cooperation network ( $N=379$ ) [19] and protein-protein interaction network in budding yeast ( $N=2224$ ). We obtain four scientific cooperation networks with average path lengths of 4.4783, 5.5061, 5.5832, and 6.0419 through the random

Figure 7. (a) (c) (e) Evolution of ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ with respect to $\stackrel{\u02dc}{N}$ for scientific cooperation network with $N=379$ and protein-protein interaction network with $N=2224$ ; (b) (d) (f) Evolution of the corresponding absolute error.

interchanging algorithm. Moreover, we get three protein-protein interaction networks with average path lengths of 4.0122, 4.2111, and 4.3763 through the random interchanging algorithm. The results are shown in Figure 7. We find that the experimental results are consistent with above conclusion. As the average path length increases, the effect of coarse-graining method is better.

5. Conclusion

This paper studies the influence of the average path length and the heterogeneity of the degree distribution on the ability of spectral coarse-graining method in keeping the network synchronizability. According to the large number of simulation experiments, the average path length and the heterogeneity of degree distribution are closely related to keep the ${\stackrel{\u02dc}{\lambda}}_{2}$ and ${\stackrel{\u02dc}{\lambda}}_{2}/{\stackrel{\u02dc}{\lambda}}_{N}$ of the initial networks when applying spectral coarse-graining method to reduce the network size. The longer the average path length, the more inhomogeneous the degree distribution, which can better maintain the network synchronization ability in the coarse-grained process.

Acknowledgements

This project is supported by National Natural Science Foundation of China (Nos.61563013, 61663006) and the Natural Science Foundation of Guangxi (No.2018GXNSFAA138095).

Conflicts of Interest

The authors declare no conflicts of interest.

[1] |
Chen, J., Lu, J.A., Lu, X.F., Chen, G.R., et al. (2013) Spectral Coarse Graining of Complex Clustered Networks. Communications in Nonlinear Science and Numerical Simulation, 18, 3036-3045. https://doi.org/10.1016/j.cnsns.2013.03.020 |

[2] |
Jost, J. and Joy, M.P. (2001) Spectral Properties and Synchronization in Coupled Map Lattices. Physical Review E, 65, Article ID: 016201.
https://doi.org/10.1103/PhysRevE.65.016201 |

[3] |
Wang, X.F. and Chen, G.R. (2002) Synchronization in Scale-Free Dynamical Networks: Robustness and Fragility. IEEE Transactions on Circuits and Systems I, 49, 54-62. https://doi.org/10.1109/81.974874 |

[4] |
Barahona, M. and Pecora, L.M. (2002) Synchronization in Small-World Systems. Physical Review Letters, 89, Article ID: 054101.
https://doi.org/10.1103/PhysRevLett.89.054101 |

[5] |
Wang, X.F. and Chen, G.R. (2002) Synchronization in Small-World Dynamical Networks. International Journal of Bifurcation and Chaos, 12, 187-192.
https://doi.org/10.1142/S0218127402004292 |

[6] |
Motter, A.E., Zhou, C.S. and Kurths, J. (2005) Network Synchronization, Diffusion, and the Paradox of Heterogeneity. Physical Review E, 71, Article ID: 016116.
https://doi.org/10.1103/PhysRevE.71.016116 |

[7] |
Nishikawa, T. and Motter, A.E. (2006) Maximum Performance at Minimum Cost in Network Synchronization. Physica D: Nonlinear Phenomena, 224, 77-89.
https://doi.org/10.1016/j.physd.2006.09.007 |

[8] |
Zhou, J., Lu, J.A. and Lü, J.H. (2006) Adaptive Synchronization of an Uncertain Complex Dynamical Network. IEEE Transactions on Automatic Control, 51, 652-656.
https://doi.org/10.1109/TAC.2006.872760 |

[9] |
Fortunato, S. (2010) Community Detection in Graphs. Physics Reports, 486, Article ID: 75174. https://doi.org/10.1016/j.physrep.2009.11.002 |

[10] |
Lafon, S. and Lee, A.B. (2006) Diffusion Maps and Coarse-Graining: A Unified Framework for Dimensionality Reduction, Graph Partitioning, and Data Set Parameterization. IEEE Transactions on Pattern Analysis and Machine Intelligence, 28, Article ID: 13931403. https://doi.org/10.1109/TPAMI.2006.184 |

[11] |
Li, T. and Vanden-Eijnden, E. (2008) Optimal Partition and Effective Dynamics of Complex Networks. Proceedings of the National Academy of Sciences of the United States of America, 105, 7907-7912. https://doi.org/10.1073/pnas.0707563105 |

[12] |
Li, T. and Liu, J. (2009) Probabilistic Framework for Network Partition. Physical Review E, 80, Article ID: 026106. https://doi.org/10.1103/PhysRevE.80.026106 |

[13] |
Liu, J. and Liu, T. (2010) Coarse-Grained Diffusion Distance for Community Structure Detection in Complex Networks. Journal of Statistical Mechanics: Theory and Experiment, 2010, Article ID: 12030.
https://doi.org/10.1088/1742-5468/2010/12/P12030 |

[14] |
Gfeller, D. and De, L.R.P. (2007) Spectral Coarse Graining of Complex Networks. Physical Review Letters, 99, Article ID: 038701.
https://doi.org/10.1103/PhysRevLett.99.038701 |

[15] |
Gfeller, D. and De, L.R.P. (2008) Spectral Coarse Graining and Synchronization in Oscillator Networks. Physical Review Letters, 100, Article ID: 174104.
https://doi.org/10.1103/PhysRevLett.100.174104 |

[16] | Zhou, J., Jia, Z. and Li, K.Z. (2017) Improved Algorithm of Spectral Coarse Graining Method of Complex Network. Acta Physica Sinica, 66, Article ID: 060502. |

[17] |
Pecora, L.M. and Carroll, T.L. (1998) Master Stability Functions for Synchronized Coupled Dynamical Systems. Physical Review Letters, 80, 2109-2012.
https://doi.org/10.1103/PhysRevLett.80.2109 |

[18] |
Zhao, M., Zhou, T., Wang, B.H., Yan, G., Yang, H.J. and Bai, W.J. (2006) Relations between Average Distance, Heterogeneity and Network Synchronizability. Physica A: Statistical Mechanics and Its Applications, 371, 773-780.
https://doi.org/10.1016/j.physa.2006.03.041 |

[19] |
Newman, M.E.J. (2003) The Structure and Function of Complex Networks. SIAM Review, 45, 167-256. https://doi.org/10.1137/S003614450342480 |

Journals Menu

Contact us

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2022 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.