abstract: Traffic congestion is one of the most notable problems arising in worldwide urban areas, importantly compromising human mobility and air quality. Current technologies to sense real-time data about cities, and its open distribution for analysis, allow the advent of new approaches for improvement and control. Here, we propose an idealized model, the Microscopic Congestion Model, based on the critical phenomena arising in complex networks, that allows to analytically predict congestion hotspots in urban environments. Results on real cities' road networks, considering, in some experiments, real-traffic data, show that the proposed model is capable of identifying susceptible junctions that might become hotspots if mobility demand increases.
author: Albert Solé-Ribalta; Sergio Gómez; Alex Arenas
title: A model to identify urban traffic congestion hotspots in complex networks

# Introduction

Urban life is characterized by a huge mobility, mainly motorized. Amidst the complex urban management problems there is a prevalent one: traffic congestion. Several approaches exist to efficiently design road networks and routing strategies, however, the establishment of collective actions, given the complex behavior of drivers, to prevent or ameliorate urban traffic congestion is still at its dawn. Usually, congestion is not homogeneously distributed around all city area but there are salient locations where congestion is settled. We call this locations congestion hotspots. These hotspots usually correspond to junctions and are problematic for the efficiency of the network as well as for the health of pedestrians and drivers. It has been shown that drivers in-queue in a traffic jam are the most affected individuals to car exhaust pollution inhalation. In addition, these hotspots are usually located in the city center, magnifying the problem. Assuming that congestion is an inevitable consequence of urban motorized areas, the challenge is to develop strategies towards a sustainable congestion regime at which delays and pollution are under control. The first step to confront congestion is the modelling and understating of the congestion phenomena.

The modeling of traffic flows it is prevalent hot topic since the late 70's when the Gipps' model appear . The Gipps' model and other car-following models have evidenced the necessity of modeling traffic flows to improve road network efficiency and also have shown how congestion severely affects the traffic flows. Since ten years ago also the complex networks' community has proposed stylized models to analyze the problem of traffic congestion in networks and design optimal topologies to avoid it. The focus of attention of the previous works was the onset of congestion, which corresponds to a critical point in a phase transition, and how it depends on the topology of the network and the routing strategies used. However, the proper analysis of the system after congestion has remained analytically slippery. It is known that when a transportation network reaches congestion, the system becomes highly non-linear, large fluctuation exists and the travel time and the amount of vehicles queued in a junction diverge. This phenomenon is equivalent to a phase transition in physics, and its modeling is challenging. Here, we propose an idealized model to predict the behaviour of transportation networks after the onset of congestion. The presented model is analytically tractable and can be iteratively solved up to convergence. To the best of our knowledge, this is the first analytical model that is able to give predictions beyond the onset of congestion. We present the model in terms of road transportation networks but it could also be applied to analyze other types of transportation networks such as: computer networks, business organizations or social networks.

# Results

To identify congestion hotspots in urban environments we propose a model based on the theory of critical (congestion) phenomena on complex networks. The model, that we call *Microscopic Congestion Model* (MCM), is a mechanistic model (yet simple) and analytically tractable. It is based on assuming that the growth of vehicles observed in each congested node of the networks is constant. This usually happens in real transportation networks at the stationary state. The assumption allows us to describe, with a set of balance equations (one for each node), the increment of vehicles in the junction queues' and the number of vehicles arriving or traversing each junction from neighboring junctions. Mathematically, the increment of the vehicles per unit time at every junction $i$ of the city, $\Delta{q}_{i}$, satisfies the following balance equation: $$\Delta{q}_{i} = g_{i} + \sigma_{i} - d_{i},
    \label{BE}$$ where $g_{i}$ is the average number of vehicles entering junction $i$ from the area surrounding $i$, $\sigma_{i}$ is the average number of vehicles that arrive to junction $i$ from the adjacent links of that junction, and $d_{i} \in [0,\tau_i]$ corresponds to the average of vehicles that actually finish in junction $i$ or traverse towards other junctions. Note that the value of $d_i$ is upper-bounded by the maximum amount of vehicles $\tau_i$ that can traverse junction $i$ in a time step. This simulates the physical constraints of the road network. A graphical explanation of the variables of the model is shown in Fig. .

The system of Eqs.  defined for every node $i$, is coupled through the incoming flux variables $\sigma_{i}$, that can be expressed as $$\sigma_{i} = \sum^{S}_{j=1} P_{ji} p_{j} d_{j},
    \label{sigma}$$ where $P_{ji}$ accounts for the routing strategy of the vehicles (probability of going from $j$ to $i$), $p_j$ stands for the probability of traversing junction $j$ but not finishing at $j$ and $S$ is the number of nodes in the network (see *Materials and Methods* for a detailed description of the MCM).

For each junction $i$, the onset of congestion is determined by $d_i=\tau_i$, meaning that the junction is behaving at its maximum capability of processing vehicles. Thus, for any flux generation ($g_i$), routing strategy ($P_{ij}$) and origin-destination probability distribution, Eqs.  can be solved using an iterative approach (see *Materials and Methods*) to predict the increase of vehicles per unit time at each junction of the network ($\Delta{q}_{i}$). The only hypothesis we use is that the system dynamics has reached a stationary state in which the growth of the queues is constant. It is worth commenting here that the MCM model considers a fixed average of new vehicles entering the system $g_i$. However, $g_i$ certainly changes during day time, with increasing values in rush hours and lower values during off-peak periods. MCM can easily consider evolving values of $g_i$ provided the time scale to reach the stationary state in the MCM (which is usually of the order of minutes in real traffic systems) is shorter than the rate of change in the evolution of $g_i$ (which is usually of the order of hours for the daily peaks).

## Validation on synthetic networks

To validate MCM we conducted experiments on several synthetic networks and with two different routing strategies: local search strategy and shortest path strategy. In both routing strategies we assume, for simplicity, that all vehicles randomly choose the starting and ending junctions of their journey uniformly within all junctions of the network. Thus, each junction generates new vehicles with the same rate $g_i = \rho$. For shortest path strategy, vehicles follow a randomly selected shortest path towards the destination. Without loss of generality we fix $\tau=1$ and analyze the performance of MCM for different values of $\rho$.

Figure  shows the accuracy on predicting the values of the order parameter $\eta=\frac{\sum \Delta q_i}{\rho S}$ and $d_i$ for shortest paths routing strategy. As in refs. , this order parameter $\eta$ corresponds to the ratio between in-transit and generated vehicles. All experiments show that the MCM achieves high accuracy in predicting the macroscopic and microscopic variables of the stylised transportation dynamics.

## Application to real scenarios

