Abstract

This paper investigates the time lags effects of innovation input on output in the national innovation systems (NIS). Firstly, we analyze the intrinsic properties of China’s NIS based on the conceptual framework and the causal loop diagram. Secondly, we construct a time lags distribution calculation model and employ the main innovation input indicators to measure the specific characteristics of time lags effects of innovation input on output for China in the period 2000–2012. The results indicate that there are considerable time lags in the four major feedback loops, and the distribution of time lags is various according to the characteristics of innovation input and influencing factors in the internal transformation. The trends of time lags from R&D personnel and industry-academia-research collaboration show steady growth, and the trends of time lags from R&D expenditure and government’s macro-control take inverted U-shaped pattern. Finally, we combine the estimation results of lags distribution calculation model with the analysis results of the causal loop diagram (including positive loops and negative loops) to provide some innovation policy suggestions. This study provides important implications for our understanding of the long-term complex lags effects of innovation input on output as well as for policy-makers designing and implementing the innovation strategies.

1. Introduction

The concept of the national innovation systems (NIS) emerged in the late 1980s [1] and further developed in the years that followed [24]. NIS is a crucial driver of economic growth and competitiveness [5] and enjoys widespread attention in both academic and policy-making contexts [6]. Although there is no single definition of NIS, a good synthetic summary of the prevailing definitions is given as follows: “… the system of interacting private and public firms (either large or small), universities and government agencies, aiming at the production of science and technology (S&T) within national borders. Interaction among those units may be technical, commercial, legal, social and financial, in as much as the goal of the interaction is the development, protection, financing or regulation of new S&T” [7].

Many studies have been carried out over the last two decades. These studies can be listed in three categories. The first category focuses on the structure of NIS [1, 2, 8]. These studies identify and describe the main components of the system and investigate their mutual interactions as well as their relationships with the social and institutional environment in which the system is embedded [9]. The aims of these studies are to recognize the main determinants influencing the production, use, and dissemination of new technologies [10]. The studies of the second category concentrate on various functions of NIS. The functions are the activities that contribute to the goal of innovation systems (both positive and negative) [11], and they are often used in relation to particular institutions [12] or to the system as a whole [2, 13, 14]. The third category of studies is related with the effectiveness of NIS which mainly cover the efficiency evaluation of input-output system [15, 16] and international comparisons of innovation capability [1719]. Much attention has so far been given to the efficiency evaluation of NIS, because this study helps both to identify the best innovation practitioners for benchmarking and to shed light on ways to improve efficiency by distinguishing weakness [16]. While a substantial amount of research has been devoted to the investigation of the innovation input-output ratio and the effect of public policy intervention on the innovation efficiency, the analysis of the dynamics of innovation input-output in NIS was neglected over time.

Some researchers [2022] have found that the research process takes time and the current research and development (R&D) may not have an effect on measured productivity until several years have elapsed. And the NIS comprises a network of interacting policies, institutions, and organizations [23] whose innovative performance depends not only on how the individual component part performs in isolation, but on the quality of interaction and cooperation between the various elements [24, 25], which is subject to dynamic processes [26]. Therefore, the transform process from innovation input to innovation output may be delayed because of the interaction and influence of various factors in NIS. The innovation input cannot generate innovation output and profits immediately or there is the multiple-period influence of one-time innovation input on output.

Because of the significance of the time lags in NIS and the lack of relevant research, this paper focuses on the time lags effects of innovation input on output in NIS. We construct a time lags distribution calculation model to investigate the structure and extent of time lags of major innovation inputs on output, and then we discuss the time lags effects regarding the national innovative factors in the long-run evolution of NIS from both quantitative and qualitative point of view.

The rest of the paper is as follows. Section 2 provides an overview of the study area: China’s national innovation system. Section 3 is devoted to describing the methodology and presenting the data for the research study of time lags. The analysis results of time lags effects of innovation input on output for China are reported in Section 4. In Section 5, policy implications drawn from empirical results are discussed. Finally, Section 6 provides a brief summary and potential future research aspects.

2. Overview of China’s NIS

China has experienced three decades of sustained and strong economic growth as its transitions from a centrally planned economy to a free market, and the associated reforms have necessitated changes in the field of S&T. Since the early 1980s, China has promulgated a series of policies for promoting R&D nationwide [27]. And the development process of China’s NIS is reflected by the six important National Science and Technology Conferences (1978, 1995, 1999, 2006, 2012, and 2016) (see Table 1).

