Abstract

High energy consumption in Ethernet Passive Optical Networks (EPONs) has caused intense research on energy-saving methods in recent years. The most common method for EPON energy saving is to cause optical network units (ONUs) to sleep which are idle or low loaded. With the popularity of multimedia applications, the ONUs carry not only unicast traffic but also more and more multicast traffic. However, existing studies mainly focus on ONUs that support only unicast traffic. The many differences between multicast traffic and unicast traffic make it necessary to design novel energy-saving methods. This paper proposes a multicast-traffic-oriented energy-saving algorithm for EPONs, called the energy-saving algorithm for multicast traffic (ESMT). To save as much energy as possible, the proposed algorithm uses a hybrid sleep mode (composed of a deep sleep state and an independent sleep mode) in which not only can ONU enter a deep sleep state, but also the ONU transmitter and receiver can sleep independently. Simulation results show that the proposed algorithm, which is oriented to multicast traffic in EPONs, is more energy efficient than other algorithms.

1. Introduction

In recent decades, the rapid development of modern industries has intensified the greenhouse effect and has threatened the survival and development of humankind. Broadband access networks, one of the last-mile solutions, consume about 75% of all network energy [1]. Therefore, decreasing the energy consumption of broadband access networks is an effective way to green all networks. Furthermore, the Ethernet Passive Optical Network (EPON) has become one of the most feasible broadband access technologies, for its low energy consumption and wide deployment [2, 3]. Consequently, energy saving is a topic of paramount importance in EPON contexts today.

The typical EPON is structured as a tree topology, where the optical line terminal (OLT) acts as the root node and optical network units (ONUs) act as the leaf nodes. The splitter bridges the OLT and the ONUs through bidirectional optical fibers. Upstream, ONUs contend for the same shared wavelength channel using time division. Downstream, the OLT can support three kinds of traffic: unicast, multicast, and broadcast. The tree topology enables the EPON to be perfectly suitable for multicast and broadcast traffic. With the popularity of multimedia applications, more and more multicast traffic and unicast traffic are sharing the same ONU. As shown in Figure 1, the OLT keeps separate buffers for multicast and unicast traffic. The same-colored ONUs belong to the same multicast group and receive the same multicast traffic. All ONUs in a multicast group must be active simultaneously to receive the traffic designated for the group. This is the major difference between an EPON oriented to unicast traffic and one oriented to multicast traffic.

Both unicast and multicast traffic are broadcast by the OLT to all ONUs. The traffic can only be received by the destination ONUs and is discarded by the other irrelevant ONUs. Studies show that more than 95% of the traffic is discarded, which means a great deal of energy is wasted on unnecessary processing by the irrelevant ONUs. Therefore, the irrelevant ONUs stay in a low-power or sleep state to save energy. However, most literature focuses on energy saving oriented to unicast traffic [418]. In this work, we propose an energy-saving algorithm for multicast traffic (ESMT) to reduce the energy consumption of EPON ONUs oriented to multicast traffic, while maintaining the traffic delay at an acceptable level.

Our ESMT considers not only downstream traffic, but also upstream traffic. ONUs in the same group must receive the downstream multicast traffic at the same time; however, the upstream transmissions of these ONUs cannot proceed simultaneously. Therefore, the ONU hybrid sleep mode presented in our previous work [5] is a better way to solve the energy-saving problem in this situation. In other words, the ONU receivers in the same group can sleep simultaneously after the reception; meanwhile, the ONU transmitters can sleep independently, according to their own idle time. If there is no upstream or downstream traffic, the ONU goes into a deep sleep state. However, due to multicast characteristics, the energy saving of ONUs oriented to multicast traffic in hybrid sleep mode is much more complicated than in [5].

To the best of our knowledge, there is no prior study concentrating on energy saving that considers both upstream and downstream multicast traffic. More importantly, because it makes the ONU transmitter and receiver sleep separately, our ESMT algorithm is more energy efficient than other algorithms.

The remainder of this paper is organized as follows. Section 2 overviews related work. In Section 3, we formulate the problem of minimizing ONU energy consumption. Section 4 illustrates the proposed ESMT algorithm. In Section 5, we detail the simulation results. Finally, Section 6 concludes the paper.

