Abstract

Software-defined networking for IoT (SDN-IoT) has become popular owing to its utility in smart applications. However, IoT devices are limited in computing resources, which makes them vulnerable to Low-rate Distributed Denial of Service (LDDoS). It is worth noting that LDDoS attacks are extremely stealthy and can evade the monitoring of traditional detection methods. Therefore, how to choose the optimal features to improve the detection performance of LDDoS attack detection methods is a key problem. In this paper, we propose DIAMOND, a structured coevolution feature optimization method for LDDoS detection in SDN-IoT. DIAMOND is consisted of a reachable count sorting clustering algorithm, a group structuring method, a comutation strategy, and a cocrossover strategy. By analysing the information of SDN-IoT network features in the solution space, the relationship between different SDN-IoT network features and the optimal solution is explored in DIAMOND. Then, the individuals with associated SDN-IoT network features are divided into different subpopulations, and a structural tree is generated. Further, multiple structural trees evolve in concert with each other. The evaluation results show that DIAMOND can effectively select optimal low-dimension feature sets and improve the performance of the LDDoS detection method, in terms of detection precision and response time.

1. Introduction

Internet of Things (IoT) ecosystem is one of the most critical aspects of human lives, which facilitates a wide variety of applications in different domains such as smart homes, agriculture, healthcare, smart cities, smart grids, industrial automation (Industry 4.0), smart driving, and elderly assistance [1]. Hence, guaranteeing communication quality for IoT devices is an issue worth exploring. In this work, WiFi is utilized as the communication protocol between IoT devices and access points (APs), which are implemented through Mininet-wifi. In fact, WiFi provides low latency and low latency jitter necessary for IoT applications such as industrial production lines, IP telephony, video conferencing, and virtual desktop infrastructures. However, the computing resources, storage resources, and network capacity of IoT devices are limited to carry high-speed data transmission.

Recently, SDN-IoT is proposed to improve transmission quality [2]. Figure 1 shows the SDN-IoT architecture. In SDN-IoT, SDN changes the limitation of the network infrastructure, gives the network more flexibility, and simplifies policy implementation and network configuration [3]. SDN-IoT is anticipated to smartly route traffic and use underutilized network resources to deliver IoT data to the Internet. However, due to the stronger flexibility that SDN gives to the network, SDN-IoT faces many security threats [4].

LDDoS (Low-rate Distributed Denial of Service) attack is a serious threat for SDN-IoT, which exploits the vulnerabilities in network protocols to launch attacks and often realizes superior attack effects at a smaller attack cost. LDDoS attacks can elude the monitor of the traditional detection approach by sending low-rate packets in the way of periodic pulse to a victim [5]. Hence, how to effectively detect LDDoS is a huge challenge [6]. Feature selection makes an outstanding contribution to overcoming the invisibility of LDDoS, exploring LDDoS-related features, and improving the efficiency of LDDoS detection.

In recent years, feature optimization algorithms based on Evolutionary Algorithms (EAs) are popular [713]. However, the existing studies have suffered from the following limitations. (1)Existing EA-based feature optimization methods usually evolve individuals in a population, leading to the loss of diversity when evolving individuals(2)The existing methods do not consider the coevolution of similar individuals during evolution, which leads to the reduction of convergence ratio

In order to overcome the shortcomings mentioned above, a structured coevolution feature optimization method for LDDoS detection in SDN-IoT (DIAMOND) is proposed. The main contributions are concluded as below: (1)By dividing the population into several subpopulations based on the proposed reachable count sorting clustering algorithm and a group structuring method and executing coevolution based on the designed comutation strategy and cocrossover strategy, DIAMOND is proposed(2)A reachable count sorting clustering algorithm (BONNET) is designed to divide the population into subpopulations with different SDN-IoT network feature information, and each subpopulation is considered as a suboptimal solution set in the solution space(3)A group structuring method is designed to further structure a subpopulation into structural trees in order to sort the individuals in the subpopulation orderly based on SDN-IoT network information, where multiple structural trees evolve in concert(4)A comutation strategy is proposed based on the optimal subpopulation guidance direction, amount of information on SDN-IoT network features, and the evolutionary trajectory of the subpopulation, to move the individual towards the optimal solution in an orderly manner(5)A cocrossover strategy is designed to facilitate information exchange between different structural trees by exchanging individuals between different structural trees and the genes of individuals in the same tree

