Theory and Practice of Quantitative Assessment of System Harmonicity: Case of Road Safety in Russia before and during the COVID-19 Epidemic

: People COVID lockdown played a positive role in this process. Harmonicity of the Russian road safety provision system, assessed by the normalized general utility GU n RSS , increased by 0.46% from 2019 to 2020.


Introduction
The year 2020 will go down in history as a very unusual year.The Great Reset is a relevant slogan, the authors of which in the context of the events of 2020 and subsequent years are K. Schwab and T. Malleret [1].Widespread lockdowns forced the transfer of activity from off-line to on-line.The forced decline in activity has changed people's lifestyles everywhere.Naturally, this caused a decrease in the transport mobility of people.
To a greater extent, this concerned citizens in the largest cities of the world whose transport traffic, according to the Tom Tom Traffic Index [2], decreased by 20-80% during the active lockdown phase (April-May 2020) [3].Rural residents were also forced to reduce transport activity, although to a lesser extent than urban residents [4].This could not but affect various aspects of the functioning of the transport sector of the economy.Primarily, this concerns road safety.
Data analysis of the Road Safety Annual Report 2020 [5] indicates two important facts.In countries with a hard lockdown in the spring of 2020, the number of deaths in road accidents significantly decreased; when comparing "April 2020 to April 2019", road traffic mortality reduced by 30-80%.Conversely, in those few countries where lockdown was not introduced, the death rate in road accidents increased; in Sweden and the Netherlands, the countries that are long-standing world leaders in road safety, the number of people killed in road accidents in April 2020 increased by 6% compared to April 2019.Of course, all this could not go unnoticed by the regional authorities, or by the global level of government.
One of the most important decisions of the 74th session (18.08.2020) of the United Nations (UN) General Assembly was resolution A/RES/74/299, "Improving global road safety" [6].This document highlights the importance of the work already conducted in this area over the past decade.It also formulates the need to achieve a new goal by 2030, which is to reduce the number of deaths on roads by 50%.It should be assumed that this ambitious goal was formulated at the time when the first results of the assessment of the COVID lockdown's impact on road traffic accidents had already been received, and the connection between the restriction in transport mobility and the number of fatalities in road accidents was comprehended.Perhaps these statistics largely served as the basis for cautious optimism in promoting the concept of Mobility-as-a-Service (MaaS), which is the ideological basis for abandoning personal transport in favor of a transportation service provider.Expectation of the upcoming mass rejection of the individual car is perhaps the main reason for optimism about potential trends in road safety.In some countries, targets have been announced for a multiple reduction in the number of people killed in road accidents over the next decade [7,8].An example of such goal-setting is the Russian Federation [9], where the task is to reach a level of human risk of four deaths in road accidents per hundred thousand people by 2024, and a level of zero mortality in road accidents by 2030 [10].For reference, in 2020, the level of human risk in Russia was equal to 10.80 people killed in road accidents per hundred thousand people [11].Meanwhile, no assessment of the real possibility of achieving this goal was carried out.
How possible is the success of the practical implementation of this target setting?How did the COVID restrictions affect the effectiveness, usefulness and harmony of state road safety systems?These are the most important issues and are discussed by the authors from the standpoint of the theory of system harmony, when comparing the situation in the field of assessing the organization of the state system for ensuring road safety in the Russian Federation during characteristic periods (between 2006 and 2020, and comparing the situation of 2019, a conditionally normal year, to 2020, the lockdown year).

Related Works 2.1. Synergetics, Orderliness and System Harmonicity
The subject of synergetics can be defined as the study of the process of system selforganization.Social synergetics studies the patterns of self-organization of society, i.e., the relationship between social order and social chaos (their mutual transition into one another and their synthesis).R. Benedict [12] is the author of the "synergy" concept and initially the authors specializing in the interaction between the individual and society were engaged in the problems of synergy [13].Today, synergetics is mainly a tool of cognition in the social sciences [14][15][16][17].
However, it is generally recognized that the main ideas of synergetics were laid down in the works of scientists of the natural science profile: I. Prigogine, I. Stengers, and G. Nicolis [18][19][20][21].
Nevertheless, H. Haken can be considered the founder of synergetics [22][23][24].He defined synergetics as the science of self-organization, the theory of "the joint action of many subsystems, that in result arises a structure and corresponding functioning at the macroscopic level" [23].This is mentioned in analytical works [25,26].Today, synergetics is recognized as a scientific direction all over the world.In Russia, the ideas of synergetics were actively promoted and developed by scientists of the scientific group of S. P. Kurdyumov [27,28].
Synergetics made a revolution by creating a new image of the world that continuously evolves according to nonlinear laws.
The basic concepts of synergetics are "order" and "chaos".The concept of "order" is usually used when describing a system of stable, repetitive in space and time relationships between elements of any nature.Conversely, the concept of "chaos" is usually used to describe a set of elements between which there are no stable repeating relationships.Dissipative systems, which include any socio-technical systems, can exist only in conditions of a constant exchange of matter, energy and information with the external environment [18].The existence of a dissipative system is based on a constant synthesis of order and chaos.It has two aspects.Firstly, its system order exists only due to the chaos injected into the environment.Secondly, the system acquires the ability to adequately respond to chaotic influences of the external environment and thereby maintain its stability due to its system order.This remark primarily concerns structural stability preservation.It is possible if the system has a high level of orderliness.
Orderliness is the property of a system to preserve the structural composition, i.e., to restrict the freedom of change of both the system's element set and the connections between the elements.Usually, in a society, orderliness is a consequence of the practical implementation of functioning of a set of attitudes, laws and prohibitions that structure the system and organize its functioning within a clear system of rules.System orderliness is connected with its harmony.
The functioning of complex systems, which include all forms of transportation systems, is based on two global laws: the first law of dialectics and the law of the generalized golden ratio (GGR).
The first law of dialectics is the law of unity and the struggle of the opposites.Its essence can be expressed in the following form: Positive + negative = universe which, after normalization, is reduced to the form (1): where a-dominant; b-subdominant.
In the structural identification of the concept of system sustainability, we will use the first law of dialectics in form (1), taking an increase in utility as a positive and a loss of sustainability as a negative: The law of the generalized golden ratio (S-ratio) was identified by E. M. Soroko [29] and it has the following analytical form (2): where s-so-called multiplicity indicator.
In complex systems management, the law of the generalized golden ratio was the predecessor of modern F-technologies.In general, both laws are tools for building system functionality and structure close to the standard.The purpose of reference transport systems functionality is to form the general utility (GU) when consuming the provided resource.Thus, it is necessary to identify the general utility model.To do this, it is necessary to identify it structurally and parametrically.
2.2.Road Safety Is a System Property.Road Safety as a Result of Functioning of a Specialized System for Preventing Conflict Road Traffic Situations A variety of authors have considered road safety from a system perspective.In the period of 1950-1970 many authors brought the consideration of road safety to the system level [30][31][32][33][34][35].In many ways, this was promoted by the general theoretical work of L. Von Bertalanffy [36].Road safety specialists relied on this work and considered the problems of road traffic accidents not as a set of special cases, but as a demonstration of a system peculiarity.
The result of this work was the concept of "Vision Zero", as Sweden's Traffic Safety Policy, declared in 1997 [37].Two and a half decades of implementation of this concept, at first in the Scandinavian countries, and then, albeit partially, in many others, have shown the high efficiency of using the system approach in road safety provision.
P. Larsson et al. [38] believe that "the so-called zero-tolerance position, or Vision Zero approach, to road safety is built around two axioms; the system must be adapted to the psychological and physical conditions and limitations of the human being and the responsibility for road safety must be shared between the road-users and the designers and professional operators of the system".This is an important thesis that requires road traffic organizers to have a professional understanding of the cause-and-effect relationships of the process of forming of road traffic accidents and active actions to minimize the negative consequences in the chain of events.
A. Szymanek [39] has a similar opinion: "The road safety management methodology should be based on a system approach.This means that the road transport must be formalized as a complex system (CS), and then safety can be interpreted as an emergent feature of such a system.Road accidents should be interpreted as "organizational accidents".B. P. Hughes et al. systematize various approaches of road safety provision in [40].The analysis of 2620 literary sources conducted by the authors allowed them to summarize the experience in the road safety sphere in various countries in the form of 121 different types of models.This diversity of models for road safety provision is explained by a wide range of different targets and restrictions.This is an illustrative example of a wide variety of approaches, among which there may be the most unexpected ones.