In existing work, ONU energy saving can be realized by making part of the ONU elements sleep by dynamic bandwidth allocation (DBA). Yin et al. [15] studied the stability of predictor-based dynamic bandwidth allocation over EPON; bandwidth utilization can be improved by adaptive control at the OLT. The ONU dozing mode, specified in the standard ITU-T G.984 for Gigabit PON (GPON), saves energy by turning off the ONU transmitter when there is no upstream data to transmit [16]. Zhang et al. [17] proposed a sleep control scheme composed of downstream traffic scheduling rules at the OLT and sleep control rules at ONUs, which reduces energy consumption and improves bandwidth utilization effectively. Dhaini et al. [8] proposed a green bandwidth allocation (GBA) framework to achieve the maximum possible energy saving. In GBA, hybrid cyclic/deep sleep is enabled. After sending the request message to and receiving a grant message from the OLT, each ONU remains idle for a period of time, during which the ONU can enter an OLT-based or ONU-based deep sleep mode. The doze mode happens after the ONU-based deep sleep mode in order to receive a grant message.

Dias and Wong [10, 11] demonstrated that energy saving can be achieved when the idle time of a vertical-cavity surface-emitting laser (VCSEL) ONU is more than the sleep-to-active transition time. The VCSEL ONU can turn into sleep or doze modes during its idle time. In particular, the VCSEL ONU will enter sleep mode when its idle time is more than 2 ms or turn into doze mode when its idle time is between 330 ns and 2 ms. Energy-saving methods in these mentioned works are based on a combination of existing sleep modes. However, these methods are not sufficiently flexible because the upstream and downstream transmissions do not finish at the same time. Thus, in [5], we presented the independent sleep mode for ONUs, in which the transmitter and receiver of the ONU can sleep independently according to upstream and downstream traffic, respectively. Hosseinabadi and Ansari [18] also considered independent sleep of the transmitter and the receiver by dividing nine states for ONU in different situations. Both upstream and downstream traffic are taken into consideration.

Furthermore, all the above works only present energy-saving strategies for unicast traffic. Because of the differences between unicast and multicast traffic, the sleep processes of ONUs oriented to multicast traffic cannot be handled by the strategies in [418]. The authors in [1921] mainly studied the bandwidth allocation problem of multicast traffic in TDM-PON (time division multiplexing-PON), WDM-PON (wavelength division multiplexing-PON), and TDM-WDM-PON, respectively. However, they did not consider the energy-saving issue. Currently, only a few previous energy-saving methods are oriented to multicast traffic in EPONs. Valcarenghi et al. [22, 23] proposed an energy-efficient multicasting solution by setting the ONUs belonging to the groups which were not receiving data into sleep mode. The ONU sleep mode is combined with optimal multicast group scheduling to reduce the wake-up overhead and energy consumption. However, the author only considered downstream traffic. In fact, putting the ONUs into sleep mode should take into account both upstream and downstream traffic to ensure the Quality of Service (QoS). This is exactly the focus of this paper.

3. System Model

In this section, we define the system model in terms of the ONU sleep mode, the energy consumption model, and the problem formulation to minimize ONU energy consumption.

3.1. Sleep Mode

Downstream, ONUs belonging to the same multicast group can receive the multicast traffic simultaneously. However, the upstream channel is shared by all ONUs, in the way of time division multiple access (TDMA); that is, at any time, only the polled ONU can transmit its upstream traffic. Therefore, the other unpolled ONUs’ transmitters are idle. In other words, if we utilize the traditional sleep mode, in which the ONUs do not sleep unless there is no upstream or downstream traffic, unnecessary energy will be wasted. Consequently, the sleep processes of ONUs oriented to multicast traffic become more complicated. However, more energy can be saved by using the hybrid sleep mode of transmitter and receiver, which was presented in our previous work [5]. The transmitter and receiver sleep independently according to their own traffic. As illustrated in Figure 2, , , and belong to the same group and their receivers sleep simultaneously after time . Upstream, we assume that the ONUs transmit their data by ascending order of their own index numbers. Therefore, transmitter begins to sleep at time . Similarly, transmitter sleeps before time and after time . transmitter sleeps before time and after time . As a result, and can enter deep sleep mode at time , and can enter deep sleep mode at time . Since the hybrid sleep mode in [5] can save more energy than traditional sleep modes, this result is also appropriate for ONUs oriented to multicast traffic. In other words, we efficiently combine the characteristics of multicast traffic and hybrid sleep mode, intending to minimize the energy consumption of ONUs with both upstream and downstream traffic (including unicast and multicast traffic).