This section presents a review of intrusion detection systems for malicious traffic detection in SDN-IoT networks as well as the deployment of malicious network attack detection technologies for feature selection tasks.

Mafarj and Mirjalili [7] proposed a new wrapper feature selection approach based on whale optimization algorithm (WOA). In this work, two binary variants of the WOA algorithm are proposed to search the optimal feature subsets for classification. Emary et al. [8] proposed a new grey wolf optimization (GWO) algorithm. In that approach, the sigmoidal function is responsible for squeezing the continuously updated position; then, these values are stochastically compared with a threshold to find the updated binary grey wolf position. Abdel-Basset et al. [9] proposed a new grey wolf optimizer algorithm integrated with a two-phase mutation to solve the feature selection for classification problems based on the wrapper methods. Mafarja and Mirjalili [10] used two hybridization models to design different feature selection techniques based on WOA. In the first model, simulated annealing (SA) algorithm is embedded in WOA, while it is used to improve the best solution found after each iteration of WOA in the second model.

Roopak et al. [11] have implemented the nondominated sorting algorithm with its adapted jumping gene operator to solve the optimization problem. The proposed method exploits an extreme learning machine as the classifier for feature selection based on six important objectives for an IoT network. Mafarja et al. [12] presented a novel wrapper feature selection approach based on augmented WOA, which was adopted in the context of IoT attacks detection. Transfer functions (sigmoid function and tanh function) are applied in the real-value step vector such that a high probability of change is given to the dimension with a large value. Haddadpajouh et al. [13] proposed a multikernel support vector machine (SVM) for IoT cloud-edge gateway malware hunting, using the GWO technique.

Table 1 describes the advantages and limitations of the related works. As shown in Table 1, the common limitations in these literatures are the high complexity of the models and the tendency of falling into local optimal solutions. This is because the existing methods do not consider the coevolution of similar individuals during evolution, leading to a complex model, which in turn leads to the decrease of the convergence rate. Likewise, existing methods lose diversity in populations during evolution, due to the fact that they usually evolve individuals in only one population, leaving the model trapped in a local optimum solution. The performance of the algorithm is limited by these problems. Therefore, DIAMOND is designed to solve these problems in this work.

3. DIAMOND: A Structured Coevolution Feature Optimization Method for LDDoS Detection in SDN-IoT

3.1. Motivation

As a result of the rapid increase in the number of mobile-connected devices such as smartphones and tablets, SDN-IoT is attracting more and more attention from industry and academia. Therefore, the security of SDN-IoT becomes a serious issue [14]. Feature selection methods for network attack detection methods have been a hot research topic among researchers in recent years. However, the existing feature selection methods have some limitations such as follows: (1)Existing EA-based feature optimization methods usually evolve individuals in a population, leading to the loss of diversity when evolving individuals(2)The existing methods do not consider the coevolution of similar individuals during evolution, which leads to the reduction of convergence ratio

Existing feature selection methods all have the common limitations of high model complexity and the tendency to fall into local optimal solutions. Therefore, in order to solve these problems, this work enables multiple subpopulations to coevolve and makes subpopulations into ordered populations. In this way, multiple subpopulations can tract each other to avoid falling into local optimal solutions. At the same time, individuals within subpopulations are clearly distinguished to reduce the consumption of computational resources.

3.2. Overall Structure

Figure 2 shows the overall structure of DIAMOND in SDN-IoT. DIAMOND is deployed on the SDN controller and collects network traffic by sending collect traffic instructions to the switch. Generally, DIAMOND includes the following steps: (i)Step 1. DIAMOND sends the instruction to collect network traffic information.(ii)Step 2. Once the SDN controller receives the instruction, it sends OFPT_STATS_REQUEST messages to the switches to get the statistics of flow entries.(iii)Step 3. The switches will reply OFPT_STATS_REPLY messages to the SDN controller. These messages contain the stream data in the flow entries of the switches.(iv)Step 4. After the controller receives the OFPT_STATS_REPLY messages, it deencapsulates these messages and sends the collected statistic data to DIAMOND.(v)Step 5. In DIAMOND, the statistic data of flow entries collected from switches are used to calculate the initial SDN-IoT network features. Subsequently, these features will be coded using the binary coding method to generate the initial population that is consisted of some individuals.(vi)Step 6. The population is divided into several subpopulations by BONNET.(vii)Step 7. The structuring factor is calculated based on the fitness factor and the membership factor, which structures the individuals in the subpopulation into a structural tree.(viii)Step 8 and Step 9. Multiple structural trees evolve collaboratively with each other through the mutation and crossover strategies.(ix)Step 10 and Step 11. After evaluation and selection, the dominant individuals are retained to the next generation.(x)Step 12, Step 13, and Step 14. After several iterations, the individual with maximum fitness will be identified as the best individual. Subsequently, the best individual will be decoded into the corresponding feature set.

