Abstract

The COVID-19 pandemic has plunged the housing affordability and living sanitary conditions of China's low-income groups into a crisis. Public rental housing has made a significant contribution to alleviating this crisis, yet its exit problem has become more acute due to the impact of the epidemic. In order to explore how to effectively play the role of multi-governance in public rental housing exit, this paper adopts game theory to analyze the evolution of the behavioral strategies of government departments and public rental housing tenants in public rental housing exit by combining four influencing parameters: the epidemic severity, the epidemic influence coefficient, the public participation degree, and the public reporting fairness. In response to exploring the conditions for the effective functioning of the multi-governance model under the impact of the epidemic, the influences of parameter change on the stability of the system are analyzed. The results show that when the COVID-19 epidemic impact coefficient is high, government departments need to reduce their own supervision costs to achieve better multi-governance effects; when public reports are distorted, the government departments’ own supervision capabilities determine the impact of public participation on governance effects.

1. Introduction

The sudden arrival of the COVID-19 epidemic has had an extremely negative impact on urban governance and public health in China. The most immediate impact has been reflected in the housing problem of low-income groups. On the one hand, the epidemic has directly led to a decline in the housing affordability of low-income groups [1]. Affected by the epidemic, urban housing costs [2] and unemployment rates have been continually on the rise. In 2020, China’s new employment decreased by 4.7 million, and the loss rate of new employment reached 29.26% [3]. Due to the overlap with the Spring Festival, it is difficult for low-income groups to return to work, which makes the structural contradiction between the supply and demand of urban and rural labor more prominent, causing the unemployment of low-income groups to reach a new high [4]. The gradual increase in unemployment has brought a greater economic shock to low-income groups, under which the budget for housing rents has exceeded their ability to pay. Housing affordability has become a major problem during the epidemic. On the other hand, the COVID-19 epidemic has threatened low-income groups’ living sanitary conditions. During the epidemic, the majority of low-income groups were renting peasant-built houses in urban-rural areas. These dwellings have neither fixed property management services nor sufficient environmental health management. After the onset of the epidemic, low-income groups either fled the city and returned to their relatively safe rural home or lived in the city in fear, increasing the risk of public health. Poor housing affordability and living sanitary conditions caused by the epidemic have adversely affected the mental health of low-income groups [5]. Cheap and affordable housing with certain supporting facilities, public health services, and government supervision has received unprecedented attention and demand. The General Office of the State Council also issued the “Opinions on Accelerating the Development of Affordable Rental Housing” in June 2021, emphasizing the expansion of the supply of affordable rental housing and the resolution to solve the prominent housing problems in big cities.

Public rental housing eases the housing and health pressure of China’s low-income groups during the COVID-19 epidemic [6]. On the one hand, the rent of public rental housing does not increase but falls, which greatly alleviates the economic pressure of low-income tenants. On the other hand, under the unified management of the government, public rental housing has basically guaranteed the tenants’ living sanitary conditions and public health management during the severe period of the epidemic. However, due to the imperfect review mechanism and exit mechanism, there is a phenomenon of “only enter but not exit” in public rental housing in many places of China. Especially under the negative impact of the COVID-19 epidemic, the problem of public rental housing exit has become more serious. Firstly, the epidemic has raised the cost of government supervision and enforcement so that the number of delay exit behaviors has increased. Furthermore, some tenants who should have quit public rental housing were affected by the epidemic, resulting in a decline in their housing affordability and choosing not to exit. As public rental housing has a social security function, if government departments fail to effectively address the problem of reneging by some tenants, it is very likely that public rental housing will lose its fairness and fall into the “public housing” dilemma again.

