Multidimensional network

In network theory, multidimensional networks, a special type of multilayer network, are networks with multiple kinds of relations.^{[1] [2][3][4][5][6][7]} Increasingly sophisticated attempts to model real-world systems as multidimensional networks have yielded valuable insight in the fields of social network analysis,^{[3][4][8][9][10][11][12]} economics, urban and international transport,^[13][14][15] ecology,^{[16][17][18][19]} psychology,^[20][21] medicine, biology,^[22] commerce, climatology, physics,^[23] computational neuroscience,^{[24][25][26][27]} operations management, and finance.

The rapid exploration of complex networks in recent years has been dogged by a lack of standardized naming conventions, as various groups use overlapping and contradictory^[28][29] terminology to describe specific network configurations (e.g., multiplex, multilayer, multilevel, multidimensional, multirelational, interconnected). To fully leverage the dataset information on the directional nature of the communications, some authors consider only direct networks without any labels on vertices, and introduce the definition of edge-labeled multigraphs which can cover many multidimensional situations.^[30] The term "fully multidimensional" has also been used to refer to a multipartite edge-labeled multigraph.^[31] Multidimensional networks have also recently been reframed as specific instances of multilayer networks.^{[1] [5][6][32]} In this case, there are as many layers as there are dimensions, and the links between nodes within each layer are simply all the links for a given dimension.

In elementary network theory, a network is represented by a graph G = ( V , E ) {\displaystyle G=(V,E)} in which V {\displaystyle V} is the set of nodes and E {\displaystyle E} the links between nodes, typically represented as a tuple of nodes u , v ∈ V {\displaystyle u,v\in V} . While this basic formalization is useful for analyzing many systems, real world networks often have added complexity in the form of multiple types of relations between system elements. An early formalization of this idea came through its application in the field of social network analysis (see, e.g.,^[33] and papers on relational algebras in social networks) in which multiple forms of social connection between people were represented by multiple types of links.^[34]

To accommodate the presence of more than one type of link, a multidimensional network is represented by a triple G = ( V , E , D ) {\displaystyle G=(V,E,D)} , where D {\displaystyle D} is a set of dimensions (or layers), each member of which is a different type of link, and E {\displaystyle E} consists of triples ( u , v , d ) {\displaystyle (u,v,d)} with u , v ∈ V {\displaystyle u,v\in V} and d ∈ D {\displaystyle d\in D} .^[6]

Note that as in all directed graphs, the links ( u , v , d ) {\displaystyle (u,v,d)} and ( v , u , d ) {\displaystyle (v,u,d)} are distinct.

By convention, the number of links between two nodes in a given dimension is either 0 or 1 in a multidimensional network. However, the total number of links between two nodes across all dimensions is less than or equal to | D | {\displaystyle |D|} .

In the case of a weighted network, this triplet is expanded to a quadruplet e = ( u , v , d , w ) {\displaystyle e=(u,v,d,w)} , where w {\displaystyle w} is the weight on the link between u {\displaystyle u} and v {\displaystyle v} in the dimension d {\displaystyle d} .

Further, as is often useful in social network analysis, link weights may take on positive or negative values. Such signed networks can better reflect relations like amity and enmity in social networks.^[31] Alternatively, link signs may be figured as dimensions themselves,^[35] e.g. G = ( V , E , D ) {\displaystyle G=(V,E,D)} where D = { − 1 , 0 , 1 } {\displaystyle D=\{-1,0,1\}} and E = { ( u , v , d ) ; u , v ∈ V , d ∈ D } {\displaystyle E=\{(u,v,d);u,v\in V,d\in D\}} This approach has particular value when considering unweighted networks.

This conception of dimensionality can be expanded should attributes in multiple dimensions need specification. In this instance, links are n-tuples e = ( u , v , d 1 … d n − 2 ) {\displaystyle e=(u,v,d_{1}\dots d_{n-2})} . Such an expanded formulation, in which links may exist within multiple dimensions, is uncommon but has been used in the study of multidimensional time-varying networks.^[36]

Whereas unidimensional networks have two-dimensional adjacency matrices of size V × V {\displaystyle V\times V} , in a multidimensional network with D {\displaystyle D} dimensions, the adjacency matrix becomes a multilayer adjacency tensor, a four-dimensional matrix of size ( V × D ) × ( V × D ) {\displaystyle (V\times D)\times (V\times D)} .^[3] By using index notation, adjacency matrices can be indicated by A j i {\displaystyle A_{j}^{i}} , to encode connections between nodes i {\displaystyle i} and j {\displaystyle j} , whereas multilayer adjacency tensors are indicated by M j β i α {\displaystyle M_{j\beta }^{i\alpha }} , to encode connections between node i {\displaystyle i} in layer α {\displaystyle \alpha } and node j {\displaystyle j} in layer β {\displaystyle \beta } . As in unidimensional matrices, directed links, signed links, and weights are all easily accommodated by this framework.