As described above, we use binary coding to translate the original SDN-IoT network feature set to the individual used in DIAMOND. We assume that the original feature set includes -dimensional gene sequences. Hence, we set -dimensional gene sequences as follows: is an individual, which presents a network feature set of SDN-IoT. When gene , it means -numbered feature is selected. Conversely, when gene , it means -numbered feature is not selected.

3.3. BONNET: Reachable Count Sorting Clustering Algorithm

In SDN-IoT, the subset SDN-IoT network features are encoded as binary individuals. The initialized individuals are defined as data points that are distributed in the solution space . When individuals share the same or similar SDN-IoT network features, the distribution of these data points in the solution space is relatively aggregated. Then, when individuals distributed in the solution space are clustered into a subpopulation, it means that SDN-IoT network features in individuals are related to the solution closest to the subpopulation. Therefore, we divide the individuals into several subpopulations in order to extract information about SDN-IoT network features in subpopulations and explore the relationship between SDN-IoT network features and optimal solutions. As a result, we design a reachable count sorting clustering algorithm (BONNET) to divide the population into several subpopulations. In BONNET, the following definition is proposed: (1)A data point is defined as a reachable point, which is within the hypersphere region with a certain data point . Here, is defined as the center of the sphere, and the length is recognized as the radius of the neighborhood(2)A subpopulation centrosome is defined as an individual located at the center of the subpopulation Ҡ

Further, the distance distribution between individuals in the population is denoted as . Subsequently, means neighborhood radius of , means the mean value of , and stands for the variance of .

We measure the distribution density of the population in terms of the number of times each individual is utilized as a reachable point. Hence, is defined as the number of times that is utilized as a reachable point. Then, means the set of descending order of , and stands for obtaining the corresponding individual. Further, the initial queue is defined as follows, where of the individuals decreases sequentially.

Subsequently, the individual is selected as the subpopulation centrosome which is removed from in turn, and the reachable points within the neighborhood radius of the individual are formed as a subpopulation {}. Then, the reachable points of are removed from . Finally, {} is joined to the subpopulation queue . As shown in Equation (4), is the number of subpopulations.

An example is shown in Figure 3, where , . When is removed from the initial queue and denoted the subpopulation centrosome , the reachable points , , and in are joined to the subpopulation . It is worth noting that there are still outliers and smaller subpopulations around the main subpopulations. Therefore, we merge these outliers and smaller subpopulations into the surrounding main subpopulations. Algorithm 1 describes the detailed steps of BONNET.

3.4. Group Structuring Method

In the subpopulations divided by BONNET, individuals are relatively close together, but

Input: Population number: N; Individuals: ; Sub-populations number:
Output: The sub-population:
,
: the add element operation.
: the reachable points of .
: the distance function.
: the function that take out first element.
: the set of descending order of .
: the function that obtain the corresponding individual.
: the function that get the largest sub-populations in .
1 for do
2 for do
3
4 end for
5 end for
6
7
8 for do
9 =0
10 for do
11 if then
12 +=1
13 end if
14 end for
15
16 end for
17 for do
18
19
20 while do
21
22
23
24
25 end while
26 while do
27
28
29
30
31 end while
32 return

it lacks group coordination during evolution. In order to distinguish the importance of different SDN-IoT network features among common individuals and enable important SDN-IoT network features to provide guidance to individuals in a subpopulation, we propose a group structuring method. During the group structuring, we combine the fitness and membership of the individuals construct the structural tree.

As shown in Equation (5), stands for the fitness of individual , and the fitness factor reflects the guidance ability of the individual in the subpopulation. When the individual is close to the optimal solution, of has a larger value. As shown in Equation (6), means the distance from the individual to the subpopulation centrosome . Then, means the number of individuals in the subpopulation, and the membership factor means the degree to which the individual belongs to the subpopulation centrosome . As shown in Figure 4, stratifies individuals in the subpopulation {} by distance, and there is some genetic similarity between individuals at the same level. However, individuals in the same layer have different values depending on whether they have the SDN-IoT network features associated with the optimal solution.