When studying the issue of urban governance, most scholars emphasized the importance of subject participation to governance effects and proposed that the construction of a multiple collaborative co-governance mechanism should be accelerated to realize the transition from government “unified” supervision to social “diversified” governance. Wood et al. (1991) [7] elaborated on the theory of collaborative governance and proposed a governance model in which all stakeholders make joint decisions in the face of public affairs. The theory of collaborative governance has played an important role in the study of various governance issues. Guoxing Zhang proposed that a collaborative governance mechanism composed of government departments, enterprises, social organizations, and other subjects should be established to better protect the environment [8]. Junqing Li aimed to solve the problem of poverty in ethnic areas by building a multi-subject collaborative governance system in the field of poverty alleviation poverty issues [9]. Yulin Zhao pointed out that the multi-coordinated governance model should be led by the government, with the participation of relevant social entities, so as to achieve the effect of co-governance in Internet governance issues [10]. The multi-governance model is of great significance to flexible governance of COVID-19 epidemic [11]. Under the multi-governance model of collaborative governance, public participation is essential. It not only enables urban planners and policy makers to consider genuine public needs [12, 13], but also supports the sustainable development of urban communities [14]. By combining the organizational structure and mechanism of public rental housing exit from the multi-governance model, it can be seen that the behavioral strategies of government departments and tenants influence each other, and there is a strong game relationship between them. At present, most scholars choose to adopt the evolutionary game method to explore effective measures to improve the governance effect when studying the regulatory issues of government departments. Different from the existing literature, we consider the roles of key stakeholders in the process of public rental housing exit, analyze the relationship between government and public rental housing tenants, explore the causes of the slow exit of tenants, and reveal the evolutionary mechanism of public rental housing exit. On the basis of theoretical analysis, the complexity scientific method is introduced and the evolutionary game model of public rental housing exit is established in order to study how to give full play to the best role of multi-governance in public rental housing exit governance. The innovation of this study mainly lies in two aspects: the first is to explore the effect of the COVID-19 epidemic and the public on public rental housing exit by considering two related parameters, the epidemic impact coefficient and the public reporting fairness; the second is to investigate the role of the government in the promotion of public rental housing exit through analyzing the dynamic decision-making process of both the government and tenants. Based on the research conclusions, this paper intends to provide reference for the government departments to give full play to the role of multi-governance in public rental housing exit. The framework of the study is shown in Figure 1.

2. Literature Review

2.1. Current Status of Public Rental Housing Exit Governance

The issue of public rental housing exit has existed for a long time in China's urban governance. Scholars have done a lot of researches on related issues, which mainly focused on two aspects: tenant exit influencing factors and exit governance mechanisms.

The behavioral decisions of tenants to exit public rental housing are affected by many factors. A large body of literature showed that family income, education level, family size, and understanding of public rental housing policies are the main internal factors [1517]. The probability of government departments’ effective review is the key external factor [18]. In addition, environmental health has gradually become an important factor influencing tenants’ decisions due to the COVID-19 epidemic [19]. However, few scholars have studied the impact of the epidemic on public rental housing exit.

Some scholars have conducted researches on typical models of public rental housing community governance in various countries, such as the American model [20], the Singapore HDB community management model, the Chongqing model [21], the Shanghai model [22] and the Ankang model [23], discussing the necessity for the establishment of dedicated agency, semiofficial people's associations, and civil charity organizations to participate in public rental community governance. The polycentric and cooperative governance model is the common point of most studies, arguing that the government should develop incentive mechanisms [24] and encourage residents to participate in the decision-making process of community management. It is generally accepted that the public participation in public rental housing exit is of great significance to exit governance. However, most studies focused on public participation in public rental housing construction [25] and ignored public’s role in monitoring and reporting when tenants exit public rental housing.

2.2. Application of Evolutionary Game in Public Rental Housing Exit

Evolutionary game has considerable unrealized potential for modeling substantive economic issues. Smith (1986) [26] described the concepts of evolutionary game theory in the context of a simple game. Friedman (1998) [27] exposited the specification of evolutionary game models and classified the possible asymptotic behavior for one- and two-dimensional models. Evolutionary game theory is the application of population dynamical methods to game theory [28], which has strong modeling capabilities for explaining past evolution and predicting the future one [29]. As an analytical tool, evolutionary game theory has been applied to various fields of research studies to explain social phenomena [3033], including construction and demolition waste supply chain [34], environmental governance [35, 36], and public health emergency management [37, 38]. In recent years, the application of the evolutionary game theory in the field of public rental housing exit has been increasing. For example, Baolong (2016) [39] pointed out the two key factors in exit mechanism of public rental housing by constructing game model of the public rental housing authorities and public rental housing tenants. Hui et al. (2017) [40] constructed an evolutionary game model between government departments and tenants based on evolutionary game theory, determined the best time and proportion for government departments to sue delay exit tenants. Liu (2020) [41] provided effective suggestions for the establishment of the flexible exit mechanism of public rental housing in Urumqi through using logistic model to analyze the game evolution of the lessee′s rent-seeking behavior.

