Distribution dynamics,regional differences,and convergence of sustainable development of cities and communities in China

China is one of the most populated and rapidly urbanizing countries worldwide and was among the earliest countries to integrate sustainable development into urban construction.To achieve high-quality development and implement the objectives of“Transforming Our World:The 2030 Agenda for Sustainable Development”,it is crucial to measure and analyze the current level of sustainable development of cities in China.Following the principles of relevance,scientific rigor,universality,reliability,and timeliness,this study constructs an assessment index system for sustainable development,covering seven themes corresponding to the UN Sustainable Development Goal 11.Through detailed calculations,we obtained sustainable development indices for 139 Chinese cities from 2016 to 2022 and analyzed them in three dimensions:distribution dynamics,regional differences,and convergence.The key findings are as follows.First,the level of sustainable development showed improvement,with the average score of included cities increasing by 11.88%from 2016 to 2022.Second,the level of sustainable development was relatively balanced,maintaining low Gini coefficients between 0.05 and 0.06.Third,a weak overallσconvergence feature existed,with increased differentiation in 2021.From a regional perspective,aσconvergence feature was observed in the northeastern but not in the western region.Fourth,both overall absoluteβ-convergence and conditionalβeffects were significant.Regional absoluteβ-convergence and conditionalβ-convergence were also significant.This study contributes to the literature by providing evidence of China’s urban sustainable development,offering policy insights for deepening the implementation of development goals in the future,and providing experiential reference for other developing countries to achieve sustainable development.

Hydrochemical Dynamic Characteristics and Evolution of Underground Brine in the Mahai Salt Lake of the Qaidam Basin Qinghai-Tibet Plateau

The mineral rock salts present in the Mahai Salt Lake of the Qaidam basin exhibit high solubilities in water. In addition, the multicomponent underground brine exhibits a high salinity and is easily precipitated. In the natural state, brine transport in the brine layer is extremely slow, and the brine is in a relatively stable chemical equilibrium state with the rock salt media. However, during mining, both the seepage and the chemical fields fluctuate significantly, thereby disrupting the equilibrium and leading to variations in the chemical composition and dynamic characteristics of the brine. Therefore, we selected underground brine from the Mahai Salt Lake, collecting a total of 183 brine samples over three stages of mining（i.e., the early stage of underground brine extraction, the initial stage of mining, and the later stage of mining）. Using a range of analytical techniques, the chemical dynamics of the underground brine water and its evolution were systematically studied. We found that evaporation and enrichment were the main mechanisms of underground brine evolution in the Mahai Salt Lake, with cation exchange and mineral dissolution/precipitation being key factors in determining the dynamic characteristics and evolution of the brine.

Dynamic evolution of shear rate-dependent behavior of rock discontinuity under shearing condition

Rock blocks sliding along discontinuities can cause serious disasters,such as landslides,earthquakes,or rock bursts.The shear rate-dependent behavior is a typical time-dependent behavior of a rock discontinuity,and it is closely related to the stability of a rock block.To further study the shear rate-dependent behavior of rock discontinuities,shear tests with alternating shear rates(SASRs)were conducted on rock discontinuities with various surface morphologies.The dynamic evolution of the shear rate dependency was studied in detail based on the shear test results,and three stages were identified with respect to the shear stress and shear deformation states.The test results revealed that dynamic changes in shear stiffness and the energy storage abilities of the rock discontinuities occurred in relation to the shear rate-dependent behavior of crack growth,which increased with an increase in normal stress and/or the joint roughness coefficient.The stage of decreasing shear stiffness corresponded to a stage of noticeable shear rate-dependency,and the shear rate was found to have no influence on the initial crack stress.

Dynamic evolution patterns of the degree of ground resistivity anisotropy and the sets mogenic process

Following a new train of thinking, this paper has explored first the potential information in the ground resistivitydata observed by the existing geoelectric observation system, investigated and proposed a new dimensionlessgeoelectric precursor factor, the degree of ground resistivity anisotropy, S, and studied the characteristics of dynamic evolution pattern of S during the seismogenic process. The results show that, during the seismogenic process, the degree of ground resistivity anisotropy (S) displays h process of 'normal' → 'abnormal strengthening(amplitude, range)' → 'abnormal weakening'→ 'earthquake occurrence'→ 'normal'. The earthquake wouldoccur at the time when the S value has entered the late stage of strengthening and turns to weaken and in the gradient belt on the margin ofS anomaly region. The dynamic evolution pattern ofS reflects the changes of the tectonicstress field during the seismogenic process. Therefore, it would be possible to trace the process of earthquake generation and occurrence from the dynamic evolution pattern ofS so as to service eaJ'thquake prediction.