Generally, the optimal solution is distributed on one side of the subpopulation rather than inside the subpopulation, and individuals close to the optimal solution are provided with greater fitness. Hence, we design fitness factor minus membership factor as a structural factor.

As shown in Figure 5, the individuals at the same level share an identical degree of membership and gene similarity. The closer an individual is to the optimal solution, the stronger the ability to guide. As a result, the structured factor integrates individual guidance ability and gene similarity. We build a structural tree based on the value of in descending order, and the individuals with the same score are sibling nodes. As shown in Equation (8), is the number of structural trees.

When the structural tree coevolves, the sibling nodes have equivalent evolutionary guidance ability and similar genes. Therefore, only partial individuals of sibling nodes are randomly selected to calculate the fitness to reduce the consumption of computing resources and accelerate the convergence rate.

Figure 5 shows the diversity of different structural trees. indicates the gradient change of fitness. indicates that the SDN-IoT network features in the individuals of information and layers are more different. Due to the width of the structural tree is greater than the width of the structural tree , there are more individuals in the same layer, indicating that the gene pool of the structural tree is rich in network feature information. In the process of the iterative evolution of the population, when the fitness of an individual in the layer is greater than the average fitness of the node in the layer, the individual will transition to the layer.

To evaluate the convergence and stability of the structural tree, we define Equation (9) for the diversity of the structural tree:

As shown in Equation (9), presents the number of nodes in the layer, means the depth of the tree, and stands for the average fitness of the nodes in the layer. When the width of the structural tree and fitness gradient change greatly, the structural tree has stronger diversity.

3.5. Comutation Strategy

In DIAMOND, the balance between global search and local optimization is extremely important. Therefore, we design a comutation strategy to guide the evolution of the population. In this strategy, we define the movement trajectory of the subpopulation centrosome , the disturbance factor , and the diversity factor of the population structural tree. It is worth noting that the subpopulation centrosome is characterized as the evolutionary benchmark of the individuals in the subpopulation . Therefore, the evolutionary trajectory and evolutionary route of the subpopulation centrosome have important contributions to the evolution of the subpopulation . The movement trajectory of the subpopulation centrosome is defined as Equation (10). In Equation (10), means the last predicted

position, and stands for the current position, and then, for , we design it as the weight . Finally, is the population generation.

In Equation (11), means that the mutation individual generated by . denotes the -th generation of the population. stands for the movement trajectory of the subpopulation centrosome . is designed as the direction of the winning subpopulation. Furthermore, presents the disturbance factor generated by the . The comutation strategy ensures global search capability while accelerating convergence.

As shown in Figure 6, the movement trajectory of the subpopulation centrosome and the winning subpopulation guides the search direction for the individual . Further, projects the direction that is the sum of direction with the subpopulation centrosome movement trajectory and the winning subpopulation . The diversity factor regulates the evolutionary step length of individual . The larger the diversity of the structural tree is, the stronger the overall optimization ability will be, and each individual in the structural tree should be given a larger explore step size. However, the weaker the diversity of the structural tree is, the stronger the local optimization ability will be, and each individual in the structural tree should be given a smaller step size. The disturbance factor produced by can prevent the population from falling into a locally optimal solution.

In summary, the mutation individual is directed to the direction of the optimal solution, which contains a large amount of network feature information related to the optimal solution. It is worth noting that the SDN-IoT network features related to the optimal solution in , which carries more information than the SDN-IoT network feature unrelated to the optimal solution. Therefore, in order to reduce redundant features, we retain features closest to the optimal solution according to the network feature information carried by the mutation individual .

As in Equation (12), refers to the average amount of information carried by each feature in the mutation individual . is composed of SDN-IoT network features whose information content is greater than , and is more closely related to the SDN-IoT network features of the optimal solution. Further, features with the largest information value in are selected through . As in Equation (14), if the network feature is in , then =1. Otherwise, =0. This method can retain important SDN-IoT network features and filter out redundant SDN-IoT network features.

3.6. Cocrossover Strategy