Notable progress has been made in the construction and the structural transformation of China’s NIS. Enterprises have been becoming the center in both the funding structure and performing structure of S&T expenditures on a state level; however, the central government remains the leading force in reforming its NIS with Chinese characteristics [28].

2.1. Conceptual Framework of China’s NIS

Figure 1 shows the conceptual framework of China’s NIS. The framework is put forward based on the current situation of China and extensive references to the literatures, such as Porter’s Diamond Model [29], FP&S Model [30], and M&J Model [31]. Our framework contains four subsystems: innovation foundation subsystem, innovation environment subsystem, innovation carrier subsystem, and innovation output subsystem.

2.1.1. Innovation Foundation Subsystem

Innovation foundation subsystem is the most basic element of China’s NIS. It represents the level of national economic and social development as well as the input and contribution of the country to innovation development.

2.1.2. Innovation Environment Subsystem

Innovation environment subsystem is the important guarantee to the effective operation of China’s NIS. Favorable innovation environment not only can effectively allocate the innovation resources but also foster a highly competitive innovation carrier. Innovation environment subsystem involves many elements, such as social environment, policy environment, business environment, and ecological sustainability.

2.1.3. Innovation Carrier Subsystem

Innovation carrier subsystem includes innovation human resources and R&D expenditure. Innovation carriers flow among the various innovation entities. The role is responsible for linking each innovation entity.

2.1.4. Innovation Output Subsystem

Innovation outputs are the results of innovative activities. The quantity and quality of innovation outputs represent the development and implementation dynamics of innovation activities.

2.1.5. Internal Links between the Subsystems

Innovation foundation and innovation environment subsystems are the foundation and guarantee of innovation carrier and innovation output subsystem, whereas innovation output can also strengthen innovation foundation, improve innovation environment, and provide more resources for innovation by promoting economic development. The growth of innovative resources drives the dynamics of scientific, technological, and innovative output. In contrast, the growth of innovation output sustains further innovative investments over time, in turn.

2.2. Causal Loop Diagram of China’s NIS

The literatures suggest system dynamic (SD) as a well-documented methodological approach for studying the dynamic behavior of innovation systems [3234]. And a causal loop diagram is a map of system with all its constituent components and their interactions. A causal loop diagram represents the major feedback mechanisms. These mechanisms are either positive feedback or negative feedback loops. The positive sign (“+”) indicates a reinforcing effect, while the negative sign (“−”) represents a balancing effect.

In order to reach a better understanding of what really takes place inside China’s NIS, we developed a causal loop diagram of China’s NIS (shown in Figure 2) based on conceptual framework. And there are four major feedback loops in the diagram.

2.2.1. Feedback Loop of R&D Personnel

The first loop (Figure 3) shows that an increase of public expenditure devoted to education may lead to the improvement of the education level. In this way, increasing public expenditure on education may stimulate the production of human capital for scientific and technological research and then promote the production of innovation achievement.

2.2.2. Feedback Loop of R&D Expenditure

The second loop (Figure 4) describes that the government fiscal policies and increasing R&D expenditure of enterprises for themselves and for research universities and research institutions effectively promote innovation outputs.

2.2.3. Feedback Loop of Government’s Macro-Control

The third loop (Figure 5) depicts that the scientific and technological advances help drive economic and social development and then strengthen the government’s macro-control. The reform of the S&T management system will be accelerated, and the innovation quantity and innovation quality will be enhanced in turn on the back of the government’s macro-control policies and measures, such as fiscal policies, monetary policies, industrial policies, and the protection and use of intellectual property rights.

2.2.4. Feedback Loop of Industry-Academia-Research Collaboration

Most major technological advances and innovations originate from interactions between industry and the scientific community [35]. Many countries are actively encouraging industry-academia-research collaboration to improve innovation efficiency and therefore enhance wealth creation. It can be seen from the fourth loop (Figure 6) that the collaboration helps promote research, development, and innovation.