A semi-analytic model for the study of 1/1 resonant dynamics of the planar elliptic restricted co-orbital problem

Mean motion resonances(MMRs)are widespread in our Solar System.Moreover,resonant dynamics has always been an essential topic in planetary research.Recently,the research about exoplanets and the potential Planet Nine with large eccentricity has given rise to our interests in the secular dynamics inside MMRs in the elliptic model.In this paper,we study the fixed points of the averaged Hamiltonian and the long-term stable regions of the 1/1 resonance(or co-orbital motion)in the elliptic restricted three-body problem(ERTBP)systematically.Numerical integrations prove those test particles trapped in"apsidal co-rotation",where both the resonant angleφ_(res)and the secular angleΔω(or apsidal longitude differences)librate simultaneously,always survive the long-term simulations.Furthermore,utilizing a semianalytical method based on the adiabatic approach,three families of long-term fixed points of the averaged Hamiltonian of the planar ERTBP inside the 1/1 resonance have been found.We call them QS-points,Hpoints,and T-points here,whose values of the(φ_(res),Δω)are(0°,180°),(180°,0°),and(±60°,±60°),respectively.All the fixed points of the averaged Hamiltonian of the co-orbital motion in the ERTBP are presented in the e-e’plane(’represents the elements of the planet in this paper).We find that QS-points and T-points always exist for the arbitrary eccentricity of a planet,while H-points only exist for the cases of low e’and very high e.Furthermore,we measure the libration width in terms of eccentricity,Δe,around these stable equilibrium points in the e-Δωphase-space portraits.The"apsidal co-rotation"around all the stable equilibrium points is presented in the e-e’plane.All these results are effectively confirmed by numerical experiments.The long-term stable zones around these periodic orbits in the e-e’plane are significant for the research of the co-orbital motion in the ERTBP.Above all,these practical approaches that we proposed can also be used to study the secular dynamics of other MMRs.

Solving chemical dynamic optimization problems with ranking-based differential evolution algorithms

Dynamic optimization problems(DOPs) described by differential equations are often encountered in chemical engineering. Deterministic techniques based on mathematic programming become invalid when the models are non-differentiable or explicit mathematical descriptions do not exist. Recently, evolutionary algorithms are gaining popularity for DOPs as they can be used as robust alternatives when the deterministic techniques are invalid. In this article, a technology named ranking-based mutation operator(RMO) is presented to enhance the previous differential evolution(DE) algorithms to solve DOPs using control vector parameterization. In the RMO, better individuals have higher probabilities to produce offspring, which is helpful for the performance enhancement of DE algorithms. Three DE-RMO algorithms are designed by incorporating the RMO. The three DE-RMO algorithms and their three original DE algorithms are applied to solve four constrained DOPs from the literature. Our simulation results indicate that DE-RMO algorithms exhibit better performance than previous non-ranking DE algorithms and other four evolutionary algorithms.

Analysis of dynamical properties for the two-site Bose Hubbard model with an algebraic method

In this work, we propose an algebraic recursion method to study the dynamical evolution of the two-site Bose- Hubbard model. We analyze its properties from the viewpoints of single partite purity, energy, and trace distance, in which the model is considered as a typical bipartite system. The analytical expressions for the quantities are derived. We show that the purity can well reflect the transition between different regimes for the system. In addition, we demonstrate that the transition from the delocalization regime to the self-trapping regime with the ratio r/increasing not only happens for an initially local state but also for any initial states. Furthermore, we confirm that the dynamics of the system presents a periodicity for η = 0 and the period is tc =π/2J when the initial state is symmetric.

Time Evolution,Dynamical Quantum Fluctuation and High-Order Squeezing Feature in Polariton System——Ⅰ

We have set up a new reduced model Hamiltonian for the polariton system, in which the nonlinear interaction contains the rotating term k l （ a ＋ b ＋ ab＋） and the attractive two-mode squeezed coupling - k2 （ a ＋ b＋ ＋ ab ） . The dynamical evolution of this system has been solved and the nonclassical features relevant to the second-order and high-order squeezing have been obtained in an analytical form. For the first time, in contrast to the existing result, we have confirmed for the phonon field that the attractive two-mode squeezed interaction will not only result in the second-order and high-order squeezing in X-component with the time evolution, but also in time average. Furthermore, the phenomena of collapse and revival of inversion will occur as well in the time evolution of the average number of photon and phonon, as also in the second-order and high-order squeezing of photon field, particularly, in the high-order squeezing of phonon field.

