Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Public Health, 09 September 2025

Sec. Public Health Policy

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1640294

This article is part of the Research TopicEnvironment and Healthcare, a two-way traffic: Challenges, Impacts, and Sustainable SolutionsView all 8 articles

Game theoretic research on strategies for medical carbon neutrality stakeholders under time delay effects


Bo Xie,Bo Xie1,2Xuyang GaoXuyang Gao1Yingying Cheng
Yingying Cheng1*Muqing NiuMuqing Niu3Shijia LiuShijia Liu4
  • 1Business School and MBA Education Center of Henan University of Science and Technology, Luoyang, Henan, China
  • 2Henley Business School, University of Reading, Reading, United Kingdom
  • 3Affiliated First Hospital, Henan University of Science and Technology, Luoyang, Henan, China
  • 4School of Business Administration, Henan University of Economics and Law, Zhengzhou, Henan, China

Introduction: Carbon emissions from the healthcare sector have drawn increasing attention as global climate change intensifies. Achieving carbon neutrality in healthcare is vital for sustainable development, yet the complexity of stakeholder interactions and regulatory mechanisms creates substantial challenges.

Methods: This study develops a dynamic evolutionary game model with time-delay effects to examine the strategic behaviors of four stakeholders: government, public hospitals, pharmaceutical enterprises, and the public. Stability analysis is applied to determine equilibrium strategies, while MATLAB simulations explore the impacts of regulatory, financial, and reputational parameters on system dynamics.

Results: The findings reveal that government regulation is constrained by high costs and limited credibility, indicating the need to reduce costs and enhance credibility through market-based carbon trading and information-driven supervision. Public hospitals' weak compliance, driven by short-term profit motives, can be improved by strengthening financial incentives and penalties. Pharmaceutical enterprises' self-discipline is promoted by raising fines, reinforcing reputational incentives, and expanding public oversight channels. Public participation exerts a significant influence on carbon neutrality outcomes, and optimized online engagement mechanisms coupled with flexible compensation ratios are found to enhance effectiveness. Time delays accelerate system convergence but do not alter the evolutionary direction.

Discussion: These results underscore the critical role of multi-stakeholder interaction in advancing healthcare carbon neutrality. The study provides both theoretical insights and practical policy implications for designing more effective carbon neutrality strategies in the healthcare sector.

1 Introduction

In the 21st-century global agenda, climate change has emerged as an urgent problem with significant implications for human health (1) and wellbeing (2). As a key component of societal infrastructure, the healthcare industry's energy usage and greenhouse gas emissions during operation increasingly draw attention. Statistically speaking, the carbon emissions of the worldwide healthcare industry account for roughly 5% of the global total emissions (3), while certain industrialized nations, such as the United States, even exceed 10% (4). These data highlight the urgency and importance of reducing carbon emissions in the healthcare industry.

Against this backdrop, the healthcare industry has adopted the goal of medical carbon neutrality, aiming for net-zero carbon emissions by reducing emissions, improving efficiency, utilizing sustainable energy, and implementing carbon-offset technologies, thus balancing greenhouse gas outputs in clinical, operational, and supply chain processes with equivalent removals to achieve a net-zero environmental impact. For instance, the UK's National Health Service (NHS) commitment to carbon neutrality, targeting net-zero emissions by 2040 (5), exemplifies the emission reduction initiatives across the worldwide healthcare industry. The team from Solihull Hospital at the University Hospital of Birmingham in the UK completed the world's first carbon-neutral surgery in May 2022, reducing carbon emissions by 80% through equipment optimization, energy substitution, and other measures (6). The Cleveland Clinic in the United States has effectively decreased energy usage and carbon emissions via the adoption of energy-efficient practices and sustainable construction requirements. It aims to attain carbon neutrality by 2030 (7). Umea University Hospital in Sweden has lowered carbon emissions by 40 to 70 times relative to in-person treatment using remote rehabilitation therapy (8), this remarkable reduction exemplifies the pivotal role of digital therapeutics in advancing medical carbon neutrality.

The aforementioned instances illustrate that technology innovation and system optimization may proficiently regulate carbon emissions within the healthcare industry. Currently, the research focal points in healthcare carbon neutrality include energy management, the establishment of sustainable supply chains, the advancement and implementation of low-carbon medical technology, patient education (9), and public engagement. These studies focus on reducing the direct carbon emissions of the healthcare industry. Concurrently, these investigations aim to enhance the efficiency and quality of medical services while indirectly mitigating carbon emissions resulting from the increase in medical demand. In addition, the carbon neutrality measures in the healthcare industry also entail various levels, including policy formulation, financial investment, technological transfer, and public collaboration. Ultimately, these examples highlight that achieving medical carbon neutrality requires a comprehensive approach—one that combines stringent policy frameworks, strategic investment, cutting-edge technological innovations, and collaborative stakeholder engagement—to redefine healthcare as a sustainable, climate-resilient ecosystem.

The attainment of carbon neutrality in the medical industry may mitigate climate change, boost the efficiency of medical services, and bolster the resilience of the public health system on a worldwide scale. It is a crucial metric for advancing global sustainable development. Accordingly, we examined the policies and practices of different nations and regions in promoting medical carbon neutrality and developed a four-party evolutionary game model. This model examines the evolutionary dynamics among four principal roles: government, public hospitals, pharmaceutical enterprises, and the public. Subsequently, we delineate the strategic roles and decision-making environments of each participant inside the four-party game structure: the government acts as a worldwide policymaker and regulator, capable of directing the healthcare industry's low-carbon growth via law and financial incentives. Public hospitals, as direct partners in carbon neutrality initiatives, have simultaneous difficulties of improving healthcare service quality while minimizing environmental effect. Pharmaceutical enterprises serve as innovators and solution providers in this process. The public's health knowledge and consumption habits significantly influence medical carbon neutrality as end-users of healthcare services.

Furthermore, under uncertainty, decision-making processes are often delayed as stakeholders weigh long-term impacts (10). Therefore, incorporating time delay into the game framework is critical for accurately simulating strategic evolution and coordination among the four parties. This model contributes to a deeper understanding of how policy design and stakeholder behavior collectively shape the pathway toward achieving carbon-neutral healthcare systems.

2 Related literature

2.1 Medical carbon neutrality

The healthcare industry plays an important role in global carbon emissions and is currently facing significant sustainable development challenges. According to data from the World Health Organization (WHO), the global healthcare industry accounts for ~5% of total greenhouse gas emissions. If viewed as a country, it would become the fifth largest source of global emissions (11). The carbon emissions in the medical industry mainly come from energy consumption, anesthesia gases, the use of disposable products, and transportation and travel. Hospitals, especially operating rooms, have become the main source of carbon emissions due to their high energy consuming equipment (12, 13). Additionally, the use of anesthetics such as desflurane (14) and disposable medical devices such as sharps containers (15) greatly increases the carbon footprint. The transportation of patients and drugs (16) are also major contributors to carbon emissions in the healthcare industry. To address the above issues and achieve medical carbon neutrality, emission reduction strategies are essential. Currently, the main global emission reduction strategies include telemedicine (17), optimized management of medical equipment and operating room energy, and the use of low-carbon energy. Telemedicine, especially in replacing patients' short distance visits, has been proven to reduce carbon emissions by 40–70 times (16). By optimizing the design of operating rooms, using low global warming potential (GWP) anesthetics (18), and reusable medical devices (19), medical institutions can significantly reduce their carbon footprint. Furthermore, promoting clean energy (such as solar energy) and energy-saving technologies, as well as strengthening the classification and recycling management of medical waste, are also important paths for medical carbon neutrality (20). However, despite the increasing number of emission reduction measures and strategies for medical carbon neutrality, its implementation still faces many challenges. Firstly, the insufficient carbon emission data in the healthcare industry limits the comprehensive evaluation of emission reduction effects, especially the lack of sufficient empirical research on the comprehensive benefits of remote healthcare and low-carbon energy use throughout the entire lifecycle (8). Secondly, the lack of cross departmental cooperation and policy support makes it difficult to effectively implement carbon reduction measures throughout the entire healthcare system (21). The promotion of carbon neutrality in the healthcare industry requires not only technological innovation, but also policy incentives (21) and the transformation of industry culture (22). Strengthening cooperation between medical institutions, governments, patients, and suppliers, and promoting carbon footprint accounting and emission reporting systems in the healthcare industry will be key to achieving carbon neutrality in the future (21). In conclusion, attaining carbon neutrality in the healthcare industry demands both a reduction in energy consumption and an improvement in the efficiency of medical facilities and equipment, along with the promotion of multi-party collaboration through social engagement, policy suggestions, and technological progress. Despite facing certain challenges, the healthcare industry is expected to achieve carbon neutrality goals in the future by continuously improving technological means, strengthening policy implementation, and promoting public participation, making positive contributions to global environmental protection and public health.

2.2 The application of game theory in the research of medical carbon neutrality

The application of game theory in medical carbon neutralization has gradually attracted attention, which is mainly reflected in the strategic design of optimizing resource allocation, promoting cooperation and reducing carbon emissions. Recently, researchers have proposed several models based on game theory to meet the challenge of carbon neutrality in the medical field. Foremost, the game model of carbon trading market provides a theoretical framework for carbon emissions trading through blockchain Technology (23). Secondly, the cooperative game in medical waste treatment provides new ideas for reducing carbon emissions. Zhao et al. (24) constructed a tripartite game model among the government, medical institutions, and disposal enterprises, which promotes the resource utilization and disposal of medical waste through reasonable incentive and punishment mechanisms and reduces carbon emissions. The game model of the medical supply chain also demonstrates the potential of game theory in resource optimization. By mobilizing the game behavior of hospitals, patients, and the government, effective allocation of resources can be achieved, reducing the overall carbon emissions of the healthcare system (25). Yu (26) optimized internal control in public hospital supply chains through a game theory combination weighting method (AHP-entropy weighting method), established a comprehensive evaluation index system, and provided a new approach for the scientific management of medical supply chains. In addition, the dynamic regulatory and incentive model involving multiple parties also provides practical guidance for medical carbon neutrality (27). Although these game models provide theoretical basis for medical carbon neutrality, they still face challenges in data quantification, mechanism design, and cross domain integration. In the future, by combining evolutionary game theory with empirical data (13), more accurate carbon game frameworks will be developed, especially the application of technologies such as blockchain, which will provide more feasible solutions for medical carbon trading.

Evolutionary game theory is based on the assumption of bounded rationality, emphasizing the group dynamic evolution of strategies rather than individual optimal decisions, which is highly consistent with the complex interactions involving multiple entities (government, medical institutions, patients, suppliers, etc.) in medical carbon neutrality (28). Its core tools include Evolutionary Stability Strategy (ESS) (29), replicative dynamic (30), and reward and punishment mechanism design (31, 32). This provides a powerful tool for analyzing the long-term stability of various subject strategies in the context of medical carbon neutrality. In the study of specific application scenarios, the design of policy incentives is the primary field of evolutionary game theory. Scholars constructed a dynamic reward and punishment model including carbon tax, carbon subsidy and purchasing restriction to analyze the evolution path of healthcare organizations' strategies under different incentive strengths. For example, the National Health Service (NHS) in the UK requires providers to have a carbon reduction plan and imposes a purchasing ban on those who fail to meet the target (33, 34); the results of the game model show that when the penalty cost is higher than the abatement cost, low-carbon cooperation gradually becomes an evolutionarily stable strategy. Guo et al. (35) constructed a four-party evolutionary game model and pointed out that carbon prices, subsidy intensity, and punishment measures significantly affect the emission reduction decisions of supply chain enterprises. Among them, suppliers are sensitive to carbon prices, while manufacturers are more concerned about subsidies. This provides a reference for the government to assist in the implementation of medical carbon neutrality. Similarly, green technology promotion (e.g., introduction of solar energy facilities) can be evaluated through an evolutionary game framework to assess the impact of initial investment sharing vs. long-term energy returns on healthcare organizations' decision-making (36, 37). Evolutionary game theory also has forward-looking significance in patient participation and behavior guidance. Zhang et al. (38) used replicated dynamic equations to simulate patients' choices between remote healthcare and offline visits, and found that when carbon footprint information is transparent and the cost difference does not exceed 10%, the proportion of patients who prefer low-carbon diagnosis and treatment significantly increases. In addition, strategic synergy between physicians and patients through shared decision-making models can reduce carbon emissions by 20% to 30% in preventive healthcare promotion. This highlights the importance of public participation in building a carbon-neutral system for healthcare (39). However, although evolutionary game theory has made many advances in the field of healthcare carbon neutrality, scholars have mostly used the two-party or three-party evolutionary game approach to study the healthcare carbon neutrality problem, while the four-party evolutionary game approach to study healthcare carbon neutrality can further enrich the application of game theory in this field.

2.3 Application of time delays in supply chain carbon neutrality

Time delay, defined as the temporal delay between a system's input and its corresponding response, are commonly modeled by delay differential equations (DDEs). Such delays may manifest as discrete delay, distributed (continuous) delays, or impulsive delays (40). At present, time delays are widely employed across various disciplines. They occur in engineering systems (41), biological models (40), neural network architectures (42, 43), and in economic as well as traffic-flow analyses (44).