China has become the world’s second largest economy in the world today, with the gross domestic product (GDP) of 56.88 trillion yuan (9.3 1 trillion U.S. dollars) in 2013, according to the National Bureau of Statistics (NBS). Nevertheless, today China also faces the challenge of making the transition from sustained to sustainable growth from social, economical, ecological, and environmental points of view. Innovation has been identified as a main engine to fuel the next phase of growth, and China’s “National Guideline on Medium- and Long-Term Program for Science and Technology Development (2006–2020)” reflects its determination to overcome growing problems through S&T and to become a world leader in innovation. Therefore investigating the operation mechanism of China’s NIS is useful for adjusting the development strategies of NIS and for enhancing S&T capabilities.

3. Methodology and Data

To study the questions raised in the earlier discussion we implement a time lags distribution calculation model. Then key indicators and data in the time lags analysis are described.

3.1. Methodology

In this section, we construct a time lags distribution calculation model based on the grey system theory. Grey system theory, which is a mathematical method to deal with the dynamic system, was formulated by Deng [36]. The grey correlation analysis (GCA) is to determine whether the relationship among a series of data is close, based on the degree of similarity among the geometric shapes of the data series’ curves. The closer the distance between the curves becomes, the greater the correlation among the relative data series is. Compared with traditional regression statistical methods, which require more data and time, GCA can determine the relationship between variables based on the input-output data, without assuming a concrete production function. There are many other advantages, such as no special requirement on the sample size and data distribution, a small amount of calculation, and ease of programming. The method has been applied in many fields [37, 38]. Measuring the time lags effects, GCA has been proven to be useful empirically [39, 40]. For example, Wang et al. [39] and Liu and Guo [40] applied GCA in an analysis of the lag effects of R&D expenditure on economic growth in China.

Accordingly, to cope with the multicollinearity of NIS, we applied GCA for estimating the time lags structure of the effect of innovation input on output [41]. The steps are as follows.

3.1.1. Sequence Determination

There are two kinds of lag sequences: reference sequence and comparative sequence. The reference sequence reflects the behavior features of the system. The data sequence consisting of factors which affect the behavior of system is called comparative sequence. The lag sequences are as follows:where is the lag step length. The minimum amount of data in the GCA is 3. In order to ensure the basic condition of the model, we set (). is the reference sequence, and is the comparative sequence.

3.1.2. Data Preprocessing

Data preprocessing should be performed before the calculation of the grey correlation coefficients. Data preprocessing converts the original sequence to a dimensionless sequence. The original sequence is normalized as follows:where and are the normalized element in the reference and comparative sequence, respectively.

3.1.3. Grey Relative Correlation Grade Calculation

The grey relative correlation grades in different lag steps are calculated. The lag step of the maximum grey relative correlation grade, which is taken as the time difference of the reference sequence caused by the comparative sequence, is the final time lag of this comparative sequence.

Using (3), we define the following calculation model of grey relative correlation grade:where is the grey relative correlation grade of ; is the maximum grey relative correlation grade; is the lag step of the maximum grade; if , the reference sequence and the comparative sequence are synchronized; otherwise, if , we take the comparative sequence as the lag indicator of the reference sequence. is the final time lag of this indicator.

3.2. Data

Based on the analysis results in Section 2.2, the four major feedback loops describe the major feedback mechanisms of R&D personnel, R&D expenditure, government’s macro-control, and industry-academia-research collaboration promoting innovation outputs. Accordingly, four input indicators representing the four feedback loops are considered: researchers in R&D (per million people) (RRD), R&D expenditure (% of GDP) (RDE), regulatory quality (RQ), and university-industry research collaboration (UIRC).

RQ is from the Worldwide Governance Indicators (WGI) which are a research dataset summarizing the views on the quality of governance provided by a large number of enterprises, citizens, and expert survey respondents in industrial and developing countries. It reflects perceptions of the ability of the government to formulate and implement sound policies and regulations that permit and promote private sector development. Estimate of governance ranges from approximately −2.5 (weak) to 2.5 (strong) governance performance [42].

UIRC is from the average answer to the survey question of the World Economic Forum: In your country, to what extent do business and universities collaborate on R&D? [1 = do not collaborate at all; 7 = collaborate extensively].

Since patents are often considered an important indicator of technological change and innovation [4346], this paper will take the patent application as the variable for the innovation output indicator. Table 2 provides the indicators, the feedback loops which belong to, their units of measurement, and the sources of data. The detailed indicator values are given in Table 3. The empirical analysis accounts for annual time series data for China from 1996 to 2012.

4. Results and Discussion