INRIX Traffic Scorecard (http://www.inrix.com/) reports the rankings of the most congested countries worldwide in 2014. US, Canada and most of the European countries are in the top 15, with averages that range from 14 to 50 hours per year wasted in congestion, with their corresponding economical and environmental negative consequences. To demonstrate that the MCM model can be applied to real scenarios to obtain real predictions, in the following we apply the MCM model to the ninth most congested cities according to the INRIX Traffic Scorecard (see Table ).

We first focus on the city of Milan, the city with largest INRIX value. To evaluate the outcome of the MCM model, we first gather data about the road network topology using Open Street Map (OSM). OSM data represents each road (or way) with an ordered list of nodes which can either be road junctions or simply changes of the direction of the road. We have obtained the required abstraction of the road network building a simplified version of the OSM data which only accounts for road junctions (nodes). Then, for each pair of adjacent junctions we have queried the real travel distance (i.e. following the road path) using the API provided by Google Maps. The resulting network corresponds to a spatial weighed directed network where the driving directions are represented and the weight of each link indicates the expected traveling time between two adjacent junctions.

Second, we build up the dynamics of the model analyzing real traffic data provided by Telecom Italia for their Big Data Challenge. The data provides, for every car entering the cordon pricing zone in Milan during November and December 2013, an encoding of the car's plate number, time and gate of entrance (a total of 9183475 records). This allows us to obtain the (hourly) average incoming and outgoing traffic flow, for each gate of the cordon taxed area.

Given the previous topology and traffic information, we generated traffic compatible with the observations, and evaluated the outcome of the MCM model. Specifically, the simulated dynamics is as follows: for each vehicle entering the Area-C we fix a randomly selected location as destination and use the shortest path route towards it. After the vehicle has arrived to its destination, it randomly chooses an exit door and travels to it also using the shortest path route. This is similar to the well-known *Home-to-Work* travel pattern (see details in *Materials and Methods*). Figures  and  show the obtained results. Figures  displays the predicted congestion hotspots on a map of Milan, panel **A** of the same figure shows a real traffic situation obtained with google maps. We see that the predicted congestion hotspots are located in the circular roads of Milan as well as on the arterial roads of the city; this agrees with the real traffic situation shown in panel **A**. Figures  shows the distribution of the mean increments each junction has to deal with. This might be a good indicator to decide about future planning actuations to improve city mobility. However, differently from what is described in , the improvement of the throughput of a single junction might not be enough to improve city mobility since this might end up with the collapse of neighbouring junctions (their incoming rate $\sigma_i$ will increase). This is situation is similar to the Braess' paradox . Figures  shows the mean increment of vehicles (in vehicles per minute) for each hour of the weekday. The figure clearly shows the morning and evening rush hours as well as the lunch time.

```latex
\begin{figure*}[!ht]\begin{tabular}{ll}
        {\bf A} & {\bf B}\\
        \includegraphics[width=0.90\columnwidth]{fig03a.png} ~~~ &~~~
        \includegraphics[width=0.90\columnwidth]{fig03b.png}
     \end{tabular}
    \caption{Congestion hotspot analysis of the city of Milan. Panel {\bf A} shows the typical situation around 9 a.m. for a week day. The image and the data has been obtained with Google Maps. Google maps displays traffic information considering historical data and real-time car velocity reported by smartphones \cite{barth2009googlemapstraffic}. Panel {\bf B} shows the prediction of the MCM model considering the real road topology obtained using Open Street Map and real traffic data provided by Telecom Italia for their Big Data Challenge. For all congestion hotspots the model has predicted, we shown its mean increment of the queue size, $\langle \Delta q_{i} \rangle$.}
    \label{fig:milano_hotspot_map}
\end{figure*}
```

For the other top nine congested cities, we do not have previous traffic information, neither about the real flux of vehicles nor about the vehicle source and destination distributions (to obtain a fair comparison between all the analysed cities we have not consider the Telecom traffic data for Milan here). Thus, for each city, we consider homogeneously distributed source and destination locations and the required road traffic to obtain an order parameter $\eta$ compatible with the congestion effects recorded by INRIX sensing of real traffic. By relating the INRIX value and $\eta$, we are assuming that there exists a relation between the fraction of global congestion and the fraction of extra time wasted reported by INRIX. The obtained results are summarized in Table , which shows that the amount of hotspots is correlated with the INRIX value. This shows evidences that the percentage increase in the average travel time to commute between to city locations is related to the number of congestion hotspots and with the excess of vehicles within the city.

```latex
\begin{threeparttable}
    \begin{tabular}{lcccc}
    City & INRIX\tnote{a} & hotspots & nodes & links\\
      \hline
      Milano & 36.2 & 108 & 6924 & 14315 \\
      London &  32.4 & 93 & 6378 & 14662 \\
      Los Angeles & 32.2 & 57 & 6799 & 19368 \\
      Brussels & 30.5 & 50 & 6645 & 15624 \\
      Antwerpen & 28.6 & 44 & 6530 & 15252 \\
      San Francisco & 27.9 & 45 & 8854 & 25530 \\
      Stuttgart & 21.9 & 34 & 8330 & 19946 \\
      Nottingham & 21.6 & 28 & 7337 & 16723 \\
      Karlsruhe & 21.3 & 19 & 4257 & 10379 \\
      \hline
    \end{tabular}
    \begin{tablenotes}
        \item[a] The INRIX index is the percentage increase in the average travel time of a commute above free-flow conditions during peak hours, e.g.\ an INRIX index of~30 indicates a 40-minute free-flow trip will take 52~minutes. Each city has been mapped to a graph with the indicated numbers of nodes and links.
        \end{tablenotes}
    \end{threeparttable}
```

# Discussion

The previous results show that the MCM (Microscopic Congestion Model) can be used to predict the local congestion before and beyond the onset of congestion of a transportation network. Up to the knowledge of the authors, this is the first analytical model that is able to give predictions beyond the onset of congestion where the system is highly non-linear, large fluctuation exists and the amount of vehicles on transit diverge with respect to time. Our model is based on assuming that the growth of vehicles observed in each congested node of the networks is constant, which allowed us to derive a set of balance equations that can accurately predict microscopic, mesoscopic and microscopic variables of the transportation network.

Traffic congestion is a common and open problem whose negative impacts range from wasted time and unpredictable travel delays to a waste of energy and an uncontrolled increase of air pollution. A first step towards the understanding and fight of congestion and its related consequences is the analytical modelling of the congestion phenomena. Here, we have shown that the MCM model is detailed enough to give real predictions considering real traffic data and topology. These results pave the way to a new generation of stilyzed physical models of traffic on networks in the congestion regime, that could be very valuable to assess and test new traffic policies on urban areas in a computer simulated scenario.

# Materials and Methods

## Microscopic Congestion Model

Let node $i$ denote a road junction, edge $a_{ij}$ the road segment between junctions $i$ and $j$, $N^{\mbox{\scriptsize in}}_{i}$ and $N^{\mbox{\scriptsize out}}_{i}$ the sets of ingoing and outgoing neighbouring junctions of junction $i$ respectively, and $S$ the number of junctions in the road network of the city. Incoming vehicles to node $i$ at each time step can be of two types: those coming from other junctions $N^{\mbox{\scriptsize in}}_{i}$ and those that start its trip with node $i$ as its first crossed junction. We consider this second type of incoming vehicles as generated in node $i$. Our Microscopic Congestion Model (MCM) describes the increment of the vehicles per unit time at every junction $i$ of the city, $\Delta{q}_{i}$, as: $$\label{oneLayerMMC}
    \Delta{q}_{i} = g_{i} + \sigma_{i} - d_{i}\,,$$ where $g_{i}$ is the average number of vehicles generated in node $i$, $\sigma_{i}$ is the average number of vehicles that arrive to node $i$ from $N^{\mbox{\scriptsize in}}_{i}$ junctions, and $d_{i} \in [0,\tau_i]$ corresponds to the average number of vehicles that actually finish in it, or traverse this junction towards neighboring nodes in $N^{\mbox{\scriptsize out}}_{i}$. Parameter $\tau_i$ represents the maximum routing rate of junction $i$. As described in the main text, we decompose the incoming flux of vehicles $\sigma_i$ to node $i$ as $$\label{sigmadef}
    \sigma_{i} = \sum_{j \in N^{\mbox{\scriptsize in}}_{i}} P_{ji} p_{j} d_{j}\,,$$ where $p_i$ is the probability that a vehicle waiting in node $i$ has not arrived to its destination (i.e., it is going to visit at least another junction in the next step) and $P_{ji}$ is the probability that a vehicle crossing node $j$ goes to node $i\in N^{\mbox{\scriptsize out}}_{j}$ in its next step.

Since vehicles just generated in a certain node are not affected by the congestion in the rest of the network, we separate their contributions in the computation of probabilities $p$ and $P$. Thus, we decompose $p_i$ as $$\begin{aligned}
\label{smallP}
    p_{i} & = & p^{\mbox{\scriptsize gen}}_{i}p^{\mbox{\scriptsize loc}}_{i} + (1-p^{\mbox{\scriptsize gen}}_{i}) p^{\mbox{\scriptsize ext}}_{i}\,,
\end{aligned}$$ where the first term accounts for vehicles generated in node $i$ ($p^{\mbox{\scriptsize gen}}_{i}$) whose destination is not $i$ ($p^{\mbox{\scriptsize loc}}_{i}$) and the second term accounts for vehicles not generated in $i$ whose destination is not $i$ ($p^{\mbox{\scriptsize ext}}_{i}$). Supposing trips consist in traveling through two or more junctions we have that $p^{\mbox{\scriptsize loc}}_{i}=1$. Probability $p^{\mbox{\scriptsize gen}}_{i}$ is equal to the fraction of vehicles generated in $i$ with respect to the total amount of incoming vehicles: $$\begin{aligned}
    p^{\mbox{\scriptsize gen}}_{i} = \frac{g_i}{g_i+\sigma_i}\,.
\end{aligned}$$ Considering the distribution of origins, destinations, the routing strategy and the congestion in the network, probability $p^{\mbox{\scriptsize ext}}_{i}$ can be expressed in terms of the effective node betweenness $\tilde{B}_i$ and the effective vehicle arrivals $\tilde{e}_i$ (the amount of vehicles with destination node $i$ that arrive to node $i$ at each time step): $$\begin{aligned}
    p^{\mbox{\scriptsize ext}}_{i}=\frac{\tilde{B}_i}{\tilde{B}_i + \tilde{e}_i} \,.
\end{aligned}$$ The effective betweenness $\tilde{B}_i$ of a node $i$ accounts for the expected amount of vehicles each node $i$ receives per unit time considering the routing algorithm and the overall congestion of the network. See *Materials and Methods* subsection *Effective betweenness in congested transportation networks* for an extended description and computation of the effective node betweenness $\tilde{B}_i$ and the effective vehicle arrivals $\tilde{e}_i$.

In the same spirit, we decompose the probability $P_{ji}$ that a vehicle waiting in node $j$ goes to node $i$ as: $$\begin{aligned}
\label{BigP}
    P_{ji} &=& p^{\mbox{\scriptsize rgen}}_{j}P^{\mbox{\scriptsize loc}}_{ji} + (1-p^{\mbox{\scriptsize rgen}}_{j})P^{\mbox{\scriptsize ext}}_{ji}\,.
\end{aligned}$$ The first term corresponds to the routed vehicles generated in node $j$ ($p^{\mbox{\scriptsize rgen}}_{j}$) that go to node $i$ ($P^{\mbox{\scriptsize loc}}_{ji}$) and the second term to the routed vehicles not generated in $j$ that go to node $i$ ($P^{\mbox{\scriptsize ext}}_{ji}$). Similarly as before, $p^{\mbox{\scriptsize rgen}}_{j}$ can be expressed as the rate between the vehicles generated in $j$ and the total amount of routed vehicles: $$\begin{aligned}
    p^{\mbox{\scriptsize rgen}}_{j} = \frac{g_j}{g_j+\sigma_j p^{\mbox{\scriptsize ext}}_{j}}\,,
\end{aligned}$$ and, $P^{\mbox{\scriptsize loc}}_{ji}$ and $P^{\mbox{\scriptsize ext}}_{ji}$ can be computed in terms of the normalized effective edge betweenness of the network: $$\begin{aligned}
    P^{\mbox{\scriptsize loc}}_{ji} &=& \frac{\tilde{E}^{\mbox{\scriptsize loc}}_{ji}}{\displaystyle\sum_{k=1}^S \tilde{E}^{\mbox{\scriptsize loc}}_{jk}} \,, \\
    P^{\mbox{\scriptsize ext}}_{ji} &=& \frac{\tilde{E}^{\mbox{\scriptsize ext}}_{ji}}{\displaystyle\sum_{k=1}^S \tilde{E}^{\mbox{\scriptsize ext}}_{jk}} \,, \label{oneLayerMMC_lastEq}
\end{aligned}$$ where the computation of $\tilde{E}^{\mbox{\scriptsize loc}}_{ji}$ only considers paths that start on node $j$ and $\tilde{E}^{\mbox{\scriptsize ext}}_{ji}$ only considers paths that do not start on node $j$. Equivalently to the effective node betweenness $\tilde{B}_i$, computation of $\tilde{E}^{\mbox{\scriptsize loc}}_{ji}$ and $\tilde{E}^{\mbox{\scriptsize ext}}_{ji}$ consider, if required, all congested junctions in the network, as described in a later section, as well as the distribution of the vehicle sources and destinations. Note that the sum of $E^{\mbox{\scriptsize loc}}_{ji}$ and $E^{\mbox{\scriptsize ext}}_{ji}$ corresponds to the classical edge betweenness. Moreover, $P_{ji}$ is an exact expression before and after the onset of congestion.

Eventually, the MCM is composed by a set of $S$ equations ($\Delta{q}_{i} = g_{i} + \sigma_{i} - d_{i}$), one for each junction, and, in principle, a set of $2S$ unknowns, $\Delta{q}_{i}$ and $d_i$ for each junction. To see that the system is indeed determined we need to note that for congested junctions $\Delta q_{i} > 0$ and, thus, after the transient state, $d_{i}=\tau_i$. For the non-congested junctions we have that $\Delta q_{i}=0$ and consequently $d_i =  g_i + \sigma_i$. In conclusion, for any node $i$, either $d_i=\tau_i$ or $d_i=g_i + \sigma_i$ which reduces the amount of unknowns to $S$.

To solve the model given a fixed generation rate $g_i$, we start by considering that no junction is congested and we solve the set of equations Eqs. – by iteration. It is possible that some nodes exceed their maximum routing rate. If this is the case, we set the node with maximum $d_i$ as congested and we solve the system again. This process is repeated until no new junction exceeds its maximum routing rate.

## Onset of congestion using the Microscopic Congestion Model

Most of the works that consider static routing strategies assume that the generation rate of vehicles is the same for all nodes, $g_i=\rho$. In that case, it is possible to compute the critical generation rate $\rho_c$ such that for any generation rate $\rho > \rho_c$ the network will not be able to route or absorb all the traffic. After this point is reached, the amount of vehicles $Q(t)$ in the network will grow proportionally with time, $Q(t) \propto t$, since some of the vehicles get stacked in the queues of the nodes. This transition to the congested state is characterized using the following order parameter: $$\label{oneLayerOrderParameter}
    \eta(\rho)= \lim_{t\rightarrow\infty}\frac{\langle{\Delta Q}\rangle}{\rho S}\,,$$ where $\langle\Delta Q\rangle$ represents the average increment of vehicles per unit of time in the stationary state. Basically, the order parameter measures the ratio between in-transit and generated vehicles.

In the non-congested phase, the amount of incoming and outgoing vehicles for each node can be computed in terms of the node's algorithmic betweenness $B_i$, see ref. . In particular, $$\label{oneLayerSigma}
    \sigma_{i} = \rho\left(\frac{B_i}{S-1} + 1\right)\,,$$ where the second term inside the parentheses accounts for the fact that, in our model, vehicles are also queued at the destination node, unlike in ref. . When no junction is congested we have that $\Delta{q}_{i} = 0$ for all nodes and consequently $$\label{oneLayerD^i}
    d_i = \rho + \sigma_i = \rho\left(\frac{B_i}{S-1} + 2\right)\,.$$ A node $i$ becomes congested when it is required to process more vehicles than its maximum processing rate, $d_i > \tau$. Thus, the critical generation rate at which the first node, and so the system, reaches congestion is: $$\label{rhocFirstNode}
    \rho_{c} = \min_i \frac{\tau\left(S-1\right)}{B_i + 2(S-1)}\,.$$

This is one of the most important analytical results on transportation networks with static routing strategies. In the following, we show that we can recover Eq. () before the onset of congestion using our MCM approach. After substitution of the expression of the probabilities in Eq. : $$\sigma_i = \sum\limits_j \frac{\rho(B_j + S-1) P^{\mbox{\scriptsize loc}}_{ji} + \sigma_j B_j P^{\mbox{\scriptsize ext}}_{ji}}{(\rho+\sigma_j)(B_j + S-1)} d_j \,,$$ and, given we do not have congestion (i.e., $d_j = \rho + \sigma_j$), it simplifies to $$\label{eq:sigmaMCM}
    \sigma_i = \sum\limits_j \frac{\rho(B_j + S-1) P^{\mbox{\scriptsize loc}}_{ji} + \sigma_j B_j P^{\mbox{\scriptsize ext}}_{ji}}{B_j + S-1}\,.$$ Equation () in matrix form becomes $$(I - M)\boldsymbol{\sigma} = \rho\boldsymbol{\pi} \,,$$ where $$\begin{aligned}
    M_{ij} &=& \frac{B_j P^{\mbox{\scriptsize ext}}_{ji}}{B_j + S - 1}\,,\\
    \pi_i &=& \sum_j P^{\mbox{\scriptsize loc}}_{ji}\,,
\end{aligned}$$ and then $$\boldsymbol{\sigma} = \rho(I-M)^{-1}\boldsymbol{\pi}\,.$$ This expression can be shown to be equivalent to Eq. () by using the following relationship between node and edge betweenness: $$\label{nodeEdgeBWRelation}
    B_i + (S-1)=  \sum_{j} \left(B_{j} P^{\mbox{\scriptsize ext}}_{ji} + (S-1)P^{\mbox{\scriptsize loc}}_{ji}\right)\,.$$ The right hand side corresponds to the accumulated fractions of paths that pass through the neighbors of node $i$ and then go to $i$. Each neighbor contributes with two terms, the paths that go through $j$ coming from other nodes, and the paths that start in $j$.

## Effective betweenness in congested transportation networks

The effective betweenness $\tilde{B}_i$ of a node $i$, as defined in ref. , accounts for the expected amount of vehicles each node $i$ receives per unit time. When the network is not congested and the vehicle generation rate $g_i$ is equal for all nodes, $g_i = \rho$, the number of vehicles each node receives can be obtained using Eq. . However, if the network is congested, the traffic dynamics becomes highly non-linear and the value of $\sigma_i$ computed in Eq.  becomes a poor approximation.

Suppose we focus on a particular congested node $j^{\ast}$ of the network. For $j^{\ast}$, being congested means that it is receiving more vehicles that the ones it can process and route. In particular, from the $\sigma_{j^{\ast}}+g_{j^{\ast}}$ vehicles that arrive to the node, only $\tau_{j^{\ast}}$ can be processed at each time step.

Therefore, the contribution to the effective betweenness $\tilde{B}_i$ of the paths from a source/destination pair, $(s,t)$, that traverse the congested node $j^{\ast}$ before reaching $i$, must be rescaled by the fraction of processable vehicles: $$\label{rescalingEffBetweenness}
    s_{j^{\ast}} = \frac{\tau_{j^{\ast}}}{\sigma_{j^{\ast}} + g_{j^{\ast}}}\,.$$ When a path traverses multiple congested nodes ${j^{\ast}},{k^{\ast}},\dots$, the remaining fraction of paths that will reach the target node will be the result of the application of the multiple re-scalings $s_{x^{\ast}}$.

The computation of $s_{j^{\ast}}$ is not straightforward. In general, $\sigma_i$ is not known after the onset of congestion and depends on the effective betweenness that requires, at the same time, to know the $s_{j^{\ast}}$ fraction for all congested nodes. Thus, an iterative calculation is needed to fit all the parameters at the same time as we do in our Microscopic Congestion Model.

The effective arrivals $\tilde{e}_i$ account for the amount of vehicles with destination node $i$ that arrive to node $i$ at each time step. This value in the non-congested phase can be obtained, considering homogeneous source and destination nodes, as $$\label{effectiveArrivals}
    e_i = \rho(S-1)\,.$$ However, congestion affects the variable $e_i$ as well, and needs to be corrected accordingly using the same procedure presented above.

## Traffic Dynamics

To simulate the traffic dynamics of the road network, we assign a first-in-first-out queue to each junction that simulates the blocking time of vehicles before they are allowed to cross it and continue their trip. We suppose these queues have infinite capacity and a maximum processing rate that simulates the physical constraints of the junction. Vehicles origins and destinations may follow any desired distribution. In this work, we have considered two distributions: a random uniform distribution for the synthetic experiments, and one obtained considering the ingoing and outgoing flux of vehicles of the city of Milan. At each time step (of 1 minute duration) vehicles are generated and arrive to their first junction. During the following time steps, vehicles navigate towards their destination following any routing strategy. Here, we have used two different routing strategies: shortest-path and random local search.

For the particular case of simulating traffic in the city of Milan we assume a traffic dynamics similar to the "Home-to-Work" travel pattern where vehicles arrive from the outskirts of the city, go to the city center and then return to the outskirts. Specifically, in our simulation, traffic is generated in the peripheral junctions of the network, goes to a randomly selected junction within the city and then returns back to a randomly selected peripheral junction. We do not consider trips with origin and destination inside the city center since public transportation systems (e.g., train or subway) usually constitute a better alternative than private vehicles for those trips.

The maximum crossing rate of each junction $\tau_i$ accounts, among others, for the existence of traffic lights governing the junction, the width of the street as well as its traffic. We have not been able to get this information for the studied cities, and consequently we cannot set to each junction its precise value. Instead, without loss of generality and for the sake of simplicity, we set to all junctions the same maximum crossing rate, $\tau_i = 15$ (an estimation of the average of their real values).

# Acknowledgements

This work has been supported by Ministerio de Economı́a y Competitividad (Grant FIS2015-71582-C2-1) and European Comission FET-Proactive Projects MULTIPLEX (Grant 317532). A.A. also acknowledges partial financial support from the ICREA Academia and the James S. McDonnell Foundation.


**Quantum information processing, science of**[^1] - The theoretical, experimental and technological areas covering the use of quantum mechanics for communication and computation. Quantum information processing includes investigations in quantum information theory, quantum communication, quantum computation, quantum algorithms and their complexity, and quantum control. The science of quantum information processing is a highly interdisciplinary field. In the context of mathematics it is stimulating research in pure mathematics (e.g. coding theory, $*$-algebras, quantum topology) as well as requiring and providing many opportunities for applied mathematics.

The science of quantum information processing emerged from the recognition that usable notions of information need to be physically implementable. In the 1960s and 1970s researchers such as R. Landauer, C. Bennett, C. Helstrom and A. Holevo realized that the laws of physics give rise to fundamental constraints on the ability to implement and manipulate information. Landauer repeatedly stated that "information is physical", providing impetus to the idea that it should be possible to found theories of information on the laws of physics. This is in contrast to the introspective approach which led to the basic definitions of computer science and information theory as formulated by A. Church, A. Turing, C. Shannon and others in the first half of the 20th century.

Early work in studying the physical foundations of information focused on the effects of energy limitations and the need for dissipating heat in computation and communication. Beginning with S. Wiesner's work on applications of quantum mechanics to cryptography in the late 1960s, it was realized that there may be intrinsic advantages to using quantum physics in information processing. Quantum cryptography and quantum communication in general were soon established as interesting and non-trivial extensions of classical communication based on bits. That quantum mechanics may be used to improve the efficiency of algorithms was first realized when attempts at simulating quantum mechanical systems resulted in exponentially complex algorithms compared to the physical resources associated with the system simulated. In the 1980s, P. Benioff and R. Feynman introduced the idea of a quantum computer for efficiently implementing quantum physics simulations. Models of quantum computers were developed by D. Deutsch, leading to the formulation of artificial problems that could be solved more efficiently by quantum than by classical computers. The advantages of quantum computers became widely recognized when P. Shor (1994) discovered that they can be used to efficiently factor large numbers — a problem believed to be hard for classical deterministic or probabilistic computation and whose difficulty underlies the security of widely used public key encryption methods. Subsequent work established principles of quantum error-correction to ensure that quantum information processing was robustly implementable. See  for introductions to quantum information processing and a quantum mechanics tutorial.

In the context of quantum information theory, information in the sense of C. Shannon is referred to as *classical* information. The fundamental unit of classical information is the *bit*, which can be understood as an ideal system in one of two states or configurations, usually denoted by $0$ and $1$. The fundamental units of quantum information are qubits (short for "quantum bits"), whose states are identified with all "unit superpositions" of the classical states. It is common practice to use the bra-ket conventions for denoting states. In these conventions, the classical configurations are denoted by $|{0}\rangle$ and $|{1}\rangle$, and superpositions are formal sums $\alpha|{0}\rangle+\beta|{1}\rangle$, where $\alpha$ and $\beta$ are complex numbers satisfying $|\alpha|^2+|\beta|^2 = 1$. The states $|{0}\rangle$ and $|{1}\rangle$ represent a standard orthonormal basis of a two-dimensional Hilbert space. Their superpositions are unit vectors in this space. The state space associated with $n>1$ qubits is formally the tensor product of the Hilbert spaces of each qubit. This state space can also be obtained as an extension of the state space of $n$ classical bits by identifying the classical configurations with a standard orthonormal basis of a $2^n$ dimensional Hilbert space.

Access to qubit states is based on the postulates of quantum mechanics with the additional restriction that they are *local* in the sense that elementary operations apply to one or two qubits at a time. Most operations can be expressed in terms of standard measurements of a qubit and two-qubit *quantum gates*. The standard qubit measurement has the effect of randomly projecting the state of the qubit onto one of its classical states; this state is an output of the measurement (accessible for use in a classical computer if desired). For example, using the tensor product representation of the state space of several qubits, a measurement of the first qubit is associated with the two projection operators $P^{(1)}_0=P_0\otimes I\otimes\ldots$ and $1-P^{(1)}_0$, where $P_0|{0}\rangle = |{0}\rangle$ and $P_0|{1}\rangle = 0$. If $\mathbf{\psi}$ is the initial state of the qubits, then the measurement outcome is $0$ with probability $p_0=\|P_0\mathbf{\psi}\|^2$, in which case the new state is $P_0\mathbf{\psi}/p_0$, and the outcome is $1$ with probability $1-p_0=\|P_1\mathbf{\psi}\|^2$ with new state $P_1\mathbf{\psi}/(1-p_0)$. This is a special case of a *von Neumann measurement*. A general two-qubit quantum gate is associated with a unitary operator $U$ acting on the state space of two qubits. Thus $U$ may be represented by a $4\times 4$ unitary matrix in the standard basis of two qubits. The quantum gate may be *applied* to any two chosen qubits. For example, if the state of $n$ qubits is $\mathbf{\psi}$ and the gate is applied to the first two qubits, then the new state is given by $(U\otimes I\otimes\ldots
)\mathbf{\psi}$. Another important operation of quantum information processing is preparation of the $|{0}\rangle$ state of a qubit, which can be implemented in terms of a measurement and subsequent applications of a gate depending on the outcome.

Most problems of theoretical quantum information processing can be cast in terms of the elementary operations above, restrictions on how they can be used and an accounting of the physical *resources* or *cost* associated with implementing the operations. Since classical information processing may be viewed as a special case of quantum information processing, problems of classical information theory and computation are generalized and greatly enriched by the availability of quantum superpositions. The two main problem areas of theoretical quantum information processing are quantum computation and quantum communication.

In studies of quantum computation (cf. **quantum computation**) one investigates how the availability of qubits can be used to improve the efficiency of algorithmic problem solving. Resources counted include the number of quantum gates applied and the number of qubits accessed. This can be done by defining and investigating various types of quantum automata, most prominently quantum Turing machines, and studying their behavior using approaches borrowed from the classical theory of automata and languages. It is convenient to combine classical and quantum automata, for example by allowing a classical computer access to qubits as defined above, and then investigating the complexity of algorithms by counting both classical and quantum resources, thus obtaining trade-offs between the two.

Most of the complexity classes for classical computation have analogues for quantum computation, and an important research area is concerned with establishing relationships between these complexity classes. Corresponding to the classical class $\mathbf{P}$ of polynomially decidable languages is the class of languages decidable in bounded error quantum polynomial time, $\mathbf{BQP}$. While it is believed that $\mathbf{P}$ is properly contained in $\mathbf{BQP}$, whether this is so is at present an open problem. $\mathbf{BQP}$ is known to be contained in the class $\mathbf{P}^{\raisebox{.5pt}{\#}\mathbf{P}}$ (languages decidable in classical polynomial time given access to an oracle for computing the permanent of $0$-$1$ matrices), but the relationship of $\mathbf{BQP}$ to the important class of nondeterministic polynomial time languages $\mathbf{NP}$ is not known.

In quantum communication one considers the situation where two or more entities with access to local qubits can make use of both classical and quantum (communication) channels for exchanging information. The basic operations now include the ability to send classical bits and the ability to send quantum bits. There are two main areas of investigation in quantum communication. The first aims at determining the advantages of quantum communication for solving classically posed communication problems with applications to cryptography and to distributed computation. The second is concerned with establishing relationships between different types of communication resources, particularly with respect to noisy quantum channels, thus generalizing classical communication theory.

Early investigations of quantum channels focused on using them for transmitting classical information by encoding a source of information (cf. **information, source of**) with uses of a quantum channel (cf. **quantum communication channel**). The central result of these investigations is A. Holevo's bound (1973) on the amount of classical information that can be conveyed through a quantum channel. Asymptotic achievability of the bound (using block coding of the information source) was shown in the closing years of the twentieth century. With some technical caveats, the bound and its achievability form a quantum information-theoretic analogue of Shannon's capacity theorem for classical communication channels.

Quantum cryptography, distributed quantum computation and quantum memory require transmitting (or storing) quantum states. As a result it is of great interest to understand how one can communicate quantum information through quantum channels. In this case, the source of information is replaced by a source of quantum states, which are to be transmitted through the channel with high *fidelity*. As in the classical case, the state is encoded before transmission and decoded afterwards. There are many measures of fidelity which may be used to evaluate the quality of the transmission protocol. They are chosen so that a good fidelity value implies that with high probability, quantum information processing tasks behave the same using the original or the transmitted states. A commonly used fidelity measure is the Bures-Uhlmann fidelity, which is an extension of the Hilbert space norm to probability distributions of states (represented by *density operators*). In most cases, asymptotic properties of quantum channels do not depend on the details of the fidelity measure adopted.

To improve the reliability of transmission over a noisy quantum channel, one uses *quantum error-correcting codes* to encode a state generated by the quantum information source with multiple uses of the channel. The theory of quantum codes can be viewed as an extension of classical coding theory. Concepts such as minimum distance and its relationship to error-correction generalize to quantum codes. Many results from the classical theory, including some linear programming upper bounds and the Gilbert-Varshamov lower bounds on the achievable rates of classical codes have their analogues for quantum codes. In the classical theory, linear codes are particularly useful and play a special role. In the quantum theory, this role is played by the *stabilizer* or *additive* quantum codes, which are in one-to-one correspondence with self-dual (with respect to a specific symplectic inner product) classical $\mathrm{GF}_2$-linear codes over $\mathrm{GF}_4$ (cf. **finite fields**).

The capacity of a quantum channel with respect to encoding with quantum codes is not as well understood as the capacity for transmission of classical information. The exact capacity is known only for a few special classes of quantum channels. Although there are information theoretic upper bounds, they depend on the number of channel instances, and whether or not they can be achieved is an open problem. A further complication is that the capacity of quantum channels depends on whether one-way or two-way classical communication may be used to restore the transmitted quantum information .

The above examples illustrate the fact that there are many different types of information utilized in quantum information theory, making it a richer subject than classical information theory. Another physical resource whose properties appear to be best described by information-theoretic means is *quantum entanglement*. A quantum state of more than one quantum system (e.g. two qubits) is said to be entangled if the state can not be factorized as a product of states of the individual quantum systems. Entanglement is believed to play a crucial role in quantum information processing, as demonstrated by its enabling role in effects such as quantum key distribution, superdense coding, quantum teleportation, and quantum error-correction. Beginning in 1995 an enormous amount of effort has been devoted to understanding the principles governing the behavior of entanglement. This has resulted in the discovery of connections between quantum entanglement and classical information theory, the theory of positive maps  and majorization .

The investigation of quantum channel capacity, entanglement, and many other areas of quantum information processing involves various quantum generalizations of the notion of entropy, most notably the von Neumann entropy. The von Neumann entropy is defined as $H(\rho)=\mbox{tr}\rho\log_2(\rho)$ for density operators $\rho$ ($\rho$ is positive Hermitian and of trace $1$). It has many (but not all) of the properties of the classical information function $H(\cdot)$ (cf. **information, amount of**). Understanding these properties has been crucial to the development of quantum information processing (see  for reviews). Probably the most powerful known result about the von Neumann entropy is the strong subadditivity inequality. Many of the bounds on quantum communication follow as easy corollaries of strong subadditivity. Whether still more powerful entropic inequalities exist is not known.

An important property of both classical and quantum information is that although it is intended to be physically realizable, it is abstractly defined and therefore independent of the details of a physical realization. It is generally believed that qubits encapsulate everything that is finitely realizable using accessible physics. This belief implies that any information processing implemented by available physical systems using resources appropriate for those systems can be implemented as efficiently (with at most polynomial overhead) using qubits. It is noteworthy that there is presently no proof that information processing based on quantum field theory (cf. **quantum field theory**) is not more efficient than information processing with qubits. Furthermore, the as-yet unresolved problem of combining quantum mechanics with general relativity in a theory of quantum gravity prevents a fully satisfactory analysis of the information processing power afforded by fundamental physical laws.

Much effort in the science of quantum information processing is being expended on developing and testing the technology required for implementing it. An important task in this direction is to establish that quantum information processing can be implemented robustly in the presence of noise. At first it was believed that this was not possible. Arguments against the robustness of quantum information were based on the apparent relationship to analogue computation (due to the continuity of the amplitudes in the superpositions of configurations) and the fact that it seemed difficult to observe quantum superpositions in nature (due to the rapid loss of phase relationships called *decoherence*). However, the work on quantum error-correcting codes rapidly led to the realization that provided the physical noise behaves locally and is not too large, it is at least in principle possible to process quantum information fault tolerantly. Research in how to process quantum information reliably continues; the main problems is improving the estimates on the maximum amount of tolerable noise for general models of quantum noise and for the types of noise expected in specific physical systems. Other issues include the need to take into consideration restrictions imposed by possible architectures and interconnection networks.

There are many physical systems that can potentially be used for quantum information processing . An active area of investigation involves determining the general mathematical features of quantum mechanics required for implementing quantum information. More closely tied to existing experimental techniques are studies of specific physical systems. In the context of communication, optical systems are likely to play an important role, while for computation there are proposals for using electrons or nuclei in solid state, ions or atoms in electromagnetic traps, excitations of superconductive devices etc. In all of these, important theoretical issues arise. These issues include how to optimally use the available means for controlling the quantum systems (*quantum control*), how to best realize quantum information (possibly indirectly), what architectures can be implemented, how to translate abstract sequences of quantum gates to physical control actions, how to interface the system with optics for communication, refining the theoretical models for how the system is affected by noise and thermodynamic effects, and how to reduce the effects of noise.

|               |
|:--------------|
| E. H. Knill   |
| M. A. Nielsen |

AMS 2000 Subject Classification: 81P68, 68Q05

[^1]: Article by E. H. Knill and M. A. Nielsen accepted for Supplement III, Encyclopaedia of Mathematics (publication expected Summer 2001). See also "http://www.wkap.nl/series.htm/ENM". E. H. Knill is with the Los Alamos National Laboratory, MS B265, Los Alamos NM 87545, USA, and M. A. Nielsen is with the Center for Quantum Computer Technology, Department of Physics, University of Queensland 4072, Australia.


abstract: Are there parallels between the furthest reaches of our universe, and the foundations of thought, awareness, perception, and emotion? What are the connections between the webs and structures that define both? What are the differences? "As Above As Below" was an exhibition that examined these questions. Conceptualized, curated, and produced by Esther Mallouh, it consisted of six artworks, each of them the product of a collaboration that included at least one artist, astrophysicist, and neuroscientist. The installations explored new parallels between intergalactic and neuronal networks through media such as digital projection, virtual reality, and interactive multimedia, and served to illustrate diverse collaboration practices and ways to communicate across very different fields.
author: Mark Neyrinck, Tamira Elul, Michael Silver, Esther Mallouh,; Miguel Aragón-Calvo, Sarah Banducci, Cory Bloyd, Thea Boodhoo, Benedikt Diemer,; Bridget Falck, Dan Feldman, Yoon Chung Han, Jeffrey Kruk, Soo Jung Kwak,; Yagiz Mungan, Miguel Novelo, Rushi Patel, Purin Phanichphant, Joel Primack,; Olaf Sporns, Forest Stearns, Anastasia Victor, David Weinberg, Natalie M. Zahr
bibliography: refs.bib; refs_res.bib
date: SciArt Magazine, February 2020 <https://www.sciartmagazine.com/collaboration-as-above-as-below.html>
title: Exploring Connections Between Cosmos & Mind
  Through Six Interactive Art Installations in
  "As Above As Below"

A short video walkthrough of the exhibition is at <https://youtu.be/RM3gDjb3phU>, and a panel discussion at the opening, featuring many of the authors, is at <https://youtu.be/qNhtyvPz8yw>. See <https://asaboveasbelow.com/> for even more information about the project.

# Introduction

The visual similarities between cosmic and neural webs are striking. Fig. A shows a single mouse hippocampal neuron with an elaborate dendritic arbor stemming from a central cell body (large green blob); Fig. B shows a cluster from a cosmological simulation, in which a whole galaxy like ours would be but a middling yellow dot.

However, many networks in the universe exhibit these kinds of visual similarities. For example, other biological systems, such as circulatory and respiratory networks in animals, and branching structures in plants and fungi, appear individually to have nearly maximal fractal dimension, and therefore look similar . These features are also evident in networks created by humans, such as traffic in road and air-travel systems, and non-spatial networks such as the Internet. See for more about the network-science properties of the cosmic web.

With so many types of similar spatial networks to pick from, why focus on cosmic and neural webs? One compelling reason is that they are each fascinating on their own: one is the network that is able to reflect on and study itself, and the scale of the other is nearly as big as the observable Universe. It is also curious that they are of similar complexity, with numbers of nodes that are within a couple of orders of magnitude of each other .

The term "cosmic web" was first coined by . It refers to the network of filaments and walls of dark matter and gas that connect nearby galaxies, in the standard cosmological paradigm. The cosmic web exists on multiple spatial scales; there are also filaments of galaxies that connect clusters of galaxies. The cosmic web, and structural trussworks like spiderwebs and trees, also share a geometry that has been a source of artistic inspiration . To anthropomorphize: gravity acts like a haunted house explorer, clearing away flexible cobwebs from cosmic voids and gathering structural elements together into thicker filaments . The cosmic web moves at an aeonically slow pace in the expanding Universe, set in motion just after the Big Bang and nearly deterministically forming since then.

The cosmic web is not a pathway of communication – the distances are far too great. But the whole point of the neural net is communication, on highways under constant construction and remodeling as the brain processes information and learns. In fact, interneuronal communication is likely to be the driving force in the evolution of neural webs. Like cosmic webs, neural webs exist on various scales: nervous systems are composed of neurons and clusters of neurons that exhibit network organization , dendrites and axons of individual neurons have complex branching structures (Fig. A), and each neuron contains a molecular cytoskeleton that is made up of a network of microfilaments.

The principles underlying the development of specific neuron connectivity patterns in the neural web, hypothesized originally by the Spanish neuroanatomist Santiago Ramon y Cajal (and still relevant for theoretical and experimental neuroscientific studies today) are maximum efficiency of communication and meeting the metabolic challenge of making and supporting new connections . This leads to neurons that develop morphologies that approximate minimal spanning trees that contain a set of nodes that are branch points and/or synapse locations (Cuntz, 2012). Neurons have also inspired many artistic interpretations and visualizations .

There are yet more differences between neural and cosmic webs. Neurons have objective physical membrane boundaries, while the boundaries of the filaments and walls of the cosmic web are subject to definition . Also, each neuron has a defined polarity (dendritic input and axonal output). In the cosmic web, galaxy clusters can locally look like neurons, but there is no analogous polarity, and no objective boundaries separating galaxy clusters. Another difference is that neurons have objective centers (somas, containing cell nuclei), but the Universe does not, nor do patches of the cosmic web. Any given region of the Universe on really large scales looks pretty much like any other large region of the Universe – the large scale architecture of the observable Universe doesn't have any special regions. However, our brains are organized very differently – vertebrate brains are networks with clearly identified specialized regions.

# Art Installations

## *Chamber of (In)finite Potential*

An interactive installation by Purin Phanichphant (artist), Benedikt Diemer (astrophysicist), and Natalie Zahr (neuroscientist).

In this piece, we consider a fundamental similarity between the cosmic web of dark matter in the Universe and the complex connections of the human brain: the Universe has the potential to form a nearly limitless number of structures (such as galaxies); brain cells have the potential to form a nearly limitless number of connections. However, both the Universe and the brain have – at this point in time – achieved only a fraction of their potential.

Furthermore, both fields of study are in the midst of lively debates as to whether the current understanding fully explains reality. While we tend to think of the Universe as infinite, the growth of structure within it is not: as dark energy accelerates cosmic expansion, gravity will eventually be unable to overcome the stretching, and structures will be frozen. Similarly, while the abundance of neurons in the brain is limited (at birth, the brain has almost all the neurons that it will ever have), the number of connections is potentially infinite; nevertheless, the number of connections reaches a peak and then begins to decay with aging.

These considerations lead to fundamental questions about the nature of the two systems. Will there be a countable number of galaxies in the Universe at the end of its life? Is the capability of our minds fundamentally limited by a finite number of neurons or their connections? These are the kinds of questions we seek to explore in our installation.

The seeming contradiction between the infinite and finite nature is expressed through a six-foot wide geodesic dome that is finite and graspable, but with reflective surfaces within that create an infinite visual field. When the viewer enters the dome, his/her perspective is from the center of this perceptual field, metaphorically at the center of the Universe and self-awareness. The artist aims to reveal the infinite potential that exists inside each and every one of us, and that our existence, both in the neural and cosmic realms, are one and the same.

## *Natural Science*

An interactive installation by Thea Boodhoo (scenery design, concept), Cory Bloyd (interactive design, concept), Yagiz Mungan (sound), David Weinberg (astrophysicist), and Dan Feldman (neuroscientist), with contributions by Gary Boodhoo (concept), James Morgan (concept, materials), and Lyn Collier (feltwork).

Natural science, or the scientific study of nature, includes ecology, biology, geology, and physics, as well as cosmology and neuroscience. This surpasses the everyday meaning of nature, going far beyond wildlife and forests to include the farthest reaches of the Universe, as well as the deepest inner workings of the mind.

It could be concluded, from this point of view, that all endeavors of natural science are in fact nature studying itself. This piece invites you to broaden your definition of nature, and what you consider natural, dramatically. Don't forget to include yourself.

The Science: What do the brain at rest and the cosmos at rest have in common? Science finds that both are surprisingly active. Brain waves cascade from sleeping brains; gravitational waves flood the quietest regions of space. Dan Feldman studies how our neurons process, and often create, what we experience. His ongoing research inspired the setup of our "experiment," as well as a question: could we someday step into another being's dream? What would we encounter? One answer comes from astrophysics. Every brain, and so every dream, is part of the cosmos – David H. Weinberg's specialty. His research asks: How do galaxies form and cluster? What happens when black holes merge? What was the Universe like before stars, and what will it be like when it finally rests?... Will it dream?

The Art: A wild fox has chosen a secluded mossy outcrop as a place to nap before getting back to fox business. Little does he know, researchers studying the inner workings of the mind are running an experiment in this part of the woods. Secret prototype devices tune into the dreams of any animal who sleeps here.

The fox's reverie is projected back as sound and vision, and its contents are quite a surprise. Explorers in these woods will hear massive celestial bodies merging and tiny neuron cells firing. If they approach quietly, they may even influence the fox's dream – for they, like dreams, are part of the cosmos.

## *I am a  Neuron*

An interactive installation by Yoon Chung Han (artist), Soo Jung Kwak (artist, sound), Rushi Patel (engineer, software), Jeffrey Kruk (astrophysicist), and Tamira Elul (neuroscientist).

This installation is inspired by the concept of interactions among neighboring galaxies over cosmic time and among neurons during development of the brain. In the forming and expanding Universe, galaxies interact through gravity, winds, and radiation, whereas in the developing brain, neurons are influenced by mechanical, molecular, and electrical cues from other neurons.

*I am a  Neuron* is an interactive WebVR artwork of creative exploration within a virtual cosmic world. Viewers create their own personalized 3D artificial neurons, interact with them, and observe their interactions. Live feedback from participants controls three factors (gravity, winds, and radiation) to change the environment and neurons. When neurons collide with other neurons, special sounds are emitted, and changes occur in the neurons such as scale, design, opacity, and colors that visually represent interactions of actual neurons in our brains.

Through this artwork, participants also explore self-representations that interact with other humans and grow, die, interact, collide, metabolize, reproduce and emit sounds. The participants choose one adjective that represents their personalities for a virtual 3D neuron. Participants relate to these neuron avatars by observing how audiovisual representations change through time. The morphing neurons mimic how we are changed through environmental factors, social groups, and human relationships. Through this artwork participants can find personal meanings into what nature is.

One factor contributing to the success of this collaboration was that the scientists trusted the artist to develop an interpretation of this concept independently, without further trying to influence her process beyond the initial concept.

## *One in the Universe*

A virtual reality experience by Anastasia Victor (artist), Mark Neyrinck (astrophysicist), Miguel Angel Aragón-Calvo (astrophysicist), and Michael Silver (neuroscientist).

Despite vast differences in spatial scale, there is a conceptual link between networks in the brain and those in the cosmic web: they both contain branching structures that connect neurons and nodes of dark matter to one another. Emergent network properties that arise from connectivity of individual elements are also evident for human beings and their relationships.

*One in the Universe* is an interactive installation that asynchronously connects people across time and space. Within the current climate of filter bubbles and social division, this piece provides moments of introspection and connection to others.

Participants in this installation navigate through a virtual reality composed of structures inspired by dark matter webs and the shapes of neurons, and their movements guide the growth of neural branches within this space. Participants also provide verbal responses to prompts about their own experiences of feeling connected, and these responses are recorded and integrated into the virtual Universe.

The virtual Universe is additive, with each new participant's movements through the space and their verbal responses contributing to the growing patterns of connectivity. Participants directly experience this connectivity by traversing the virtual Universe and hearing verbal responses of previous participants at specific locations in the network. Other visitors to the exhibit observe participants' interactions with the virtual Universe through a projected image that mirrors the visual experience in the virtual reality headset.

A scientific project grew out of this piece, as well: we found preliminarily, and cannot yet explain, that the scaling of the total wiring length of the cosmic web with the number of its nodes is similar to the scaling of the wiring length of a neuron with its branch points. They scale roughly with a power law obeyed by a minimal spanning tree, found for neurons by . We came upon this relationship through using the Cuntz et al. algorithm to grow branches of the neuron based on input from the participant's hands.

## *voi!drum!emory*

An interactive installation by Miguel Novelo (artist), Bridget Falck (astrophysicist), and Sarah Banducci (neuroscientist).

A void is not merely an absence, a loss, an emptiness, or a lack, though it can be these things. By calling something a void, we give it a life of its own. It becomes a negative space that defines the positive. In the Universe, voids push matter away, and where they collide with each other, they create the cosmic web, and galaxies. In brains, the web of interconnected neurons deteriorates due to aging and disease; empty space takes over as the brain atrophies. Both types of voids are not static: they reflect the growing, aging, and dying of the Universe and the people that live in it. Our team of a neuroscientist (Sarah Banducci), a cosmologist (Bridget Falck), and an artist (Miguel Novelo) explored this connection between brain and cosmic voids over a year of emails and video calls.

The collaboration kicked off with Falck and Banducci, two strangers on opposite sides of the country, sharing their love of science rooted in two very different fields. Pleasantries aside, both scientists dove past the superficial and into the meat of their respective interests. The first brainstorm ended with a list of four potential ideas for overlap between cosmology and neuroscience: density, connectivity, folds, and voids. By the end of the second call two weeks later, it became clear that voids would provide the greatest opportunity for exploration between the two fields.

For Falck, voids are exciting because they are the largest structures in the cosmic web – except they aren't structures at all. Voids are the emptiest regions of space, and they are growing. As the Universe ages, gravity causes matter to collapse in on itself and pushes matter away from voids, making relatively empty regions of space become even emptier. The aging of the Universe thus grows empty space.

On a much smaller scale, Banducci views voids as a byproduct of human aging. In contrast to the cosmic web, where the force of pulling matter together forms emptiness, in the brain, voids are formed by matter deteriorating. This deterioration of neurons and connections between them gives the appearance of a shrinking brain. Evidence of this atrophy looks different for each individual but translates to characteristic cognitive and motor impairments in older adults.

As an artist, Novelo's initial thoughts were to materialize the void with objects that might have a void or vacuum. From an empty cylinder, his sketch led to a percussive instrument. Once a drum is hit on the top membrane, air moves inside and sound waves reverberate in concave space. This sound illustrates a wholeness of the object, the echo, and decay symbolizing the void in a sonic way. Taking the concept a step further, Novelo explored the idea of experiencing a void by sensing the end of reverberation, paying attention to the time between echoes and the sound fading out. This became the physical representation of the void.

To add an emotional component to this experience, the team considered, "What is left after the void?" Discussing memory in this context, Novelo, Falck, and Banducci came up with poetic voids – remembrance of what is not there anymore. In the cosmic web, what is left could be the photons finally reaching us from an event that happened millions of years ago; in the brain, the rest of our neurons or gray matter are still intact; and in our memory, previous experiences, nostalgia, and trauma. The art object we conceptualized and created has the characteristic of activating sound and visual cues that symbolize time prior to the void, rapid decay, and lastly a void surrounded by memory.

The final expression of this project should leave the viewer with the assurance that even when there is a void, some piece of the experience will continue to linger.

## *The Undulating Architecture Illuminating the Individual Sciences*

A video installation by Forest Stearns (artist), Joel Primack (astrophysicist), and Olaf Sporns (neuroscientist).

As a collective, our team of three spent a good amount of time investigating the superficial similarities between dark matter of the universe and the neural networks of our brains. After much discussion, we collectively agreed that the differences greatly superseded the similarities.

This being the case, the artist focused on the architectures of both sciences. Using soap bubbles as a medium, the artist found that the undulating bubbles with their walls, filaments, and nodes illuminated the architecture of the universe. The relationships between cosmic webs and soap bubble foams has also been explored by . The soap running down filaments towards central nodes and into a bigger pool is similar to how impulses move across the neural net. The physical mass of bubbles were an abstracted screen, projected onto videos that graphically explain the specific sciences. The camera collected the photons rippling across the bubble field.

# Acknowledgments

MCN is grateful for funding from Basque Government grant IT956-16.

The first three authors were primarily responsible for this article as it appears, but all others contributed ideas and/or text. Curator Esther Mallouh originated the project.