As can be seen from the abovementioned literature review, the evolutionary game model has been applied to the study of public rental housing exit, but the related research mainly focused on the government supervision and paid less attention to the epidemic impact and public participation. Prior studies mainly considered the relationship between government and tenants. However, few scholars analyzed public rental housing exit from the perspective of multi-stakeholder relationship. It is difficult to see the full application of evolutionary game in public rental housing exit.

3. Model Construction and Analysis

3.1. Parameter Assumptions and Modeling

This paper models the issue of multi-governance for public rental housing exit under the impact of the COVID-19 epidemic, which involves two participants: government departments and the tenants who should have quit public rental housing but did not. During the epidemic period, the government’s strategic actions include “positive supervision” and “negative supervision.” The strategic behaviors of the tenants include “active exit” and “delay exit.”

In order to clarify the relationship between the two participants in the game, this paper makes the following rational assumptions on the premise of fully integrating the public rental housing reality(1)In the process of multi-governance during the epidemic, the probability of government agencies adopting positive supervision is x (0 < x ≤ 1). At this time, the government will spend a fixed supervision cost and incur additional supervision costs caused by the epidemic; i.e., . Simultaneously, the government can obtain benefits , which means the improvement of the total social welfare and government credibility.(2)The probability that the government adopts negative supervision is . At this point, although the government does not have to pay supervision cost, it will suffer a loss if the tenants’ irregularities are exposed, indicating the loss of credibility, suffered as a result of government inaction.(3)In the process of multi-governance, the probability of the tenants taking active exit strategy is y (0 < y ≤ 1). At this point, the government will give a certain housing subsidy to alleviate the financial pressure of the tenants, while the active exit of the tenants will raise the total social welfare.(4)The tenants’ delayed exit probability is 1-y. At this time, it will save a certain amount of housing cost but reduce the total social welfare. Moreover, when government departments investigate and deal with delay exit tenants, they will be fined.(5)When the public perceives that the tenants have violated the regulations, they can report this to the government departments, at which point the probability that the tenants will be investigated and punished by the government is . Nevertheless, considering that the information reported by the public may not be completely accurate, government departments will suffer losses and tenants will suffer losses when the public incorrectly reports tenants for violations.(6)The parameters of the game model: Under each game strategy combination, according to the cost, income, and loss of the government and tenants, the parameter settings are shown in Table 1.

3.2. Model Analysis

According to the model assumptions and the game parameters shown in Table 1, the benefits and payments of the two-party game subjects, the government and the tenant, are shown in Table 2.

Through the above game payment matrix, it can be concluded that the expected benefits of government departments under different behavior strategies are as follows:

The average expected revenue of the government is as follows:

In the same way, the expected benefits of tenants under different behavior strategies are as follows:

The average expected revenue of the tenants is:

It can be obtained that the dynamic equations of replication between government departments and tenants are as follows:

Let dx/dt = 0 and ; the five local equilibrium points of the evolutionary game system can be obtained as (0, 0), (0, 1), (1, 0), (1, 1), (, ).

3.3. Equilibrium and Stability Analysis

In this paper, the local stability of the equilibrium point of the evolutionary game system is judged by the Jacobian matrix composed of the partial derivatives of the dynamic equations.

Through calculation, the value of the game system at the five local equilibrium points can be obtained, as shown in Table 3.

If the Jacobian matrix determinant value is greater than zero and the determinant trace value is less than zero, then the evolutionary game system has an evolutionary stability strategy (ESS). Since it is necessary to comprehensively consider the influence of the four parameters p, , q, and on the stable evolution state of the system and it is difficult to determine the threshold space corresponding to each parameter, without affecting the purpose of this study, we select two reference values and as the judgment intermediary for analyzing the impact of parameter value range change on system stability, which can be divided into the following five situations:(1)When , the ESS is (0, 0). This means that after continuous selection and evolution, the system finally forms a stable state in which both the government's supervision probability and the tenant's exit probability tend to 0, which indicates that low subsidy and penalties are unfavorable for tenants to exit public rental housing.(2)When , the ESS is (0, 1). It can be seen that after continuous selection and evolution, the system finally forms a stable state in which the government's supervision probability tends to 0 and the tenant's exit probability tends to 1, which indicates that high epidemic impact is unfavorable for the government to supervise.(3)When , the ESS is (1, 0). We can see that after continuous selection and evolution, the system finally forms a stable state in which the government's supervision probability tends to 1 and the tenant's exit probability tends to 0, which indicates that low epidemic impact and public participation are favorable for government to supervise.(4)When , the ESS is (1, 1). This means that after continuous selection and evolution, the system finally forms a stable state in which both the government's supervision probability and the tenant's exit probability tend to 1, which indicates that low epidemic impact and high penalties are favorable for tenants to exit public rental housing.(5)When , there is no ESS.