General Trends of Road Safety Changes
At various times, scientists from all over the world were interested in the issues of road safety provision [41][42][43][44][45][46][47][48][49][50][51].The analysis of the content of most of these studies allows us to conclude about the positive global changes in trends in road safety.The statistics of such organizations as the World Health Organization [52], IRTAD Group and permanent working group on road safety of the International Transport Forum [53] confirm this fact.R. Elvik and R. Goel [54] state that "recent studies find a stronger tendency towards safety-in-numbers than older studies".This remark applies to the majority of the economically and socially developed countries of the world [55], but not to all of them.The economically and socially underdeveloped countries of Africa and Southeast Asia are characterized not only by a relatively high level of road traffic accidents, but also by trends in its increase.This spatial heterogeneity of road safety trends on a global scale is explained by the fact that different paradigms of road safety provision are currently relevant in different countries [56].This means that some developing countries now have a level of road safety provision at approximately the same level as it was 30-70 years ago, from 1950 to 1990 in currently leading countries in the field of road safety.
While we consider the general trends in road safety, we should understand that the attention of states to the issues of reducing road traffic accidents and deaths in road accidents is largely related to the socio-economic aspects of life, particularly to one of the most important-the damage caused by the decrease in the quality of citizens' human potential.Thus, it is shown in [57] that the total costs of road crashes in 31 European countries are equivalent to 0.4-4.1% of Gross Domestic Product.Cost per fatality varies in different European countries in the range from €0.7 million to €3.0 million.
According to L. J. Blincoe et al. [58], the economic cost of US motor vehicle crashes that occurred in 2010 totaled $242 billion.This is equivalent to 1.6% of the US Gross Domestic Product.
These facts indicate that the problem of road traffic accident rate is highly significant both from the social and economic sides for most countries of the world.This is the reason for a serious increase in attention to the issue of road safety that has been paid in recent years in almost all countries of the world.As a result, in the vast majority of countries of the world, the general trends in changes in road safety are quite encouraging.

The Road User Factor
Long-term experience  allowed us to establish that the contribution of the human factor to the road traffic accident rate is at least 90%.That is why autonomous car driving technologies have been rapidly developing in recent years [59][60][61][62], and the possibility of removing a person from the process of driving a vehicle has been actively studied [63].However, in recent decades, there has been a huge improvement in the road transport infrastructure in various countries [64][65][66].Many organizational and technical solutions are devoted to the minimization of the possibility of direct contact between a pedestrian and a driver [67,68].Another direction of road safety improvement involves administrative actions, such as reducing the speed of vehicles in cities [69] and monitoring drivers' behavior [70,71].In this regard, it is important to analyze the main approaches to assessing the driving style of drivers.
In 2004, O. Taubman-Ben-Ari et al. [72] published an article that proposed to classify the style of driving according to danger to others.The authors proposed to select eight typical driving styles.In the same article, the results of statistical research were presented, which aimed to identify a factor relationship between belonging to a particular driving style and the characteristic features of Israeli drivers who were representatives of different styles.
As indicated in the article by O. Taubman-Ben-Ari et al., the ideas of classifying drivers by driving styles are based on earlier works [73,74].
A little later, the method by O. Taubman-Ben-Ari et al. was modernized by Spanish [75] and Romanian authors [76].
In 2019, similar studies were conducted by Bulgarian researchers Z. Totkova and R. Racheva, who modernized the method and presented their version of factor analyses of the multidimensional driving style inventory (MDSI-BG) in [77].
The above studies were highly useful for understanding the psycho-physiological differences between different drivers and their different contributions to the formation of the actual road safety level.
A lot of works are also devoted to the issues of the impact of road users' sex and age on the road traffic accident rate [78][79][80].Almost all of these studies indicate a greater tendency to unjustified risk, which is a risk factor for accidents, in young inexperienced male drivers with a low level of education and internal culture [81][82][83].Conversely, women are more attentive and careful on the road than men [84].It is also established that more educated professionals with a reliable workplace perform significantly less hazardous driving techniques than uneducated unemployed people [85].There is also evidence of an inverse positive relationship between road traffic accidents and the general employment of people [86,87].

The Factor of Road Transport Infrastructure and Traffic Control Systems
Many works are devoted to the issues of the connection between road transport accidents and the quality of road transport infrastructure [88,89].Numerous studies have been devoted to this topic in the context of the Third Paradigm of road safety [56].Today, most developed countries make the highest requirements to the quality of the roadway, as well as to the engineering equipment of highways and the road network of cities.Nevertheless, even today, research on the further improvement of engineering systems for road safety provision is being conducted.One of the most important areas of such research is the identification of potentially dangerous locations for road users in the city.On the basis of these locations, new engineering solutions contributing to greater security will be developed and implemented [90][91][92][93][94][95][96][97][98].
The analysis of these works shows that creation of a high-quality road network is an expensive and resource-intensive task.Its solution has many limitations.In this regard, in the current situation, one of the most efficient ways to increase traffic safety is to identify dangerous locations, notify drivers about their existence and increase control of driving behavior on these sections of the road network.

The Factors of the Technical Level and Quality of the Vehicle
Researchers have been paying attention to car safety issues for more than 100 years.In this direction, especially effective research was carried out in the period from 1960 to 2010.During this time, cars have become much safer and much less demanding on the driver's qualifications.The level of active, passive and post-accident safety of individual cars has significantly increased.Duplication of safety systems, decrease in the probability of an accident and decline of the injury risks in an accident are the results of many years of work of designers, production technologists and researchers in the field of road safety.According to [99,100], the poor technical condition of vehicles and related technical failures were the cause of only 1-1.5% of accidents in EU countries (but 5-7% of accidents were caused by a combination of the technical failure of the car and incompetent actions of the driver in this critical situation).

The Factors of the External Environment That Negatively Affect Traffic Conditions
Many studies have been devoted to the influence of negative manifestations of the external environment on road safety [101][102][103].In general, summarizing the results of these studies, it can be concluded that in adverse weather and weather-related road conditions, the accident rate increases significantly.That is why in countries with a high level of socio-economic development, active road signs and information boards are used on federal highways and in cities to regulate traffic flow modes taking into account environmental conditions [104][105][106].This significantly reduces the risk of accidents.Unfortunately, such advanced technologies are used in a small number of countries.In most countries, the driver chooses the speed mode independently, without external informational support [107].