In the case of multiplex networks, which are special types of multilayer networks where nodes can not be interconnected with other nodes in other layers, a three-dimensional matrix of size ( V × V ) × D {\displaystyle (V\times V)\times D} with entries A i j α {\displaystyle A_{ij}^{\alpha }} is enough to represent the structure of the system^[8][37] by encoding connections between nodes i {\displaystyle i} and j {\displaystyle j} in layer α {\displaystyle \alpha } .

In a multidimensional network, the neighbors of some node v {\displaystyle v} are all nodes connected to v {\displaystyle v} across dimensions.

A path between two nodes in a multidimensional network can be represented by a vector r = ( r 1 , … r | D | ) {\displaystyle =(r_{1},\dots r_{|D|})} in which the i {\displaystyle i} th entry in r is the number of links traversed in the i {\displaystyle i} th dimension of G {\displaystyle G} .^[38] As with overlapping degree, the sum of these elements can be taken as a rough measure of a path length between two nodes.

The existence of multiple layers (or dimensions) allows to introduce the new concept of network of layers,^[3] peculiar of multilayer networks. In fact, layers might be interconnected in such a way that their structure can be described by a network, as shown in the figure.

The network of layers is usually weighted (and might be directed), although, in general, the weights depends on the application of interest. A simple approach is, for each pair of layers, to sum all of the weights in the connections between their nodes to obtain edge weights that can be encoded into a matrix q α β {\displaystyle q_{\alpha \beta }} . The rank-2 adjacency tensor, representing the underlying network of layers in the space R L × L {\displaystyle \mathbb {R} ^{L\times L}} is given by

Ψ δ γ = ∑ α , β = 1 L q α β E δ γ ( α β ) {\displaystyle \Psi _{\delta }^{\gamma }=\sum \limits _{\alpha ,\beta =1}^{L}q_{\alpha \beta }E_{\delta }^{\gamma }(\alpha \beta )}

where E δ γ ( α β ) {\displaystyle E_{\delta }^{\gamma }(\alpha \beta )} is the canonical matrix with all components equal to zero except for the entry corresponding to row α {\displaystyle \alpha } and column β {\displaystyle \beta } , that is equal to one. Using the tensorial notation, it is possible to obtain the (weighted) network of layers from the multilayer adjacency tensor as Ψ δ γ = M j δ i γ U i j {\displaystyle \Psi _{\delta }^{\gamma }=M_{j\delta }^{i\gamma }U_{i}^{j}} .^[3]

In a non-interconnected multidimensional network, where interlayer links are absent, the degree of a node is represented by a vector of length | D | : k = ( k i 1 , … k i | D | ) {\displaystyle |D|:\mathbf {k} =(k_{i}^{1},\dots k_{i}^{|D|})} . Here | D | {\displaystyle |D|} is an alternative way to denote the number of layers L {\displaystyle L} in multilayer networks. However, for some computations it may be more useful to simply sum the number of links adjacent to a node across all dimensions.^[3][39] This is the overlapping degree:^[4] ∑ α = 1 | D | k i α {\displaystyle \sum _{\alpha =1}^{|D|}k_{i}^{\alpha }} . As with unidimensional networks, distinction may similarly be drawn between incoming links and outgoing links. If interlayer links are present, the above definition must be adapted to account for them, and the multilayer degree is given by

k i = M j β i α U α β u j = ∑ α , β = 1 L ∑ j = 1 N M j β i α {\displaystyle k^{i}=M_{j\beta }^{i\alpha }U_{\alpha }^{\beta }u^{j}=\sum _{\alpha ,\beta =1}^{L}\sum _{j=1}^{N}M_{j\beta }^{i\alpha }}

where the tensors U α β {\displaystyle U_{\alpha }^{\beta }} and u j {\displaystyle u^{j}} have all components equal to 1. The heterogeneity in the number of connections of a node across the different layers can be taken into account through the participation coefficient.^[4]