3.4. Model Construction and Analysis

According to the above stability analysis, it is observed that the evolution of the stable state of the system depends on the relative magnitude relationship between the influence parameters and the reference values, where represents the difference between the benefits gained and the fixed cost required when the government adopts an positive supervision behavior strategy, i.e., the maximum benefits gained without considering the epidemic and public participation in the governance process; represents the additional benefits that the tenants can obtain if they adopt the delayed exit strategy, i.e., the opportunity benefits that can be obtained without considering the epidemic and public participation. From this, it can be seen that these two values are stable in the actual environment of public rental housing exit governance, independently of the epidemic and public participation. Therefore, the variation of the system stability equilibrium point within the range of different parameter values can be derived from a comprehensive analysis of the epidemic severity p, the epidemic impact coefficient ω, the public participation degree q, and the public reporting fairness a.

3.5. Special Case Analysis

Case 1. , .
The epidemic impact coefficient  = 1 indicates that the epidemic has a great restriction on the supervision ability of the local government, causing a higher cost when the government adopts positive supervision strategy. The public reporting fairness  = 1 indicates that the public has high information recognition ability; that is, they can better identify whether the tenants delay exit so as to avoid additional losses to the government and tenants.
According to the calculation results of the model, in order to make the evolution state of the system toward (1, 1) stable, the relevant influencing parameters need to meet the following conditions:This shows that when the negative effect of the epidemic is great and the public reporting information is totally accurate, lowering epidemic severity and increasing public participation level can improve the effectiveness of the multi-governance model for public rental housing exit.

Case 2. , .
The epidemic impact coefficient  = 0 means that the epidemic will not affect the government's supervision costs. At this time, in order to make the evolution state of the system toward (1, 1) stable, the relevant influencing parameters need to meet the following conditions:This shows that when the negative effect of the epidemic is minimal, the severity of the epidemic has no impact on the effect of government governance. By comparing the parameter value range of Case 2 with that of Case 1, it can be found that the parameter value range of Case 2 is easier to meet in the context of the system evolution toward point (1, 1); i.e., in order to achieve better governance effect of multi-governance model for public rental housing exit, the parameter setting of Case 2 has more advantages. It shows that the government needs to reduce their own supervision cost to achieve better multi-governance effect when the epidemic has a great impact.

Case 3. , .
The public reporting fairness  = 0 means that public’s reporting information is distorted; i.e., the public will consider the active exit strategy adopted by the tenants as delayed exit strategy and report that to the government departments, which brings additional losses to government departments and tenants.
In order to make the evolution state of the system toward (1, 1) stable, the relevant influencing parameters need to meet the following conditions:In this case, the influence of public participation on the governance effect of government is not stable. Therefore, the discussion is divided into two different situations according to the government departments’ own supervision ability.(1)When , the above conditions are converted to and . At this time, means that the government’s supervision ability is low and the parameter conditions for the system to evolve toward the point (1, 1) need to satisfy ; i.e., in order to achieve better governance effect, it is essential to increase public participation enthusiasm.(2)When , the above conditions are converted to and . At this time, means that the government’s supervision ability is high and the parameter conditions for the system to evolve toward the point (1, 1) need to satisfy ; i.e., in order to achieve better governance effect, it is essential to reduce public participation enthusiasm.Through the above analysis, it is found that the difference in the government's own supervision capacity will affect the role of public participation in the governance process. In the case where public reporting is not completely impartial, government departments relying solely on public reporting information can hardly achieve the optimal effect of multi-governance, which indicates that the government's own supervision capacity plays a key role in strategy formulation and governance effectiveness.
By comparing the parameter value range of Case 3 with that of Case 1, it can be found that since , in practice, in order to make the system evolve toward point (1, 1), the parameter value range of Case 1 is easier to meet, showing the key role of the public reporting fairness in the governance process.