We propose a cocrossover strategy to facilitate information exchange between different structural trees. The cocrossover strategy includes a crossover strategy between multiple subpopulations and a crossover strategy for structural trees within a subpopulation. As shown in Equations (15) and (16), is a crossover intermediate set that contains all individuals generated by the crossover operation. (1)As shown in Equation (15), individuals at the same level in multiple structural trees are randomly selected for exchange. For instance, individual and individual are randomly selected from the -th level individuals of structural tree and for exchange, which can increase the information interaction between subpopulations and avoid falling into local optimum(2)Individuals at the layer and layer of the structural tree are randomly selected for single-point crossover. For instance, the -th gene of the individual is exchanged with the allele of the individual

3.7. Detailed Design of DIAMOND

The IoT devices may include cameras, smartwatches, smartphones, various activity sensors, smart speakers, and so on. These IoT devices are connected to SDN switches, which are the endpoints of the service provider’s network [15].

In SDN-IoT, the LDDoS detection method needs to highly adaptively, fast, and accurately detect LDDoS. Hence, in DIAMOND, the AUC of detector performance, detection time (time), accuracy (ACC), and feature dimensions (DIM) are considered as the optimization objective, which is defined as follows:

Each individual within the population has its binary coded genotype. In Equation (18), is the number of iterations. is the population number. is the feature dimension.

After encoding the individuals, DIAMOND executes BONNET and group structuring method using Equations (2)–(8). Then, the comutation strategy is executed through Equations (11)–(14), and the cocrossover strategy is executed through Equations (15) and (16). In addition, by comparing the corresponding advantages and disadvantages of the new generation of the individual with the contemporary population, individuals who are more suitable for survival are selected as the new generations. As shown in Algorithm 2, individuals are encoded in lines 1~4. BONNET and group structuring in lines 6 and 7. The comutation strategy in lines 7~11. The cocrossover strategy and optimal individuals selecting strategy are implemented in lines 12~20.

Input: Network feature: ; Population Number: N; Feature Dimension: d; Iterations Number: G;
Output: The best individual and best fitness: ,
1 for do
2
3 Calculate
4 end for
5
6 Execute BONNET
7 Execute group structuring
6 for do
7 Update using Eq. 9
8 Update using Eq. 10
9 for do
10 Execute mutation using Eq. 11–14
11 end for
12 for do
13 Execute crossover using Eq. 15–16
14 if then
15
16 if then
17
18 end if
19 end if
20 end for
21
22 end for
23 return ,

The flow features are shown in Table 2. These features are extracted and calculated for each flow. Here, a flow is defined as traffic with the same source/destination MAC/IP address and the same source/destination port number.

4. Experimental Design and Result Analysis

4.1. Experiment Environment

To verify the logic of the program and evaluate its performance, DIAMOND is implemented on the RYU controller. In this work, Mininet-wifi implements the communication between IoT devices and access points (APs), using WiFi as the communication protocol between them. As shown in Figure 7, Mininet-wifi is utilized to simulate an SDN-IoT network that consists of eighty IoT endnotes, eight access points, and eight links, as previous work in [16]. The background traffic is to simulate the traffic environment of the IoT in the campus, which includes 40% HTTPS, 10% HTTP, 10% TCP, 10% SSDP, 5% DNS, and 5% NTP [17]. Specifically, the traffic is generated by applications on some IoT devices that require latency and low latency jitter, such as IP telephony and video conferencing. Further, the attack traffic consists of traffic generated by 20 IoT endnotes to launch LDDoS attacks against IoT servers. The attackers, victims, and service requests are all randomly generated. At the same time, the attack traffic accounts for 20% of normal traffic.

Further, in order to verify the detection ability of the optimized feature subset, eight different classifiers including support vector machine (SVM), K-Nearest Neighbor (KNN), Naïve Bayes (NB), Logistic Regress (LR), decision tree (DT), C4.5, Random Forest (RF), and AdaBoost (AB) have been applied to evaluate the feature optimization methods. In addition, DIAMOND is compared with other related state-of-the-art algorithms, which are BDIE [18] and BGDE3 [19]. As shown in Table 3, it describes the initial parameter settings in the DIAMOND.

4.2. Fitness Curve

Figure 8 shows the fitness curve of DIAMOND in the first 50 iterations. DIAMOND_FA, DIAMOND_FB, and DIAMOND_FL, respectively, represent the average fitness value, the best fitness value, and the lowest fitness value of DIAMOND. According to Figure 8, it can be