Regional difference and dynamic evolution of development quality of power industry in China

To achieve the goals of carbon peaking and carbon neutrality and maintain high-quality economic growth,China is currently striving to improve the quality of development of its power sector.In this regard,revealing the regional differences and evolutionary trends in the development quality of China’power sector has a high value to inspire the next improvement direction toward how to integrate regional power recourses to an overall optimization level.Motived by this purpose,this paper uses the entropy method to evaluate the com‐prehensive and subsystem indices of the development quality of the power industry,and reveals their re‐gional differences and evolutionary trends with the help of the Dagum Gini coefficient and Kernel density es‐timation methods.The findings show that:There are obvious regional differences in the development quality of China’s power industry,and the differences are steadily declining in all regions except the West.Regional differences are mainly derived from inter-regional differences,with the largest inter-regional differences in the East-Northeast region.Intra-regional differences show a distribution pattern of East>West>North‐east>Center.

A Dynamic Evolution Model of Airline Networks

Empirical data show that most of the degree distribution of airline networks assume a double power law. In this work, firstly, we assume cities as sites, flight between two cities as an edge between two sites, and build a dynamic evolution model for airline networks by improving the BA model, in which the conception of attractiveness plays a decisive role in the course of evolution of the networks. To this end, we discuss whether the attractiveness depends on the site label s or not separately, finally we obtain analytic degree distribution. As a result, if the attractiveness of a site is independent of the degree distribution of sites, which will follow the double power law, otherwise, it will be scale-free. Moreover, degree distribution depends on the parameters of the models, and some parameters aye more sensitive than others.

Features and physical proces of the dynamic evolution pattern of ground resistivity precursor front