Case 4. , .
In order to make the evolution state of the system toward (1, 1) stable, the relevant influencing parameters need to meet the following conditions:In this case, the influence of public participation enthusiasm on the governance effect of government departments is also unstable, and it is affected by the government's own supervision ability as well. At this time, the influence of the government's own supervision ability on the enthusiasm of public participation is similar to that of Case 3.
By comparing and analyzing the parameter values of the system evolution toward the point (1, 1) under the four cases, we can clearly see that the government departments' multi-governance of public rental housing exit is influenced by many factors, among which the epidemic influence coefficient, the public reporting fairness, and the government's own supervision ability play a key role.

3.6. General Situation Analysis

In order to more clearly understand the impact of both the epidemic and the public related parameters on the effectiveness of governance, this paper analyzes the general condition for the evolution of the system toward the (1, 1) point, i.e., analyzes the situation when each parameter of the system satisfies the following condition:

On the premise of not affecting the conclusion of the analysis, the following part will also take and as reference values to analyze how to set parameters to make the limitation more easily satisfied, i.e., by solving the system of (13) to discuss the effect of parameter changes on the function values, which leads to more general conclusions.

3.6.1. Analysis of Epidemic-Related Parameters

Analysis of the epidemic severity p.

Let the function be derived with respect to the epidemic severity p.

According to (15), it can be seen that the function is an increasing function of the epidemic severity p; i.e., alleviating the epidemic severity p reduces the value of the function . Therefore, to better satisfy the range of parameters for the evolution of the system toward the (1, 1) point, the epidemic spread needs to be suppressed.

Analysis of the epidemic impact coefficient ω.

Let the function be derived with respect to the epidemic impact coefficient p.

According to (16), it can be seen that the function is an increasing function of the epidemic impact ; i.e., reducing the epidemic impact factor reduces the value of the function . Therefore, to better satisfy the range of parameters for the evolution of the system toward the (1, 1) point, it is necessary to reduce the impact of the epidemic on government surveillance capacity.

3.6.2. Analysis of Public-Related Parameter

Analysis of the public participation degree q.

Let the function derive from the level of public participation degree q.

According to (17), it can be seen that the plus or minus characteristic of the equation is unknown, and it is related to the government department's own supervision capacity. Therefore, this paper discusses two different cases according to the government department's own supervision capacity.

When , then . The function is positively correlated with the public participation degree q; increasing the public participation degree q can make the function increase. Therefore, in order to better meet the range of parameters for the evolution of the system toward the point (1, 1), the public participation level needs to be increased.

When , then . The function is negatively correlated with the public participation degree q; decreasing the public participation degree q can make the function increase. Therefore, in order to better meet the range of parameters for the evolution of the system toward the point (1, 1), the public participation level needs to be reduced.

Analysis of the fairness of public reporting.

Let the function derive from the public reporting fairness .

According to (18), it can be seen that the function is an increasing function of the public reporting fairness . Increasing the public reporting fairness can make the function increase. Therefore, in order to better satisfy the range of parameters for the evolution of the system toward the (1, 1) point, the public reporting fairness needs to be increased.

4. Numerical Simulation

By solving and analyzing the model, we can see that the severity and impact coefficients of the epidemic, the level of public participation, and the degree of public report fairness play a complex role in the stable evolutionary equilibrium of the game system. In order to explore the continuous trend of parameters on the evolution of the system toward the (1, 1) point, this section conducts numerical simulations of the outcome conditions of the parameters in the general case of the system and analyzes the impact trends of the parameters on the system by the simulation images.

4.1. Analysis of Strategic Options for Government Action

To investigate the effect of the epidemic severity p and the epidemic impact coefficient on the government strategy selection function , let k = 5, while and are substituted into MATLAB software for numerical simulation; the resulting images are shown in Figure 2.

As can be seen from Figure 1, the function that influences the choice of government strategy is an increasing function of the epidemic severity p and the impact coefficient . The result of simulation coincides with the conclusion that the COVID-19 epidemic expansion will increase homelessness and government costs [42]. Homelessness and the many vulnerabilities that inevitably accompany it are driving COVID-19 outbreaks in shelters and other congregate living situations. Public rental housing is an investment that will curb the spread of COVID-19 and help protect most of low-income groups from future pandemics. Therefore, in order for the system to evolve toward the (1, 1) point, the values of the epidemic severity p and the impact coefficient need to be reduced; i.e., in order to optimize the effect of multi-governance by government departments on public rental housing exit during the epidemic, government departments should strengthen their own supervision capacity to weaken the impact of the epidemic.