When extended to interconnected multilayer networks, i.e. those systems where nodes are connected across layers, the concept of centrality is better understood in terms of versatility.^[10] Nodes that are not central in each layer might be the most important for the multilayer systems in certain scenarios. For instance, this is the case where two layers encode different networks with only one node in common: it is very likely that such a node will have the highest centrality score because it is responsible for the information flow across layers.

As for unidimensional networks, eigenvector versatility can be defined as the solution of the eigenvalue problem given by M j β i α Θ i α = λ 1 Θ j β {\displaystyle M_{j\beta }^{i\alpha }\Theta _{i\alpha }=\lambda _{1}\Theta _{j\beta }} , where Einstein summation convention is used for sake of simplicity. Here, Θ j β = λ 1 − 1 M j β i α Θ i α {\displaystyle \Theta _{j\beta }=\lambda _{1}^{-1}M_{j\beta }^{i\alpha }\Theta _{i\alpha }} gives the multilayer generalization of Bonacich's eigenvector centrality per node per layer. The overall eigenvector versatility is simply obtained by summing up the scores across layers as θ i = Θ i α u α {\displaystyle \theta _{i}=\Theta _{i\alpha }u^{\alpha }} .^[3][10]

As for its unidimensional counterpart, the Katz versatility is obtained as the solution Φ j β = [ ( δ − a M ) − 1 ] j β i α U i α {\displaystyle \Phi _{j\beta }=[(\delta -aM)^{-1}]_{j\beta }^{i\alpha }U_{i\alpha }} of the tensorial equation Φ j β = a M j β i α Φ i α + b u j β {\displaystyle \Phi _{j\beta }=aM_{j\beta }^{i\alpha }\Phi _{i\alpha }+bu_{j\beta }} , where δ j β i α = δ j i δ β α {\displaystyle \delta _{j\beta }^{i\alpha }=\delta _{j}^{i}\delta _{\beta }^{\alpha }} , a {\displaystyle a} is a constant smaller than the largest eigenvalue and b {\displaystyle b} is another constant generally equal to 1. The overall Katz versatility is simply obtained by summing up the scores across layers as ϕ i = Φ i α u α {\displaystyle \phi _{i}=\Phi _{i\alpha }u^{\alpha }} .^[10]

For unidimensional networks, the HITS algorithm has been originally introduced by Jon Kleinberg to rate Web Pages. The basic assumption of the algorithm is that relevant pages, named authorities, are pointed by special Web pages, named hubs. This mechanism can be mathematically described by two coupled equations which reduce to two eigenvalue problems. When the network is undirected, Authority and Hub centrality are equivalent to eigenvector centrality. These properties are preserved by the natural extension of the equations proposed by Kleinberg to the case of interconnected multilayer networks, given by ( M M t ) j β i α Γ i α = λ 1 Γ j β {\displaystyle (MM^{t})_{j\beta }^{i\alpha }\Gamma _{i\alpha }=\lambda _{1}\Gamma _{j\beta }} and ( M t M ) j β i α Υ i α = λ 1 Υ j β {\displaystyle (M^{t}M)_{j\beta }^{i\alpha }\Upsilon _{i\alpha }=\lambda _{1}\Upsilon _{j\beta }} , where t {\displaystyle t} indicates the transpose operator, Γ i α {\displaystyle \Gamma _{i\alpha }} and Υ i α {\displaystyle \Upsilon _{i\alpha }} indicate hub and authority centrality, respectively. By contracting the hub and authority tensors, one obtains the overall versatilities as γ i = Γ i α u α {\displaystyle \gamma _{i}=\Gamma _{i\alpha }u^{\alpha }} and υ i = Υ i α u α {\displaystyle \upsilon _{i}=\Upsilon _{i\alpha }u^{\alpha }} , respectively.^[10]

PageRank, originally introduced to rank web pages, can also be considered as a measure of centrality for interconnected multilayer networks.

It is worth remarking that PageRank can be seen as the steady-state solution of a special Markov process on the top of the network. Random walkers explore the network according to a special transition matrix and their dynamics is governed by a random walk master equation. It is easy to show that the solution of this equation is equivalent to the leading eigenvector of the transition matrix.

Random walks have been defined also in the case of interconnected multilayer networks^[15] and edge-colored multigraphs (also known as multiplex networks).^[40] For interconnected multilayer networks, the transition tensor governing the dynamics of the random walkers within and across layers is given by R j β i α = r T j β i α + ( 1 − r ) N L u j β i α , {\displaystyle R_{j\beta }^{i\alpha }=rT_{j\beta }^{i\alpha }+{\frac {(1-r)}{NL}}u_{j\beta }^{i\alpha },} , where r {\displaystyle r} is a constant, generally set to 0.85, N {\displaystyle N} is the number of nodes and L {\displaystyle L} is the number of layers or dimensions. Here, R j β i α {\displaystyle R_{j\beta }^{i\alpha }} might be named Google tensor and u j β i α {\displaystyle u_{j\beta }^{i\alpha }} is the rank-4 tensor with all components equal to 1.

