Rodrick Wallace
email: wallace@nyspi.columbia.edu

Canonical instabilities of autonomous vehicle systems

Keywords: control theory, driving experience, groupoid, information theory, turbulence, vehicle density ...[T]he psychological actions of drivers make traffic different

Formal argument suggests that command, communication and control systems can remain stable in the sense of the Data Rate Theorem that mandates the minimum rate of control information required to stabilize inherently unstable 'plants', but may nonetheless, under fog-of-war demands, collapse into dysfunctional modes at variance with their fundamental mission. We apply the theory to autonomous ground vehicles under intelligent traffic control in which swarms of interacting, self-driving devices are inherently unstable as a consequence of the basic irregularity of the road network. It appears that such 'V2V/V2I' systems will experience large-scale failures analogous to the vast propagating fronts of power network blackouts, and possibly less benign, but more subtle patterns of 'psychopathology' at various scales.

Introduction

We use a variety of mathematical models to explore the dynamics of rapidacting, inherently unstable command, communication and control systems (C 3 ) that are cognitive in the sense that they must, in an appropriate 'real time', evaluate a large number possible actions and choose a small subset for implementation. Such choice decreases uncertainty, in a precise formal manner, and reduction in uncertainty implies the existence of an information source. See [START_REF] Wallace | Consciousness, crosstalk, and the mereological fallacy: an evolutionary perspective[END_REF]Wallace ( , 2015a) ) for more details. We first examine dynamics related to V2V/V2I automated ground vehicle systems from the perspective of the Data Rate Theorem (e.g., [START_REF] Nair | Feedback control under data rate constraints: an overview[END_REF], extend the argument to more general phase transition analogs for C 3 structures, study instability as onset of 'turbulent' modes, and then use an information bottleneck method to derive similar results, developing statistical tools useful at different scales and levels of organization. V2V/V2I atuomated ground vehicle systems operate along geodesics in a densely convoluted 'map quotient space' -more fully described later -that is in contrast to the problem of air traffic control, where locally stable vehicle paths are seen as thick braid geodesics in a simpler Euclidean quotient space [START_REF] Hu | Hybrid geodesics as optimal solutions to the collision-free motion planning problem[END_REF]. Such geodesics are generalizations of the streamline characteristics of hydrodynamic flow (e.g., [START_REF] Landau | Fluid Mechanics[END_REF].

Data Rate Theorem

Unlike aircraft, which can be constructed to be inherently stable in linear flight by placing the aerodynamic center of pressure sufficiently behind the mechanical center of gravity, the complex nature of road geometry and density of dense vehicular traffic ensures that V2V/V2I systems will be inherently unstable, requiring constant input of control information to prevent, at the least, traffic jams and tie-ups. For a linear 'plant' described by an n-dimensional state vector x t at time t, the dynamics are given by the vector equation x t+1 = Ax t + Stochastic and Control terms

(1

)
where A is an appropriate n × n system matrix.

Let H be the control information rate needed to stabilize a rapidly-responding but inherently unstable C 3 system. By the Data Rate Theorem (DRT), H must be greater than the rate at which 'topological information' is generated. For the linear system of Eq.( 1), that rate is

α ≡ log[| det(A u )|]
where det is the determinant and A u is the decoupled component of the matrix A having eigenvalues ≥ 1 [START_REF] Nair | Feedback control under data rate constraints: an overview[END_REF].

For a fixed road network, it is evident that the central topological variate is the vehicle density that defines system dynamics under the 'fundamental diagram' relating vehicular flow per unit time to vehicle density per unit length. See figure (1), showing, for a Rome street, the number of vehicles/hour as a function of vehicles/mile. Behavior shifts from regular to 'uncontrolled' at about 40 v/mi.

Under such circumstances, vehicle density ρ is the only parameter defining the rate of topological information generation, and we can write the stability relation of the DRT as

H > α = f (ρ)α 0 (2)
Figure 1: Vehicles per hour as a function of vehicle density per mile for a street in Rome [START_REF] Blandin | A general phase transition model for vehicular traffic[END_REF]. Both streamline geodesic flow and the phase transition to 'crystallized' turbulent flow at critical traffic density are evident. Some of the states may be 'supercooled', i.e., delayed 'crystallization' in spite of high traffic density. 'Fine structure' can be expected within both geodesic and turbulent modes.

where α 0 is a road network constant and f (ρ) is a positive, monotonically increasing function. The Mathematical Appendix uses a classic Black-Scholes argument to approximate the 'cost' of H as a function of the 'investment' ρ as a linear function. Similarly, we can, again in first order, approximate f (ρ) as a linear function, giving the limit condition for stability as

κ 1 ρ + κ 2 κ 2 ρ + κ 4 > α 0 (3) For ρ = 0 the stability condition is κ 2 /κ 4 > α 0 . At large ρ this becomes κ 1 /κ 2 > α 0 . If κ 2 /κ 4
κ 1 /κ 2 , the stability condition may be violated at high traffic density, and instability becomes manifest, as at higher densities in figure (1).

A more complex model, based explicitly on the cognitive nature of control systems, provides deeper insight.

Phase transition

Again, the essence of cognition is active choice. As described above, a cognitive system, confronted with uncertainty, must choose a small set of actual responses from a larger set of those available to it. Such choice reduces uncertainty, and the reduction of uncertainty implies the existence of an information source. The argument can be made quite formally [START_REF] Wallace | Consciousness, crosstalk, and the mereological fallacy: an evolutionary perspective[END_REF](Wallace , 2015a)).

Given an information source associated with an inherently unstable, rapidfire cognitive C 3 system -said to be 'dual' to it -an equivalence class algebra emerges by choosing different system origin states a 0 and defining the equivalence of two subsequent states at times m, n > 0, written as a m , a n , by the existence of high-probability meaningful paths connecting them to the same origin point. Disjoint partition by equivalence class, essentially similar to orbit equivalence classes in dynamical systems, defines a symmetry groupoid associated with the cognitive process. Groupoids represent generalizations of the group concept in which there is not necessarily a product defined for each possible element pair [START_REF] Weinstein | Groupoids: unifying internal and external symmetry[END_REF]. The simplest example would be a disjoint union of groups.

The equivalence classes define a set of cognitive dual information sources available to the inherently unstable C 3 system, creating a large groupoid, with each orbit corresponding to a transitive groupoid whose disjoint union is the full groupoid. Each subgroupoid is associated with its own dual information source, and larger groupoids will have richer dual information sources than smaller.

Let X Gi be the C 3 dual information source associated with the groupoid element G i , and let Y be the information source associated with incoming environmental information, in a large sense See [START_REF] Wallace | Consciousness, crosstalk, and the mereological fallacy: an evolutionary perspective[END_REF]Wallace ( , 2015a) ) for details of how environmental regularities imply the existence of an environmental information source.