3.2. Energy Consumption Model

Since we utilize the hybrid sleep mode, we can divide the ONU into three parts: the transmitter, the receiver, and the common modules [24]. The common modules refer to the parts that always keep working, even if the ONU goes into deep sleep mode. Certainly, these parts also work when the transmitter sends upstream traffic or the receiver receives downstream traffic. The notations related to ONU energy consumption are shown in Notations.

The energy consumption of each ONU in Figure 2 is represented in Figure 3. Since Figure 2 only shows the sleep process in one cycle, we are not sure whether to wake up the ONU at the beginning of the next cycle. Hence, the overhead time is drawn by a dotted line, and is not marked in Figure 3. The ONU energy consumption can be composed of four parts, namely, energy consumption of the transmitter, energy consumption of the receiver, energy consumption of the common modules, and energy consumption produced by the ONU wake-up process.

We use to denote the energy consumption of . Therefore, can be formulated as follows: where , , and denote the energy consumed by the transmitter and the receiver of during their active time and the common modules of during the running time, respectively. denotes the energy consumption of during its wake-up process from deep sleep mode. The wake-up process consists of clock recovery and synchronization. is large enough that it cannot be omitted. Hence, researchers are engaged in reducing the time and energy consumption of clock recovery and synchronization [3]. Frequent wake-ups significantly increase the ONU energy consumption. Hence, it is necessary to reduce the ONUs’ wake-up times.

3.3. Problem Formulation

Based on the hybrid sleep mode, we attempt to minimize the total energy consumption of all ONUs. Therefore, the objective of this paper can be defined as the minimum energy consumption (MEC) problem. For simplicity in problem formulation, we introduce the following notations:: the set of all ONUs and multicast groups, respectively.: index number of ONUs, .: index number of multicast groups, .: index number of upstream transmission slots, .: index number of downstream transmission slots, .: matrix with rows and columns.: binary variable, taking 1 if the receiver of is active in downstream transmission slot and 0 otherwise.: binary variable, taking 1 if the transmitter of is active in upstream transmission slot and 0 otherwise.: binary variable, taking 1 if needs to wake up at the beginning of downstream transmission slot and 0 otherwise.: binary variable, taking 1 if is subscribed to multicast group and 0 otherwise.

We assume that there are upstream transmission slots. Each slot is assigned to an ONU, and the length of the slot is the transmission time of this ONU’s upstream traffic. We also assume that there are downstream transmission slots, and each downstream transmission slot is assigned to a multicast group. In other words, downstream transmission slot is assigned to multicast group , where ; and the length of downstream transmission slot is the transmission time of the downstream traffic of multicast group . It is worth noting that unicast traffic can be seen as multicast traffic of a multicast group with only one member. Hence, the energy consumption of in (1) can also be represented as follows:where denotes the working period of downstream transmission slot of the receiver. As a result, the MEC problem can be formulated as follows:

The objective formulated in (7) is to minimize the energy consumption of . The constraint in (8) ensures that all ONU receivers belonging to the same multicast group must stay active during the transmission of multicast traffic for the multicast group. The constraint in (9) determines whether to wake up the ONU at the beginning of downstream transmission slot. The constraint in (10) ensures that each ONU can stretch across one, and only one, upstream transmission slot. The constraint in (11) enables each upstream transmission slot to be assigned to only one ONU. To solve the MEC problem, we propose the ESMT algorithm described in the next section.

4. Energy-Saving Algorithm for Multicast Traffic (ESMT)

Due to the transmission characteristics of upstream and downstream traffic, this paper introduces the hybrid sleep mode—that is, the combination of independent sleep mode and deep sleep mode mentioned in Section 3.1—to solve the sleep problem of ONUs oriented to multicast traffic in EPONs. Based on the hybrid sleep mode, we propose the ESMT algorithm.

In this paper, multicast group is denoted by , where denotes the number of ONUs contained in multicast group and is used to store the index numbers of the ONUs in multicast group. For simplicity, we use to represent . Each multicast group stretches across only one slot, and the slot length depends on the traffic transmission time of the multicast group that stretches across this slot. Furthermore, each slot can only carry the downstream traffic of one multicast group. In addition, we consider downstream unicast traffic as a particular kind of multicast traffic where .