COVID Lockdown and Its Influence on Road Safety
In 2020 and 2021, a lot of works were published on the study of the impact of the COVID lockdown on changing people's transport behavior and, as a result, on changing the characteristics of road traffic accidents in various countries of the world [108][109][110][111][112][113].
Studies by [108,109] are devoted to changes in the patterns of using public and individual transport in the conditions of COVID-19.It was established that in 2020, people were afraid of the risks of infection and significantly reduced the frequency of using transport services.The same conclusion was made with regard to the use of individual cars and the mobility of citizens in general.
A study by Indian authors [110] shows that during the spring lockdown period (from 24.03.2020 to 31.06.2020), in comparison with the same period of 2019, India recorded a decrease in the number of road accidents by 68% and in the number of the dead and injured in road accidents by 62% and 72%, respectively.
In [112], the Greek authors, on the basis of studies conducted in Greece and Saudi Arabia in March-April 2020, showed that "reduced traffic volumes due to lockdown, led to a slight increase in speeds by 6-11%, but more importantly to more frequent harsh acceleration and harsh braking events (up to 12% increase) as well mobile phone use (up to 42% increase) during March and April 2020, which were the months where COVID-19 spread was at its peak.On the bright side, accidents in Greece were reduced by 41% during the first month of COVID-19-induced measures and driving in the early morning hours (00:00-05:00), which are considered dangerous, dropped by up to 81%".
Summarizing the results of the above studies, it can be concluded that road safety during the lockdown restrictions in various countries of the world was largely determined by the peculiarities of the manifestation of transport activity.Based on these observations, L. Budd and S. Ison made conclusions about the need to develop new concepts and strategies of transport policy [113].The main idea of these authors is that it is necessary to increase awareness of the choice of the method of mobility from the perspective of health safety in the future.

Identification of the Concept of "Reference Systems"
Reference systems are open, interacting with the external environment, multicomponent systems with a number of specific features such as:

•
they work in the mode of generalized golden ration (GGR) [29]; • their normalized utility function is described by the Equation (3): where GUthe normalized general utility; x-the share of consumed resource; g = 1 + s-Q-factor-system quality indicator; s-the multiplicity parameter, independent of x.
• they have real (independent of argument x) values of the quality factor g.

Functionality of Reference System
The essence of the functionality of reference systems of urban mobility management is to deliver the general utility of GU in the consumption of the provided resource.The macro-model of the normalized general utility GU n is represented as (3).A significant advantage of this macro-model type is that it has fractal properties.Suppose that the argument x is an analogous function of the new parameter T n , i.e., then normalized general utility GU n can be represented as: Obviously, the new result ( 3) is again a fractal that has normalized time T n = t/T as an argument (here T is a lifecycle time).The creation of the new model (3) allows us to proceed to the analysis of its dynamic characteristics.

Analysis of the Dynamic Characteristics of Reference Systems' Functionality
For the development of an algorithmic provision for road safety management systems it is necessary to have models of the basic dynamic characteristics of the information process.Article [114] is devoted to the topic of the development of such provision.
The functional scheme of the process (Figure 1) identifies dynamic models of the following process characteristics: speed of the change of process utility SGU (speed of general utility): process utility GU (general utility): vital power (performance) of process VP (vital power): reserve of process vital power S (stamina): Mathematics 2021, 9, x FOR PEER REVIEW 8 of 35 then normalized general utility GUn can be represented as: Obviously, the new result ( 3) is again a fractal that has normalized time as an argument (here T is a lifecycle time).The creation of the new model ( 3) allows us to proceed to the analysis of its dynamic characteristics.

Analysis of the Dynamic Characteristics of Reference Systems' Functionality
For the development of an algorithmic provision for road safety management systems it is necessary to have models of the basic dynamic characteristics of the information process.Article [114] is devoted to the topic of the development of such provision.
The functional scheme of the process (Figure 1) identifies dynamic models of the following process characteristics: speed of the change of process utility SGU (speed of general utility): process utility GU (general utility): vital power (performance) of process VP (vital power): reserve of process vital power S (stamina): Let us introduce the normalization of the above mentioned dynamic characteristics (Table 1).Table 1.Dynamic models of system characteristics.

Model
Determination Algorithm Normalized speed of the change of process utility SGUn ( ) Normalized process q Figure 1.Functional scheme of the road safety management process (system).
Let us introduce the normalization of the above mentioned dynamic characteristics (Table 1).
Table 1.Dynamic models of system characteristics.

Model Determination Algorithm
Normalized speed of the change of process utility Normalized process general utility Normalized reserve of process vital power Let us explain the models (10)(11)(12)(13).
From the perspective of the quality of system changes management process it is preferable to have not only a high level of technology perfection (determined by the level of general utility GU(T n )), but also a high growth rate SGU p (T n ) = dGU(T n )/dT n .
In this regard, we can take as a criterion of effectiveness (vital power) of the system management process the next value (14): where From the engineering point of view the question about the extremum of vital power VP(T n ) is interesting.Its maximum is reached under the condition (15): The solution of the Equation ( 15) is ( 16): The maximum possible value of the process vital power VP is reached under the condition (17): The normalized value of the process vital power VP n is determined by the algorithm (18): In Table 1 model ( 18) is denoted as (12).Argument T n opt , corresponding to the maximum of effectiveness of the controlled process, is determined by (19): The reserve of vital power of population mobility management process can be interpreted as cumulative sum of vital powers on the interval T n (20): The level of reserve of vital powers at the end of the lifecycle is (21): Introduction of the normalization generates model ( 22): In Table 1 model ( 22) is denoted as (13).

The Research of the Structural Features of Utility Function of System Management
When analyzing the structure of general utility GU in multicomponent systems we should take into account two main points: • the general utility GU of a complex system is an additive composition of individual functions of utility of system components; • the general utility GU of urban mobility management is a consequence of interaction between chaotic state and orderliness of processes in this sphere.
When assessing the general utility GU of a complex system we need to identify the contribution of each component to the total success.A Pareto chart (PC) is a traditional tool for this type of analysis.It represents a cumulative function of the contribution of individual components into a total balance.As a rule, PC is approximated by a function (23): where x-normalized rank; x = r i /r max ; r i and r max -correspondingly current and maximal ranks; g = 1 + s-Q-factor-system quality indicator; g and s-indicators.
Since PC is the increasing sum of "weights" of system components, it allows us to analytically deduce "weight" coefficients as (24): Essentially, the normalized Pareto chart also represents the function of general utility GU, where x is a share of the consumed resource.A comparison of models ( 3) and (23) indicates their functional equivalence PC(x) = GU(x).
Found "weight" coefficients ω i are used to solve two significant problems: the estimation of relative structural entropy H n [29,115,116], characterizing the level of orderliness (and structural perfection) of a process or system (25): representation of utility function GU(x) in additive form (26): where GU i (x)-an individual utility function (the share of contribution of i-component to the final system result).
Let us examine the second approach in detail.Its specifics are generated by the wide usage in modern computer-aided manufacturing systems of F-technologies (F-Fibonacci).F-technologies exploit Fibonacci algorithms to build highly efficient and harmonic technical, socio-economic and organizational systems.One of the leading experts in this sphere is E. M. Soroko [29,115,116].He found that in the majority of situations the principle of generalized golden ratio (GGR) works when the proportion between parts of the unit interval f r and f s (meeting the condition of normalization f r + f s = 1) is right (27): From ( 27) it follows, that (28): where f r -dominant; 1 − f r = f s -subdominant; g = 1 + s-Q-factor-system quality indicator g and s-indicators.
The generalized golden ratio regime (GGR regime) directly relates to the functionality of reference management systems, presented as a macromodel of normalized general utility GU n (3).Its graphic representation as Pareto chart PC(x) (taking into account equivalency PC(x) = GU(x)) is shown in Figure 2.
From ( 27) it follows, that (28):   The GGR regime is presented in the figure as a diagonal y = (1 -x), where every point meets the condition y + x = 1.The point of intersection of the diagonal and Pareto chart (with coordinates f s and f r ) is called the first Pareto point.We want to highlight its properties.Firstly, it complies with the GGR principle.Secondly, it is related to indicator g (system quality factor) by the correlation (29): Thirdly, it is predetermined by proportion «order/chaos» in the management system.The meaning of the dominant f r and f s the subdominant is fundamentally important.The approach, evaluating levels of orderliness and chaos in an informational metric, where the maximal relative entropy of combination is taken as a unit, is foundational.The use of such approach allowed the creation of the theory of systems harmony [29].However, at the same time, there is some non-transparency of the harmony concept in terms of engineering of efficient systems of urban mobility management.Let us consider the methodology of the optimal synthesis of harmonic systems in the examined domain sphere.
If we assume that the target management function is the retention of a process at the first Pareto point, then we need to prove that it provides an optimum in terms of predetermined parameter; in our case by default it is harmony.However, the harmonicity concept is not clearly formalized yet, and the question of an acceptable metric remains open.It is obvious that the total process utility in the first Pareto point is higher, the larger value of parameter g.However, the choice of the optimal value of the indicator g is not simple.With a decline in the value of indicator g, system losses should increase due to the poor quality of the system.With the growth of the indicator g the system efficiency should increase and, as a consequence, the additional expenses on the implementation of a system's functioning would grow.This problem belongs to the optimization category, where we need to establish which conditions provide minimal total expenses Z and what is the value of the reached extremum Z min .
The solution to this problem is given in the work [117].A short explanation is given below.
In the paper it is supposed that the total expenses Z are caused by the simultaneous effect of two mechanisms: • on the one hand, a low value of the indicator g indicates a low system efficiency and, as a consequence, it is related to unwanted losses Z 1 ; • on the other hand, a growth of the value g is impossible withoud additional expenses Z 2 .
Equation (28), written in the form 1 − f = f g (f -is a dominant of GGR regime), reaches a compromise.The left part characterizes the negative of decline of the f dominant, while the right part characterizes the positive, provided by its growth.
The negative (1 − f ) leads to proportianal losses Therefore total expenses are (28): It is important that the price C 1 can have a broad meaning.Price C 1 can be understood as money, energy, information, time and other resource types.The essence of article [117] consists in the fact that the harmonic regime is the regime of provision of the minimal total resource expenses in the first Pareto point.
It is established in [117] that the level of total expenses Z is (31): where m = 1 .
The normalized value of expenses equals to Z n (32): The minimum of expenses is reached under the condition dZ n /d f r = 0. Value of extremum Z n min equals to (33): and the level of minimal total expenses equals to (34):