The results presented in this section are based on the time lags distribution calculation model (see Section 3.1). Table 4 and Figures 710 provide the estimated lag structures and trends of the key innovation input indicators (RRD, RDE, RQ, and UIRC) on the innovation output indicator (PAR) over time. We will discuss them in greater detail combined with the four major feedback loops of China’s NIS, in which both quantitative and qualitative research strategies can be used to complement each other.

Before the computation of time lags, we test all variables for unit roots. Augmented Dickey Fuller (ADF) [47] test which is one of the most popular tests for a unit root in a time series sample is applied to check the stationary of data. The tests are computed with automatic lag length selection relying on Schwarz’ information criterion (SIC) and with a maximum lag length of 3. The results show that RRD, RDE, RQ, and UIRC are stationary at first difference level and PAR is stationary at second difference level. All computations are performed with the help of EViews.

The estimation results of GCA are presented in Table 4.

4.1. Time Lags from R&D Personnel

The results for RRD are shown in Table 4 (column (I)) and Figure 7. , and the maximum lag step . The stable effects dominate during the period between the first year after start-up and the thirteenth year, suggesting a relatively long innovation accumulation period. The positive impact of RRD on PAR reaches peak around 14 years after entry.

Promoting the education sector normally entails increasing public expenditure on education [48], as well as further integration of educational resources. Effective integration of educational resources is a long-term and systematic project. Personnel training process takes time and graduates may not contribute to innovation in short time. Researchers in R&D need to undergo a process of germination, growth, and maturity. That is the time lag from the process of training R&D personnel.

From the negative causal effect (in Figure 3) point of view, the growth of TNSTP may lead to reduction of the demand for TNSTP, decrease the number of personnel training and the personnel introduction, and then slow down the pace of innovation output.

4.2. Time Lags from R&D Expenditure

For RDE, , and the maximum lag step (Table 4, column (II), and Figure 8). The peak of the positive impact of RDE on PAR is reached about 12 years after entry and then fades away.

The main role of R&D expenditure is to provide cash flow for the whole innovation process. R&D expenditure provides an important guarantee for the successful operation of the innovation process, and in the meantime, the innovation process influences circulation of R&D expenditure in turn. It can be said that the time lags from innovation process can cause a certain delay in the use of R&D expenditure. The time lags in the innovation process mainly come from two aspects. On the one hand, the innovation process, which integrates the existing technology, R&D expenditure, personnel, and so forth into a holistic dynamic consideration, needs to go through a period of time. This period is called the normal time lags. On the other hand, whether or not such an innovation process can succeed is uncertain because there are a variety of risk factors, including decision-making risk, personnel risk, financial risk, technology risk, and industry risk. Any risks are likely to have a negative impact on the development process, to extend the innovation process, and even to ruin the entire research project. The delay caused by the risk factors is called the nonnormal time lags. It is the two aspects of delay that explain the time lags of the innovation process from input to output as well as the use of R&D expenditure.

4.3. Time Lags from Government’s Macro-Control

Table 4 (column (III)) and Figure 9 present the estimated lag structures of RQ, , and the maximum lag step . According to the results, the magnitude of such effects over time takes an inverted U-shaped pattern with a peak for the start-up activity from ten years earlier, while no significant decrease is identified after the peak.

There are still the time lags in the process of government policy formulation, development, and implementation. The reason is twofold. Firstly, a fundamental role of the governments is the configuration of regulations and legislation via which a nation’s innovative activity is influenced [29]. The formulation of government innovation policy, which takes a number of steps, may require a lot of time. Secondly, there may be considerable time lags involved when governments launch numerous innovation policies to improve the environment for innovation. The transformation needed in the lagging innovation environments is hampered by the limitations and disadvantages of conventional and traditional system of organization.

As can be seen from the negative causal effect in Figure 5, irrational government regulations may enhance the market risk of enterprises innovation. And the end result may be the slowdown of innovation output.

4.4. Time Lags from Industry-Academia-Research Collaboration

As shown in Table 4 (column (IV)) and Figure 10, for UIRC, , and the maximum lag step . The positive effect of UIRC on PAR has steadily increased since the start-up activity, tending to peak around the fourteenth year. With the accumulation of data, UIRC will lead to more PAR in the longer run.