As its unidimensional counterpart, PageRank versatility consists of two contributions: one encoding a classical random walk with rate r {\displaystyle r} and one encoding teleportation across nodes and layers with rate 1 − r {\displaystyle 1-r} .

If we indicate by Ω i α {\displaystyle \Omega _{i\alpha }} the eigentensor of the Google tensor R j β i α {\displaystyle R_{j\beta }^{i\alpha }} , denoting the steady-state probability to find the walker in node i {\displaystyle i} and layer α {\displaystyle \alpha } , the multilayer PageRank is obtained by summing up over layers the eigentensor: ω i = Ω i α u α {\displaystyle \omega _{i}=\Omega _{i\alpha }u^{\alpha }} ^[10]

Like many other network statistics, the meaning of a clustering coefficient becomes ambiguous in multidimensional networks, due to the fact that triples may be closed in different dimensions than they originated.^[4][41][42] Several attempts have been made to define local clustering coefficients, but these attempts have highlighted the fact that the concept must be fundamentally different in higher dimensions: some groups have based their work off of non-standard definitions,^[42] while others have experimented with different definitions of random walks and 3-cycles in multidimensional networks.^[4][41]

While cross-dimensional structures have been studied previously,^[43][44] they fail to detect more subtle associations found in some networks. Taking a slightly different take on the definition of "community" in the case of multidimensional networks allows for reliable identification of communities without the requirement that nodes be in direct contact with each other.^{[3][8][9][45]} For instance, two people who never communicate directly yet still browse many of the same websites would be viable candidates for this sort of algorithm.

A generalization of the well-known modularity maximization method for community discovery has been originally proposed by Mucha et al.^[8] This multiresolution method assumes a three-dimensional tensor representation of the network connectivity within layers, as for edge-colored multigraphs, and a three-dimensional tensor representation of the network connectivity across layers. It depends on the resolution parameter γ {\displaystyle \gamma } and the weight ω {\displaystyle \omega } of interlayer connections. In a more compact notation, making use of the tensorial notation, modularity can be written as Q ∝ S i α a B j β i α S a j β {\displaystyle Q\propto S_{i\alpha }^{a}B_{j\beta }^{i\alpha }S_{a}^{j\beta }} , where B j β i α = M j β i α − P j β i α {\displaystyle B_{j\beta }^{i\alpha }=M_{j\beta }^{i\alpha }-P_{j\beta }^{i\alpha }} , M j β i α {\displaystyle M_{j\beta }^{i\alpha }} is the multilayer adjacency tensor, P j β i α {\displaystyle P_{j\beta }^{i\alpha }} is the tensor encoding the null model and the value of components of S a i α {\displaystyle S_{a}^{i\alpha }} is defined to be 1 when a node i {\displaystyle i} in layer α {\displaystyle \alpha } belongs to a particular community, labeled by index a {\displaystyle a} , and 0 when it does not.^[3]

Non-negative matrix factorization has been proposed to extract the community-activity structure of temporal networks.^[46] The multilayer network is represented by a three-dimensional tensor T i j τ {\displaystyle T_{ij}^{\tau }} , like an edge-colored multigraph, where the order of layers encode the arrow of time. Tensor factorization by means of Kruskal decomposition is thus applied to T i j τ {\displaystyle T_{ij}^{\tau }} to assign each node to a community across time.

Methods based on statistical inference, generalizing existing approaches introduced for unidimensional networks, have been proposed. Stochastic block model is the most used generative model, appropriately generalized to the case of multilayer networks.^[47][48]

As for unidimensional networks, principled methods like minimum description length can be used for model selection in community detection methods based on information flow.^[9]

Given the higher complexity of multilayer networks with respect to unidimensional networks, an active field of research is devoted to simplify the structure of such systems by employing some kind of dimensionality reduction.^[22][49]