Establishment of the Relation between Functionality and Structure
Associating the system functionality with the function of general utility GU, it was detected that the key characteristic of GU is the indicator g, which determines the level of success (the positive) in the consumption of the available resource and the structure of "weight" coefficients (relation ( 23)).However, the "weight" structure in turn predetermines the level of relative entropy H n (25).Therefore, the relation between g and H n is quite expected (Figure 3).
Associating the system functionality with the function of general utility GU, it was detected that the key characteristic of GU is the indicator g, which determines the level of success (the positive) in the consumption of the available resource and the structure of "weight" coefficients (relation ( 23)).However, the "weight" structure in turn predetermines the level of relative entropy Hn (25).Therefore, the relation between g and Hn is quite expected (Figure 3).

Identification of Model of Evaluation of Reference Systems' Efficiency
The successful solution of management problems is based on competent goal-setting, the quantitative measure of which is an efficiency.Formally, efficiency is a relation of performance to the resource expenses on results' achievement.
Let us consider the interpretation of the concept of "performance".
Traditionally, success is associated with the dominant fr; the higher the dominant, the more positive the result.Since fr is a measure of the total utility, as a consequence, it can be a characteristic of the performance of the examined process.Model GU = f(g) is a useful tool for success evaluation.It indicates that the quantitative characteristic of management performance can be presented both as general utility GU and as a system (process) characteristic g, defined as a quality factor or Q-factor.

Identification of Model of Evaluation of Reference Systems' Efficiency
The successful solution of management problems is based on competent goal-setting, the quantitative measure of which is an efficiency.Formally, efficiency is a relation of performance to the resource expenses on results' achievement.
Let us consider the interpretation of the concept of "performance".
Traditionally, success is associated with the dominant f r ; the higher the dominant, the more positive the result.Since f r is a measure of the total utility, as a consequence, it can be a characteristic of the performance of the examined process.Model GU = f(g) is a useful tool for success evaluation.It indicates that the quantitative characteristic of management performance can be presented both as general utility GU and as a system (process) characteristic g, defined as a quality factor or Q-factor.
Therefore, there are two possible options for performance evaluation.Meanwhile, it is important to correctly evaluate resources expenses on the result achievement.
Analysis of efficiency of reference systems.If we define efficiency as a relation of performance to resource expenses, then we can consider two alternative options: The analysis showed that the second option, i.e., analysis of the system efficiency with regard to the criterion 2 (36), is more preferable than the first, because it provides higher elasticity EF relative to the total utility GU.The relation EF 2m /EF 1m = g/ f r in the mode of the golden ratio (when g = 2; f r = 0.618) is equal to 3.23.
The plot of dependency of efficiency EF 2 on dominant f r for reference systems (with different levels of quality factor Q-factor estimated g) is shown in Figure 4.
In the analysis of the model of process (system) efficiency, presented in Figure 4, it may seem as if the growth of the efficiency of reference systems is not limited at all.In fact, with the growth in the quality factor Q-factor estimated g, the system cost is increasing and according to the paper [118], this cost is proportional to g. Due to the high elasticity of the process efficiency relative to Q-factor estimated g, the system sustainability stock is decreasing and the possibility of bifurcation is increasing.The issue of the boundary of the efficient systems' sustainability is quite important but it requires individual consideration.
The analysis showed that the second option, i.e., analysis of the system efficiency with regard to the criterion 2 (36), is more preferable than the first, because it provides higher elasticity EF relative to the total utility GU.The relation in the mode of the golden ratio (when ) is equal to 3.23.The plot of dependency of efficiency EF2 on dominant fr for reference systems (with different levels of quality factor Q-factor estimated g) is shown in Figure 4.In the analysis of the model of process (system) efficiency, presented in Figure 4, it may seem as if the growth of the efficiency of reference systems is not limited at all.In fact, with the growth in the quality factor Q-factor estimated g, the system cost is increasing and according to the paper [118], this cost is proportional to g. Due to the high elasticity of the process efficiency relative to Q-factor estimated g, the system sustainability stock is decreasing and the possibility of bifurcation is increasing.The issue of the boundary of the efficient systems' sustainability is quite important but it requires individual consideration.

Conceptual Model of Reference System Analysis
Key properties of any rational systems are integrity and goal-orientation.System goals can differ but their main humanistic destination is a creation of some positive (utility).Management of complex systems is directed towards the growth of the positive.However, there are always forces that prevent the growth.Counteraction of multidirectional vectors forms dialectical results of complex systems' functioning.The antagonistic force of the utility is a system's sustainability.
With the growth of the consumed resource, the utility should increase, but at the same time the remaining resource is declining and this fact determines system sustainability and its ability to stably function.Therefore, the utility growth is not unlimited.In this situation the suggestion about some optimal (harmonious) system working mode