Figure 4 shows the traffic scheduling in one cycle. The first multicast group contains two ONUs, that is, and . Hence, we can obtain the first multicast group , where . Similarly, we can obtain the second multicast group , where , the third multicast group , where , and the fourth multicast group , where . We can see that the OLT first sends a GATE control message to each ONU, in which the upstream transmission slot, sleep start time, and wake-up time for each ONU are included. Therefore, before sending the GATE control message, the OLT should obtain the information about transmission slots for both upstream and downstream traffic. Especially downstream, a certain ONU may receive downstream traffic in successive slots, when sending a GATE control message makes the computation for the ONU’s sleep and wake-up processes much more complicated.

To determine the sleep processes of the ONUs’ transmitters, we sort all ONUs in a multicast group by ascending order of their own index numbers. Therefore, we can obtain , , , and from all the sequential multicast groups. Then, removing the nonoverlapping indices in sequence from the four number sets, we get the upstream transmission sort as .

Since the wake-up processes consume considerable energy that cannot be neglected, we must be aware of the wake-up times for each ONU. First, the start time and the end time of both upstream and downstream transmissions should be determined. Then, the ONUs’ wake-up times can be counted. The entire process is described in the following steps.

Step 1. Compute the start time and the end time of the upstream transmission. Let denote the index number of ONU that is sending upstream traffic; then, we need to compute the start time and the end time for each as and , where denotes the allocated bandwidth of , is the guard time between two successive transmissions, and is the transmission rate. Thus, we can obtain the start time and end time of the upstream transmission of (). Based on these two time points, we can determine the sleep time of the transmitter and compute the sleep start and end time of the transmitter of [5].

Step 2. Compute the start time and the end time for each downstream transmission slot as and , where denotes the traffic load of downstream transmission slot. If belongs to () multicast group, we use and to denote the start time and the end time of in downstream transmission slot. Based on these two time points, we can determine the receiver’s sleep time and compute its sleep start and end time [5].

Step 3. At the beginning of each polling cycle, the OLT will send a GATE control message to each , according to the order obtained above. Hence, the wake-up times of each are at least equal to 1. In addition to the wake-up process at the beginning of each polling cycle, whether wakes up at the beginning of the other downstream transmission slots within the cycle is determined by the following situations.

The end time of the upstream transmission is during downstream transmission slot, and is in a deep sleep state between and . Take in Figure 4 as an example and let . As the receiver has entered the sleep state before , only the transmitter and the common modules consume energy during . Assuming that goes into deep sleep state during , its saved energy during can be computed as . Therefore, if and only if both and (the energy consumption constraint) are met, can enter deep sleep state during and wake up before the transmission of downstream data , and the wake-up times increase by 1. Consequently, we need to insert a GATE control message into the downstream queue to enable to receive it before and update the start time and end time of all the downstream transmissions after .

As is shown in Figure 5, the start time of ’s upstream transmission is between two downstream transmissions; let , where is the time used to abstract the useful information from the GATE control message. As the transmitter has entered the sleep state between and , only the receiver and the common modules consume energy during . Assuming that goes into deep sleep state during , the saved energy of during can be computed as . Therefore, if and only if both and (i.e., the energy consumption constraint) are met, can enter deep sleep state during and wake up before the upstream transmission, and the wake-up times increase by 1. Consequently, we need to insert a GATE control message into the downstream queue to enable to receive it before the upstream transmission, which leads to , where denotes the transmission time of a GATE control message. Then, update the start and end time of all the downstream transmissions after .

As is shown in Figure 6, both the start time and end time of the upstream transmission are between two downstream transmissions. We address the case before following situation and the case after following situation . If both cases need to wake up , we insert a GATE control message for each case; the locations of the GATEs are shown in Figure 6. After inserting GATE control messages, we let and update the start time and end time of all downstream transmissions after .

As is shown in Figure 7, there is no upstream transmission between two downstream transmissions. We denote the end time of the first downstream transmission by and the start time of the second downstream transmission by and let . As the transmitter has entered the sleep state before , only the receiver and the common modules consume energy during . Assuming that goes into deep sleep state during , the saved energy of during can be computed as . Therefore, if and only if both and (i.e., the energy consumption constraint) are met, can enter the deep sleep state during and wake up before , and the wake-up times increase by 1. Consequently, we must insert a GATE control message into the downstream queue to enable to receive it before , which leads to . Then, update the start and end time of all downstream transmissions after .