The analyses of the time lags for industry-academia-research collaboration are presented on the two aspects: process and factors (barriers). On the one hand, the process of industry-academia-research collaboration is a knowledge transfer process in a certain sense. Szulanski [49] conceptualized the phases of knowledge transfer using the milestones. Nonaka and Takeuchi [50] proposed the SECI (Socialization, Externalization, Combination, Internalization) model to present the spiraling knowledge processes of interaction between explicit knowledge and tacit knowledge. The knowledge transfer process of industry-academia-research collaboration is nonlinear paths of knowledge acquisition, knowledge sharing, and knowledge application by and among the various partners in the innovation process. In addition, prior to setting up a research collaboration, there is a “step by step” process of how to include a consultancy agreement [51]. In conclusion, the whole process of industry-academia-research collaboration may take a while or a long time. Thus, there may be a long delay at this stage. On the other hand, a variety of barriers and constraints may adversely affect industry-academia-research relationships and reduce the collaboration efficiency. Industry and universities might be well motivated to collaborate, but there are many barriers to such involvement and the considerable potential benefits are often not realized in practice [52]. Published literature suggests that the most significant barriers are in general related to differences or conflicts in the missions and objectives, institutional differences (different cultures and structures), different research interests of companies and universities, ineffective communication, restrictions in information dissemination due to confidentiality issues, and conflicts over intellectual property rights [53, 54].

As shown in Figure 6, the negative causal effect is that unreasonable innovation services policy may increase the costs of innovation collaboration, reduce the intensity of industry-academia-research cooperation, and thereby slow down the innovation output.

5. Implications for Policy

The time lags effects regarding the national innovative factors in the long-run evolution of China’s NIS have been identified and analyzed. The analyses and discussions presented above suggest important implications for policy-makers to ensue.

5.1. Personnel Policy

The results shown in Figure 7 indicate that researchers need to go through the technical accumulation cycle of 14 years in order to maximize their usefulness to innovation output. How can the accumulation cycle be shortened? These are two important aspects to improve efficiency of input-output.

From the point of view of the positive causal effect: First, the basic problem is to increase the quality and quantity of R&D personnel by implementing higher educational reforms. Second, make the training mechanism of R&D personnel better and strive to shorten the period from novice to expert. The training mechanism reform contains the following: Develop strong quality standards and procedures for training the new R&D personnel; propose a personalized training system (including the training content and the training means) for different groups; carry out training evaluation to achieve the maximum training effectiveness; pay attention to the training of the organization culture and core values. Third, the direct introduction of foreign, highly skilled “talent” will play a significant role in shortening the technical accumulation cycle in China. Creating the appropriate working and living environments is an effective means of attracting and retaining foreign talent. And another means is to establish China as an education hub, attracting foreign talent that may remain in the country, bringing in foreign revenue, and helping to draw more multinational R&D companies to China.

When considering the negative causal effect shown in Figure 3, the government should balance the demand for TNSTP to achieve a balanced development of training and introduction of S&T personnel.

5.2. R&D Funding Policy

From Figure 8, the positive effect of RDE on PAR gradually increases within a 12-year period that continues for a relatively long time. There is a close connection between R&D process and efficiency of R&D funding utilization (see Section 4.2). Therefore, optimizing R&D process as well as improving R&D performance is an effective method of reducing the normal time lags, decreasing the nonnormal time lags, and then shortening the time lags of R&D funding input on output. These methods include the following: Eliminate the boundaries between departments in order to achieve the goal-oriented team without boundary; build and optimize market-oriented R&D process; focus on risk management, including risk identification, risk analysis, risk planning, risk tracking, and risk response. In addition, increase the quantity of innovation input, at the same time, paying more attention to the sustainability of innovation input. A certain amount of innovation input may have no significant effect on innovative performance growth immediately. Safeguard the persistence of innovative investment, especially in the rising phase of the contribution rate of RDE on PAR.

5.3. Government Innovation Policy

The results (Figure 9) may be attributed to the relatively long time lags that are required for the national macro-control policies to become visible, the durability of public policy, and the enhancement of the market risk faced by enterprises innovation.