A popular method is based on the calculation of the quantum Jensen-Shannon divergence between all pairs of layers, which is then exploited for its metric properties to build a distance matrix and hierarchically cluster the layers. Layers are successively aggregated according to the resulting hierarchical tree and the aggregation procedure is stopped when the objective function, based on the entropy of the network, gets a global maximum. This greedy approach is necessary because the underlying problem would require to verify all possible layer groups of any size, requiring a huge number of possible combinations (which is given by the Bell number and scales super-exponentially with the number of units). Nevertheless, for multilayer systems with a small number of layers, it has been shown that the method performs optimally in the majority of cases.^[22]

The question of degree correlations in unidimensional networks is fairly straightforward: do nodes of similar degree tend to connect to each other? In multidimensional networks, what this question means becomes less clear. When we refer to a node's degree, are we referring to its degree in one dimension, or collapsed over all? When we seek to probe connectivity between nodes, are we comparing the same nodes across dimensions, or different nodes within dimensions, or a combination?^[6] What are the consequences of variations in each of these statistics on other network properties? In one study, assortativity was found to decrease robustness in a duplex network.^[50]

Given two multidimensional paths, r and s, we say that r dominates s if and only if: ∀ d ∈ ⟨ 1 , | D | ⟩ , r l ≤ s l {\displaystyle \forall d\in \langle 1,|D|\rangle ,r_{l}\leq s_{l}} and ∃ i {\displaystyle \exists i} such that r l < s l {\displaystyle r_{l}<s_{l}} .^[38]

Among other network statistics, many centrality measures rely on the ability to assess shortest paths from node to node. Extending these analyses to a multidimensional network requires incorporating additional connections between nodes into the algorithms currently used (e.g., Dijkstra's). Current approaches include collapsing multi-link connections between nodes in a preprocessing step before performing variations on a breadth-first search of the network.^[28]

One way to assess the distance between two nodes in a multidimensional network is by comparing all the multidimensional paths between them and choosing the subset that we define as shortest via path dominance: let M P ( u , v ) {\displaystyle MP(u,v)} be the set of all paths between u {\displaystyle u} and v {\displaystyle v} . Then the distance between u {\displaystyle u} and v {\displaystyle v} is a set of paths P ⊆ M P {\displaystyle P\subseteq MP} such that ∀ p ∈ P , ∄ p ′ ∈ M P {\displaystyle \forall p\in P,\nexists p'\in MP} such that p ′ {\displaystyle p'} dominates p {\displaystyle p} . The length of the elements in the set of shortest paths between two nodes is therefore defined as the multidimensional distance.^[38]

In a multidimensional network G = ( V , E , D ) {\displaystyle G=(V,E,D)} , the relevance of a given dimension (or set of dimensions) D ′ {\displaystyle D'} for one node can be assessed by the ratio: Neighbors ( v , D ′ ) Neighbors ( v , D ) {\displaystyle {\frac {{\text{Neighbors}}(v,D')}{{\text{Neighbors}}(v,D)}}} .^[39]

In a multidimensional network in which different dimensions of connection have different real-world values, statistics characterizing the distribution of links to the various classes are of interest. Thus it is useful to consider two metrics that assess this: dimension connectivity and edge-exclusive dimension connectivity. The former is simply the ratio of the total number of links in a given dimension to the total number of links in every dimension: | { ( u , v , d ) ∈ E | u , v ∈ V } | | E | {\displaystyle {\frac {|\{(u,v,d)\in E|u,v\in V\}|}{|E|}}} . The latter assesses, for a given dimension, the number of pairs of nodes connected only by a link in that dimension: | { ( u , v , d ) ∈ E | u , v ∈ V ∧ ∀ j ∈ D , j ≠ d : ( u , v , j ) ∉ E } | | { ( u , v , d ) ∈ E | u , v ∈ V } | {\displaystyle {\frac {|\{(u,v,d)\in E|u,v\in V\wedge \forall j\in D,j\neq d:(u,v,j)\notin E\}|}{|\{(u,v,d)\in E|u,v\in V\}|}}} .^[39]

Burstiness is a well-known phenomenon in many real-world networks, e.g. email or other human communication networks. Additional dimensions of communication provide a more faithful representation of reality and may highlight these patterns or diminish them. Therefore, it is of critical importance that our methods for detecting bursty behavior in networks accommodate multidimensional networks.^[51]

Diffusion processes are widely used in physics to explore physical systems, as well as in other disciplines as social sciences, neuroscience, urban and international transportation or finance. Recently, simple and more complex diffusive processes have been generalized to multilayer networks.^[23][52] One result common to many studies is that diffusion in multiplex networks, a special type of multilayer system, exhibits two regimes: 1) the weight of inter-layer links, connecting layers each other, is not high enough and the multiplex system behaves like two (or more) uncoupled networks; 2) the weight of inter-layer links is high enough that layers are coupled each other, raising unexpected physical phenomena.^[23] It has been shown that there is an abrupt transition between these two regimes.^[53]