As is shown in Figure 8, there is no downstream transmission between two upstream transmissions. We denote the end time of the first upstream transmission by and the start time of the second upstream transmission by and let . As the receiver has entered sleep state at , only the transmitter and the common modules consume energy during . Assuming that goes into deep sleep state during , the saved energy of during can be computed as . Therefore, if and only if and (i.e., the energy consumption constraint), can enter deep sleep state during and wake up before , and the wake-up times increase by 1. Consequently, we must insert a GATE control message into the downstream queue to enable to receive it before . Then, update the start and end time of all downstream transmissions after .

In our proposed algorithm, the OLT takes charge of the upstream and downstream slot allocations and determines the wake-up times of each ONU. In addition to the downstream data, the OLT also sends GATE control messages to the ONUs, which contain the slot allocations and the information about sleep and wake-up. However, the ONU functions are much simpler than the OLT. The ONUs transmit upstream data and a REPORT control message with a bandwidth request to the OLT. In addition, the ONUs comply with the GATE control message to sleep or wake up.

5. Simulation Results

For simulation, we name the proposed ESMT algorithm ESMT-N if no energy constraints are involved (here, N denotes no energy consumption constraints). That is, ESMT-N determines that ONU enters a deep sleep state only when . In addition, we take the previous algorithm in [18, 22] as comparison algorithm which focuses on energy saving oriented to multicast traffic in an EPON without upstream traffic. In simulation, we investigate the performance of our proposed algorithm and the above two algorithms on an EPON platform implemented by the OPNET simulator. On the EPON platform, the network topology is composed of an OLT, a splitter, and four ONUs about 25 km away from the OLT. In each ONU, there are two queues to buffer the upstream real-time and non-real-time packets, respectively. Similarly, in the OLT, there are two queues to buffer the downstream real-time and non-real-time packets for each ONU. We assume that these queues have infinite length and the packets queued in them are accessed by First In First Out (FIFO). The packet generation for both upstream and downstream cases follows a Poisson distribution. The average packet size is 500 bytes. For simplicity, we assume that each multicast group contains two ONUs at most, and the number of ONUs is generated randomly. is computed as 0.48125 mW [4]. For other parameters, refer to [3].

5.1. Comparison of ESMT and ESMT-N [18, 22], Given an Upstream Packet Arrival Rate

Given an upstream packet arrival rate (10 packets/ms), we compare the four algorithms in terms of the average number of wake-up times (ANWT), average energy consumption (AEC) of each ONU, and the average packet delay of each downstream packet.

In Table 1, we show the four algorithms with an increasing downstream packet arrival rate. We can see that the ANWT of [18] is the smallest due to listen state, followed by ESMT and ESMT-N, and that of [22] is the largest. The result demonstrates the effectiveness of our proposed algorithm in terms of the reduction of wake-up times for each ONU. ANWT in [22] is constant because of the fixed sleep duration. Furthermore, the ANWT of ESMT is smaller than ESMT-N due to the additional constraints allowing ONUs to enter deep sleep state.

To highlight the advantage of ESMT over other algorithms, we show the AEC of the four algorithms in Figure 9. It is obvious that ESMT consumes the least energy. This is because the ANWT contributes a large part to the ONU energy consumption for the four algorithms. In addition to the ANWT, the AEC is also more dependent on the energy consumption of the ONUs’ transmitters and receivers. Despite the presence of larger ANWT, ESMT do not include listen state which is more energy saving than deep sleep state, so ESMT is more energy efficient.

In Figure 10, we show the average packet delay of [18] and real-time delay (RD) and non-real-time delay (NRD) of ESMT, ESMT-N, and [22], as [18] does not consider the real-time performance of packets. Due to independent sleep mode in our proposed algorithm and [18], the packet delay outperforms that in [22].

From Table 1 and Figures 9 and 10, we can see that the proposed ESMT and ESMT-N algorithms outperform that in [22] in terms of wake-up times, energy consumption, and traffic delay. However, the cost is increased network overhead due to more GATE control messages. Nevertheless, the proposed algorithm is still more energy efficient. In other words, the effect produced by the increased network overhead is so little that it can be ignored in this paper. On the other hand, the proposed algorithm consumes less energy than that in [18] at the cost of packet delay, while they still meet QoS of the real-time packets.