First, improve the efficiency of formulation, development, and implementation of government innovation policy; strengthen the macro-coordination of innovation policies and regulations across China and deal with the central-local government relations, internal relations of science departments, and science-nonscience department relations well to avoid repeated policy deployments. Second, get rid of the ideological obstacles and institutional barriers and dredge the transmission channel of innovation policy, for speeding up playing the effect of innovation policies and regulations. Third, play the continuing effect of policies and regulations on innovation output. Coordinate the relationship between policies and regulations for innovation and other types of policies and laws, seek common ground while reserving differences, and then explore the breadth and depth of the coordination. Furthermore, perfect the relevant laws and regulations for keeping up with the constant change of the objective conditions and legal environment. Finally, in order to accelerate innovation output, the government should make more reasonable innovation-related regulations, which help to improve the market sophistication and reduce the market risk of enterprises innovation.

5.4. Industry-Academia-Research Collaboration Policy

There is a 14-year maximum lag structure and contribution rate of UIRC on PAR varies at a slow speed (Figure 10), which indicates that the run-in period of collaboration persists longer.

In order to shorten the run-in period, we can implement the following strategies: Firstly, improve the efficiency of knowledge transfer from the knowledge acquisition, knowledge sharing, and knowledge application [49]; secondly, establish consulting agency [51], which provides high-standard consultancy services (e.g., consultancy agreement, technology transfer, intellectual property, and investment and financing decisions); thirdly, as shown in Figure 10, the gradual increasing effect of UIRC on PAR suggests that, after a period of run-in, both parties may easily boost their credibility, visibility, and reputation; in other words, they can work well with each other. Therefore, an option to help accelerate the development of business academia interactions is to plan more long-term value-added collaborations based on the existing collaborations. Strive to develop deep-level relationship and invent the novel collaboration pattern models; finally, reduce the negative causal effect of unreasonable industry-academia-research collaboration policy on innovation output (Figure 6). Government should implement innovation services policy in line with their actual situation based on the in-depth investigation. Practicable innovation services policy can contribute to the costs reduction of innovation collaboration and the enhancement of industry-academia-research cooperation intensity.

6. Conclusion

The existing literatures about the NIS have made considerable theoretical progress, but they fail to empirically investigate the time lags of internal transformation process of innovation input-output and the multiple-period influence of one-time innovation input on output based on systems thinking.

In order to explore this new direction of research, firstly, our study has presented an overview of China’s NIS based on the conceptual framework and the causal loop diagram. Secondly, we have made use of GCA method measuring the structure and extent of time lags of four major innovation inputs (RRD, RDE, RQ, and UIRC), which represent the four major feedback loops (R&D personnel, R&D expenditure, government’s macro-control, and industry-academia-research collaboration), respectively, on PAR for China in the period 1996–2012. The distribution of time lags of different innovation input on output varies according to the characteristics, process, and factors. Thirdly, we have provided some policy suggestions taking into account simultaneously shortening the significant time delays which result in accelerating the exponential growth, as well as the side effects of interventions in shortening the delays on dominating negative loops.

On the whole, this study contributes to theory-building and policy-making in the following ways. For theory-building, the investigation of time lags properties of the long-run evolution of NIS can complement and inform the NIS research and reveal the long-term complex relationships between the innovation input and output to further understand the operation mechanism of NIS. So it will be an important research direction. For policy-making, identifying the time lags effects of NIS helps policy-makers better understand why an innovation input policy may have no significant effect on innovative performance growth immediately. Therefore, public policy should safeguard the persistence of innovative investment and avoid everything that could disturb the persistent transformation process.

The limitations of our study are particularly interesting and worthy of further exploration in future research. First, NIS is a dynamic system such that its evolution is driven by a complex set of factors. Any given change in one of the factors composing the NIS has a set of direct and indirect effects on the process of transforming innovation input into output. Although this study explicitly analyzes the time lags effects of the four major input-output innovation processes, it is still a simplified structure within which some practical complexities are not taken into account. In our future research, we will attempt to investigate the detailed structure for a better simulation of NIS. Second, this paper takes China as an example to investigate the time lags; this could actually be useful for many countries that are in similar transitions from command economy to market economy. But whether our models are appropriate for all other countries needs further verification. Hence, a direction for our future research is to analyze the time lags of more countries by using our hybrid analytical method and to conduct comparative studies between the selected countries.

Competing Interests

The authors confirm that this paper content has no conflict of interests.

Acknowledgments

This research is supported by Special Study Project on G20 Innovation Topics, The Ministry of Science and Technology, China; the Fundamental Research Funds 576 for the Central Universities (Grant no. 2572015BX12); and the Educational Science and Planning Project of Heilongjiang Province (Grant no. GJD1214015).