Conceptual Model of Reference System Analysis
Key properties of any rational systems are integrity and goal-orientation.System goals can differ but their main humanistic destination is a creation of some positive (utility).Management of complex systems is directed towards the growth of the positive.However, there are always forces that prevent the growth.Counteraction of multidirectional vectors forms dialectical results of complex systems' functioning.The antagonistic force of the utility is a system's sustainability.
With the growth of the consumed resource, the utility should increase, but at the same time the remaining resource is declining and this fact determines system sustainability and its ability to stably function.Therefore, the utility growth is not unlimited.In this situation the suggestion about some optimal (harmonious) system working mode occurs.The traditional approach supposes hidden (latent) work of some global mechanism, called the golden ratio or golden mean.
Lately the principle of the general golden ratio (GGR), established by Belarusian scientist E. M. Soroko [29] started to be used for analysis.Researchers of the V.A. Trapeznikov Institute of Management Problems of the Russian Academy of Science [119,120] created new scientific field: F-technologies (F-Fibonacci), successfully solving complex application problems of efficient management in the systems of various purposes.
The conceptual model (Figure 5) that we use also relies on the GGR principle and the first dialectic law.This model defines two new entities: normalized total utility and normalized sustainability.Lately the principle of the general golden ratio (GGR), established by Belarusian scientist E. M. Soroko [29] started to be used for analysis.Researchers of the V.A. Trapeznikov Institute of Management Problems of the Russian Academy of Science [119,120] created new scientific field: F-technologies (F-Fibonacci), successfully solving complex application problems of efficient management in the systems of various purposes.
The conceptual model (Figure 5) that we use also relies on the GGR principle and the first dialectic law.This model defines two new entities: normalized total utility and normalized sustainability.The research into these characteristics allows us to solve two important problems: the problem of the optimization of the system utility GU and the problem about the logical consequences generated by this optimization (casual, structural and functional system harmonies).

Definition of the System Utility Function
The system utility is characterized by the level of system functionality.Analytically it is set by the model of total utility GU.
In the structural identification of this model, we will assume that the first dialectic law works: Positive + negative = universum.
We interpret the positive as a saved resource and the negative as an under-received utility.After normalization we have a + b = 1 where a-dominant, a = 1 − x c ; x c -the share of resource consumed on goal achievement; b-subdominant; b = 1 − GU; GU-the share of received utility relatively to maximal possible.
The GGR principle works as follows: which gives (37): i.e., where g-Q-factor-system quality indicator.
The plots of total utility functions for different values of the quality control g are presented on Figure 6.
In Figure 6 the point C, corresponding to the classical golden ratio (under the condition g = 2), is highlighted.All sustainable solutions, conformed to GGR, are placed on the diagonal The analysis of extremely possible utility PU (37) showed that GU meets the necessary requirements, i.e., it declines with the growth of consumed resource x c (39): where g-Q-factor-system quality indicator.
The plots of total utility functions for different values of the quality control g are presented on Figure 6.

Figure 6. Plots of total utility functions GU(xc).
In Figure 6 the point C, corresponding to the classical golden ratio (under the condition g = 2), is highlighted.All sustainable solutions, conformed to GGR, are placed on the diagonal The analysis of extremely possible utility PU (37) showed that GU meets the necessary requirements, i.e., it declines with the growth of consumed resource c x (39):

Optimization of the Total Utility (Criteria of Optimality, System Harmonization)
Taking into account the essence of the first dialectic law, we should expect a confrontation between the system positive and the system negative.The sustainable balance between them forms the system balance.Many authoritative scientists [121][122][123] suppose that system balance is identical to the concept of system harmony.In that point (denoted as opt x ) two conditions should simultaneously meet:

Optimization of the Total Utility (Criteria of Optimality, System Harmonization)
Taking into account the essence of the first dialectic law, we should expect a confrontation between the system positive and the system negative.The sustainable balance between them forms the system balance.Many authoritative scientists [121][122][123] suppose that system balance is identical to the concept of system harmony.In that point (denoted as x opt ) two conditions should simultaneously meet:

•
with regard to the first dialectic law, a + b = 1, where a = 1 − x c ; b = 1 − GU), the next condition should be satisfied • with regard to GGR, GU = 1 − (1 − x c ) g .Therefore, in the optimum point (where which implies that the quality factor equals to (40): In the regime of classical golden ration, when g = 2, x opt = 0.382, substitution (40) into (37) gives an important result (41): Therefore, It means that: optimal solutions are always placed on the diagonal (42); -system quality factor in the optimum point equals to g = ln(x opt ) ln(1−x opt ) ; -work state in optimum should be considered as harmonious system state; -"weights" ω i of system components, as shown in work [124], are predetermined by the system quality factor g (43): where n-the number of components of the system under study; i-the rank of the component.

Evaluation of the Complex Systems' Sustainability
In recent times, in strategic management publications [125] the concept of sustainability is broadly used.Different authors variously interpret this concept [126].The majority define it only from the economic or ecological point of view, but in fact this concept has a wider meaning.There is no unified point of the view on this issue yet.At the same time the sustainability concept is significant in system analysis.It predetermines the successfulness of functioning of an analyzed complex system.
Sustainability can be evaluated by different metrics.We will use the metric of the generalized golden ratio (GGR).As was shown earlier in Section 2.1: The positive (P) assumes the growth of system utility, while the negative (N) is defined by the loss of system sustainability.
In the growth of the system utility the first summand-the Positive (P)-is dominant and after the normalization it will take the form of GU opt .GU opt is total utility in the operating point.The second summand-the Negative (N)-is a subdominant and it defines the share of the lost sustainability (1 -ST opt.); here ST opt. is a normalized sustainability in the operating point.
We will use GGR in the form 1 GU opt s = GU opt 1−ST opt , and as a result we will have (45): (45) where g-Q-factor-system quality indicator; g = s + 1.
Analysis of the formula (45) allows us to make a significant conclusion: The growth of the utility is not unlimited: it is restrained by the loss of the sustainability.Sustainability (as an ability to execute prescribed functionality) is defined by the remaining resource and with its decline it also declines.

Establishment of the Relation between System Sustainability and Its Utility
The relation between the sustainability and the utility in the operating point (i.e., x c = x c opt ) is characterized by the relation (45).In points x c the structure of this relation is inherited and the sustainability model can be represented in the form (46): Plots of the functions ST(x c ) and GU(x c ) are presented in Figure 7.
) is characterized by the relation (45).In points xc the structure of this relation is inherited and the sustainability model can be represented in the form (46): Plots of the functions ST(xc) and GU(xc) are presented in Figure 7.

Figure 7. Plots of functions GU(xc) and ST(xc).
Operating points, corresponding to values of the quality factor g = 2, 3, 4 and 5 are highlighted.The pair GU2 and GT2 belongs to the golden ratio regime (g = 2).In the operating point the condition GU2(xc) = STc(xc) = 0.618 is met.Operating points, corresponding to values of the quality factor g = 2, 3, 4 and 5 are highlighted.The pair GU 2 and GT 2 belongs to the golden ratio regime (g = 2).In the operating point the condition

Analysis of the Possibilities of System Sustainability Application
There is a series of spheres in which arises a need to evaluate the sustainability level.From the positions of the PEST analysis [125] there are several types of sustainability: In cases of sustainable urban mobility development, we should use vector (generalized) sustainability evaluation (47): where ST i -sustainability in thei-sphere; ω i -the "weight" coefficients, meeting the normalization condition

Harmonization of Complex Systems
Usually, harmony is considered from three points of view [29].Firstly, it is considered from the positions of the character and mechanism of harmonious impact of one system part on another (casual harmony).Secondly, it is considered as proportionality of the parts in the whole (structural harmony).Thirdly, it is considered as the mutual consistency and synchronicity of system components functioning (functional harmony).The authors agree with E. M. Soroko's point of view [29,[115][116][117], consisting of the observation that system harmony is observed in the regime of the classical golden ratio, when the system quality factor equals g = 2.