4.2. Analysis of the Behavioral Strategy Choices of Tenants

In order to investigate the effect of the public participation degree q and the reporting fairness on the strategy selection function , let = 1,  = 0.5,  = 0.9, and F = 2, while and are substituted into MATLAB software for numerical simulation; the resulting images are shown in Figure 3.

From Figure 3, it can be seen that there are two cases that affect the tenant strategy selection function . When the public reporting fairness is low, the function is negatively correlated with the public participation degree q and the reporting fairness ; when the public reporting fairness is high, the function is positively correlated with the public participation degree q and the reporting fairness . The result shows that the effect of multi-governance is not directly related to the amount of public participation. In fact, since public reporting is more based on interests than authenticity, the lack of authenticity of public reporting often becomes a major obstacle to cooperative governance [43]. Therefore, government departments should focus on the fairness of public reporting in public rental housing exit governance during the epidemic to achieve a better effect of multi-governance.

4.3. The Influence of Subsidy on Strategy Choices of Tenants

To derive the impact of the subsidy on evolutionary results, we fix other parameters and conduct experiments when = 2, 3, respectively. The dynamic evolutionary processes of the tenants in different situation are demonstrated in Figure 4.

From Figure 4, it can be seen that the subsidy can effectively promote the tenants to select the “active exit” strategy. Considering the maximization of the payoff, the willingness of tenants to choose active exit strategy would increase with the increase in the subsidy. However, the high subsidy brings a heavy burden on the local government, so the willingness of the government to adopt strict governance would decline with the increase in the subsidy [44]. Thus, appropriate subsidy is a key factor that the government should consider in the multi-governance of public rental housing exit during the epidemic.

5. Conclusions

The government departments’ implementation of multi-governance for public rental housing exit requires considering the participation of multiple factors, such as epidemic and public, in the governance process and exploring the conditions under which each acts effectively. This paper adopts game theory to construct an evolutionary game model of the multi-governance for public rental housing exit, analyzes the influence of the epidemic and public parameters on the stable evolutionary state of the game system, and explores the conditions for the effective performance of the multi-governance model. The results show the following:(1)The epidemic impact coefficient, the public reporting fairness, and the government supervision capacity are the critical factors affecting the multi-governance of public rental housing exit.(2)Reducing the impact of the epidemic, improving the public reporting fairness, and increasing the housing subsidy are the principal factors affecting tenants’ adoption of active exit strategy.(3)The impact of the public participation degree on the effectiveness of governance is related to the level of the government departments' own supervision capacity, which decreases with the increase of supervision capacity and increases with the decrease of supervision capacity.

This study contributes to the present literature on collaborative governance theory in two areas. First, evolutionary game theory is utilized to study the multi-governance of public rental housing exit based on multiple stakeholders. This study compensates for the current lack of application of collaborative theory in public housing exit. Second, the effects of public participation on public rental housing exit are studied, which is of great value in the formulation and implementation of exit policies.

This study also has policy implications. First of all, government departments should pay full attention to the impact of the epidemic and try their best to suppress the epidemic transmission to promote tenants’ active exit. Second, government departments need to enhance their supervision capability to reduce the losses caused by information misreported by the public. At last, the power of the public and other participating parties needs to be utilized in the multi-governance model in order to better maintain the healthy development of public rental housing in China.

However, this manuscript still has some limitations. In order to simplify the model, this manuscript assumes that the relationship between epidemic impact and government supervision cost is linear, but the actual situation may be more complicated. Moreover, we have only analyzed the general trend of the influence of each factor on the evolution results but not solved the critical value. Accordingly, the solution of critical value and inflection point and the validation of numerical theoretical models with real case data will be the main direction of future research.

Data Availability

The 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 are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was funded by the Fundamental Research Funds for the Central Universities, China (Project no. 2019CDJSK03XK20), and jointly supported by the Fundamental Research Funds for the Central Universities (Project nos. 2017CDJXY and 2020CDJSK03PT08).