5.2. Comparison of ESMT and ESMT-N, Given Two Different Upstream Packet Arrival Rates

In this subsection, we compare ESMT, ESMT-N, and [18] given two different upstream packet arrival rates (i.e., 10 packets/ms and 100 packets/ms) to study the influence of the upstream traffic load on the three algorithms. Since upstream traffic is not involved in [22], we do not consider it in this subsection. In addition to the ANWT and AEC, we also compare the AEC of the transmitter (AEC-T) and of the receiver (AEC-R) of the ONUs among ESMT, ESMT-N, and [18]. For simplicity, we define scenario 1 as having an upstream packet arrival rate of 10 packets/ms and scenario 2 as having an upstream packet arrival rate of 100 packets/ms.

Figure 11 shows the ANWT of the three algorithms under two different upstream packet arrival rates. First, it is obvious that the ANWT of ESMT is smaller than ESMT-N, and the reason is the same as in Table 1. In addition, when the downstream packet arrival rate is 20–80 packets/ms, the ANWTs of ESMT and ESMT-N for scenario 2 are smaller than those for scenario 1. Thus, the increase of upstream traffic load does not increase the ONU wake-up times; however, the ANWT decreases conversely. The reason is that the increased upstream traffic load increases the probability of an upstream transmission stretching across two downstream transmission slots, where the ONU does not need to enter a deep sleep state. On the other hand, when the downstream packet arrival rate is 100 packets/ms, the ANWTs of ESMT and ESMT-N in scenario 2 are larger than those in scenario 1. Here, the downstream traffic load is so large that it greatly reduces the probability of an upstream transmission stretching across two downstream transmission slots; therefore, the ONU wake-up times increase accordingly. In addition, compared with [18], the ANWTs of our proposed algorithm are much larger and the result is the same as in Table 1.

Figures 12 and 13, respectively, show the comparison of the AEC-T and AEC-R of each ONU under the two different upstream packet arrival rates. We can see that the changing trend of the AEC-T and AEC-R under the two scenarios is just the opposite of that of the ANWT in Figure 11. This is because the smaller the ONU wake-up times are, the longer the ONUs remain out of the deep sleep state where the transmitter (in Figure 12) and receiver (in Figure 13) are turned off. Therefore, the transmitter/receiver consumes more energy. The results also demonstrate that more energy is consumed in [18] in comparison to our proposed algorithm. It is worth noting that the energy consumption of transmitter and receiver increases with the traffic load in [18].

Figure 14 compares the AEC of the ONUs for ESMT, ESMT-N, and [18] under the two scenarios. The AEC of the ONUs is composed of AEC-T, AEC-R, and the energy consumed by the common modules and the wake-up processes. As is shown in Figure 14, under the same scenario, ESMT is more energy efficient than ESMT-N. Although the upstream traffic load of scenario 2 is larger than that of scenario 1, the AEC of both ESMT and ESMT-N does not always increase. Thus, it can be seen that increasing the upstream traffic load has a certain influence on the energy consumption of ESMT and ESMT-N but does not always increase them. The AEC of [18] is still larger than our proposed algorithm, which is in accord with the results in Section 5.1, and the increasing trend is in keeping with the energy consumption of the transmitter and receiver.

6. Conclusion

This paper proposed an ESMT algorithm to optimize the energy consumption of EPONs oriented to multicast traffic. In ESMT, the independent sleep mode and the deep sleep mode were combined with the characteristics of multicast traffic. Both upstream and downstream traffic loads were considered when designing the energy-saving algorithm. In addition, the ONU wake-up times were also considered. The simulation was run from three different perspectives, and it was demonstrated that our proposed ESMT is effective in terms of both energy consumption and packet delay.

Notations

The power consumed by the transmitter of each ONU
The power consumed by the receiver of each ONU
The power consumed by the common modules of each ONU
The power consumed by the wake-up process of each ONU
The transmission period of ’s upstream traffic
The receiving period of ’s downstream traffic
The deep sleep period of ’s downstream traffic
The number of wake-up times of
The running time of .

Conflict of Interests

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China (61501308), Liaoning BaiQianWan Talents Program (no. 2014921042), Key Laboratory Open Foundation of Shenyang Ligong University, and the Special Foundation of Doctoral Program of Shenyang Ligong University.