Time delays are also gradually receiving attention in the research of medical carbon neutrality technology. Biobased medical materials, such as surgical instrument packaging, can delay the release of CO2 beyond the target lifespan (such as 100 years) by extending their carbon storage time, and their emission reduction value is similar to permanent storage (45). From the perspective of technology promotion, the large-scale application of low-carbon drugs and devices that replace high GWP inhalers (such as salbutamol MDI) faces multiple delays in research and development verification, production transformation, and clinical acceptance, and usually takes several years to achieve significant emission reduction and efficiency improvement (46). Digitization and telemedicine are also expected to take 5–10 years to achieve significant results in patient transportation carbon emissions due to changes in population behavior habits and lagging infrastructure construction (47). At the policy and governance level, Portnoy et al. (48) draw on Norway's phased implementation of cervical cancer prevention policies to demonstrate how multi-stage evaluation processes can delay emissions reductions. Specifically, the approval and standard-setting for novel environmental technologies—such as low-carbon anesthesia protocols—are often prolonged by these successive assessments, thereby slowing the pace of decarbonization. Tennison et al. found that, between 1990 and 2019, the UK NHS achieved a 25% reduction in emissions by phasing out CFC-based inhalers, decarbonizing the electricity grid, and optimizing its supply chain; however, Scope 3 emissions (from pharmaceuticals and devices) remain difficult to curb due to delays in multinational production and procurement (49). Current literature reveals that research on time delays in medical supply chains and their carbon neutrality implications remains relatively scarce. Therefore, in this study, we integrated the time delay model into the four-party evolutionary game framework, enriching the research on time delay in the context of medical carbon neutrality.

Although some scholars have used game theory to study the issue of medical carbon neutrality, they still face challenges in mechanism design and cross disciplinary integration. Therefore, we use evolutionary game theory to study the issue of medical carbon neutrality and construct a more accurate framework for medical carbon neutrality games. Meanwhile, most existing literature has used two-party and tripartite evolutionary games to study medical carbon neutrality. Hence, we adopt the four-party evolutionary game method to investigate the issue of medical carbon neutrality. Finally, we incorporated time delays into the four-party evolutionary game model and further discussed the impact of time delays on medical carbon neutrality.

Table 1 summarizes the differences between our study and the existing literature reviewed previously. Specifically, most scholars use a tripartite evolutionary game model to study the carbon neutrality issue in the healthcare industry. Our research, however, has developed a four-party evolutionary game model that incorporates time delay factors, offering a new perspective for implementing carbon neutrality strategies in the medical industry. It further strengthens the literature on healthcare carbon neutrality by adjusting incentive and punishment policy mechanisms for stakeholders within a broader environmental sustainability framework.

Table 1
www.frontiersin.org

Table 1. Comparison between our study and correlated literature.

3 Problem description and model assumptions

3.1 Problem description

The issues of energy consumption and carbon emissions in the healthcare industry have become a focus of social concerns. Certain pharmaceutical enterprises prioritize economic gains above their environmental obligations, leading to inadequate investment in emission reduction initiatives. The government, as policymakers and regulators, may direct and oversee public hospitals and pharmaceutical enterprises by establishing pertinent regulations to advance carbon neutrality in the healthcare industry. Public hospitals, as healthcare providers, ought to lead in sustainable healthcare practices by minimizing carbon emissions through enhanced energy management and service quality, while collaboratively assuming responsibility with governments to advance the green transformation of the healthcare industry. However, conflicts of interest between public hospitals and pharmaceutical enterprises, along with insufficient governmental oversight, have emerged as ambiguous elements in attaining medical carbon neutrality and protecting public health rights. The public, as the primary consumers of healthcare services, is essential in advancing the sustainable growth of the medical business. But the fact is that insufficient public understanding of and participation in green healthcare services has become a major barrier to achieving carbon neutrality in healthcare.

Collaboration and games among the government, public hospitals, pharmaceutical enterprises, and the public are essential for achieving carbon neutrality in the healthcare industry. The interrelated strategies of participants create a dynamic environment where each player's choice affects and is affected by others. Figure 1 illustrates the logical links among the game elements of the four players.

Figure 1
Flowchart illustrating the interactions between government, public hospitals, pharmaceutical enterprises, and the public. It shows regulatory and enforcement processes, including strict or loose regulation by the government, policy enforcement by hospitals, enterprise promotion, and public participation. Arrows indicate various enforcement strengths and mutual responsibilities such as compensations, regulations, and participation in carbon neutrality.

Figure 1. Relationship diagram of four players in the game of medical carbon neutrality.

3.2 Model assumption

3.2.1 Hypothesis 1

Government, public hospitals, pharmaceutical enterprises, and the public are selected as the subjects of the four-party game. The probability of government enforcing strict regulation or loose regulation is (x, 1−x). The probability of public hospitals performing strong enforcement or weak enforcement is (y, 1−y). The probability of pharmaceutical enterprises adopting self-discipline or not self-discipline is (z, 1−z). The probability that the public chooses participation or non-participation is (w, 1−w). Here, x, y, z, w∈[ 0, 1].

3.2.2 Hypothesis 2

We assume that the strict regulation cost of the government is Cg. If the public hospitals carry out strong enforcement, pharmaceutical enterprises practice self-discipline, which can benefit the government, and we assume that this benefit is Rg. Additionally, under public participation, this scenario can enhance the government's social credibility, represented by kg. Conversely, the government will incur losses Lg if either condition is met: (1) public hospitals adopt weak enforcement strategies, or (2) pharmaceutical enterprises engage in not self-discipline practices. For the public, the cost of participation is defined as Cp. Meanwhile, when public hospitals choose strong enforcement and pharmaceutical enterprises choose self-discipline, the public will gain benefits Sp. While when the above conditions are not met, the public will suffer losses Lp.

3.2.3 Hypothesis 3

Under strict government regulation, if public hospitals strongly enforce medical carbon neutrality policies and pharmaceutical enterprises adopt self-discipline strategies, they can receive government rewards, denoted as Rh and Rc, respectively. Simultaneously, with public participation, if public hospitals strongly enforce medical carbon neutrality policies and pharmaceutical enterprises adopt self-discipline strategies, the credibility of both will be enhanced, denoted as kh and kc, respectively. Assume the cost of strong enforcement for public hospitals is Ch, and the cost of self-discipline for pharmaceutical enterprises is Cc. If public hospitals choose weak enforcement, they can obtain short-term benefits Sh. However, if the government imposes strict regulations at this time, the public hospitals will be required to pay a fine Dh to the government. Similarly, when pharmaceutical enterprises choose not self-discipline, they can obtain short-term income, which is Sc. If the government carries out strict regulations in this scenario, pharmaceutical enterprises need to pay a fine Dc to the government.

3.2.4 Hypothesis 4

The fines paid by public hospitals and pharmaceutical enterprises can serve as government revenue. Furthermore, when the public participates, these fines can be used to compensate the public. Suppose that α indicates the proportion of penalties obtained by the government during strict regulation, and 1−α expresses the proportion allocated to public compensation obtained by the public when the public participates, where α ∈ (0,1). It is essential that the public can only receive compensation if they participate when the government enforces strict supervision; otherwise, no compensation will be provided.

The definitions of the main parameter symbols are summarized in Table 2.

Table 2
www.frontiersin.org

Table 2. Parameter definitions.

4 Stability analysis of players' strategy choices

Based on Hypotheses 1–4 in Section 3, a game payment matrix of the government, public hospitals, pharmaceutical enterprises, and the public are constructed under different strategy choices, which are shown in Table 3.

Table 3
www.frontiersin.org

Table 3. Payoff matrix for the quadripartite game.

4.1 Stability analysis for government

From the game payment matrix, it can be seen that the expected benefits for a government with strict regulation is Ex1, the expected benefits for a government with loose regulation is Ex2, and the average expected income for the government is E¯x:

Ex1=yzw(Rg-Cg-Rh-Rc+kg)+yz(1-w)   (Rg-Cg-Rh-Rc)+y(1-z)w(αDc-Cg-Lg)   +y(1-z)(1-w)(Dc-Cg-Lg)+(1-y)   zw(αDh-Cg-Lg)+(1-y)z(1-w)(Dh-Cg-Lg)   +(1-y)(1-z)w[α(Dc+Dh)-Cg-Lg]   +(1-y)(1-z)(1-w)(Dh+Dc-Cg-Lg)     (1)
Ex2=yzw(Rg-Cg)+yz(1-w)Rg+y(1-z)   w(-Lg-kg)+y(1-z)(1-w)(-Lg)+(1-y)   zw(-Lg-kg)+(1-y)z(1-w)(-Lg)+(1-y)(1-z)   w(-Lg-kg)+(1-y)(1-z)(1-w)(-Lg)     (2)
E¯x=xEx1+(1-x)Ex2     (3)

The replication of dynamic (50) is the core of evolutionary game theory, which describes the evolution process of the frequency of each player's strategy. In this paper, the replication dynamic equation of the government can help us better understand the impact of different strategies of the government as a subject on healthcare carbon neutrality. Based on this, we calculate the government's replication dynamic equation based on its expected and average benefits as follows:

F(x)=dxdt=x(Ex1Ex2)=x(1x)(Ex1E¯x)=x(1x){wkgCgy[(αww+1)Dh+z(Rc+RhwCg)](wαw1)(Dc+Dh)z[(α1)w+1]Dc}    (4)

The stability analysis of replicating dynamic equations is a touchstone for testing the feasibility of economic equilibrium evolution. This can provide dynamic optimization directions for institutional design. Therefore, we perform first-order differentiation on the government's replication dynamic equation for stability analysis. The first derivative of F(x) is as follows:

F(x)=(12x){wkgCgy[(αww+1)Dh            +z(Rc+RhwCg)](wα1)(Dc+Dh)            z[(α1)w+1]Dc}    (5)

It can be seen that whether the government chooses strict regulation depends on the cost and benefit of its different regulatory strategies and the probability of the other three parties' decision-making choices. According to the stability theorem of differential equations, the government's choice of strict supervision in a steady state must satisfy conditions F(x) = 0 and F′(x) < 0.

4.1.1 Proposition 1

When y>y1, z>z1, w<w1, the government's stable strategy is loose regulation. When y<y1, z<z1, w>w1, the government's stable strategy is strict regulation. When y = y1, z = z1, w = w1, the government's stable strategy cannot be determined. The threshold value is as follows:

{y1=wkgCg(wαw1)(Dc+Dh)z[(α1)w+1]Dc(αww+1)Dh+z(Rc+RhwCg)z1=wkgCgy(αww+1)Dh(wαw1)Dcy(Rc+RhwCg)+[(α1)w+1]Dcw1=(z1)Dc+(y1)Dh+zy(Rc+Rh)+Cgwkgy[(1α)Dh+zCg]+(α1)(Dc+Dh)z(α1)Dc    (6)

Among them, y1 represents the threshold of public hospitals under the evolution of government strategies, which is the boundary point between weak or strong enforcement strategies adopted by public hospitals. z1 represents the threshold for pharmaceutical enterprises under the evolution of government strategies, which is the boundary between self-discipline or not self-discipline strategies adopted by pharmaceutical enterprises. w1 represents the threshold of the public under the evolution of government strategies, which is the dividing point between the public's adoption of participation or non-participation strategies. We have taken the threshold y1 of public hospitals as an example to explain how to calculate the threshold. The specific calculation process is as follows:

Construct function G(y, z, w) = wkgCgy[(αww+1)Dh+z(Rc+RhwCg)]−(w−αw−1)(Dc+Dh)−z[(α−1)w+1]Dc, as G(y, z, w) = 0, where y is an unknown variable and z and w are definite, solve for the threshold y1=wkg-Cg-(w-αw-1)(Dc+Dh)-z[(α-1)w+1]Dc(αw-w+1)Dh+z(Rc+Rh-wCg) of y. Similarly, the process of solving the threshold z1 for pharmaceutical companies and the threshold w1 for the public is also the same.

4.1.2 Proof of proposition 1

As can be seen from the previous text, G(y, z, w) = wkgCgy[(αww+1)Dh+z(Rc+RhwCg)]−(w−αw−1)(Dc+Dh)−z[(α−1)w+1]Dc, as G(y,z,w)y<0, G(y, z, w) is strictly decreasing function in terms of y. When y>y1, G(y, z, w) < 0, F(x)|x = 0 = 0 and F(x)|x = 0 < 0, the equilibrium point x = 0 is asymptotically stable; when y<y1, G(y, z, w)>0,F(x)|x = 1 = 0 and F(x)|x = 0 < 1, then x = 1 is asymptotically stable; when y = y1, F(x) = 0 and F′(x) = 0, the governments stabilization strategy cannot be determined. By condition G(y,z,w)z<0 and G(y,z,w)w>0, the influence of the critical thresholds z1 and w1 on the strategy's stability can be similarly demonstrated. The proof of Proposition 1 is complete.

Proposition 1 shows that in the process of medical carbon neutral policy enforcement, if the probability of strong enforcement of public hospitals and the probability of self-discipline of pharmaceutical enterprises increase, and the probability of public participation decreases, the government will choose the loose regulation strategy; on the contrary, the government will shift from loose regulation to strict regulation.

According to Proposition 1, the phase diagram of government regulation strategy selection is drawn in Figure 2.

Figure 2
Three diagrams of three-dimensional graphs with axes labeled x, y, and z. (A) shows a cube with a diagonal plane at y = y1. (B) shows a partial cube with a diagonal plane at y < y1. (C) shows a cube with a diagonal plane at y ≥ y1 and marked Vx1. Arrows indicate directions along the axes.

Figure 2. Phase diagram of government strategy selection. (A) y = y1. (B) y < y1. (C) y > y1.

As shown in Figure 2, the volumes of regions Vx0 and Vx0 represent the probabilities of strict and loose government regulation, respectively. Since we cannot determine the specific values of the points on the y-axis in Figure 2, we set z = 0 and w = 0; simultaneously calculate y=Dc-yDh-Cgy(Rc+Rh)+Dc based on y1. Finally, to better calculate the volumes of Vx0 and Vx0, we construct the equation v=Dc-yDh-Cgy(Rc+Rh)+Dc, the calculation yields are as follows:

Vx0=01z0vd1ydx       =wkgCg(wαw1)(Dc+Dh)Rc+RhCg        ln[1+v(Rc+RhwCg)(αww+1)Dh][(α1)w+1]Dc        ln{v(Rc+RhwCg)(αww+1)Dh        ln[1+v(Rc+RhwCg)(αww+1)Dh]}     (7)
Vx1=101zv1d1ydx        =1wkgCg(wαw1)(Dc+Dh)Rc+RhCg        ln[1+v(Rc+RhwCg)(αww+1)Dh]+[(α1)w+1]Dc        ln{v(Rc+RhwCg)(αww+1)Dh        ln[1+v(Rc+RhwCg)(αww+1)Dh]}     (8)

4.1.3 Corollary 1

Government choices regarding strict regulation strategies are affected by a variety of factors. These factors include the improvement of government credibility; the increase in fines paid by public hospitals adopting weak enforcement strategies; the increase in penalties paid by pharmaceutical enterprises adopting not self-discipline strategies; the reduction in government regulatory costs; and the reduction in the reward scale when public hospitals adopt strong enforcement strategies and pharmaceutical enterprises maintain self-discipline. If these conditions are met, the likelihood of the government implementing strict regulation will significantly increase. Conversely, governments tend to prefer loose regulation strategies if these conditions are reversed.

4.1.4 Proof of corollary 1

By solving the first-order partial derivatives of the government's choice of the strict regulation probability Vx0, with respect to kg, Dc, Dh, Cg, Rc, and Rh, we can obtain: Vx0kg>0, Vx0Dc>0, Vx0Dh>0, Vx0Cg<0, Vx0Rc<0,Vx0Rh<0, certificate completed.

4.2 Stability analysis for public hospitals

From the game payment matrix, it can be seen that the expected benefits for the public hospitals with strong enforcement is Ey1, the expected benefits for the public hospital with weak enforcement is Ey2, and the average expected income for the public hospitals is E¯y:

Ey1=xzw(Rh-Ch+kh)+xz(1-w)(Rh-Ch)+x(1-z)   w(-Ch+kh)+x(1-z)(1-w)(-Ch)+(1-x)   zw(-Ch+kh)+(1-x)z(1-w)(-Ch)+(1-x)   (1-z)w(-Ch+kh)+(1-x)(1-z)(1-w)(-Ch)     (9)
Ey2=xzw(Sh-Dh-kh)+xz(1-w)(Sh-Dh)+x(1-z)   w(Sh-Dh-kh)+x(1-z)(1-w)(Sh-Dh)   +(1-x)zw(Sh-kh)+(1-x)z(1-w)Sh+(1-x)(1-z)   w(Sh-kh)+(1-x)(1-z)(1-w)Sh     (10)
E¯y=yEy1+(1-y)Ey2     (11)

The replication dynamic equation of public hospitals can help us better understand the impact of different strategies of public hospitals on healthcare carbon neutrality. Based on this, we calculated the public hospitals' replication dynamic equation based on the expectations and average benefits of public hospitals. The equation for the replication dynamics of public hospitals is as follows:

F(y)=dydt=y(Ey1-Ey2)   =y(1-y)(Ey1-E¯y)   =y(1-y)[x(Dh+zRh)+2wkh-Sh-Ch]     (12)

To conduct stability analysis on the replication dynamic equation of public hospitals and identify the impact of the other three parties on public hospital strategies, we performed first-order differentiation on the replication dynamic equation of public hospitals. The first derivative of F(y) is as follows:

F(y)=(1-2y)[x(Dh+zRh)+2wkh-Sh-Ch]     (13)

From Equations 12, 13, it can be seen that the main factors affecting the strategy selection of public hospitals include when public hospitals enforce strong enforcement strategy, the government's strict regulation brings rewards Rh to public hospitals; the cost Ch of public hospitals adopting a strong enforcement strategy; the fine Dh that public hospitals have to pay when adopting a weak enforcement strategy under strict government regulation; the reputation kh of public hospitals; and the short-term benefits Sh of public hospitals adopting a weak enforcement strategy. According to the stability theorem of differential equations, the steady-state conditions for selecting strong enforcement in public hospitals should be satisfied: F(y) = 0 and F′(y) < 0.

4.2.1 Proposition 2

When x<x1, z<z2, w<w2, the stable strategy of public hospitals is to choose weak enforcement; When x>x1, z>z2, w>w2, the stable strategy of public hospitals is to choose strong enforcement; When x = x1, z = z2, w = w2, the public hospitals' stable strategy cannot be determined. The threshold value is as follows:

{x1=Sh+Ch2wkhDh+zRhz2=Sh+ChxDh2wkhxRhw2=Sh+ChxDhxzRh2kh    (14)

Among them, x1 represents the threshold of the government's strategy evolution in public hospitals, which is the boundary point between the government's adoption of strict regulatory strategies or loose regulatory strategies. z2 represents the threshold for pharmaceutical enterprises in the evolution of public hospital strategies, which is the boundary between the adoption of self-discipline strategies or non-self-discipline strategies by pharmaceutical enterprises. w2 represents the threshold of the public's participation in the evolution of public hospital strategies, which is the dividing point between the public's adoption of participation strategies or non-participation strategies. The calculation process has been explained in Section 4.1, so it will not be elaborated here.

4.2.2 Proof of proposition 2

Construct the function H(x, z, w) = x(Dh+zRh)+2wkhShCh. Since H(x,z,w)x>0, H(x, z, w) is an increasing function with regard to x. When x<x1, we have that H(x, z, w) < 0, F(y)|y = 0 = 0, F(y)|y = 0 < 0, and then y = 0 has stability. When x>x1, H(x, z, w)>0, F(y)|y = 1 = 0, F(y)|y = 0 < 0, and then y = 1 has stability. When x = x1, we can obtain that H(x, z, w) = 0, and then F(y) = 0 and F′(y) = 0,. Now, no stable strategy can be identified. By condition H(x,z,w)z<0 and H(x,z,w)w<0, the influence of the critical thresholds z2 and w2 on the strategy's stability can be similarly demonstrated. The proof of Proposition 2 is complete.

Proposition 2 indicates that during the enforcement of the medical carbon neutrality strategy, if the government chooses strict regulation, pharmaceutical enterprises choose self-discipline, and the probability of public participation increases, the probability of public hospitals choosing strong enforcement will also increase. Conversely, public hospitals will shift from strong enforcement to weak enforcement.

According to Proposition 2, Figure 3 shows the phase diagram for the strategy selection of public hospitals.

Figure 3
Three diagrams of a cube with x, y, and w axes. (A) shows x = x1 with vertical y-direction arrows. (B) shows x < x1 with Vy0 direction arrow. (C) shows x > x1 with Vy1 direction arrow.

Figure 3. Phase diagram of public hospitals' strategy selection. (A) x = x1. (B) x < x1. (C) x > x1.

As can be seen from Figure 3, the volumes of the Vy0 and Vy1 sections represent the probability of weak and strong enforcement by public hospitals, respectively, which are calculated as follows:

Vy0=01x01d1ydw =Sh+ChkhDh+zRh    (15)
Vy1=101x01d1ydw=1Sh+ChkhDh+zRh    (16)

4.2.3 Corollary 2

Many factors influence public hospitals to adopt weak enforcement strategies. The factors include the increase in short-term incomes of public hospitals; public hospitals' costs increase when they adopt a strong enforcement strategy; the credibility of public hospitals declines; the reward from the government to the public hospitals is reduced when the public hospitals adopt a strong enforcement strategy; the government reduces penalties for public hospitals implementing weak enforcement strategies. If these conditions are met, the likelihood of the public hospitals implementing weak enforcement will significantly increase. Conversely, when these conditions reverse, public hospitals tend to shift toward strong enforcement strategies. The proof process is similar to the analysis of government strategy stability.

4.3 Stability analysis for pharmaceutical enterprises

From the game payment matrix, it can be seen that the expected benefits for the pharmaceutical enterprises with self-discipline is Ez1, the expected benefits for the pharmaceutical enterprises with not self-discipline is Ez2, and the average expected income for the pharmaceutical enterprises is E¯z:

Ez1=xyw(Rc-Cc+kc)+xy(1-w)(Rc-Cc)+x(1-y)   w(-Cc+kc)+x(1-y)(1-w)(-Cc)+(1-x)   yw(-Cc+kc)+(1-x)y(1-w)(-Cc)+(1-x)(1-y)   w(-Cc+kc)+(1-x)(1-y)(1-w)(-Cc)     (17)
Ez2=xyw(Sc-Dc-kc)+xy(1-w)(Sc-Dc)   +x(1-y)w(Sc-Dc-kc)   +x(1-y)(1-w)(Sc-Dc)   +(1-x)yw(Sc-kc)+(1-x)y(1-w)Sc   +(1-x)(1-y)w(Sc-kc)+(1-x)(1-y)(1-w)Sc     (18)
E¯z=zEz1+(1-z)Ez2     (19)

The replication dynamic equation of pharmaceutical enterprises can help us better understand the impact of different strategies of pharmaceutical enterprises on healthcare carbon neutrality. Based on this, we calculated the replication dynamic equation for pharmaceutical enterprises based on their expectations and average benefits. The replication dynamic equation for pharmaceutical enterprises is as follows:

F(z)=dzdt=z(Ez1E¯z)=z(1z)[x(Dc+yRc)          +2wkcScCc]     (20)

To conduct stability analysis on the replication dynamic equation of pharmaceutical enterprises and determine the impact of the other three parties on the strategy of pharmaceutical enterprises, we performed first-order differentiation on the replication dynamic equation of pharmaceutical enterprises. The first derivative of F(z) is as follows:

F(z)=(1-2z)[x(Dc+yRc)+2wkc-Sc-Cc]     (21)

It can be seen from Equations 20, 21 that the main factors influencing the strategic choices of pharmaceutical enterprises include: under strict government regulation, the government rewards self-disciplined pharmaceutical enterprises with Rc; the cost of self-discipline Cc for pharmaceutical enterprises; the reputation enhancement degree kc; the short-term benefits Sc from not self-discipline; and when the government strictly regulates, the pharmaceutical enterprises do not self-discipline, pay the fine Dc. According to the stability theorem of differential equations, the steady-state conditions for selecting self-discipline in pharmaceutical enterprises should be satisfied: F(z) = 0 and F′(z) < 0.

4.3.1 Proposition 3

When x<x2, y<y2, w<w3, the pharmaceutical enterprises' stable strategy is not self-discipline. When x>x2, y>y2, w>w3, the pharmaceutical enterprises' stable strategy is self-discipline; When x = x2, y = y2, w = w3, the pharmaceutical enterprises' stability strategy cannot be determined. The threshold value is as follows:

{x1=Sc+Cc2wkcDc+zRcz2=Sc+CcxDc2wkcxRcw2=Sc+CcxDcxzRc2kc    (22)

Among them, x2 represents the threshold of the government's strategic evolution in the pharmaceutical enterprises, which is the boundary point between the government's adoption of strict regulatory strategies or loose regulatory strategies. y2 represents the threshold of public hospitals under the evolution of pharmaceutical enterprises strategies, which is the boundary point for public hospitals to adopt strong enforcement strategies or weak enforcement strategies are implemented. w3 represents the threshold of the public in the evolution of pharmaceutical enterprises strategies, which is the dividing point between the public adopting participation strategies or not participating strategies. The specific calculation process has been explained in Section 4.1.

4.3.2 Proof of proposition 3

Construct the function P(x, y, w) = x(Dc+yRc)+2wkcScCc. Because P(x,y,w)x>0, P(x, y, w) is an increasing function with regard to x. When x<x2, we have that P(x, y, w) < 0, F(z)|z = 0 = 0, F(z)|z = 0 < 0, and then z = 0 has stability. When x>x2, P(x, y, w)>0, F(z)|z = 1 = 0, F(z)|z = 1 < 0, and then z = 1 has stability. When x = x2, we can obtain that P(x, y, w) = 0, and then F(z) = 0 and F′(z) = 0,. Now, no stable strategy can be identified. By condition P(x,y,w)y<0 and P(x,y,w)w<0, the influence of the critical thresholds y2 and w3 on the strategy's stability can be similarly demonstrated. The proof of Proposition 3 is complete.

Proposition 3 implies that during the implementation of medical carbon neutrality strategies, when the government enforces strict regulation, public hospitals demonstrate strong enforcement, and the probability of public participation increases, the probability of pharmaceutical enterprises practicing self-discipline also rises; Conversely, pharmaceutical enterprises will shift from self-discipline to not self-discipline.

According to Proposition 3, the phase diagram for the strategy selection of a pharmaceutical enterprise is drawn, as shown in Figure 4.

Figure 4
Three diagrams labeled A, B, and C show different 3D plots with x, y, z, and w axes. Diagram A depicts a cube when w equals w3. Diagram B shows a cube with a diagonal when w is less than w3, labeled Vz0. Diagram C illustrates a cube with an upward diagonal when w is greater than w3, labeled Vz1.

Figure 4. Phase diagram of pharmaceutical enterprises' strategy selection. (A) w = w3. (B) w < w3 (C) w > w3.

As shown in Figure 4, the volumes corresponding to regions Vz0 and Vz1 represent the probabilities of pharmaceutical enterprises' not self-discipline and self-discipline, respectively. The calculation yields:

Vz0=01w01d3xdz=Sc+CckcDc+ykc    (23)
Vz1=101w01d3xdz=1Sc+CckcDc+ykc    (24)

4.3.3 Corollary 3

Pharmaceutical enterprises' choices regarding self-discipline strategies are affected by various factors. These factors include the increase of short-term profits of pharmaceutical enterprises; increased cost of self-discipline in pharmaceutical enterprises; when the government strictly regulates, the fines that pharmaceutical enterprises have to pay for no self-discipline are reduced; pharmaceutical enterprises' credibility declines. If these conditions are met, the probability of the pharmaceutical enterprises adopting not self-discipline will significantly increase. Conversely, if these conditions reverse, pharmaceutical enterprises tend to prefer self-discipline strategies.

4.4 Stability analysis for public

From the game payment matrix, it can be seen that the expected benefits for the public with participation is Ew1, the expected benefits for the public with non-participation is Ew2, and the average expected income for the public is E¯w:

Ew1=xyz(Sp-Cp)+xy(1-z)[(1-α)Dc-Cp-Lp]   +x(1-y)z[(1-α)Dh-Cp-Lp]+x(1-y)(1-z)   [(1-α)(Dh+Dc)-Lp-Cp]+(1-x)yz(Sp-Cp)   +(1-x)y(1-z)(-Cp-Lp)+(1-x)(1-y)z(-Lp-Cp)   +(1-x)(1-y)(1-z)(-Lp-Cp)    (25)
Ew2=xyzSp+xy(1-z)(-Lp)   +x(1-y)z(-Lp)+x(1-y)(1-z)(-Lp)   +(1-x)yzSp+(1-x)y(1-z)(-Lp)+(1-x)   (1-y)z(-Lp)+(1-x)(1-y)(1-z)(-Lp)     (26)
E¯w=wEw1+(1-w)Ew2     (27)

The public's replication of dynamic equations can help us better understand the impact of different public strategies on healthcare carbon neutrality. Based on this, we calculated the public's replication dynamics equation based on public expectations and average returns. The replication dynamic equation for pharmaceutical companies is as follows:

F(w)=dwdt=w(Ew1E¯w)=w(1w){x(1α)(Dc+Dh)            x[(1α)yDh+(1α)zDc]Cp}    (28)

To conduct stability analysis on the public replication dynamic equation and determine the impact of the other three parties on public strategy, we performed first-order differentiation on the public replication dynamic equation. The first derivative of F(w) is as follows:

F(w)=(12w){x(1α)(Dc+Dh)               x[(1α)yDh+(1α)zDc]Cp}    (29)

It can be seen from Equations 28, 29 that the main factors affecting the public's strategy choice include the fine Dh that public hospitals have to pay when they adopt a weak enforcement strategy under strict government regulation, the fine Dc that pharmaceutical enterprises have to pay when they adopt a not self-discipline strategy, and the cost Cp when the public participates.

4.4.1 Proposition 4

When x<x3, y>y3, z>z3, the public's stable strategy is non-participation; When x>x3, y<y3, z<z3, the public's stable strategy is participation; When x = x3, y = y3, z = z3, the public stability strategy cannot be determined. The threshold value is:

{x3=Cp(1α)(Dc+Dh)(1α)yDh(1a)zDcy3=x(1α)(Dc+Dh)xz(1α)DcCpx(1α)Dhz3=x(1α)(Dc+Dh)xy(1α)Dhx(1α)Dh    (30)

Among them, x3 represents the threshold of the government under the evolution of public strategies, which is the boundary point between the government adopting strict regulatory strategies or loose regulatory strategies. y3 represents the threshold of public hospitals under the evolution of public strategies, which is the boundary point for public hospitals to adopt strong or weak enforcement strategies. z3 represents the threshold for pharmaceutical enterprises to adopt self-discipline strategies or not self-discipline strategies under the evolution of public policies. The specific calculation process has been explained in Section 4.1 and will not be further elaborated here.

4.4.2 Proof of proposition 4

Construct the function Q(x, y, z) = x(1−α)(Dc+Dh)−x[(1−α)yDh+(1−α)zDc]−Cp. Since Q(x,y,z)x>0, Q(x, y, z) is an increasing function about x. When x<x3, we have that Q(x, y, z) < 0, F(w)|w = 0 = 0, F(w)|w = 0 < 0, and then w = 0 has stability. When x>x3, Q(x, y, z)>0, F(w)|w = 1 = 0, F(w)|w = 1 < 0, and then w = 1 has stability. When x = x3, we can obtain that Q(x, y, z) = 0, and then F(w) = 0 and F′(w) = 0,. Now, no stable strategy can be identified. By condition Q(x,y,z)y>0 and Q(x,y,z)w>0, the influence of the critical thresholds y3 and z3 on the strategy's stability can be similarly demonstrated. The proof of Proposition 4 is complete.

Proposition 4 suggests that during the implementation of medical carbon neutrality strategies, if the government adopts loose regulation, public hospitals adopt strong enforcement, and the probability of pharmaceutical enterprises adopting self-discipline increases, the public will opt for non-participation; Conversely, the public will shift from non-participation to participation.

According to Proposition 4, the phase diagram for the public's strategy selection is drawn, as shown in Figure 5.

Figure 5
Three diagrams labeled (A), (B), and (C) depict cubes with axes labeled x, y, and w. Diagram (A) shows arrows on the w axis, and both (B) and (C) feature curved surfaces. Specific conditions for x, such as x = x3, x < x3, and x > x3, are noted below each diagram.

Figure 5. Phase diagram of public strategy selection. (A) x = x3. (B) x < x3. (C) x > x3.

As shown in Figure 5, the volumes corresponding to regions Vw0 and Vw1 represent the probabilities of public' participation and non-participation, respectively. The calculation yields:

Vw0=01x01d3ydw=Cp(1α)Dhln[(1α)(Dc+Dh)         (1α)yDh(1α)zDc]     (31)
Vw1=101d01ydw=1+Cp(1α)Dhln[(1α)(Dc+Dh)        (1α)yDh(1 α)zDc]    (32)

4.4.3 Corollary 4

When the public adopts participation strategies, the increased fines imposed on public hospitals for choosing weak enforcement strategies, the increased fines imposed on pharmaceutical enterprises for choosing not self-discipline strategies, and the reduction in public participation costs will all increase the probability of the public choosing participation strategies. Conversely, the public will opt for non-participation strategies. The proof process is similar to the stability analysis of government strategies.

5 Stability analysis of strategy combinations

In a four-party evolutionary game replication dynamic system, there are 16 potential strategy combination scenarios, each of which may become the equilibrium point of the group evolutionary game—conducting stability analysis on the equilibrium points of the replication dynamic system, which helps to grasp the emergence of different evolutionary scenarios. Furthermore, adjusting key parameters promotes the evolution of the game system toward Pareto optimal scenarios and avoids the emergence of poor game equilibria. The stability of the game players' strategy combination can be judged according to Lyapunov's first rule (51). If all eigenvalues of the Jacobian matrix (52) have negative genuine parts, the equilibrium point is an asymptotically evolutionarily stable strategy (ESS) (53). The equilibrium point is unstable if at least one eigenvalue of the Jacobian matrix has a positive real part. If all eigenvalues of the Jacobian matrix, except for zero eigenvalues, have negative real parts, the equilibrium point is in a critical state with uncertain stability. Based on the replication dynamics equations of each game player, the Jacobian matrix of the replication dynamic system is:

J=[F(x)/xF(x)/yF(x)/zF(x)wF(y)/xF(y)/yF(y)/zF(y)/wF(z)/xF(z)/yF(z)/zF(z)/wF(w)/xF(w)/yF(w)/zF(w)/w]     (33)

As a supervisory body, the government is responsible and obligated to implementing medical carbon neutrality strategies. To study the impact of government regulation strategies on system evolution through comparative analysis, we divided the 16 pure strategy equilibrium points into two scenarios: loose regulation and strict regulation. This can focus on the evolution and stability of public hospitals and pharmaceutical enterprises under regulatory strategies and provide a theoretical basis for building a more effective government regulatory system.

5.1 Asymptotic stability analysis of equilibrium points in replication dynamic systems under loose government regulation

During the implementation of medical carbon neutrality strategies, when the government's stable strategy is loose regulation, that is, when the following condition is satisfied: wkgCgy[(αww+1)Dh+z(Rc+RhwCg)]−(w−αw−1)(Dc+Dh)−z[(α−1)w+1]Dc < 0, the stability analysis of the equilibrium points in the replication dynamic system is summarized in Table 4.

Table 4
www.frontiersin.org

Table 4. Analysis of the asymptotic stability of the equilibrium point of replication dynamic systems under loose government regulation.

As shown in Table 4, in the case of loose regulation by the government, there is no stable strategy that can effectively achieve medical carbon neutrality. This is because the government's loose regulation strategy cannot effectively restrain the improper behavior of public hospitals and pharmaceutical enterprises. To avoid this, the government should increase the punishment for public hospitals that adopt weak enforcement strategies and pharmaceutical enterprises that choose not self-discipline strategies. Simultaneously, the government should also reduce its regulatory costs through relevant measures. Therefore, government must increase penalties for public hospitals adopting weak enforcement strategies and pharmaceutical enterprises employing not self-discipline strategies, while reducing their regulatory costs, to achieve medical carbon neutrality.

5.2 Asymptotic stability analysis of the equilibrium points of replication dynamic systems under strict government regulation

During the implementation of medical carbon neutrality strategies, when the government's stable strategy is strict regulation, that is, when the following condition is satisfied: kgCgy[(αww+1)Dh+z(Rc+RhwCg)]−(w−αw−1)(Dc+Dh)−z[(α−1)w+1]Dc>0, the stability analysis of the equilibrium points in the replication dynamic system is summarized inTable 5.

Table 5
www.frontiersin.org

Table 5. Analysis of the asymptotic stability of the equilibrium point of replicating dynamic systems under strict government regulation.

From Table 5, it can be seen that under the strict regulation of government, there are 6 possible stable strategy combinations, namely (1, 0, 0, 0), (1, 0, 0, 1), (1, 0, 1, 0), (1, 0, 1, 1), (1, 1, 0, 0), (1, 1, 0, 1); There are 3 stable strategy combinations of public participation (1, 0, 0, 1), (1, 0, 1, 1), (1, 1, 0, 1). When the public participates and the government chooses strict regulation, pharmaceutical enterprises will adopt self-discipline strategies to avoid excessive losses; public hospitals will tend to adopt a weak enforcement strategy due to the strategic choice of the other three players in the game. In this scenario, the strategy combination (1, 0, 1, 1) becomes the evolutionarily stable policy combination of the system. However, even when government strictly regulates, the public participates, and public hospitals enforce rigorously, some pharmaceutical enterprises may still choose not self-discipline strategies in pursuit of excessive profits, which is the underlying reason why the strategy combination (1, 1, 0, 1) can become a stable strategy. To prevent (1, 1, 0, 1) from stabilizing, the government must significantly increase penalties for pharmaceutical enterprises adopting not self-discipline strategies. Simultaneously, consumers should also actively participate in public opinion supervision of pharmaceutical enterprises to ensure the smooth realization of carbon neutrality in the medical industry. We can get some enlightenment from the observation of these stable equilibrium points: to safeguard public health rights and interests for a long time and promote the green transformation of the medical industry, government must enforce effective supervision over public hospitals and pharmaceutical enterprises to ensure their proactive engagement in medical carbon neutrality initiatives. The government must invest corresponding regulatory resources to monitor the progress of these key subjects in reducing carbon emissions and improving energy efficiency. Furthermore, the government's strict regulation of public hospitals and pharmaceutical enterprises can effectively curb inefficient or dishonest emission reduction behaviors and prevent them from becoming an asymptotic stability strategy. In the numerical simulations section, further in-depth analysis will be conducted on how government regulation affects the decision-making of public hospitals and pharmaceutical enterprises, and how these behaviors jointly promote the development of the medical industry toward high-quality and low-carbon emissions, to ensure the smooth realization of medical carbon neutrality.

6 Numerical simulations

To effectively illustrate the influence of critical parameters in the replicator dynamic system on the evolutionary processes and results of multi-agent games, we use MATLAB 2021b to conduct numerical simulations of each agent's evolutionary paths.

Considering the research hypothesis of this article, based on evolutionary stable strategies (strict government regulation, weak enforcement in public hospitals, self-regulation of pharmaceutical companies, and public participation) and convergence conditions (Cg−αDhkg < 0,DhChSh+Rh+2kh < 0,−Dc+Cc+Sc−2kc < 0,(α−1)Dh+Cp < 0). And combined with the real situation, we have made initial settings for the model parameters, as follows:

The government regulatory constraint conditions Cg−αDhkg < 0 indicates that in the context of public participation, the increase in government credibility kg and the sum of fines αDh paid by public hospitals under non-self regulatory strategies are greater than the cost Cg paid by the government in strict regulation. In other words, the net benefit of strict regulation is higher than that of loose regulation, so the government is more inclined to adopt strict regulation. Based on this, we set kg = 5,Cg = 6.

The constraint condition DhChSh+Rh+2kh < 0 for the implementation strategy of public hospitals indicates that under strict government supervision, if public hospitals adopt weak implementation of medical carbon neutrality policies, their short-term benefits Sh minus the fines paid Dh are greater than the benefits obtained when adopting strong implementation strategies Rh plus the credibility gained through public participation 2kh, minus the cost of strong implementation strategies Cg. At this point, the net profit from weak execution is higher than that from strong execution, and hospitals tend to lean toward weak execution. According to constraints Cg−αDhkg < 0 and DhChSh+Rh+2kh < 0, we set Dh = 6,Ch = 4,Sh = 10,Rh = 2,kh = 2.

The constraint condition −Dc+Cc+Sc−2kc < 0 for the self-discipline strategy of pharmaceutical enterprises indicates that under strict government regulation, if pharmaceutical enterprises adopt non-self-discipline strategies, their short-term benefits Sc minus fines Dc are less than the credibility obtained under self-discipline strategies 2kh minus self-discipline costs Cc. Therefore, the net benefits of self-discipline are higher, and companies tend to lean toward self-discipline. Based on this, we set Dc = 6, Cc = 3, Sc = 5, kc = 6.

The constraint condition for public participation (α−1)Dh+Cp < 0 indicates that under strict government regulation, the cost of public participation Cg is less than the compensation received, resulting in a higher net benefit of participation than non-participation. Based on this and combining Cg−αDhkg < 0 with this condition, we set Cp = 1,α = 0.5. In addition, according to reference (54), in the field of environmental protection, the government's benefits Rg under reasonable regulation are greater than its regulatory costs Cg, and therefore Rg = 10 is set. The conclusion of reference (55) is that when corporate self-regulation is combined with strict government regulation, the benefits Rc obtained are greater than the self-regulation costs Cc. Based on this, Rc = 4 is set. The above settings are shown in Table 6.

Table 6
www.frontiersin.org

Table 6. Parameter value.

Finally, in order to observe how the dynamic process of evolutionary game deviates from equilibrium due to parameter influence, we set the initial strategy probabilities of each player to x = 0.5, y = 0.5, z = 0.5, w = 0.5; And numerical simulation was conducted through 100 iterations.

6.1 Impact of the cost of strict government supervision

The cost of strict government supervision (Cg) was set to 6, 9, and 12 to examine the sensitivity of stakeholders' strategic evolution to changes in regulatory expenditure. These values were chosen to span low, medium, and high-cost scenarios, ensuring comparability with other parameters of similar scale as in the initial settings. Such a range enables the identification of threshold effects predicted by the theoretical analysis, particularly the transition between strict and loose regulatory strategies. The selection also facilitates intuitive visualization of evolutionary trajectories and supports reproducibility of results under varying cost conditions. Figure 6 shows the evolution process and results of the strategies of the main players in the four-party game.

Figure 6
Three line graphs labeled (A), (B), and (C), showing probabilities over time (t) for different organizations: Government (red), Public Hospitals (green), Pharmaceutical Enterprises (blue), and Public (cyan). Graph (A) has stable high probabilities. Graph (B) shows oscillating probabilities. Graph (C) has probabilities quickly dropping then leveling. Each graph title varies: Cy 6, Cy = 9, Cy = 12.

Figure 6. The impact of strict supervisory costs Cg on the evolution of strategies of various parties. (A) Cg = 6. (B) Cg = 9. (C) Cg = 12.

Figure 6 shows that as the cost of strict government regulation increases, the probability of the government adopting a strict regulation approach exhibits a consistently declining trend, ultimately stabilizing at 0. This phenomenon illustrates that, due to severe cost restrictions, the government's capacity for strict regulation diminishes, making it difficult to maintain consistent and successful regulatory initiatives. Simultaneously, as the likelihood of government regulation diminishes, the probability of public engagement in its capacity as a “watchdog” exhibits an increasing tendency, nearing 1 during the stabilization period. This implies that when the public perceives that the government is becoming less inclined toward strict regulation, they will compensate for the regulatory deficit by increasing their participation. Thus, establishing collective supervision of the medical carbon neutrality process. However, when the cost of government regulation rises to the point where the public believes that the government can no longer adopt strict regulatory strategies, the motivation for active public participation will sharply decrease and eventually stabilize at 0. The “synergistic failure” impact indicates that the implementation of carbon neutrality in the healthcare industry will be significantly obstructed if both governmental regulatory willingness and public participation are minimal. Consequently, the government ought to concentrate on refining the design and execution of the regulatory framework, while minimizing regulatory costs through enhanced efficiency, judicious resource allocation, and the incorporation of technological solutions. Simultaneously, it should establish a multi-faceted incentive structure and fortify institutionalized avenues for public engagement and feedback loops, thereby augmenting public motivation for ongoing oversight. Only with the dual drive to reduce regulatory costs and activate public participation can a stable synergistic evolution between the government and the public be formed. This, in turn, enables the successful realization of carbon neutrality in healthcare and provides a replicable empirical model for public health governance.

6.2 The impact of government social credibility on the strategic choices of various parties

To investigate how government credibility affects strategic evolution, we set the added value of public credit kg to 1, 3, and 5 based on initial parameters. When kg = 1, strict regulation brings only a small credibility gain. This tests the stability of strategies under weak public endorsement. When kg = 3, the credibility gain is moderate. This reflects common real-world cases and captures the transition from low to high credibility. When kg = 5, strict regulation greatly boosts public trust. This represents a strong public support scenario. These settings allow an analysis of how different credibility levels shape the evolutionary paths of the government, public hospitals, pharmaceutical enterprises, and the public in achieving medical carbon neutrality. Figure 7 shows the evolution process and results of the strategies of the main players in the four-party game.

Figure 7
Three line graphs (A, B, C) show probability trends over time (labeled as t) for different institutions: Government, Public Hospitals, Pharmaceutical Enterprises, and Public. Graph A (ks = 1) shows fluctuating trends; Graph B (ks = 3) shows convergence after t=3; Graph C (ks = 5) shows rapid convergence by t=2. Each graph includes a legend for the represented institutions.

Figure 7. The impact of government social credibility kg on the strategic choices of various parties. (A) kg = 1. (B) kg = 3. (C) kg = 5.

Figure 7 illustrates that as public trust in the government's strict regulatory approach increases during the medical carbon neutralization process, the system's evolutionary dynamics exhibit a pronounced “two-way enhancement” effect: firstly, the likelihood of the government adopting a strict regulatory strategy rises consistently with increasing social credibility, ultimately approaching 1; secondly, the probability of pharmaceutical enterprises implementing self-discipline strategies also escalates rapidly in the high-pressure regulatory environment, driven by a substantial rise in not self-discipline costs, eventually stabilizing at 1. This is because as public trust increases, the marginal benefit of the government's choice of strict regulation continuously exceeds that of loose regulation, which not only effectively safeguards the public interest, but also enhances the legitimacy of regulation with the support of public opinion. Conversely, when strict government regulation becomes the dominant strategy, pharmaceutical enterprises that choose the not self-regulation strategy will face higher risks of penalties and reputational losses, and their expected returns are significantly lower than those of the self-discipline strategy, thus tilting toward self-discipline in their strategy evolution. Consequently, improving the public's social credibility in government regulators is a crucial impetus for encouraging the government to proactively implement a strict regulatory strategy. Additionally, it is essential to enhance information transparency and societal regulation of the regulatory process, refine the public participation evaluation system, and utilize public opinion to inform and incentivize mechanisms that promote self-discipline among pharmaceutical enterprises. This ensures that the goal of healthcare carbon neutrality is efficiently achieved through the evolution of the four-party interaction between the government, public hospitals, pharmaceutical companies and the public.

6.3 The impact of short-term financial profit on the strategic choices of various parties

In our model, the tendency of public hospitals to choose weak enforcement is governed by DhChSh+Rh+2kh < 0 and Cg−αDhkg < 0 The former condition reflects that when the short-term revenue Sh of public hospitals is high, after deducting the fine Dh and the cost Ch under the strong enforcement strategy, their net revenue may still exceed the sum of the direct revenue Rh and reputation gain kh brought by strong enforcement, thereby prompting hospitals to lean toward weak enforcement; The latter condition ensures that the government still chooses strict regulation when the regulatory benefits outweigh the costs. Based on this, we set Sh to 4, 7, and 10, representing low, medium, and high short-term return scenarios, respectively, covering three typical states of not meeting, approaching, and significantly meeting weak enforcement conditions. This facilitates sensitivity analysis and critical effect identification, to better analyze the evolutionary path of the four party game subjects in the process of implementing medical carbon neutrality. Figure 8 shows the evolution process and results of the strategies of the main players in the four-party game.

Figure 8
Three line graphs labeled (A), (B), and (C) show probability trends over time for

Figure 8. The impact of short-term financial profit Sh on the strategic choices of various parties. (A) Sh = 4. (B) Sh = 7. (C) Sh = 10.

The evolutionary trajectory depicted in Figure 8 illustrates that as the short-term benefits of public hospitals adopting weak enforcement strategies for medical carbon neutrality policies rise the probability of their inclination toward strong enforcement strategies diminishes. This ultimately results in evolutionary stability with weak enforcement strategies. During this process, a dual compensation mechanism is established within the system: “weak enforcement by public hospitals necessitates strict government regulation, and strict government regulation fosters heightened public participation.” On one hand, the rising likelihood of weak enforcement by public hospitals jeopardizes the public interest in healthcare carbon emission management, prompting the government to adopt a strict regulatory strategy, ultimately stabilizing the probability of its regulation at 1. On the other hand, the public, after observing the increased willingness of the government to regulate, correspondingly increased their motivation to participate in social supervision, and their probability of participation eventually stabilized at 1 as well. At the same time, under the dual pressure of government regulation and public opinion, pharmaceutical enterprises are experiencing a reversal in the cost-benefit structure of self-discipline enforcement, and the probability of adopting self-discipline strategies is gradually approaching 1. The outcomes of the aforementioned evolutionary game suggest that the government should implement precise incentive and constraint policies. These policies aim to tighten the short-term revenue space of public hospitals' weak performance. This also can be achieved through financial incentives, performance evaluations, and other measures. Simultaneously, the government and the public should strengthen their joint oversight of pharmaceutical enterprises downstream in the healthcare supply chain. This can be achieved through information disclosure, third-party evaluation, and other means. Only under the triple driving force of weakening the short-term interests of hospitals, strengthening the linkage between government and the public, and guiding the self-discipline of enterprises can the synergistic evolution system of carbon neutrality in the healthcare industry be firmly built. This, in turn, ensures that the goal of carbon neutrality is achieved on schedule.

6.4 The impact of penalty fines on the strategic choices of various parties when the pharmaceutical enterprise selects no self-discipline

Based on the constraint −Dc+Cc+Sc−2kc < 0, the fines Dc imposed on pharmaceutical enterprises for failing to self-discipline are set at 2, 4, and 6. Different fine levels allow analysis of how penalty intensity influences the evolution of pharmaceutical strategies. Low fines may fail to incentivize self-discipline, while moderate adjustments help identify an effective range that promotes self-regulation, ensuring reasonable and meaningful parameter settings. This approach enables the model to better reflect real regulatory environments and enhances the practical relevance of the results. The evolution process and results of the strategies of the main players in the four-party game are shown in Figure 9.

Figure 9
Three line charts labeled A, B, and C show probability curves for different entities: Government (x), Public Hospitals (y), Pharmaceutical Enterprises (z), and Public (w). The x-axis represents the evolution time, and the y-axis shows probability from 0 to 1. Each chart has a different parameter D0 value: 2, 4, and 6 respectively. All entities show a similar distribution pattern with varying probabilities across f-values in each subplot.

Figure 9. The impact of penalty fine Dc on the strategic choices of various parties when the pharmaceutical enterprise selects no self-discipline. (A) Dc = 2. (B) Dc = 4. (C) Dc = 6.

Figure 9 illustrates that following the government's augmentation of fines for pharmaceutical enterprises employing not self-discipline strategies, the probability of self-discipline among these enterprises will markedly rise, ultimately nearing 1. Consequently, when pharmaceutical enterprises implement not self-discipline tactics, the government can significantly diminish the likelihood of such behavior by increasing fines against pharmaceutical enterprises. This initiative can facilitate the attainment of medical carbon neutrality.

6.5 The impact of pharmaceutical enterprises' short-term profit on the strategic choices of various parties

Based on the constraint −Dc+Cc+Sc−2kc < 0, the short-term benefit Sc together with fines Dc, costs Cc, and other regulatory parameters determine the net payoff for pharmaceutical enterprises when they fail to self-discipline. To simulate firms' behavioral responses under different economic incentives, the short-term benefits Sc for non-self-discipline are set at 5, 9, and 13. These values cover a range from low to high economic gains, reflecting possible benefit levels under varying policy environments. A lower Sc indicates limited gains from non-compliance, which may not outweigh corresponding penalties and costs, whereas a higher Sc represents greater short-term incentives for firms to choose non-self-regulation. Setting these three levels allows systematic analysis of pharmaceutical enterprises' strategic evolution under different benefit scenarios, revealing the dynamic interplay between policy pressures and firm behavior, and providing theoretical support and strategic guidance for the implementation of carbon neutrality policies. Figure 10 shows the evolution process and results of the strategies of the main players in the four-party game.

Figure 10
Three line graphs labeled (A), (B), and (C) show probability distributions for different entities: government, public hospitals, pharmaceutical enterprises, and the public. The x-axis represents the evolution time, and the y-axis shows probability from zero to one. As Sc increases from five to thirteen, the curves for government and pharmaceutical enterprises approach a probability of one more rapidly, while public hospitals and the public maintain low probabilities.

Figure 10. The impact of pharmaceutical enterprises' short-term profit Sc on the strategic choices of various parties. (A) Sc = 5. (B) Sc = 9. (C) Sc = 13.

As shown in Figure 10, with the increase of short-term benefits of not self-discipline in pharmaceutical enterprises, the benefits of not self-discipline gradually exceed those of self-discipline. This makes pharmaceutical enterprises more interested in pursuing short-term profits, resulting in the probability of pharmaceutical enterprises taking self-discipline strategies gradually decreasing and eventually stabilizing at 0. However, the government maintained strict regulation during this process, and the public participated. Yet this still leads pharmaceutical enterprises to tend toward not self-discipline strategies. The reason is that the government and the public impose insufficient punishment and fail to cause a significant enough loss of reputation when pharmaceutical enterprises adopt no self-discipline. Therefore, the government should take stronger measures to weaken the short-term benefits obtained by pharmaceutical enterprises when they lack self-discipline. Similarly, the public should participate in the supervision process through more channels, strengthening the public opinion pressure on not self-disciplined pharmaceutical enterprises. This will impel pharmaceutical enterprises to adopt self-discipline strategies and participate in implementing medical carbon neutrality strategies.

6.6 The impact of the reputation of pharmaceutical enterprises on the strategic choices of various parties

Based on the constraint−Dc+Cc+Sc−2kc < 0, the credibility enhancement factor kc plays a critical role in balancing the net payoff for pharmaceutical enterprises. Setting kc at 2, 4, and 6 allows the model to capture varying levels of credibility incentives that influence firms' strategic behavior. These values represent a gradually increasing scale of credibility reinforcement, which is essential to analyze how changes in the perceived reputation benefits affect the decision-making process of pharmaceutical enterprises. A lower kc implies weaker credibility incentives, making self-regulation less attractive, while higher values reflect stronger incentives that encourage compliance. By selecting these three representative values, the model can systematically explore the threshold effects of credibility incentives on evolutionary dynamics, providing theoretical support and decision-making references for the implementation of medical carbon neutrality policies. Figure 11 shows the evolution process and results of the strategies of the main players in the four-party game.

Figure 11
Three plots labeled A, B, and C show probability versus variable f for different actors in a network: Government, Public Hospitals, Pharmaceutical Enterprises, and Public. Each plot shows varying degrees of probability change as f increases from zero to ten. Labels indicate node counts: kc = 2 for A, kc = 4 for B, and kc = 6 for C, affecting the curve shapes and steepness.

Figure 11. The impact of the reputation of pharmaceutical enterprises kc on the strategic choices of various parties. (A) kc = 2. (B) kc = 4. (C) kc = 6.

As shown in Figure 11, the probability of pharmaceutical enterprises choosing self-discipline increases significantly when public participation confers higher social credibility on these enterprises. This higher credibility encourages them to adopt a self-discipline strategy. Eventually, the probability stabilizes toward 1. Specifically, within the evolutionary game framework, social credibility functions as an external incentive for pharmaceutical enterprises. It drives an increase in the marginal gains they achieve by adopting a self-discipline strategy beyond what they would earn not self-discipline. As an increasing number of pharmaceutical enterprises implement self-discipline tactics, public acknowledgment and social reputation concurrently enhance. This growth in reputation has had a positive feedback effect, further reinforcing the superiority and attractiveness of self-discipline strategies adopted by pharmaceutical enterprises. Based on this, it is advised that the government and industry associations enhance public knowledge regarding the sustainable development accomplishments of pharmaceutical enterprises. And through scientific assessment, self-regulatory commitment disclosure and third-party certification, the public's attention and trust in the green governance practices of pharmaceutical enterprises will be enhanced. On this basis, pharmaceutical enterprises can be guided to transform their accumulated social reputation capital into long-term competitive advantages. This, in turn, enables them to more actively adopt self-discipline strategies in the healthcare carbon-neutralization process and achieve a win-win outcome for both environmental and economic benefits.

6.7 The impact of compensation proportion on the strategic choices of various parties

The proportion of public compensation 1−α reflects the government's support for public participation in healthcare carbon neutrality and the role of compensation mechanisms in motivating engagement. Based on the constraint (α−1)Dh+Cp < 0, the compensation ratio must ensure that the public receives sufficient compensation to offset participation costs Cp and maintain motivation, while avoiding excessive compensation that could lead to resource inefficiency. Therefore, the public compensation ratio is set at 0.1, 0.15, and 0.2, covering a range from low to moderate compensation levels. This range guarantees reasonable returns for the public to encourage active participation in carbon neutrality efforts while maintaining economic sustainability. By establishing these three representative values, the model can effectively reveal how varying compensation levels influence the evolutionary strategies of the public, thereby offering theoretical foundations and policy recommendations for designing compensation mechanisms in healthcare carbon neutrality policies. Figure 12 shows the evolution process and results of the strategies of the main players in the four-party game.

Figure 12
Three line graphs labeled (A), (B), and (C) compare probabilities over 100 time units with parameters 1 minus alpha equals 0.1, 0.15, and 0.2. Each graph includes four lines for Government, Public Hospitals, Pharmaceutical Enterprises, and Public, showing different oscillation patterns and stabilization trends across the graphs.

Figure 12. The impact of compensation ratios on the strategic choices of various parties. (A) 1 − α= 0.1. (B) 1 − α= 0.15. (C) 1 − α= 0.2.

Figure 12 shows that as the compensation proportion for public participation increases, people obtain higher marginal benefits from participating in medical carbon neutrality. Consequently, their participation probability rises significantly and eventually stabilizes at 1. Under the combined effects of increased public participation and strong government regulation, the net benefits pharmaceutical enterprises obtain from a self-discipline strategy begin to exceed those from not self-discipline. This change in benefits compels pharmaceutical enterprises to progressively implement self-discipline measures, ultimately stability at 1. Based on the above analyses, it is suggested that the government should take multi-dimensional measures: on the one hand, it should guide the public to participate in monitoring in a continuous and rational manner by dynamically adjusting the compensation ratio and setting an incentive cap; on the other hand, a third-party assessment organization can be introduced. This organization would publicize the performance of pharmaceutical enterprises in carbon emission reduction and green governance. It could also certify their efforts, thereby strengthening the social credibility effect. In addition, policy tool combinations should be optimized by integrating financial subsidies, tax incentives, and credit endorsements to form a robust benefit–linkage mechanism. This approach will maximize the benefits of public participation and the self-discipline motivation of pharmaceutical enterprises throughout the healthcare carbon-neutralization process, thereby facilitating the achievement of healthcare carbon neutrality.

6.8 The impact of public participation costs on the strategic choices of various parties

Based on the constraint (α−1)Dh+Cp < 0, there is a trade-off relationship between public participation costs Cp and the public compensation ratio 1−α. To ensure that the public receives adequate compensation to maintain their motivation despite bearing participation costs, these costs need to be set within a reasonable range. Therefore, public participation costs are set at 0, 3, and 6, representing scenarios of no cost, moderate burden, and higher burden, respectively. These values cover various participation barriers that the public may realistically face, facilitating analysis of how changes in costs affect public willingness to participate and their strategic evolution. This layered cost setup, consistent with the constraint, enables the model to scientifically reveal the balance mechanism between compensation and cost, thereby providing theoretical support for designing public participation incentives in healthcare carbon neutrality policies. Figure 13 shows the evolution process and results of the strategies of the main players in the quadripartite game.

Figure 13
Three line graphs (A, B, C) show probability over time with varying Ct values: Ct = 0, Ct = 3, and Ct = 6. Each graph includes lines for Government, Public Hospitals, Pharmaceutical Enterprises, and Public, showing similar patterns. All lines mostly start near 1.0 and quickly drop to stabilize near 0.0 by t = 10. The probability remains constant thereafter until t = 100.

Figure 13. The impact of public participation costs Cp on the strategic choices of various parties. (A) Cp = 0. (B) Cp = 3. (C) Cp = 6.

Figure 13 illustrates that the probability of public participation in healthcare carbon neutrality markedly decreases as costs associated with participation increase. In this context, due to the lack of public participation, the regulatory responsibility is left only to the government, and the regulatory pressure is insufficient, resulting in a gradual decrease in the probability of pharmaceutical enterprises adopting self-discipline strategies, ultimately stabilizing at adopting not self-discipline strategies. To reverse this trend, the government should enhance public participation by increasing participation compensation, establishing dynamic incentive mechanisms, optimizing participation processes, and reducing time and economic costs. At the same time, third-party review and public opinion monitoring are being introduced to boost social exposure and public pressure on pharmaceutical enterprises that do not self-regulate. The aforementioned actions can improve public oversight and pharmaceutical enterprises' self-discipline, create a multi-party cooperative regulatory force, and contribute in the achievement of medical carbon neutrality.

7 Model extension

7.1 Analysis of time-delay effect in public hospitals

In this section, we explored the effect of time delays in public hospitals on their evolutionary strategies under evolutionary game theory. We analyzed the influence of time delays on the stability speed of the optimal strategy. To better illustrate the conclusions of this subsection, we introduce parameter τh into the differential equation to represent the time delay effect in the game of public hospitals. The specific model is as follows:

dx(t)dt=x(t)(1x(t)){w(t)kgCgy(tτh)               [(αw(t)w(t)+1)Dh+z(t)(Rc+Rhw(t)Cg)]       (w(t)αw(t)1)(Dc+Dh)z(t)[(α1)w(t)+1]Dc}    (34)
dy(t)dt=y(t-τh)[1-y(t-τh)]   [x(t)(Dh+z(t)Rh)+2w(t)kh-Sh-Ch]     (35)
dz(t)dt=z(t)[1-z(t)]   [x(t)(Dc+y(t-τh)Rc)+2w(t)kc-Sc-Cc]     (36)
dw(t)dt=w(t)[1w(t)]               {x(t)(1α)(Dc+Dh)x(t)[(1α)y(tτh)Dh               +(1α)z(t)Dc]Cp}    (37)

As illustrated in the system of equations, the probability y of public hospitals adopting a strong enforcement strategy becomes a function of t−τh, where τh represents the time delays parameter. To validate this model, six comparative experiments were conducted, with results summarized in Figures 14AF.

Figure 14
Three graphs labeled A, B, and C depicting the decrease of spreaders over time with different initial conditions. Each graph has four colored lines representing different initial spreads: red (n0=0.06), orange (n0=0.04), green (n0=0.02), and blue (n0=0). The y-axis is labeled

Figure 14. The impact of time delays on public hospitals. (A) (1,0,0,0) Stable point. (B) (1,0,0,1) Stable point. (C) (1,0,1,0) Stable point. (D) (1,0,1,1) Stable point. (E) (1,1,0,0) Stable point. (F) (1,1,0,1) Stable point.

As seen from the above figures, the time delays do not change the strategic evolution direction of public hospitals at different equilibrium points. This indicates that, in the evolutionary game process, time delays do not affect the final strategy selection of public hospitals. Interestingly, as the delay time increases, the speed at which public hospitals converge to the optimal strategy selection becomes faster. This is because, under time delays, the game process among players transforms from a simultaneous static game into a dynamic game with temporal lags. In this context, public hospitals under different equilibrium conditions can observe the strategic selections and stabilization trends of the other three game entities, thereby gaining an informational advantage to adopt strategies that maximize their payoffs. Compared with delay-free scenarios, public hospitals thus stabilize their strategies more rapidly. Notably, with longer delay times, the evolutionary trends of other players become more pronounced and predictable, enabling public hospitals to accelerate their convergence toward the optimal strategy.

Following the time-delay model results, it is notable that such delays—whether arising from policy implementation lags or behavioral inertia—do not destabilize the system. A plausible explanation is that the inherent strategic incentives embedded in the payoff structure are sufficiently robust to preserve equilibrium stability, even when decision-making is temporally misaligned. However, the influence of delays is not necessarily uniform across all stakeholders. For example, government regulators, pharmaceutical enterprises, and the public may respond differently to administrative delays, cognitive delays, or resource allocation delays, leading to varied adjustment speeds and convergence paths. Future research could systematically classify these delay types and investigate their heterogeneous impacts on decision-making dynamics.

From a policy and management perspective, mitigating the adverse effects of delays and enhancing system performance requires a multi-level, cross-sectoral, and integrated coordination framework. First, an adaptive incentive adjustment mechanism based on real-time behavioral feedback should be established. By regularly collecting and analyzing behavioral data from all stakeholders in strategic choices, policy makers can dynamically modify reward and punishment levels to correct deviations in policy implementation. This approach minimizes the risk of incentive obsolescence caused by information lags and ensures that the policy continues to exert a sustained influence on stakeholder behavior. Second, administrative procedures should be streamlined to reduce institutional bottlenecks. This includes clarifying roles and responsibilities at the policy formulation stage, simplifying approval processes, and accelerating cross-departmental information sharing to lower the time cost of policy implementation. Finally, a multi-stakeholder engagement platform should be established, involving government bodies, public hospitals, pharmaceutical enterprises, public representatives, and independent professional institutions. Through regular consultations, transparent information disclosure, and joint decision-making, such a platform can foster early consensus on policy goals and action pathways, thereby reducing behavioral inertia and resistance. These integrated measures can ensure that not only public hospitals but all stakeholders converge toward socially optimal strategies within a reasonable time frame, ultimately enhancing the timeliness, coordination, and sustainability of low-carbon public health policy formulation and implementation.

7.2 Analysis of the time-delay effect in pharmaceutical enterprises

In this section, we discuss the impact of time delays on pharmaceutical enterprises' evolutionary strategies in evolutionary games and analyze the influence of time delays on the stability and speed of the optimal strategy. To better illustrate the conclusions of this subsection, we introduce parameter τc into the differential equation to represent the effect of time delays in the game of pharmaceutical enterprises. The specific model is as follows:

dx(t)dt=x(t)[1x(t)]{w(t)kgCgy(t)[(αw(t)w(t)+1)             Dh+z(tτc)(Rc+Rhw(t)Cg)](w(t)αw(t)1)                 (Dc+Dh)z(tτc)[(α1)w(t)+1]Dc}     (38)
dy(t)dt=y(t)[1y(t)][x(t)            (Dh+z(tτc)Rh)+2w(t)khShCh]     (39)
dz(t)dt=z(tτc)[1z(tτc)][x(t)(Dc+y(t)Rc)            +2w(t)kcScCc]     (40)
dw(t)dt=(t)[1w(t)]{x(t)(1α)(Dc+Dh)x(t)              [(1α)y(t)Dh+(1α)z(tτc)Dc]Cp}     (41)

From the above equations, it can be seen that at the not self-discipline probability z of pharmaceutical enterprises becomes a function of t−τc. In order to better illustrate the above model conclusions, six groups of comparative experiments were carried out, as shown in Figures 15AF.

Figure 15
Six line graphs labeled (A) to (F), all showing

Figure 15. The impact of time delays on pharmaceutical enterprises. (A) (1,0,0,0) Stable point. (B) (1,0,0,1) Stable point. (C) (1,0,1,0) Stable point. (D) (1,0,1,1) Stable point. (E) (1,1,0,0) Stable point. (F) (1,1,0,1) Stable point.

Figure 15 analysis reveals that time delays do not alter the strategic evolution direction of pharmaceutical enterprises across equilibrium points, indicating that such delays do not affect their final strategy selection in the evolutionary game. Compared to delay-free scenarios, the speed at which pharmaceutical enterprises converge to the optimal strategy increases as delay time lengthens. This acceleration, consistent with the mechanism explained in Section 7.1 for public hospitals, arises from the ability of enterprises to observe and anticipate the stabilization trends of other players before finalizing their own strategy choices.

Notably, although the delay effect does not destabilize the equilibrium, its impact on pharmaceutical enterprises is shaped by industry-specific factors such as R&D investment cycles, regulatory approval timelines, and production capacity adjustments. Different types of delays—such as technology development lags, supply chain disruptions, or market demand recognition delays—may alter the responsiveness of enterprises and, consequently, the system's adaptation speed.

To enhance strategic responsiveness under such delays, pharmaceutical enterprises could focus on strengthening internal early-warning systems to detect market and policy signals in advance, diversifying supply chain structures to reduce vulnerability to disruptions, and employing flexible production planning to shorten adjustment lead times. These enterprise-level strategies complement but differ from system-wide policy measures, targeting the unique operational rhythms and risk structures of the pharmaceutical industry.

7.3 Analysis of the time-delay effect in public

In this section, we discuss the impact of time delays on the public's evolutionary strategies in evolutionary games and analyze the influence of time delays on the stability and speed of the optimal strategy. To better illustrate the conclusions of this subsection, we introduce the parameter τp into the differential equation to represent the effect of time delays in the game of the public. The specific model is as follows:

dx(t)dt=x(t)(1x(t)){w(tτp)kgCgy(t)[(αw(tτp)             w(tτp)+1)Dhz(t)(Rc+Rhw(tτp)Cg)]             [w(tτp)αw(tτp)1)(Dc+Dh)z(t)[(α1)              w(tτp)+1]Dc}     (42)
dy(t)dt=y(t)[1y(t)][x(t)(Dh+z(t)Rh)            +2w(tτp)khShCh]     (43)
dz(t)dt=z(t)[1z(t)][x(t)(Dc+y(t)Rc)              +2w(tτp)kcScCc]     (44)
dw(t)dt=w(tτ)[1w(tτp)]{x(t)(1α)(Dc+Dh)             x(t)[(1α)y(t)Dh+(1α)z(t)Dc]Cp}     (45)

From the above system of equations, it can be observed that the participation probability w chosen by public now becomes a function of t−τp. To better illustrate the above model conclusions, six groups of comparative experiments were carried out, as shown in Figures 16AF.

Figure 16
Graphs A to F show the evolution of a spread factor over time with varying parameters τ0. Each chart compares results for τ0 = 0.06, τ0 = 0.04, τ0 = 0.02, and τ0 = 0 using distinct colors. Graphs A, C, and E display decreasing trends, while B, D, and F depict increasing trends. Graph F includes a zoomed section highlighting detail around time 2.2 to 2.35.

Figure 16. The impact of time delay on the public. (A) (1,0,0,0) Stable point. (B) (1,0,0,1) Stable point. (C) (1,0,1,0) Stable point. (D) (1,0,1,1) Stable point. (E) (1,1,0,0) Stable point. (F) (1,1,0,1) Stable point.

The phase diagram analysis indicates that time delays do not alter the evolutionary direction of public strategies across equilibrium points, suggesting that such delays do not influence the ultimate strategic choice of the public in the evolutionary game. However, compared with the delay-free scenario, the public exhibits a faster convergence to the optimal strategy when time delays are present, and this acceleration becomes more pronounced as the delay period increases. This effect may be attributed to the public's enhanced ability to interpret and internalize information signals from other stakeholders during the delayed decision process, thereby reducing hesitation in adopting beneficial strategies. Different forms of delay—such as the dissemination lag of policy information and the adaptation lag in behavioral patterns—may have varying effects on decision-making speed and stability. Addressing these challenges could involve strengthening transparent information channels, improving the timeliness of public engagement campaigns, and facilitating gradual behavioral transitions through staged incentives, so that the public's strategic alignment remains both rapid and robust without destabilizing the system.

8 Conclusions

This study systematically analyzed the strategic evolution behavior of the government, public hospitals, pharmaceutical enterprises, and the public in the process of medical carbon neutrality by constructing a four-party evolutionary game model. Numerical simulations were performed using MATLAB 2021b to analyze critical parameters such regulatory costs and credibility, incentives and penalties, public compensation ratios, and participation costs. Time delay effects were incorporated into the model to elucidate the dynamic influence of these factors on the gains and losses of all stakeholders concerned. Considering the aforementioned analysis, we have summarized four critical conclusions that serve as significant references for optimizing carbon neutrality policies within the healthcare sector and improving the efficacy of governmental regulation and public engagement.

(1) The probability of strict government regulation is both constrained by regulatory costs and closely related to government credibility. To promote the effective implementation of carbon neutrality in the healthcare industry, the government should establish and strengthen the regulatory mechanism for energy conservation and emission reduction, promote the adoption of emission reduction strategies in public hospitals, and regulate the behavior of pharmaceutical enterprises. On the one hand, market-oriented means such as carbon trading can be used to reduce the cost of direct government regulation, allowing public hospitals and pharmaceutical enterprises to participate in carbon emission management independently through market trading. On the other hand, the disclosure process of medical carbon neutrality information should be simplified and optimized to improve process efficiency, ensure timely and accurate transmission of supervisory information to the public, and thereby enhance the government's image and credibility. Through the dual drive of marketization and informatization, continuously improving the regulatory system will help promote the smooth implementation of medical carbon neutrality.

(2) The decision-making motives of public hospitals in adopting a weak enforcement strategy are mainly rooted in short-term benefits maximization. Similarly, when choosing a strategy of self-regulation or not self-regulation, pharmaceutical enterprises will weigh not only their short-term gains from not self-regulation, but also the cost of fines that they may face under strict government regulation, as well as the reputational risks associated with public engagement. To this end, the government needs to improve the disciplinary mechanism while building a regulatory system for energy conservation and emission reduction in the healthcare industry. Fines and notification of criticism should be applied to public hospitals that do not execute the healthcare carbon neutralization policy with vigor. For pharmaceutical enterprises that adopt not self-discipline strategies, their institutional costs should be increased by means of revoking their business licenses and incorporating their bad behavior into public credit information platforms. Through a multi-dimensional incentive-constraint mechanism, public hospitals and pharmaceutical enterprises will be promoted to actively fulfill their energy-saving and emission reduction responsibilities. This will ensure that the healthcare carbon-neutrality target is achieved.

(3) Public participation is of great significance in the carbon neutrality process of the healthcare industry and is influenced by both compensation ratios and participation costs. On the one hand, active public participation enhances the green growth of the healthcare industry and increases the relevance and efficacy of government energy conservation and emission reduction programs. On the other hand, the public plays a key role in supervising the behavior of public hospitals and pharmaceutical companies, and their public opinion and reputation leverage can motivate industry entities to continue implementing energy-saving and emission reduction measures. In order to facilitate the green transformation of the medical and health sector, the government should moderately raise the compensation ratio and make full use of digital tools to expand low-cost participation channels like online consultation, electronic questionnaires, and online forums.

(4) Even though public hospitals, pharmaceutical enterprises, and the public may exhibit differentiated time delays in policy acceptance and execution due to their inherent operational characteristics, these factors will not impede their optimal strategy selection.

Our study establishes an evolutionary game-theoretic model based on perfect rationality and expected utility theory, incorporating time delay factors to investigate strategic interactions among game entities within the medical industry's carbon neutrality context. Our analysis examines explicitly how time delays influence the strategic decision-making processes of various stakeholders. However, real-world scenarios suggest that, beyond time delays, subjective factors such as emotional states and risk attitudes significantly impact agents' strategic choices. These psychological dimensions can induce stochastic mutations during evolutionary processes, potentially accelerating convergence to optimal strategies. Consequently, future research could extend this model by integrating stochastic perturbation terms and implementing multi-stage repeated game analyses.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

BX: Conceptualization, Data curation, Methodology, Validation, Writing – original draft. XG: Formal analysis, Software, Validation, Writing – review & editing. YC: Supervision, Visualization, Writing – review & editing. MN: Investigation, Writing – review & editing. SL: Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by Soft science project of Henan Province (grant No.:252400411280 and 252400410046), Philosophy and Social Science Planning Project of Henan Province (grant No.: 2022CZH015), Philosophy and social science education strong provincial project of Henan Province (grant No.:2025JYQS1079), Henan Provincial Postgraduate Course Ideological and Political Demonstration Course Project (grant No.: YJS2024SZ14), National Social Science Foundation of China (grant No.:20BJL135), and 2024 funding project for international training of high-level talents in Henan Province: research on the path and countermeasures of supply chain digital transformation under the background of big data.

Acknowledgments

The authors are grateful to the editor and reviewers for their critical suggestions of this manuscript Meanwhile, we would like to thank GhatGPT4.0 for improving the language of this article in the final stage. After using the tool, the authors reviewed and edited the content as needed and took full responsibility for the content of the publication.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that Gen AI was used in the creation of this manuscript. During the preparation of this work, the author(s) used ChatGPT 4.0 in the late stage of the process to improve the language of the article. After using the tool, the author(s) reviewed and edited the content as needed and took full responsibility for the content of the publication.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Rocque RJ, Beaudoin C, Ndjaboue R, Cameron L, Poirier-Bergeron L, Poulin-Rheault RA, et al. Health effects of climate change: an overview of systematic reviews. BMJ Open. (2021) 11:e046333. doi: 10.1136/bmjopen-2020-046333

PubMed Abstract | Crossref Full Text | Google Scholar

2. Kjellstrom T, McMichael AJ. Climate change threats to population health and well-being: the imperative of protective solutions that will last. Glob Health Action. (2013) 6:20816. doi: 10.3402/gha.v6i0.20816

PubMed Abstract | Crossref Full Text | Google Scholar

3. Pichler PP, Jaccard IS, Weisz U, Weisz H. International comparison of health care carbon footprints. Environ Res Lett. (2019) 14:064004. doi: 10.1088/1748-9326/ab19e1

Crossref Full Text | Google Scholar

4. Eckelman MJ, Sherman J. Environmental impacts of the US health care system and effects on public health. PLoS ONE. (2016) 11:e0157014. doi: 10.1371/journal.pone.0157014

PubMed Abstract | Crossref Full Text | Google Scholar

5. Torjesen I. NHS aims to become world's first “net zero” health service by 2040. BMJ. (2020) 371:m3856. doi: 10.1136/bmj.m3856

PubMed Abstract | Crossref Full Text | Google Scholar

6. University Hospitals Birmingham: A World First in Carbon Net Zero Surgery. (2022). Available online at: http://www.england.nhs.uk (Accessed October 16, 2024).

Google Scholar

7. Gammie AJ, Lopez JB, Scott S. Imperative: reducing the environmental impact of clinical laboratories. Clin Chem Lab Med. (2023) 61:634–7. doi: 10.1515/cclm-2022-1052

PubMed Abstract | Crossref Full Text | Google Scholar

8. Holmner Å, Ebi KL, Lazuardi L, Nilsson M. Carbon footprint of telemedicine solutions-unexplored opportunity for reducing carbon emissions in the health sector. PLoS ONE. (2014) 9:e105040. doi: 10.1371/journal.pone.0105040

PubMed Abstract | Crossref Full Text | Google Scholar

9. Pejić Bach M, Mihajlović I, Stanković M, Khawaja S, Qureshi FH. Determinants of intention to use of hospital information systems among healthcare professionals. Systems. (2024) 12:235. doi: 10.3390/systems12070235

Crossref Full Text | Google Scholar

10. Cecchini G, DePass M, Baspinar E, Andujar M, Ramawat S, Pani P, et al. Cognitive mechanisms of learning in sequential decision-making under uncertainty: an experimental and theoretical approach. Front Behav Neurosci. (2024) 18:1399394. doi: 10.3389/fnbeh.2024.1399394

PubMed Abstract | Crossref Full Text | Google Scholar

11. Global healthcare leaders accelerate industry decarbonization ahead of COP28. (2023). Available online at: https://www.businesswire.com/news/home/20231123023210/zh-CN?utm_campaign=shareaholic&utm_medium=copy_link&utm_source=bookmark (Accessed October 15, 2024).

Google Scholar

12. Keil M, Frehse L, Hagemeister M, Knieß M, Lange O, Kronenberg T, et al. Carbon footprint of healthcare systems: a systematic review of evidence and methods. BMJ Open. (2024) 14:e078464. doi: 10.1136/bmjopen-2023-078464

PubMed Abstract | Crossref Full Text | Google Scholar

13. Booth A. Carbon footprint modelling of national health systems: opportunities, challenges and recommendations. Int J Health Plann Manage. (2022) 37:1885–93. doi: 10.1002/hpm.3447

PubMed Abstract | Crossref Full Text | Google Scholar

14. Rizan C, Lillywhite R, Reed M, Bhutta MF. The carbon footprint of products used in five common surgical operations: identifying contributing products and processes. J R Soc Med. (2023) 116:199–213. doi: 10.1177/01410768231166135

PubMed Abstract | Crossref Full Text | Google Scholar

15. Rizan C, Reed M, Mortimer F, Jones A, Stancliffe R, Bhutta MF. Using surgical sustainability principles to improve planetary health and optimise surgical services following the COVID-19 pandemic. Bull R Coll Surg Engl. (2020) 102:177–81. doi: 10.1308/rcsbull.2020.148

Crossref Full Text | Google Scholar

16. Holmner Å, Rocklöv J, Ng N, Nilsson M. Climate change and eHealth: a promising strategy for health sector mitigation and adaptation. Glob Health Action. (2012) 5:18428. doi: 10.3402/gha.v5i0.18428

PubMed Abstract | Crossref Full Text | Google Scholar

17. Hou Y, Sun M, Huang X, Nan J, Gao J, Zhu N, et al. Autonomy support in telehealth: an evolutionary concept analysis. Front Public Health. (2025) 13:1544840. doi: 10.3389/fpubh.2025.1544840

PubMed Abstract | Crossref Full Text | Google Scholar

18. Bouley T, Roschnik S, Karliner J, Wilburn S, Slotterback S, Guenther R, et al. Climate-Smart Healthcare: Low-Carbon and Resilience Strategies for the Health Sector, Vol. 1. World Bank Group, United States of America. Available online at: https://coilink.org/20.500.12592/dnvt1s (Accessed August 24, 2025).

Google Scholar

19. Grimmond TR, Bright A, Cadman J, Dixon J, Ludditt S, Robinson C, et al. Before/after intervention study to determine impact on life-cycle carbon footprint of converting from single-use to reusable sharps containers in 40 UK NHS trusts. BMJ Open. (2021) 11:e046200. doi: 10.1136/bmjopen-2020-046200

PubMed Abstract | Crossref Full Text | Google Scholar

20. Fathy R, Nelson CA, Barbieri JS. Combating climate change in the clinic: cost-effective strategies to decrease the carbon footprint of outpatient dermatologic practice. Int J Womens Dermatol. (2021) 7:107–11. doi: 10.1016/j.ijwd.2020.05.015

PubMed Abstract | Crossref Full Text | Google Scholar

21. Mominkhan D, Alamri FA, Balharith M, Alharbi M, Alshebli A, Alshareef SY, et al. The current state and potential evolution of carbon emissions in the healthcare sector: a narrative review article. Front Sustain Energy Policy. (2023) 2:1230253. doi: 10.3389/fsuep.2023.1230253

Crossref Full Text | Google Scholar

22. Duindam D. Transitioning to sustainable healthcare: decarbonising healthcare clinics, a literature review. Challenges. (2022) 13:68. doi: 10.3390/challe13020068

Crossref Full Text | Google Scholar

23. Zhang G, Chen SCI, Yue X. Blockchain technology in carbon trading markets: impacts, benefits, and challenges—a case study of the shanghai environment and energy exchange. Energies. (2024) 17:3296. doi: 10.3390/en17133296

Crossref Full Text | Google Scholar

24. Zhao S, Ma G, Ding J. Symbiotic mechanism of multiple subjects for the resource-based disposal of medical waste in China in the post-pandemic context. Sustainability. (2023) 15:805. doi: 10.3390/su15010805

Crossref Full Text | Google Scholar

25. Gao L, Wang X. Healthcare supply chain network coordination through medical insurance strategies with reference price effect. Int J Environ Res Public Health. (2019) 16:3479. doi: 10.3390/ijerph16183479

PubMed Abstract | Crossref Full Text | Google Scholar

26. Yu, Z. Optimizing internal control in public hospital supply chain: a game theory-based approach. Front Public Health. (2023) 11:1240757. doi: 10.3389/fpubh.2023.1240757

PubMed Abstract | Crossref Full Text | Google Scholar

27. Gong H, Wang X, Zhang T, Li J, Chen B. How can China's new health care reform promote the balance of interest game?–based on game evolution and simulation analysis. Risk Manag Healthc Policy. (2023) 16:1435–54. doi: 10.2147/RMHP.S422296

PubMed Abstract | Crossref Full Text | Google Scholar

28. Nowak MA, Sigmund K. Evolutionary dynamics of biological games. Science. (2004) 303:793–9. doi: 10.1126/science.1093411

PubMed Abstract | Crossref Full Text | Google Scholar

29. Iqbal A. Studies in the theory of quantum games. arXiv preprint quant-ph/0503176. (2005). doi: 10.48550/arXiv.quant-ph/0503176

Crossref Full Text | Google Scholar

30. Gokhale CS, Traulsen A. Evolutionary multiplayer games. Dyn Games Appl. (2014) 4:468–88. doi: 10.1007/s13235-014-0106-2

Crossref Full Text | Google Scholar

31. Yang X, He G, Zhu Z, Zhao S, Zhang S. Evolutionary game analysis and efficiency test of water pollution control driven by emission trading: evidence from Zhejiang Province, China. Heliyon. (2024) 10:e36289. doi: 10.1016/j.heliyon.2024.e36289

PubMed Abstract | Crossref Full Text | Google Scholar

32. Lian X, Wang W, Zhang J. How to optimize dust pollution control in opencast coal mines: analysis of a joint social regulation model based on evolutionary game theory. PLoS ONE. (2023) 18:e0289164. doi: 10.1371/journal.pone.0289164

PubMed Abstract | Crossref Full Text | Google Scholar

33. Lee VS, Gerwig K, Hough E, Mate K, Biggio R, Kaplan RS. Decarbonizing health care: engaging leaders in change. NEJM Catal Innov Care Deliv. (2023) 4:CAT-22. doi: 10.1056/CAT.22.0433

Crossref Full Text | Google Scholar

34. Kmietowicz Z. Healthcare and tech companies pledge to decarbonise NHS supply chain by 2045. BMJ. (2021) 375:n2696. doi: 10.1136/bmj.n2696

PubMed Abstract | Crossref Full Text | Google Scholar

35. Guo W, Chen S, Lei M. Evolutionary game and strategy analysis of carbon emission reduction in supply chain based on system dynamic model. Sustainability. (2023) 15:8933. doi: 10.3390/su15118933

Crossref Full Text | Google Scholar

36. Leal Filho W, Luetz JM, Thanekar UD, Dinis MAP, Forrester M. Climate-friendly healthcare: reducing the impacts of the healthcare sector on the world's climate. Sustain Sci. (2024) 19:1103–9. doi: 10.1007/s11625-024-01487-5

Crossref Full Text | Google Scholar

37. Malhi JK, McEvoy JW, Blumenthal RS, Jacobsen AP. Climate change and cardiovascular health: recent updates and actions for healthcare. Am Heart J Plus: Cardiol Res Pract. (2024) 45:100443. doi: 10.1016/j.ahjo.2024.100443

PubMed Abstract | Crossref Full Text | Google Scholar

38. Zhang L, Wang X, Xiao H, Ma C, Li X, Dai G, et al. Governance mechanisms for chronic disease diagnosis and treatment systems in the post-pandemic era. Front Public Health. (2022) 10:1023022. doi: 10.3389/fpubh.2022.1023022

PubMed Abstract | Crossref Full Text | Google Scholar

39. Braithwaite J, Pichumani A, Crowley P. Tackling climate change: the pivotal role of clinicians. BMJ. (2023) 382:76963. doi: 10.1136/bmj-2023-076963

PubMed Abstract | Crossref Full Text | Google Scholar

40. Boubaker O, Balas VE, Benzaouia A, Chaabane M, Mahmoud MS, Zhu Q. Time-delay systems: modeling, analysis, estimation, control, and synchronization. Math Probl Eng. (2017) 2017:1398904. doi: 10.1155/2017/1398904

PubMed Abstract | Crossref Full Text | Google Scholar

41. Zhang L, Chen M, Wu Q, Wu B. Stabilization of networked control systems with induced delays and actuator saturation. Math Probl Eng. (2016) 2016:4208712. doi: 10.1155/2016/4208712

PubMed Abstract | Crossref Full Text | Google Scholar

42. Wang L, Lu W, Chen T. Multistability and new attraction basins of almost-periodic solutions of delayed neural networks. IEEE Trans Neural Netw. (2009) 20:1581–93. doi: 10.1109/TNN.2009.2027121

PubMed Abstract | Crossref Full Text | Google Scholar

43. Song ZG, Xu J. Stability switches and double Hopf bifurcation in a two-neural network system with multiple delays. Cogn Neurodyn. (2013) 7:505–21. doi: 10.1007/s11571-013-9254-0

PubMed Abstract | Crossref Full Text | Google Scholar

44. Guo H, Qian Y, Wang B, Qiu Y. The application of time-delay-dependent H∞ control in the transmission effect of monetary policy on real estate market. Math Probl Eng. (2014) 2014:421510. doi: 10.1155/2014/421510

Crossref Full Text | Google Scholar

45. Jørgensen SV, Hauschild MZ, Nielsen PH. The potential contribution to climate change mitigation from temporary carbon storage in biomaterials. Int J Life Cycle Assess. (2015) 20:451–62. doi: 10.1007/s11367-015-0845-3

Crossref Full Text | Google Scholar

46. Usmani O, Winders T. Prioritising patient outcomes and reducing environmental burden: how both are achievable in respiratory care. EMJ Respir. (2023) 11:29–40. doi: 10.33590/emjrespir/10305679

Crossref Full Text | Google Scholar

47. Lakatos K, Teherani A, Thottathil SE, Gandhi S, Weiser SD, Brindis CD. A race to net zero—early lessons from healthcare's decarbonization marathon. Health Aff Sch. (2023) 1:qxad006. doi: 10.1093/haschl/qxad006

PubMed Abstract | Crossref Full Text | Google Scholar

48. Portnoy A, Nygård M, Trogstad L, Kim JJ, Burger EA. Impact of delaying effective and cost-effective policy decisions: an example from cervical cancer prevention in Norway. MDM Policy Practice. (2022) 7:23814683211071093. doi: 10.1177/23814683211071093

PubMed Abstract | Crossref Full Text | Google Scholar

49. Tennison I, Roschnik S, Ashby B, Boyd R, Hamilton I, Oreszczyn T, et al. Health care's response to climate change: a carbon footprint assessment of the NHS in England. Lancet Planetary Health. (2021) 5:e84–92. doi: 10.1016/S2542-5196(20)30271-0

PubMed Abstract | Crossref Full Text | Google Scholar

50. Taylor PD, Jonker LB. Evolutionary stable strategies and game dynamics. Math Biosci. (1978) 40:145–56. doi: 10.1016/0025-5564(78)90077-9

Crossref Full Text | Google Scholar

51. Lyapunov AM. The general problem of the stability of motion. Int J Control. (1992) 55:531–4. doi: 10.1080/00207179208934253

Crossref Full Text | Google Scholar

52. Friedman D. A simple testable model of double auction markets. J Econ Behav Organ. (1991) 15:47–70. doi: 10.1016/0167-2681(91)90004-H

Crossref Full Text | Google Scholar

53. Ritzberger K, Weibull JW. Evolutionary selection in normal-form games. Econometrica. (1995) 63, 1371–99. doi: 10.2307/2171774

PubMed Abstract | Crossref Full Text | Google Scholar

54. Gray, W. B., Shadbegian, R. Environmental Regulations and Business Decisions. Bonn: IZA World of Labor (2025). p. 187. doi: 10.15185/izawol.187.v2

Crossref Full Text | Google Scholar

55. Fleckinger P, Glachant M. Negotiating a voluntary agreement when firms self-regulate. J Environ Econ Manage. (2011) 62:41–52. doi: 10.1016/j.jeem.2011.03.002

Crossref Full Text | Google Scholar

Keywords: medical carbon neutrality, medical supply chain, public health policy, four-party evolutionary game, time delay

Citation: Xie B, Gao X, Cheng Y, Niu M and Liu S (2025) Game theoretic research on strategies for medical carbon neutrality stakeholders under time delay effects. Front. Public Health 13:1640294. doi: 10.3389/fpubh.2025.1640294

Received: 03 June 2025; Accepted: 12 August 2025;
Published: 09 September 2025.

Edited by:

Beatriz S. Lima, Research Institute for Medicines (iMed.ULisboa), Portugal

Reviewed by:

Ricardo De Moraes E. Soares, Instituto Politecnico de Setubal (IPS), Portugal
Khadija Sajid, Government College University, Lahore, Pakistan
Ahmad Zaini Miftah, Universitas Padjadjaran, Indonesia

Copyright © 2025 Xie, Gao, Cheng, Niu and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yingying Cheng, OTkwNjE4MEBoYXVzdC5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.