In order to quantitatively describe the geoelectric precursor anomaly in the short-impending process of earthquakes, a new geoelectric precursor index — (monthly) relative change rate of ground resistivity, R ρ (t) , is designed. Using this index and choosing the internationally accepted ground resistivity data before the Tangshan M =7.8 earthquake of July 28, 1976, the features of dynamic evolution pattern of R ρ(t) are studied. The results show that: ① about 10～9 months before earthquake, the ground resistivity in a certain range around the epicentral region begins to display the anomaly of accelerating descent, and the rate of descent is higher in the epicentral region than in surrounding areas; ② with the shortening of countdown before earthquake, the R ρ(t) value in epicentral region increases gradually (ground resistivity value decreases at an increasing rate); ③ the R ρ(t) value has the epicentral area as a center and its contour lines propagate towards surrounding areas with the shortening of countdown before earthquake; ④after the R ρ(t) value in epicentral region has descended at increasing rate to reach an extremity [ R ρ(t) = (7.0], it turns to descend at decreasing rate (2～3 months) and earthquake occurs when it accelerates again. Meanwhile, earthquake occurs when the contour lines of R ρ(t) stop propagating towards surrounding areas and turn to shrink back (2～3 months later). Its physical process can be explained by the″ swollen hypothesis″ of Prof. Fu and the theory of ″Slip-weakening and rockmass instability″ of Mei, Niu, et al ..

Dynamic aggregation evolution of competitive societies of cooperative and noncooperative agents

We propose an evolution model of cooperative agent and noncooperative agent aggregates to investigate the dynamic evolution behaviors of the system and the effects of the competing microscopic reactions on the dynamic evolution. In this model, each cooperative agent and noncooperative agent are endowed with integer values of cooperative spirits and nonco- operative spirits, respectively. The cooperative spirits of a cooperative agent aggregate and the noncooperative spirits of a noncooperative agent aggregate change via four competing microscopic reaction schemes： the win-win reaction between two cooperative agents, the lose-lose reaction between two noncooperative agents, the win-lose reaction between a coop- erative agent and a noncooperative agent （equivalent to the migration of spirits from cooperative agents to noncooperative agents）, and the cooperative agent catalyzed decline of noncooperative spirits. Based on the generalized Smoluchowski＇s rate equation approach, we investigate the dynamic evolution behaviors such as the total cooperative spirits of all coop- erative agents and the total noncooperative spirits of all noncooperative agents. The effects of the three main groups of competition on the dynamic evolution are revealed. These include： （i） the competition between the lose-lose reaction and the win-lose reaction, which gives rise to respectively the decrease and increase in the noncooperative agent spirits; （ii） the competition between the win-win reaction and the win-lose reaction, which gives rise to respectively the increase and decrease in the cooperative agent spirits; （iii） the competition between the win-lose reaction and the catalyzed-decline reaction, which gives rise to respectively the increase and decrease in the noncooperative agent spirits.

Dynamical Evolution of an Effective Two-Level System with PT Symmetry

We investigate the dynamics of parity-and time-reversal（PT） symmetric two-energy-level atoms in the presence of two optical and one radio-frequency fields. The strength and relative phase of fields can drive the system from the unbroken to the broken PT symmetric regions. Compared with the Hermitian model, Rabi-type oscillation is still observed, and the oscillation characteristics are also adjusted by the strength and relative phase in the region of the unbroken symmetry. At the exception point, the oscillation breaks down. To better understand the underlying properties we study the effective Bloch dynamics and find that the emergence of the PT components of the fixed points is the feature of the PT symmetry breaking and the projections in the x–y plane can be controlled with high flexibility compared with the standard two-level system with the PT symmetry. It helps to study the dynamic behavior of the complex PT symmetric model.

The Investigation of the Dynamical Evolution of Extrasolar Three-planetary System GJ 3138

This article is devoted to studying the dynamical evolution and orbital stability of compact extrasolar threeplanetary system GJ 3138. In this system, all semimajor axes are less than 0.7 au. The modeling of planetary motion is performed using the averaged semi-analytical motion theory of the second order in planetary masses,which the authors construct. Unknown and known with errors orbital elements vary in allowable limits to obtain a set of initial conditions. Each of these initial conditions is applied for the modeling of planetary motion. The assumption about the stability of observed planetary systems allows to eliminate the initial conditions leading to excessive growth of the orbital eccentricities and inclinations and to identify those under which these orbital elements conserve moderate values over the whole modeling interval. Thus, it becomes possible to limit the range of possible values of unknown orbital elements and determine their most probable values in terms of stability.

Collaboration strategy for software dynamic evolution of multi-agent system

As the ability of a single agent is limited while information and resources in multi-agent systems are distributed, cooperation is necessary for agents to accomplish a complex task. In the open and changeable environment on the Internet, it is of great significance to research a system flexible and capable in dynamic evolution that can find a collaboration method for agents which can be used in dynamic evolution process. With such a method, agents accomplish tasks for an overall target and at the same time, the collaborative relationship of agents can be adjusted with the change of environment. A method of task decomposition and collaboration of agents by improved contract net protocol is introduced. Finally, analysis on the result of the experiments is performed to verify the improved contract net protocol can greatly increase the efficiency of communication and collaboration in multi-agent system.

Capture Efficiency Analysis in the Circular Restricted Three-body Problem

Temporary capture efficiency is studied in the framework of the circular restricted three-body problem in two steps.First, a non-uniform distribution of test particles around the secondary's orbit is obtained by fully accounting the secondary's gravitational influence. Second, the capture efficiency is computed based on the non-uniform distribution. Several factors influencing the result are discussed. By studying the capture efficiency in the circular restricted three-body problem of different mass ratios, a power-law relation between the capture efficiency(p) and the mass ratio(μ) is established, which is given by p ≈ 0.27 × μ^(0.53), within the range of 3.0035 ×10^(-6)≤ μ ≤ 3.0034 × 10^(-5). Taking the Sun–Earth system as an example, the influence from the orbit eccentricity of the secondary on the non-uniform distribution and the capture efficiency is studied. Our studies find that the secondary's orbit eccentricity has a negative influence on the capture efficiency.

A Novel Ephemeris Model for Martian Moons Incorporating Their Free Rotation

High-precision ephemerides not only support space missions,but can also be used to study the origin and future of celestial bodies.In this paper,a coupled orbit-rotation dynamics model that fully takes into account the rotation of the Martian moons is developed.Phobos and Deimos’rotations are first described by Eulerian rotational equations,and integrated simultaneously with the orbital motion equations.Orbital and orientational parameters of Mars satellites were simultaneously obtained by numerical integration for the first time.In order to compare the differences between our newly developed model and the one now used in the ephemerides,we first reproduced and simulated the current model using our own parameters,and then fit it to the Institut de Mécanique Céleste et de Calcul deséphémérides ephemerides using least-square procedures.The adjustment test simulations show Phobos and Deimos’orbital differences between the refined model and the current model are no more than 300 m and125 m,respectively.The orientation parameters are confirmed and the results are in good agreement with the International Astronomical Union results.Moreover,we simulated two perturbations(main asteroids and mutual torques)which were not included in our refined model,and find that their effects on the orbits are completely negligible.As for the effect on rotation,we propose to take care of the role of mutual attraction in future models.

Transformation power or development pressure:economic growth targets and urban carbon productivity

Maintaining moderate economic growth targets(EGTs)is the key for local governments to effectively implement the“carbon peak and carbon neutrality”goals under the refreshed development pattern.Utilizing panel data of 276 prefecture-level cities in China's Mainland from 2010 to 2020,and employing methods such as intermediary and threshold models,this study empirically analyzes the internal mechanism of EGT’s impact on urban carbon productivity(UCP).Our findings demonstrate that:①The overall EGT during the analyzed period is not conducive to improving UCP.This conclusion remains valid after a series of robustness tests.This effect is more pronounced in the central region and resource-based cities than in the east-west region and non resource-based cities.②EGT not only directly suppresses UCP but also exerts indirect negative impacts on UCP from three aspects:delaying the digital economy(DE),constraining financial expansion(FE),and hindering green technology innovation(GTI).This negative indirect effect is similar to or even surpasses the direct effect,suggesting that the internal relationship between EGT and“dual-carbon”goals should be re-evaluated from a new compound perspective.③EGT not only has a simple linear impact on UCP but also significantly exhibits a dynamic evolution pattern in inverted“U”shape.That is,as EGT continuously upgrades,a nonlinear impact on UCP emerges in the form of“promoting first,suppressing later”.This indicates that surpassing the“degree”limit for EGT will be detrimental to the improvement of UCP.This study broadens the scope of carbon productivity analysis by introducing a new perspective centered on EGT.The insights gleaned from this research offer valuable guidance for local governments to effectively manage economic growth expectations and promote the synchronized achievement of dual-carbon objectives.

Scale effect of coastal landscape pattern stability and driving forces:a case study of Guangdong Province,China

The long-term dynamic evolution and underlying mechanisms of coastal landscape pattern stability,driven by strong anthropogenic interference and consequently climate change,are topics of major interest in national and international scientific research.Guangdong Province,located in southeastern China,has been undergoing rapid urbanization over several decades.In this study,we quantitatively determined the scale threshold characteristics of coastal landscape pattern stability in Guangdong Province,from the dual perspective of spatial heterogeneity and spatial autocorrelation.An analysis of the spatiotemporal evolution of the coastal landscape was conducted after the optical scale was determined.Then,we applied the geodetector statistical method to quantitatively explore the mechanisms underlying coastal landscape pattern stability.Based on the inflection point of landscape metrics and the maximum value of the MoranⅠindex,the optimal scale for analyzing coastal landscape pattern stability in Guangdong Province was 240 m×240 m.Within the past several decades,coastal landscape pattern stability increased slightly and then decreased,with a turning point around 2005.The most significant variations in coastal landscape pattern stability were observed in the transition zone of rural-urban expansion.A q-statistics analysis showed that the explanatory power of paired factors was greater than that of a single driving factor;the paired factors with the greatest impact on coastal landscape pattern stability in Guangdong Province were the change in gross industrial output and change in average annual precipitation from 2010 to 2015,based on a q value of 0.604.These results will contribute to future efforts to achieve sustainable coastal development and provide a scientific basis and technical support for the rational planning and utilization of resources in large estuarine areas,including marine disaster prevention and seawall ecological restoration.

Triple the steady-state reaction rate by decorating the In_(2)O_(3)surface with SiO_(x)for CO_(2)hydrogenation

Indium oxide(In_(2)O_(3)),as a promising candidate for CO_(2)hydrogenation to C_(1) products,often suffers from sintering and activity decline,closely related to the undesirable structural evolution under reaction conditions.Based on the comprehension of the dynamic evolution,this study presents an efficient strategy to alleviate the agglomeration of In_(2)O_(3)nanoparticles by the surface decoration with highly dispersed silica species(SiO_(x)).Various structural characterizations combined with density functional theory calculations demonstrated that the sintering resulted from the over-reduction,while the enhanced stability originated from the anchoring effect of highly stable In-OSi bonds,which hinders the substantial formation of metallic In(In^(0))and the subsequent agglomeration.0.6Si/In_(2)O_(3)exhibited CO_(2)conversion rate of10.0 mmol g^(-1)h^(-1)at steady state vs.3.5 mmol g^(-1)h^(-1)on In_(2)O_(3)in CO_(2)hydrogenation.Enhanced steady-state activity was also achieved on Pd-modified catalysts.Compared to the traditional Pd/In_(2)O_(3)catalyst,the methanol production rate of Pd catalyst supported on 0.6Si/In_(2)O_(3)was enhanced by 23%,showing the potential of In_(2)O_(3)modified by SiO_(x)in serving as a platform material.This work provides a promising method to design new In_(2)O_(3)-based catalysts with improved activity and stability in CO_(2)hydrogenation.