In fact, all network descriptors depending on some diffusive process, from centrality measures to community detection, are affected by the layer-layer coupling. For instance, in the case of community detection, low coupling (where information from each layer separately is more relevant than the overall structure) favors clusters within layers, whereas high coupling (where information from all layer simultaneously is more relevant than the each layer separately) favors cross-layer clusters.^[8][9]

As for unidimensional networks, it is possible to define random walks on the top of multilayer systems. However, given the underlying multilayer structure, random walkers are not limited to move from one node to another within the same layer (jump), but are also allowed to move across layers (switch).^[15]

Random walks can be used to explore a multilayer system with the ultimate goal to unravel its mesoscale organization, i.e. to partition it in communities,^[8][9] and have been recently used to better understand navigability of multilayer networks and their resilience to random failures,^[15] as well as for exploring efficiently this type of topologies.^[54]

In the case of interconnected multilayer systems, the probability to move from a node i {\displaystyle i} in layer α {\displaystyle \alpha } to node j {\displaystyle j} in layer β {\displaystyle \beta } can be encoded into the rank-4 transition tensor T j β i α {\displaystyle T_{j\beta }^{i\alpha }} and the discrete-time walk can be described by the master equation

p j β ( t + 1 ) = ∑ α = 1 L ∑ i = 1 N T j β i α p i α ( t ) = ∑ α = 1 L ∑ i = 1 N ( T t ) j β i α p i α ( 0 ) {\displaystyle p_{j\beta }(t+1)=\sum _{\alpha =1}^{L}\sum _{i=1}^{N}T_{j\beta }^{i\alpha }p_{i\alpha }(t)=\sum _{\alpha =1}^{L}\sum _{i=1}^{N}(T^{t})_{j\beta }^{i\alpha }p_{i\alpha }(0)}

where p i α ( t ) {\displaystyle p_{i\alpha }(t)} indicates the probability of finding the walker in node i {\displaystyle i} in layer α {\displaystyle \alpha } at time t {\displaystyle t} .^[3][15]

There are many different types of walks that can be encoded into the transition tensor T j β i α {\displaystyle T_{j\beta }^{i\alpha }} , depending on how the walkers are allowed to jump and switch. For instance, the walker might either jump or switch in a single time step without distinguishing between inter- and intra-layer links (classical random walk), or it can choose either to stay in the current layer and jump, or to switch layer and then jump to another node in the same time step (physical random walk). More complicated rules, corresponding to specific problems to solve, can be found in the literature.^[23] In some cases, it is possible to find, analytically, the stationary solution of the master equation.^[15][54]

The problem of classical diffusion in complex networks is to understand how a quantity will flow through the system and how much time it will take to reach the stationary state. Classical diffusion in multiplex networks has been recently studied by introducing the concept of supra-adjacency matrix,^[55] later recognized as a special flattening of the multilayer adjacency tensor.^[3] In tensorial notation, the diffusion equation on the top of a general multilayer system can be written, concisely, as

d X j β ( t ) d t = − L j β i α X i α ( t ) {\displaystyle {\frac {dX_{j\beta }(t)}{dt}}=-L_{j\beta }^{i\alpha }X_{i\alpha }(t)}

where X i α ( t ) {\displaystyle X_{i\alpha }(t)} is the amount of diffusing quantity at time t {\displaystyle t} in node i {\displaystyle i} in layer α {\displaystyle \alpha } . The rank-4 tensor governing the equation is the Laplacian tensor, generalizing the combinatorial Laplacian matrix of unidimensional networks. It is worth remarking that in non-tensorial notation, the equation takes a more complicated form.

Many of the properties of this diffusion process are completely understood in terms of the second smallest eigenvalue of the Laplacian tensor. It is interesting that diffusion in a multiplex system can be faster than diffusion in each layer separately, or in their aggregation, provided that certain spectral properties are satisfied.^[55]

Recently, how information (or diseases) spread through a multilayer system has been the subject of intense research.^{[56][1] [57][58][59]}

Several software programs focusing on the analysis and visualization of multilayer networks have been introduced. Some popular solutions include multinet (C++ / Python / R), MuxViz (R), Pymnet (Python), with each software typically specializing in different analytical functions.^[60] However, most software currently face issues such as processing very large multilayer networks, while the interoperability between software also needs improvement.