We construct a Morse Function [START_REF] Pettini | Geometry and Topology in Hamiltonian Systems[END_REF] as follows.

Let H(X Gi , Y ) ≡ H Gi be the joint uncertainty of the two information sources. Define a Boltzmann-like pseudoprobability as

P [H Gi ] = exp[-H Gi /κH] j exp[-H Gj /κH] (4)
where κ is an appropriate constant depending on the particular system and its linkages to control signals, and the sum is over the different possible cognitive modes of the full system. A 'free energy' Morse Function F can be defined as

exp[-F/κH] ≡ j exp[-H Gj /κH] (5)
Given the inherent groupoid structure, it is possible to apply an extension of Landau's picture of phase transition [START_REF] Pettini | Geometry and Topology in Hamiltonian Systems[END_REF]. In Landau's formulation of spontaneous symmetry breaking, phase transitions driven by temperature changes occur as alteration of system symmetry, with higher energies at higher temperatures being more symmetric. The shift between symmetries is highly punctuated in the temperature index, here the minimum necessary control information rate H under the Data Rate Theorem for unstable control systems. Typically, such arguments involve only a very limited number of possible phases.

In this context, Birkhoff's (1960 p.146) perspective on the central role of groups in fluid mechanics is of considerable interest:

[Group symmetry] underlies the entire theories of dimensional analysis and modeling. In the form of 'inspectional analysis' it greatly generalizes these theories... [R]ecognition of groups... often makes possible reductions in the number of independent variables involved in partial differential equations... [E]ven after the number of independent variables is reduced to one... the resulting system of ordinary differential equations can often be integrated most easily by the use of group-theoretic considerations.

We will argue here that, for 'cognitive fluids' like vehicle traffic flows, groupoid generalizations of group theory become central.

Decline in the richness of control information H, or in the ability of that information to influence the system, characterized by κ, can lead to punctuated decline in the complexity of cognitive process possible within the C 3 system, driving it into a ground state collapse that may not be actual 'instability' but rather a kind of dead zone in which, using the armed drone example, 'all possible targets are enemies'. This condition represents a dysfunctionally simple cognitive groupoid structure roughly akin to certain individual human psychopathologies (Wallace 2015a).

Below, we will argue that, for large-scale autonomous vehicle/intelligent infrastructure systems, the ground state dead zone involves massive, propagating tie-ups that far more resemble power network blackouts than traditional traffic jams. The essential feature is the role of vehicle road density ρ. Most of the topology of the inherently unstable vehicle/road system will be 'factored out' via the construction of geodesics in a topological quotient space, so that ρ remains the only possible index of the rate of topological information generation for the DRT. Thus, in Eqs(4) and ( 5), H is replaced by the ratio H/f (ρ), where f is a dimensionless monotonic increasing positive function.

For a fixed H, increasing ρ is then equivalent to lowering the 'temperature', and the system passes from high symmetry 'free flow' to different forms of 'crystalline' structure -broken symmetries representing platoons, shock fronts, traffic jams, and the like.

Again, making an exactly-solvable Black-Scholes approximation, the Mathematical Appendix shows that the 'cost' of the control information H can, in first order, be expressed in terms of ρ as

H ≈ κ 1 ρ + κ 2 (6)
where κ i ≥ 0. Again, in first order, taking f (ρ) ≈ κ 3 ρ + κ 4 > 0, we obtain an effective 'temperature' as

κH/f (ρ) ≈ κ 1 ρ + κ 2 κ 3 ρ + κ 4 (7)
with limits κ 2 /κ 4 , ρ → 0 and κ 1 /κ 3 , ρ → ∞. Again, assuming κ 2 /κ 4 κ 1 /κ 3 , increase in traffic density will quickly bring the effective 'temperature' below critical values, triggering collapse into a dysfunctional ground state, even though, in this case, the system might remain formally 'stable' under the Data Rate Theorem at the limit κ 1 /κ 3 > 0. Thus we begin to explore variations in stability beyond the DRT itself.

The underlying dynamic can be treated in finer detail by viewing the initial phase transition as the first-order onset of a kind of 'turbulence', a transition from free flow to 'flock' structures like those studied in 'active matter' physics. Indeed, the traffic engineering perspective is quite precisely the inverse of mainstream active matter studies, which Ramaswamy (2010) describes as follows:

It is natural for a condensed matter physicist to regard a coherently moving flock of birds, beasts, or bacteria as an orientationally ordered phase of living matter. ...[M]odels showed a nonequilibrium phase transition from a disordered state to a flock with long-range order... in the particle velocities as the noise strength was decreased or the concentration of particles was raised.

In traffic engineering, the appearance of such 'long range order' is the first stage of a traffic jam (e.g., [START_REF] Kerner | Phase transitions in traffic flow on multilane roads[END_REF], a relation made explicit by Helbing (2001 Section VI) in his comprehensive review of traffic and related self-driven many-particle systems.

While flocking and schooling have obvious survival value against predation for animals in three-dimensional venues, long-range order -aggregation -among blood cells flowing along arteries is a blood clot and can be rapidly fatal.

Turbulence

The 'free energy' function F in equation ( 5) can be used to explore dynamics within a particular system phase defined by the associated groupoid.

Given a vector of system parameters K, in standard manner it is possible to define an 'entropy' from F as the Legendre transform

S ≡ F (K) -K • ∇ K F (8)
and a nonequilibrium Onsager stochastic differential equation for dynamics in terms of the gradients in S (de Groot and Mazur 1984), which can be written in one dimension as

dK t = -[µ∂S/∂K t ]dt + σK t dB t ( 9
)
where µ is a diffusion coefficient. The last term represents a macroscopic volatility -proportional to the parameter K -in which dB t is a noise term that may not be white, i.e., the quadratic variation [B t , B t ] may not be proportional to t. While details will depend on the particular circumstances, such systems are subject to a distressingly rich spectrum of possible instabilities (e.g., Khasminskii 2012). The full set of equations would involve properly indexed sums across the parameters making up the vector K.

Allowing discontinuous Levy-like stochastic jumps, equation ( 9) can be solved explicitly, provided K can be factored out and the equation set reexpressed as

dK t = K t-dY t (10)
where Y t is a stochastic process, and t-indicates left-continuous. Taking ∆Y t = Y t -Y t-as representing the jump process, the solution is in terms of the Doleans-Dade exponential [START_REF] Protter | Stochastic Integration and Differential Equations[END_REF])

K t = exp(Y t - 1 2 [Y t , Y t ] C t )Π s≤t (1 + ∆Y s ) exp(-∆Y s ) (11)
where [Y t , Y t ] C t is the path-by-path continuous part of the quadratic variation of Y t . This can be expressed

[Y t , Y t ] C t = [Y t , Y t ] - 0≤s≤t (∆Y s ) 2 (12)
The product term in equation ( 11), with jump processes having nonzero ∆, converges.

Invoking the mean value theorem in equation ( 11), if, heuristically,

dY t < 1/2d[Y, Y ] C
t , then the expectation of K, E(K), converges to zero. Otherwise, small perturbations will grow exponentially in expectation. This is an essential component of modern theories of turbulence (e.g., [START_REF] Ruelle | Turbulent dynamical systems[END_REF]) and will prove central in understanding traffic flow instabilities. The phenomenon is roughly analogous to models of 'first order' phase transition, for example the sudden crystallization of a supercooled liquid when its container is tapped or a seed crystal is introduced. Higher order phase transitions involve discontinuities in higher order derivatives of system parameters and phase transitions may involve discontinuities at multiple 'derivative' scales.

A simple example. If a system following equation ( 10) has been initially placed in a characteristic eigenmode -e.g., the smooth part of a 'fundamental diagram' flow on some traffic network -then the dynamic equation for deviations in some parameter K(t) from that mode can be written, in first order, as

dK t ≈ aK t dt + σK t dW t ( 13
)
where dW t represents white noise having uniform spectrum. Then, using the Ito chain rule,

d log[K] t ≈ (a -σ 2 /2)dt + σdW t ( 14
)
The expectation is then

E[K] t ∝ exp[(a -σ 2 /2)t] (15) so that, if a < σ 2 /2, E[K] → 0.
σ 2 then -quite counterintuitively as described in Wallace ( 2016) -is a kind of control information in the sense of the Data Rate Theorem that serves to stabilize system dynamics. For an inherently unstable traffic flow system at low traffic density σ 2 = ωH, where H is the minimum necessary control information to keep the vehicle on the road, and ω is the 'cost' of translation. Thus H represents the degree of independent control a driver or autonomous vehicle computer can exercise within local road flow, e.g., changing lanes, accelerating, slowing down, changing headway, predicting other vehicle maneuvers, detecting emerging bottlenecks, swerving, taking a different route, and so on.

Anticipating the argument, for traffic models below, we will argue that σ 2 must again be replaced by H(ρ)/g(ρ) where g(ρ) is positive monontonic increasing. Then, antiparalleling the arguments of Belletti et al. (Bellitti et al. 2015, Section 2.3), for this simple example a 'traffic Froude number' (TFN) F that defines regimes of free and turbulent flow can be defined as

F ≡ 1 -[a - 1 2 (H(ρ)/g(ρ))] ( 16
)
where H(ρ)/g(ρ) is clearly a variant of Eq.( 7) and a is 2α 0 in Eq.( 3). When F > 1, the system is in 'laminar' free-flow, and becomes 'turbulent' when F < 1.

A more precise characterization, from this perspective, is that H(ρ)/g(ρ) represents a kind of viscosity index so that F is more akin to a Reynolds number than to a classical Froude number.