The System Orderliness and Its Relation to the System Quality Factor
In Section 2.1 it was shown that the base of the system orderliness is the structural sustainability.Maintenance of the system structural sustainability in the conditions of the negative impact of chaos is a consequence of the high level of the system orderliness.The system orderliness is identified by the features of the "weight" distribution Pareto charts or the system/process elements' significance in the formation of the system functioning result.It is quantitatively evaluated by the level of the curvature of the Lorenz curve of the Pareto chart.
The normalized Pareto chart is presented as a function of the normalized generalized utility GU n , where x is the share of the used resource (Equation ( 1)).As it was shown above in Section 3.4, the Pareto chart is an increscent sum of "weights" ω i of the components of the considered system.As it was shown in Equation ( 22), the "weight" coefficients ω i are defined by the equation: In addition, earlier in Section 2.1, in Equation ( 23) it was presented that the found "weight" coefficients ω i relative structural entropy H n , characterizing the level of orderliness (or structural perfection) of the process or the system: The utility function GU(x) in the additive form can be presented as the sum of particular utility functions (or the sum of the shares of the contribution of i-components in the system result) (48): In this equation, x is a share of the used resource and g identifies the system quality factor.The determination of the relative structural entropy H n is important from the positions of the identification of the level of the system orderliness and the system quality factor g.

Method of Definition of Relative Entropy of Road Safety Provision System
Formula (41) connects "weights" ω i of researched system components and the quality factor g. Consequently, this raises the question of the system itself, its components, their "weights" and the relative importance of each component to the formation of the system functioning result.
To consider the road safety provision system we need to define the size of the system (federal, regional, district, urban or microdistrict) and the cause-and-effect relationship of road accident rate formation.

The Cause-and-Effect Model of Road Accident Rate Formation
The cause-and-effect model for a road safety provision system is shown in Figure 8.The goal of the analysis of the cause-and-effect model (Figure 8) is to establish the relationship between coefficients of the informational transfer of elements of the road accident rate formation process.
For the research into the examined process structure, we will use cybernetic modeling.The concept of cybernetic modeling is defined by V.I.Krutov et al. [127] in the following way: "cybernetic models are based on the receiving of relations between input and output function for some black box or grey box, representing examined an phenomenon, without disclosure of its inner structure".Formula (41) connects "weights" i  of researched system components and the qual- ity factor g. Consequently, this raises the question of the system itself, its components, their "weights" and the relative importance of each component to the formation of the system functioning result.
To consider the road safety provision system we need to define the size of the system (federal, regional, district, urban or microdistrict) and the cause-and-effect relationship of road accident rate formation.

The cause-and-effect model of road accident rate formation
The cause-and-effect model for a road safety provision system is shown in Figure 8.The goal of the analysis of the cause-and-effect model (Figure 8) is to establish the relationship between coefficients of the informational transfer of elements of the road accident rate formation process.
For the research into the examined process structure, we will use cybernetic modeling.The concept of cybernetic modeling is defined by V.I.Krutov et al. [127] in the following way: "cybernetic models are based on the receiving of relations between input and output function for some black box or grey box, representing examined an phenomenon, without disclosure of its inner structure".
We should pay attention not to the inner mechanism of each component of the examined process chain but only to coefficients Ki of the informational transfer between elements of the cause-and-effect model.Ki plays a special role in the analysis of the process model.It comprises the key characteristics of the elements of the complex system, the model of which is entirely defined by the behavior of its elements.The coefficient of the transfer of the i-element of the cause-and-effect model is a relation between its output (Aout) and input (A in): We should pay attention not to the inner mechanism of each component of the examined process chain but only to coefficients K i of the informational transfer between elements of the cause-and-effect model.K i plays a special role in the analysis of the process model.It comprises the key characteristics of the elements of the complex system, the model of which is entirely defined by the behavior of its elements.The coefficient of the transfer of the i-element of the cause-and-effect model is a relation between its output (A out ) and input (A in ): With regard to the structure of the road accident rate mechanism (Figure 8) we distinguish 3 subprocesses:

•
the formation of the vehicle fleet, determining the average annual intensity of road traffic (with the transfer coefficient K A = N/P); • the formation of the road accidents (with the transfer coefficient K RA = RA/N); • the formation of the deaths rate in road accidents, the number of the deceased in road accidents (with the transfer coefficient K F = F/RA);

Determination of the Priorities of the Process Elements (ABC Analysis)
The ABC analysis, based on the Pareto chart, is a tool of determination of the priorities of the examined process elements.The aim of the analysis is to identify main priorities.The essence of the ABC analysis consists of the search for the first and the second Pareto points (determining the position of the boundaries between A, B and C areas).
To hold this analysis, it is necessary to preliminarily make a decision on 3 issues: decide what to define as the positive; choose necessary metric; prepare data.
For the first problem, it is obvious that the lower the transfer coefficient K of the element, the more positive the result.Generally, for the whole cause-and-effect process chain we can say that the lower the value of the through transfer coefficient K i , the more positive the system result.
In the next step we need to learn how to measure the positive, i.e., introduce the necessary metric.It is reasonable to use an indicator of the process Q positive as a measure (49): Received results allow us to find "weak" elements in the chain of the cause-and-effect relations in the road accident rate formation process and set priorities in the road safety management sphere.

Evaluation of the "Weight" Coefficients ω i
Knowing the values of the positive of process elements Q i , we determine the total process positive Q (50): The "weight" coefficients ω i show the share of each summand Q i in the total sum of the process positive Q (51): The structure of the "weight" coefficients ω i of the road safety provision process defines the level of the orderliness.The numerical characteristic of the orderliness level is an entropy H (52): where n-the number of system elements (in our case n = 3); ω i -the "weight" coefficients, meeting the normalization condition Relative entropy H n is more visual and convenient from the positions of the application characteristic of the evaluation of the road safety provision system's orderliness ( 53): where n-the number of system elements (in our case n = 3); From now on we will use this characteristic to evaluate the orderliness of the road safety provision systems.

The Method of the Determination of the Quality Factor of the Road Safety Provision System
As was shown in Sections 2.4 and 4.1, value g is an exponent of the exponential equation, identifying the curvature of the change in the Pareto chart cumulate.Equation (27), presented in Section 3.4, illustrates the meaning of the system quality factor g: By knowing the coordinates f r and f s of the curve, which identify the utility function GU(x), we can determine the value of the quality factor of the system.

Example of the Theoretical Solution of the Problem of the Evaluation of the System Harmonicity in Relation to the Road Safety Provision System
As an example, we will provide results of the harmonic analysis of the complex systems in the road safety sphere, one of the key spheres of the urban mobility implementation.In this analysis the results of the authors' research relied on the database of the State Inspectorate for Road Traffic Safety [11].
The three "weight" coefficients of the rank Pareto chart are the subject of the analysis: "weight", characterizing the level of automobilization (ω A ); "weight", characterizing the level of road accidents (ω RA ) and "weight", characterizing the death rate (ω F ).The sum of "weights" meets the normalization condition 3 ∑ i=1 ω i = 1

Structural Harmony
In reference systems, the "weight" coefficients are always determined (see relation (43), Figure 9).

Casual Harmony
In the analysis it was taken into account that the level of automobilization A is the root cause, while the number of road accidents and the number of deceased in road accidents are consequences of the automobilization A. Therefore, it is reasonable to present the casual harmony in the form of the dependency of the "weights" RA  and F  on the "weight" A  (Figure 10).In the golden ratio regime, when g = 2, components' "weights" (Figure 9) correspond to the proportions (44) and when n = 3 they are correspondingly equal: ω A = 0.111; ω RA = 0.555; ω F = 0.333.

Casual Harmony
In the analysis it was taken into account that the level of automobilization A is the root cause, while the number of road accidents and the number of deceased in road accidents are consequences of the automobilization A. Therefore, it is reasonable to present the casual harmony in the form of the dependency of the "weights" ω RA and ω F on the "weight" ω A (Figure 10).