seen that after several iterations, the average fitness DIAMOND_FA steadily increases and reaches high values in all cases. At the early stage of evolution, the difference between the values of DIAMOND_FA, DIAMOND_FB, and DIAMOND_FL is small because they have not yet found the optimal solution set. In the middle of evolution, the difference between the values of DIAMOND_FB and DIAMOND_FL becomes larger because they have calculated different sets of solutions. At the later stage, almost all the values of DIAMOND_FB and DIAMOND_FL converge again, indicating that they have found out the optimal solution set.

4.3. Convergence Performance

Figure 9 exposes the convergence curve of DIAMOND. DIAMOND_1-3, respectively, stand for the convergence value of DIAMOND in 3 runs. DIAMOND_FA is the average convergence value of DIAMOND in 3 runs. DIAMOND_FA is the average convergence value. The convergence values of DIAMOND converge rapidly in the first half but gradually stabilize in the second half. The convergence value of DIAMOND decreases rapidly before 20 generations, while it stabilizes and decreases slowly after 20 generations. This indicates that the strong global search ability of DIAMOND leads to the exploration of an excellent solution set before 20 generations. Therefore, the population can converge to the outstanding individual quickly. Then, after 20 generations of iteration, DIAMOND performs an exact local exploration on the outstanding set of solutions that have been explored.

4.4. Dimension Reduction Ratio

Figure 10 shows the dimension reduction ratio of DIAMOND using eight classifiers, respectively. From Figure 10, it can be seen that the highest feature reduction ratio was 94.78%. Further, we compared the average feature reduction ratio of DIAMOND with BDIE and BGDE3. As shown in Figure 11, the average feature reduction ratio of BDIE, BGDE3, and DIAMOND are 91.67%, 48.72%, and 56.09%, respectively. The dimension reduction ratio indicates that applying DIAMOND in the actual SDN-IoT environment can effectively reduce the computational resource consumption of SDN controllers subject to LDDoS attacks.

4.5. Detection Performance

Table 4 shows the detection performance of different classifiers when using DIAMOND, BDIE, and BGDE3. AVG and STD represent the average detection performance and standard deviation of the algorithm, respectively. If the AVG values of the algorithms are equal, the algorithm with a smaller STD value is better. Besides, bold values in the table indicate the best value. The main purpose of Table 4 is to verify the robustness of the proposed algorithm on different classifiers. Obviously, DIAMOND performs better than BDIE and BGDE3 in most cases. It means that the features selected by the DIAMOND contain the key classification information. For instance, the average accuracy of DIAMOND is 92.30%, which is higher than 89.11% of BEID and 85.88% of BGDE3. In addition, all STD of DIAMOND is quite better than BDIE and BGDE3, which means that DIAMOND is more stable.

4.6. Detection Time

To evaluate the performance of DIAMOND in terms of detection time, the detection time is calculated in this evaluation. Table 5 describes the detection time of eight classifiers when using the BDIE, BGDE3, and DIAMOND. The detection efficiency of DIAMOND is quite better than that of BDIE and BGDE3. For instance, when using SVM, the detection time of DIAMOND is 106.72 ms, which is lower than 174.50 ms of BDIE and 164.59 ms of BGDE3. The detection time of DIAMOND is decreased by 39%, 35% compared to BDIE and BGDE3.

5. Conclusions

In this paper, we proposed DIAMOND, a structured coevolution feature optimization method for LDDoS detection in SDN-IoT. More specifically, a reachable count sorting clustering algorithm, a group structuring method, a comutation strategy, and a cocrossover strategy are proposed. The evaluation demonstrates that DIAMOND can effectively improve the detection accuracy, reduce the size of feature subsets, and achieve shortened detection time compared with the baseline methods.

Data Availability

The experimental data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Key Research and Development Program of China under Grant No. 2019YFB1803500, the National Natural Science Foundation of China (NSFC) under the Grant Nos. 61902085 and 62102111, the Guizhou Provincial Science and Technology Foundation under the Grant No. [2020]1Y267, the Scientific Research Foundation for Introduced Talents of Guizhou University under the Grant No. (2019)52, and the Development Project of Ordinary University Young Scientific and Technological Talents of Guizhou Provincial Science under the Grant No. Qian Jiao He KY Zi[2021]136.