A difference between our approach and that of Bellitti et al. lies in the central object-of-interest. They invoke a hydrodynamic perspective involving the 'flow' of individual vehicles in a channel that finds 'instability' to be associated with unconstrained travel speed. The focus here is on the stability of geodesics in a complex topological quotient space M 2n /W (r) that will be more precisely defined below. This is, in a sense, the inverse of their problem.

As argued at the end of Section 3, raising ρ is equivalent to 'freezing' the system from 'liquid flow' to 'crystallized' broken symmetries -platoons, shock fronts, jams, and myriad other 'snowflake' fine structures.

[Y j t , Y j t ] C t may be further parameterized and, using the methods of Dzhaparidze and Spreij (1994), for colored noise, can be estimated from a time series data periodogram, as described in [START_REF] Wallace | Subtle noise structures as control signals in high-order biocognition[END_REF]. Increasing complexity in spectral structure is another marker of turbulence onset [START_REF] Ruelle | Turbulent dynamical systems[END_REF].

More specifically, for a stochastic process X t and a finite stopping time T , for each real number λ, the periodogram of X at T is defined as

I T (X; λ) ≡ | T 0 exp[iλt]dX t | 2 (17)
Take as a real random variate that has a density function r assumed to be symmetric around zero and examine, for any positive real number L,