Casual Harmony
In the analysis it was taken into account that the level of automobilization A is the root cause, while the number of road accidents and the number of deceased in road accidents are consequences of the automobilization A. Therefore, it is reasonable to present the casual harmony in the form of the dependency of the "weights"   The relation of the "weights" is of practical interest.It is presented in Figure 11.

Functional Harmony
According to Socrates [29], we connect functional harmony with the utility principle.The optimal level of the normalized utility is observed in the golden ratio conditions (Figure 12) and it is equal to 0.618. for the road safety provision system case.

Functional Harmony
According to Socrates [29], we connect functional harmony with the utility principle.The optimal level of the normalized utility is observed in the golden ratio conditions (Figure 12) and it is equal to 0.618.The data source was the official statistics of the State Inspectorate for Road Traffic Safety of the Russian Federation [11].Over the past 15 years (2006-2020), the indicators of road safety in Russia have been constantly improving.A number of factors have contributed to this, such as a general increase in the transport culture of the country's population, a significant increase in the quality of the country's transport fleet and significant qualitative growth of the country's road transport infrastructure.It is also important that since 2006, a lot of attention has been paid to the issues of road safety provision by the Russian government [9,10]. for the road safety provision system case.

Functional Harmony
According to Socrates [29], we connect functional harmony with the utility principle.The optimal level of the normalized utility is observed in the golden ratio conditions (Figure 12) and it is equal to 0.618.The data source was the official statistics of the State Inspectorate for Road Traffic Safety of the Russian Federation [11].Over the past 15 years (2006-2020), the indicators of road safety in Russia have been constantly improving.A number of factors have contributed to this, such as a general increase in the transport culture of the country's population, a significant increase in the quality of the country's transport fleet and significant qualitative growth of the country's road transport infrastructure.It is also important that since 2006, a lot of attention has been paid to the issues of road safety provision by the Russian government [9,10].
The result of the improvement in the overall situation in the transport sphere was a decrease in the annual number of accidents and deaths in road accidents in Russia over the past 15 years (Table 2).The natural result of a long-term decrease in road traffic accident rate in Russia was a significant decrease in the relative entropy H n of the road safety system and an increase in its Q-factor g (Table 3).The dynamics of the change in Q-factor g of the road safety system of the Russian Federation from 2006 to 2020 is very impressive.
At the same time, the ratio of "weight" coefficients ω i changes over time.Based on these dynamics, it is possible to estimate changes in all three types of system harmony: structural, causal and functional.
Table 4 presents data on the change in the numerical values of the "weight" coefficients ω i from 2006 to 2020. Figure 13 shows the graphs of the dynamics of the "weight" coefficients of the road safety system in Russia from 2006 to 2020.
Figure 14 shows a comparison of the theoretical and actual dynamics of changes in the "weight" coefficients ω i of the road safety system in Russia over a 15-year period.
The analysis of the dynamics of the values of the weight coefficients ω i allows us to conclude that the value of "weight" ω A has decreased almost by half over the past 15 years, the value of "weight" ω RA has increased by 12%, and the value of "weight" ω F has increased by 6% from 2016 to 2018.However, in 2020, the "weight" ω F significantly dropped again.This happened due to a significant increase in the severity of road accidents in 2020.The reason is the increase in the speed of vehicle movement in Russia in 2020 in the conditions of a relative rarefaction of the traffic flow associated with pandemic restrictions on the population mobility. has increased by 6% from 2016 to 2018.However, in 2020, the "weight" F  significantly dropped again.This happened due to a significant increase in the severity of road accidents in 2020.The reason is the increase in the speed of vehicle movement in Russia in 2020 in the conditions of a relative rarefaction of the traffic flow associated with pandemic restrictions on the population mobility.This is also evidenced by the statistics of penalties for violating the requirements of traffic rules and the share of those that fall on violation of the requirements of the speed limits (Table 5). has increased by 12%, and the value of "weight" F  has increased by 6% from 2016 to 2018.However, in 2020, the "weight" F  significantly dropped again.This happened due to a significant increase in the severity of road accidents in 2020.The reason is the increase in the speed of vehicle movement in Russia in 2020 in the conditions of a relative rarefaction of the traffic flow associated with pandemic restrictions on the population mobility.This is also evidenced by the statistics of penalties for violating the requirements of traffic rules and the share of those that fall on violation of the requirements of the speed limits (Table 5).This is also evidenced by the statistics of penalties for violating the requirements of traffic rules and the share of those that fall on violation of the requirements of the speed limits (Table 5).The data analysis in Table 5 shows that over the past 5 years in Russia, not only the number of penalties for violating traffic rules, but also the share of speed violations has increased.This indirectly indicates an increase in the quality of Russian roads and an improvement in the system of monitoring the speed limits of traffic flows.It is known that at the end of 2020, about 19,000 automated complexes for photo-video recording of traffic violations by drivers were used in Russia.
Figure 15 shows a comparison of the theoretical and actual dynamics of causal harmony of the road safety provision system in Russia over a 15-year period.The data analysis in Table 5 shows that over the past 5 years in Russia, not only the number of penalties for violating traffic rules, but also the share of speed violations has increased.This indirectly indicates an increase in the quality of Russian roads and an improvement in the system of monitoring the speed limits of traffic flows.It is known that at the end of 2020, about 19,000 automated complexes for photo-video recording of traffic violations by drivers were used in Russia.
Figure 15 shows a comparison of the theoretical and actual dynamics of causal harmony of the road safety provision system in Russia over a 15-year period.The functionality of the road safety system is estimated by the ratio of the actual level of normalized functional utility GU n and the share of consumed resource X c , and is actually identified by the coordinate of the first Pareto point on the graph GU n = f (X c ). Table 6 shows the data on changes in the numerical values of the actual level of the normalized functional utility GU n RSS of the Russian road safety provision system.Comparing the data in Figure 16 with the known optimal value of the normalized utility GU opt = 0.618, it is easy to conclude about the current state of the Russian road safety provision system harmonicity.COVID-restrictions on mobility of all types (pedestrian, transport and general) have very seriously affected the change in all characteristics that define the process of road traffic accident formation in Russia.
Table 7 presents data describing the process of road accident rate formation in Russia in 2019 and 2020.Tables 8 and 9 present the intermediate (according to the methodology of Section 4.2 of the article) and final results of calculating the relative entropy of the Russian road safety system H n RSS RF in 2019 and 2020.On the basis of the numerical values of H n RSS RF-2019 = 0.7726 and H n RSS RF-2020 = 0.7553 we determined the appropriate Q-factor levels of the road safety system for the versions of 2019 and 2020.
Q-factor of the road safety system in the Russian Federation in 2019 g RSS RF-2019 = 2.3754, and its analogue for 2020 g RSS RF-2020 = 2.4646.
The calculations show that the changes in the Russian road safety system that occurred during COVID restrictions (2020) ensured a decrease in its relative entropy H n by 2.24%, from H n RSS 2019 = 0.7726 to H n RSS 2020 = 0.7553.The orderliness of the road safety system in Russia has improved.
The Q-factor of the Russian road safety system g RSS increased by 3.75% during COVID restrictions (2020).
Normalized general utility GU n RSS increased by 0.46% from 2019 to 2020 (GU n RSS 2019 = 0.649 and GU n RSS 2019 = 0.652).