E [I T (X; L ] = +∞ -∞ I T (X; Ls)r(s)ds (18)
Some work shows that, for L → ∞,

E [I T (X; L )] → [X T , X T ] (19)
Thus the quadratic variation can be estimated from observed time series data, as routinely done in financial engineering. Presumably, both high and low frequency limits could be explored in this method, parallel to what was done in [START_REF] Belletti | Prediction of traffic convective instability with spectral analysis of the Aw-Rascle-Zhang model[END_REF].

Information bottleneck

Another modeling approach is via the information bottleneck method of [START_REF] Tishby | The information bottleneck method[END_REF]. The essential idea is that the control information needed to stabilize an inherently unstable system, which we write as H, can be used to define an average distortion measure in a rate distortion calculation. This involves an iterated application of the Rate Distortion Theorem (Cover and Thomas 2006) to a control system in which a series of 'orders' y n = y 1 , ..., y n , having probability p(y n ), is sent through and the outcomes monitored as ŷn = ŷ1 , ..., ŷn . The distortion measure is now the minimum necessary control information for system stability, H(y n , ŷn ). We can thus define an average 'distortion' Ĥ as

Ĥ ≡ y n p(y n )H(y n , ŷn ) ≥ 0 ( 20
)
It then is possible to define a new, iterated, Rate Distortion Function R( Ĥ).

For simplicity, we take R to be a Gaussian RDF in Ĥ,

R( Ĥ) = 1/2 log[σ 2 / Ĥ] Ĥ < σ 2 R( Ĥ) = 0 Ĥ ≥ σ 2 (21)
Again, following Feynman (2000), information must be recognized as a form of free energy and a Rate Distortion Function can be used to define an 'entropy' as the Legendre transform

S = R( Ĥ) -ĤdR/d Ĥ (22)
Taking Onsager's nonequilibrium thermodynamics perspective, the dynamics can, in first order, be characterized in terms of the gradients of S, and we invoke an extended analog using the stochastic differential equation

d Ĥt = [-µdS/d Ĥt -F (ρ)]dt + β Ĥt dW t = [ µ 2 Ĥt -F (ρ)]dt + β Ĥt dW t ( 23
)
where dW t is standard white noise, and F (ρ) is a function of traffic density ρ, the only possible determinant of the rate of generation of system 'topological information', given the extreme topological factoring associated with travel along a network. β represents the magnitude of a 'volatility' noise term independent of σ 2 in the definition of R: higher H, higher stochastic jitter.

Applying the Ito chain rule to the expectation of Ĥ2 t , it becomes possible to explore the second moment stability of the system (Khashminskii 2012). A simple calculation finds that the expectation for Ĥ2 cannot be a real number unless the discriminant of a quadratic equation is nonnegative, giving the necessary condition

F (ρ) ≥ β √ µ (24)
We force 'closure' to the model by taking F (ρ) as given by Eq.( 3), so that, again,

κ 1 ρ + κ 2 κ 3 ρ + κ 4 ≥ β √ µ ≡ α 0 (25)
with similar restrictions on the constants κ i for stability.

Other channel forms will have analogous limits on traffic density as a consequence of the convexity of the RDF. The interested reader might carry the calculation through for the 'real' channel, having R( Ĥ) = σ 2 / Ĥ.

6 Reconsidering V2V/V2I systems [START_REF] Kerner | Empirical features of spontaneous and induced traffic breakdowns in free flow at highway bottlenecks[END_REF] explicitly apply insights from statistical physics to traffic flow, writing In many equilibrium... and dissipative metastable systems of natural science... there can be a spontaneous phase transition from one metastable phase to another metastable phase of a system. Such spontaneous phase transition occurs when a nucleus for the transition appears randomly in an initial metastable phase of the system: The growth of the nucleus leads to the phase transition. The nucleus can be a fluctuation within the initial system phase whose amplitude is equal or larger than an amplitude of a critical nucleus required for spontaneous phase transition. Nuclei for such spontaneous phase transitions can be observed in empirical and experimental studies of many equilibrium and dissipative metastable systems... There can also be another source for the occurrence of a nucleus, rather than fluctuations: A nucleus can be induced by an external disturbance applied to the initial phase. In this case, the phase transition is called an induced phase transition... A Data Rate Theorem approach to stability and flow of autonomous vehicle/traffic control systems, via spontaneous symmetry breaking in cognitive groupoids, generalizes and extends these insights, implying a far more complex picture of control requirements for inherently unstable systems than is suggested by the Theorem itself, or by 'physics' models of phase transition. That is, 'higher order' instabilities can appear. Such systems can require inordinate levels of control information. Here, we find that C 3 systems may remain 'stable' in the strict sense of the DRT, but can collapse into a ground state analogous to certain psychopathologies, or, following the arguments above, into even more complicated pathological dynamics. In biological circumstances, such failures can be associated with the onset of senescence [START_REF] Wallace | Extending Swerdlow's hypothesis: statistical models of mitochondrial deterioration and aging[END_REF](Wallace , 2015b)). Apparently, rapidly responding, and hence almost certainly inherently unstable, C 3 systems can display recognizable analogs to senility under fog-of-war demands.

Using these ideas, it becomes possible to formally represent the interaction of cognitive ground state collapse in autonomous vehicle/intelligent road systems with critical transitions in traffic flow.

Recall that, defining 'stability' as the ability to return, after perturbation, to the streamline geodesic trajectory of the embedding, topologically complex, road network, it is clear that individual autonomous vehicles are inherently unstable and require a constant flow of control information for safe operation, unlike aircraft that can, in fact, be made inherently stable by placing the center of pressure well behind the center of gravity. There is no such configuration possible for ground-based vehicles following sinuous road geometries, particularly in heavy traffic.

The argument can be made more precise using the approach of [START_REF] Hu | Hybrid geodesics as optimal solutions to the collision-free motion planning problem[END_REF] who show that, in the context of air traffic control, finding collision-free maneuvers for multiple agents on a Euclidean plane surface R 2 is the same as finding the shortest geodesic in a particular manifold with nonsmooth boundary. Given n vehicles, the geodesic is calculated for the quotient space R 2n /W (r), where W (r) is defined by the requirement that no vehicles are closer together than some critical Euclidean distance r. For autonomous ground vehicles, R 2 must be replaced by a far more topologically complex roadmap space M 2 subject to traffic jams and other 'snowflake' condensation geometries in real time. Geodesics for n vehicles are then in a quotient space M 2n /W (r) whose dynamics are subject to phase transitions in vehicle density ρ [START_REF] Kerner | Phase transitions in traffic flow on multilane roads[END_REF][START_REF] Jin | Spontaneous phase transition from free flow to synchronized flow in traffic on a single-lane highway[END_REF] that represent cognitive groupoid symmetry breaking. Recall figure 1. Again, the vertical axis shows the number of vehicles per hour, the horizontal, the density of vehicles per mile. The streamline geodesic flow, and deviations from it at critical vehicle density, are evident. Some of the phases may be 'supercooled' -fast-flowing 'liquid' at higher-than-critical densities. Additional 'fine structure' should be expected within both geodesic and turbulent modes.

Again, given the factoring out of most of the topological structure by the construction of geodesics in the quotient space M 2n /W (r), the only parameter available to represent the rate of generation of topological information in the inherently unstable traffic flow system is the vehicle density ρ.

Classic traffic flow models based on extensions of hydrodynamic perspectives involving hyperbolic partial differential equations (HPDE's) can be analogously factored using the methods of characteristic curves and Riemann invariants -streamlines (e.g., [START_REF] Landau | Fluid Mechanics[END_REF]. Along characteristic curves, HPDE's are projected down to ordinary differential equations (ODE's) that are usually far easier to solve. The ODE solution or solutions can then be projected upward as solutions to the HPDE's. Here, reduction involves expressing complex dynamics in terms of relatively simple stochastic differential equations and their stability properties. Those stability properties, marking the onset of 'turbulence', will be of central interest.

Taking a somewhat larger view, cognitive phase transitions in autonomous vehicle systems, in particular ground state collapse to some equivalent of 'all possible targets are enemies', should become synergistic with traffic flow phase transitions to produce truly monumental traffic jams, and it is possible to model this phenomenon, to first order, in terms of spontaneous symmetry breaking on groupoids.

Consider a random network of roads between nodal points -intersections. If the average probability of passage falls below a critical value, the Erdos/Renyi 'giant component' that connects across the full network breaks into a set of disjoint connected equivalence class subcomponents, with 'bottlenecks' at which traffic jams occur marking corridors between them. [START_REF] Li | Percolation transition in dynamical traffic network with evolving critical bottlenecks[END_REF], in fact, explicitly apply a similar percolation model to explain this effect for road congestion in a district of Beijing. The underlying road network is shown in figure 2, and in figure 3 a cross section taken during rush hour showing disjoint sections when regions with average velocity below 40% of observed maximum for the road link have been removed.

Such equivalence classes define a groupoid [START_REF] Weinstein | Groupoids: unifying internal and external symmetry[END_REF]. Above, we defined the cognitive groupoid to be associated with a C 3 structure, here a system of autonomous vehicles linked together in a V2V 'swarm intelligence' embedded in a larger vehicle to infrastructure (V2I) traffic management system. Individ- ual vehicle spacings, speed, acceleration, lane-change, and so on are determined by this encompassing distributed cognitive machine that attempts to optimize traffic flow and safety. The associated individual groupoids are the basic transitive groupoids that build a larger composite groupoid. Thus, under declining probability of passage, related to traffic congestion and viewed here as a temperature analog, this 'vehicle/road' groupoid undergoes a symmetry-breaking transition into a combined cognitive ground state collapse and traffic jam mode -essentially a transition from 'laminar' geodesic to 'turbulent' or 'crystallized' flow. Autonomous vehicle systems that become senile under fog-of-war demands will likely trigger traffic jams that are far different from those associated with human-controlled vehicles. There is no reason to believe that such differences will be benign.

A more precise-seeming model can be built using the turbulence limit of equations ( 9) and ( 10), under conditions of white noise, so that [Y, Y ] C ∝ t, where t is the time.

Following figure 1, suppose the system has been placed in a geodesic for the map quotient space M 2n /W (r), and is free-flowing 'laminar' at some vehicle density ρ, but may be inherently unstable. As described above, for the simplest kinds of stochastic differential equation models involving white noise -equation ( 14) -heuristically,

Y t → at (26) [Y t , Y t ] → (H(ρ)/g(ρ))t
in the exponential of equation ( 11), where α > 0 and H represents the degree of 'driver' control information necessary to keep a vehicle on the road according to the DRT. As argued, this may include lane change, speed variation, spacing change, and other maneuvers. Again, in first order, if a -1/2(H(ρ)/g(ρ)) < 0 then any perturbation K dies out in expectation -the driver/autonomous vehicle can exercise sufficient initiative to damp out occasional glitches -and the system has sufficient symmetry so that it can return to streamline geodesic flow.

If ρ increases beyond a critical value, then too limited spacing constrains the possibility of cognitive vehicle initiative, the 'temperature' falls below criticality, and perturbations grow exponentially in time so that instability causes 'turbulent' crystal formation -'fine structure' traffic jams of one form or another. As argued above, this can happen even though the system remains fully 'stable' under the Data Rate Theorem: the individual vehicle control information rate H(ρ) remains sufficient so that no crash occurs. A central concept of modern theories of turbulence is the onset of exponential growth in small perturbations [START_REF] Ruelle | Turbulent dynamical systems[END_REF].

The analysis of traffic flow on a network is, conceptually, somewhat similar to characterizing the propagation of a 'traffic jam signal' via the Markov 'network dynamics' formalism of [START_REF] Wallace | Subtle noise structures as control signals in high-order biocognition[END_REF][START_REF] Gould | Spatial structures and scientific paradoxex in the AIDS pandemic[END_REF], a method that might be used to empirically identify geodesic eigenmodes of real road network systems under different conditions, as opposed to individual vehicle dynamics or flow on a single road. Abducting the approach of [START_REF] Gould | Spatial structures and scientific paradoxex in the AIDS pandemic[END_REF], the spread of a 'signal' on a particular network of interacting sitesbetween and within -is described at nonequilibrium steady state in terms of an equilibrium distribution i 'per unit area' A i of a Markov process, where A scales with the different 'size' of each node, taken as distinguishable by a scale variable A (for example number of entering streets or average total traffic flow) as well as by its 'position' i or the associated probability-of-contact matrix (POCM). The POCM is normalized to a stochastic matrix Q having unit row sums, and the vector calculated as = Q There is a vector set of dimensionless network flows X i t , i = 1, ..., n at time t. These are each determined by some relation

X i t = f (t, i /A i ) (27)
Here, i is the index of the node of interest, X i t is the corresponding dimensionless scaled i-th signal, t the time, and f an appropriate function. Again, i is defined by the relation = Q for a stochastic matrix Q, calculated as the network probability-of-contact matrix between regions, normalized to unit row sums. Using Q, we have broken out the underlying network topology, a fixed between-and-within travel configuration weighted by usage that is assumed to change relatively slowly on the timescale of observation compared to the time needed to approach the nonequilibrium steady state distribution.

Since the X are expressed in dimensionless form, f, t, and A must be rewritten as dimensionless as well giving, for the monotonic increasing (or thresholdtriggered) function

F X i τ = F [τ, i A i × A τ ] (28)
where A τ is the value of a 'characteristic area' variate that represents the spread of the perturbation signal -evolving into a traffic jam under worst-case conditions -at (dimensionless) characteristic time τ = t/T 0 . F may be quite complicated, including dimensionless 'structural' variates for each individual geographic node i. The idea is that the characteristic 'area' A τ grows according to a stochastic process, even though F may be a deterministic mixmaster driven by systematic local probability-of-contact or flow patterns. Then the appropriate model for A τ of a spreading traffic jam becomes something like equations ( 9) or (10), with K replaced by A and t by τ . Thus, for the network, the signal Y τ must again have a 'noise'/vehicle density threshold condition like Eq.( 16) for large-scale propagation of a traffic jam across the full network -something that would look very similar to the spread of a power blackout. [START_REF] Zhang | Realistic, efficient and secure geographic routing in vehicular networks[END_REF] uses a similar Markov method to examine taxicab GPS data for transit within and between 12 empirically-identified 'hot zones' in Shanghai, determining the POCM and its equilibrium distribution. This approach is something in the spirit of a long line of work summarized by [START_REF] Cassidy | Macroscopic fundamental diagrams for freeway networks[END_REF] that attempts to extend the idea of a fundamental diagram for a single road to a full transport network. As they put it, Macroscopic fundamental diagrams (MFDs)... relate the total time spent to the total distance traveled... It is proposed that these macrolevel relations should be observed if the data come from periods when all lanes on all links throughout the network are in either the congested or the uncontested regime... Following our arguments here, such conditions might apply when A τ → 0, or when it encompasses the entire network domain. Indeed, figure 3 suggests why MFDs cannot be constructed in general: congested and free flowing sections of traffic networks will often, and perhaps usually, coexist in an essentially random manner depending on local traffic densities. Figure 4, adapted from [START_REF] Geroliminis | Properties of a well-defined macroscopic fundamental diagram for urban traffic[END_REF], shows the limitations of the MFD approach. It examines the flow, in vehicles/5min intervals, vs. percent occupancy over a three day period for the Minnesota Twin Cities freeway network that connects St. Paul and Minneapolis. See figure 1 of their paper for details of the road and sensor spacing. Evidently, while the unconstrained region of occupancy permits characterization of a geodesic mode, both strong hysteresis and phase transition effects are evident after about 8% occupancy, analogous to the 'nucleation' dyanamics of figure 1 at high traffic density. Again, as in figure 1 'fine structure' should be expected within both geodesic and turbulent modes, depending on local parameters. Daganzo et al. (2010) further find that MFD flow, when it can be characterized at all, will become unstable if the average network traffic density is sufficiently high. They find that, for certain network configurations, the stable congested state ...is one of complete gridlock with zero flow. It is therefore important to ensure that in real-world applications that a network's [traffic] density never be allowed to approach this critical value. [START_REF] Daqing | Spatial correlation analysis of cascading failures: congestions and blackouts[END_REF] examine the dynamic spread of traffic congestion on the Beijing central road network. They characterize the failure of a road segment to be a traffic velocity less than 20km/hr and use observational data to define a spatial correlation length in terms of the Euclidean distances between failed nodes. Our equivalent might be something like √ A τ . Adapting their results, figure 5 shows the daily pattern of the correlation length of cascading traffic jams over a 9 day period. The two commuting maxima are evident, and greatest correlation lengths reach the diameter of the main part of the city. Even at rush hour, no MFD can be defined, as, according to figure 3, the network will be a dynamic patchwork of free and congested components.

A next step would be to allow ρ vary in space and time, i.e., to parameterize the model using the moments of vehicle density.

Figure 6, adapted from Rand (1979, figure 6.4), provides a disturbing counterexample to these careful empirical and theoretical results on network traffic flow, one with unfortunate results. Summarizing observations carried out by the Rand Fire Project, it represents a repeated sampling of 'travel time vs. distance' for the full Trenton NJ road network in 1975 under varying conditions of time-of-day, day-of-week, weather, and so on, by fire companies responding to calls for service. This was an attempt to create a Macroscopic Fundamental Diagram in the sense used above, but without any reference at all to traffic density.

Indeed, fire service responses are a traffic flow 'best case' as fire units are permitted to bypass one-way restrictions, traffic lights, and so on, and usually able to surmount even the worst weather conditions. In spite of best-case circumstances, the scatterplot evidently samples whole-network turbulent flow, not unlike that to the right of the local geodesic in figures 1 and 4, part of a single street and a highway network, respectively, and consistent with the assertions of [START_REF] Cassidy | Macroscopic fundamental diagrams for freeway networks[END_REF] that MFD relations can only be defined under very restrictive conditions, i.e., either complete free flow or full network congestion.

The Rand Fire Project, when confronted with intractable whole-network traffic turbulence, simply collapsed the data onto a 'square root-linear' relation, as indicated on the figure. The computer models resulting from this gross oversimplification were used to determine fire service deployment strategies for high fire incidence, overcrowded neighborhoods in a number of US cities, with literally devastating results and consequent massive impacts on public health and public order. [START_REF] Wallace | A Plague on Your Houses[END_REF], produced under an Investigator Award in Health Policy Research from the Robert Wood Johnson Foundation, documents the New York City case history. The Rand models are still in use by the New York City Fire Department, for political purposes outlined in that analysis.

7 Discussion and conclusions [START_REF] Ruelle | Turbulent dynamical systems[END_REF], in his elegant keynote address on turbulent dynamics, raises a red flag for any traffic flow studies: ...[A] deductive theory of developed turbulence does not exist, and a mathematical basis for the important theoretical literature on the subject is still lacking... A purely deductive analysis starting with the Navier-Stokes equation... does not appear feasible... and might be inappropriate because of the approximate nature of the... equation.

Or, as the mathematician Garrett Birkoff (1960 p.5) put it, "...[V]ery few of the deductions of rational hydrodynamics can be established rigorously".

Similar problems afflict the exactly solvable but highly approximate Black-Scholes models of financial engineering, and institutions that rely heavily on them have gone bankrupt in the face of market turbulence [START_REF] Wallace | An Ecosystem Approach to Economic Stabilization: Escaping the neoliberal wilderness[END_REF].

Turbulence in traffic flow does not represent simple drift from steady linear or even parallel travel trajectories. Traffic turbulence involves the exponential (1979). Relation between fire company travel time and response distance for the full Trenton, NJ road network, 1975. The Rand Fire Project collapsed evident large-scale traffic turbulence into a simple 'square root-linear' model used to design fire service deployment policies in high fire incidence, high population density neighborhoods of many US cities, including the infamous South Bronx. The impacts were literally devastating [START_REF] Wallace | A Plague on Your Houses[END_REF].

amplification of small perturbations into large-scale deviations from complicated streamline geodesics in a topologically complex map quotient space. This is the mechanism of groupoid 'symmetry breaking' by which the system undergoes a phase transition from 'liquid' geodesic flow to 'crystalline' phases of shock fronts, platoons, and outright jams. Under such circumstances, cognitive vehicle initiative -the stabilizing control information we have called H -serves as a mechanism for returning to local geodesic flow. Inhibition of cognitive initiative occurs when vehicle density exceeds a critical limit, triggering complex dynamic condensation patterns and, for autonomous vehicle systems, perhaps even more disruptive behaviors.

It is, then, not enough to envision atomistic autonomous ground vehicles as having only local dynamics in an embedding traffic stream, as seems the current American and European practice. Traffic light strategies, a frequentlyshifting road map space, and the dynamic composition of the traffic stream, all create the synergistic context in which individual vehicles operate and which constitutes the individual 'driving experience'. It is necessary to understand the dynamics of that full system, not simply the behavior of a vehicle atom within it, and the properties of that system will be both overtly and subtly emergent, as will, we assert, the responses of cognitive vehicles enmeshed in context, whether controlled by humans or machines.

One inference from this analysis is that failure modes afflicting large-scale V2V/V2I systems are likely to be more akin to power blackouts than to traffic jams as we know them, and the description by [START_REF] Kinney | Modeling cascading failures in the North American power grid[END_REF] is of interest: Today the North American power grid is one of the most complex and interconnected systems of our time, and about one half of all domestic generation is sold over ever-increasing distances on the wholesale market before it is delivered to customers... Unfortunately the same capabilities that allow power to be transferred over hundreds of miles also enable the propagation of local failures into grid-wide events... It is increasingly recognized that understanding the complex emergent behaviors of the power grid can only be understood from a systems perspective, taking advantage of the recent advances in complex network theory... [START_REF] Dobson | Where is the edge for cascading failure?: challenges and opportunities for quantifying blackout risk[END_REF] puts it as follows:

[P]robabalistic models of cascading failure and power system simulations suggest that there is a critical loading at which expected blackout size sharply increases and there is a power law in the distribution of blackout size... There are two attributes of the critical loading: 1. A sharp change in gradient of some quantity such as expected blackout size as one passes through the critical loading. 2. A power law region in probability distribution of blackout size at the critical loading. We use the terminology 'critical' because this behavior is analogous to a critical phase transition in statistical physics. correction factor, a powerful regularity assumption that makes this an exactly solvable approximate model.

The conventional Black-Scholes calculation takes ∆L/∆T ∝ L. Here, at nonequilibrium steady state, we assume ∆L/∆t = ∂H/∂t = 0, so that

- 1 2 b 2 ρ 2 ∂ 2 H/∂ρ 2 = 0 (33)
By inspection,

H = κ 1 ρ + κ 2 (34)
where the κ i are nonnegative constants.

Figure 2 :

 2 Figure 2: Adapted from Li et al. (2015). Full road network near central Beijing.

Figure 3 :

 3 Figure 3: Adapted from Li et al., (2015). Disconnected subcomponents of Beijing central road network at rush hour. Sections with average vehicle velocity less than 40% of maximum observed have been removed.

Figure 4 :

 4 Figure 4: Adapted from Geroliminis and Sun (2011). Breakdown of the macroscopic fundamental diagram for the freeway network connecting St. Paul and Minneapolis at high vehicle densities. Both nucleation and hysteresis effects are evident, showing the fine structure within the turbulent mode.

Figure 5 :

 5 Figure 5: Adapted from Daqing et al. (2014). Daily cycle of traffic jam correlation length over a 9 day period in central Beijing. The maxima cover most of the central city. Even for rush hour, however, no macroscopic fundamental diagram can be defined since the region is characterized by a patchwork of free and congested parts, as shown in figure 3.

Figure 6 :

 6 Figure 6: Adapted from figure 6.4 ofRand (1979). Relation between fire company travel time and response distance for the full Trenton, NJ road network, 1975. The Rand Fire Project collapsed evident large-scale traffic turbulence into a simple 'square root-linear' model used to design fire service deployment policies in high fire incidence, high population density neighborhoods of many US cities, including the infamous South Bronx. The impacts were literally devastating[START_REF] Wallace | A Plague on Your Houses[END_REF].

Acknowledgments

The author thanks Dr. D.N. Wallace for useful discussions.

 [START_REF] Daqing | Spatial correlation analysis of cascading failures: congestions and blackouts[END_REF], in fact, explicitly link traffic jams and power failures: Cascading failures have become major threats to network robustness due to their potential catastrophic consequences, where local perturbations can induce global propagation of failures... [that] propagate through collective interactions among system components....

[W]e find by analyzing our collected data that jams in city traffic and faults in power grid are spatially long-range correlated with correlations decaying slowly with distance. Moreover, we find in the daily traffic, that the correlation length increases dramatically and reaches maximum, when morning or evening rush hour is approaching... While clever V2V/V2I management strategies might keep rush hour traffic streams in supercooled high-flow mode beyond critical densities, such a state is notoriously unstable. More subtle patterns of autonomous vehicle 'psychopathology' may be even less benign.

Mathematical Appendix

Suppose the system is stable, in the control theory sense, and take H as the 'cost' of stability at the traffic density ρ. What is the form of H(ρ) under conditions of volatility i.e., variability in ρ proportional to it? Let

where dW t is taken as white noise.

Assuming stability, let H(ρ t , t) be the minimum needed incoming rate of control information under the Data Rate Theorem, and expand in ρ using the Ito chain rule

We define a quantity L as a Legendre transform of the rate H, by convention having the form

Since H is an information index, it is a kind of free energy in the sense of Feynman ( 2000) and L is a classic entropy measure.

Heuristically, replacing dX with ∆X in these expressions and applying Eq.( 30) gives

As in the classical Black-Scholes model (Black and Scholes 1973), the terms in g and dW t cancel out, and the effects of noise are subsumed into the Ito