Conclusions
During the research it was established that the Q-factor g of the Russian road safety provision system changed from g 2006 = 1.9565 to g 2020 = 2.4646.This promotes the decline of the relative entropy of the Russian road safety provision system from H n RSS 2006 = 0.8623 to H n RSS 2020 = 0.7553.The deep reason for that change was the modification of the relation between "weights" or the significance of the contribution of different elements of the cause-and-effect chain in the formation of the factual level of the road accident rate in Russia in the last 15 years.
The main conclusion is that the harmonicity of the Russian road safety provision system, assessed by the normalized functional general utility GU n , has been increased and it already exceeded the level of harmonious reference systems GU n = 0.618.
Actually, the normalized functional general utility GU n RSS of the Russian road safety provision system increased from GU n RSS 2006 = 0.615 to GU n RSS 2020 = 0.652 (by 6.0%) from 2006 to 2020.Simultaneously, the share of the normalized used resource X n declined.That allowed us to make a conclusion about a significant improvement in the balance "efficiency-quality" of the Russian road safety provision system.The COVID lockdown played a positive role in this process.The harmonicity of the Russian road safety provision system, assessed by the normalized general utility GU n RSS , increased by 0.46% from 2019 to 2020.
It is important to note that the methodology of the quantitative assessment of the system harmonicity can be applied to a variety of spheres of human activity.The road safety system is just one example.Exactly the same calculations can be carried out to assess the harmonicity of the functioning of a variety of socio-technical and socio-cultural systems.

Figure 1 .
Figure 1.Functional scheme of the road safety management process (system).
quality indicator g and s-indicators.The generalized golden ratio regime (GGR regime) directly relates to the functionality of reference management systems, presented as a macromodel of normalized general utility GUn (3).Its graphic representation as Pareto chart PC(x) (taking into account equivalency PC(x) = GU(x)) is shown in Figure2.

Figure 2 .
Figure 2. Macromodel of normalized general utility GUn in the form of Pareto chart.The GGR regime is presented in the figure as a diagonal y = (1 -x), where every point meets the condition y + x = 1.The point of intersection of the diagonal and Pareto chart (with coordinates fs and fr) is called the first Pareto point.We want to highlight its

Figure 2 .
Figure 2. Macromodel of normalized general utility GU n in the form of Pareto chart.

Mathematics 2021, 9 ,
x FOR PEER REVIEW 15 of 35 occurs.The traditional approach supposes hidden (latent) work of some global mechanism, called the golden ratio or golden mean.

Figure 5 .Figure 5 .
Figure 5. Conceptual model of complex system harmonycity formation.The research into these characteristics allows us to solve two important problems: the problem of the optimization of the system utility GU and the problem about the log-

Figure 7 .
Figure 7. Plots of functions GU(x c ) and ST(x c ).

Figure 8 .
Figure 8.The cause-and-effect model of road accident rate formation.

Figure 8 .
Figure 8.The cause-and-effect model of road accident rate formation.

. 4 .
Evaluation of the Entropy as A Characteristic of the Orderliness of the Road Safety Provision System

Figure 9 .
Figure 9. Structural harmony of the "weights" coefficients i  .

Figure 10 .
Figure 10.Casual harmony of the "weights" coefficients i  .

Figure 11 .
Figure 11.The relation of the "weights" coefficients i

Figure 12 .
Figure 12.Functional harmony.6. Quantitative Assessment of Changes in the Road Safety System Harmonicity in the Russian Federation 6.1.Assessment of the General Trend in Changes in the Road Safety System Harmonicity in Russia from 2006 to 2020

Figure 11 .
Figure 11.The relation of the "weights" coefficients ω i for the road safety provision system case.

Figure 11 .
Figure 11.The relation of the "weights" coefficients i

Figure 12 .
Figure 12.Functional harmony.6.Quantitative Assessment of Changes in the Road Safety System Harmonicity in the Russian Federation6.1.Assessment of the General Trend in Changes in the Road Safety System Harmonicity in Russia from 2006 to 2020The data source was the official statistics of the State Inspectorate for Road Traffic Safety of the Russian Federation[11].Over the past 15 years (2006-2020), the indicators of road safety in Russia have been constantly improving.A number of factors have contributed to this, such as a general increase in the transport culture of the country's population, a significant increase in the quality of the country's transport fleet and significant qualitative growth of the country's road transport infrastructure.It is also important that since 2006, a lot of attention has been paid to the issues of road safety provision by the Russian government[9,10].

Figure 12 .
Figure 12.Functional harmony.6. Quantitative Assessment of Changes in the Road Safety System Harmonicity in the Russian Federation 6.1.Assessment of the General Trend in Changes in the Road Safety System Harmonicity in Russia from 2006 to 2020

Figure 13 .
Figure 13.The actual dynamics of the "weights" coefficients of the road safety system in Russia.

Figure 14 Figure 14 .
Figure14shows a comparison of the theoretical and actual dynamics of changes in the "weight" coefficients i  of the road safety system in Russia over a 15-year period.

Figure 13 . 35 Figure 13 .
Figure 13.The actual dynamics of the "weights" coefficients of the road safety system in Russia.

Figure 14 Figure 14 .
Figure14shows a comparison of the theoretical and actual dynamics of changes in the "weight" coefficients i  of the road safety system in Russia over a 15-year period.

Figure 14 .
Figure 14.Comparison of the theoretical and actual dynamics of changes in the "weights" coefficients of the road safety system in Russia over a 15-year period (2006-2020).

Figure 15 .
Figure 15.Comparison of the theoretical and actual dynamics of causal harmony of the road safety system in Russia over a 15-year period (2006-2020).The functionality of the road safety system is estimated by the ratio of the actual level of normalized functional utility n GU and the share of consumed resource

Figure 16
Figure 16 shows the dynamics of the actual numerical values of the normalized general utility RSS n GU

Figure 15 .
Figure 15.Comparison of the theoretical and actual dynamics of causal harmony of the road safety system in Russia over a 15-year period (2006-2020).

Figure 16 35 Figure 16 .
Figure 16 shows the dynamics of the actual numerical values of the normalized general utility GU n RSS of the Russian road safety provision system from 2006 to 2020.Mathematics 2021, 9, x FOR PEER REVIEW 28 of 35

Figure 16 .
Figure 16.The actual dynamics of the values of the normalized general utility GU n .

6. 2 .
Comparison of the Road Safety System Harmonicity in Russia during the COVID Restriction Period (2020) and during the Pre-COVID Period (2019)

Table 2 .
Changes in the main road safety indicators in the Russian Federation from 2006 to 2020.

Table 3 .
Changes in relative entropy H n RSS and Q-factor g values of the road safety system in Russia.

Table 4 .
Changes in the numerical values of the "weight" coefficients ω i from 2006 to 2020.

Table 5 .
General statistics of penalties for violators of traffic rules in Russia (2016-2020).
The share of speed limit violations in the total number of

Table 6 .
Changes in the numerical values of the normalized general utility

Table 6 .
Changes in the numerical values of the normalized general utility GU n RSS of the Russian RSS from 2006 to 2020.

Table 7 .
Numerical values of the characteristics of the causal path of road traffic accident formation in Russia before (2019) and during (2020) the COVID-19 pandemic.

Table 8 .
Intermediate and final results of calculating the relative entropy H n RS SRF of the Russian road safety system in 2019.

Values of the Elements of the Three-Link Mechanism of Informational Transformation of the Road Safety System in the Russian Federation, 2019
Numerical values of the positiveQ i = ln(1/K i ) Numerical value of absolute entropy H RSS RF = 0.8487 Numerical value of relative entropy H n RSS RF = 0.7726

Table 9 .
Intermediate and final results of calculating the relative entropy H n RSS RF of the Russian road safety system in 2020.

Values of the Elements of the Three-Link Mechanism of Informational Transformation of the Road Safety System in the Russian Federation, 2020
Numerical values of the positiveQ i = ln(1/K i ) Numerical values of relative "weight" of the positive ω i (Q i ) Numerical values (ω i • ln ω i )Numerical value of absolute entropy H RSS RF = 0.8298 Numerical value of relative entropy H n RSS RF = 0.7553