1. Introduction
Symmetry is a principle which has served as a guide for the spectacular advances that have been made in modern science, especially physics. For example, the continuous translational symmetry of ordinary space and time guarantees the invariance of the laws of physics under such translations. Thus, any mathematical expression describing a physical system, whether subatomic or macroscopic, must be invariant under space and time translations.
Group theory is the mathematical language used to describe symmetry and its associated invariant properties (recall that an abstract group is a set of elements together with a law of composition such that for (i) ; (ii) ; (iii) there is an identity element such that ; and (iv) for there is an inverse such that ). As a simple example, the set of and rotations in the plane of a square about its fixed center under “composition of rotations” form a symmetry group for the square ( is the identity element and the inverse of an rotation is a rotation). Each of these rotations is a symmetry which brings the square into coincidence with itself, i.e., they preserve the invariant shape of the square. A much more complicated example are the so called gauge symmetries of the standard model of physics which classify and describe three fundamental forces of nature (i.e., the electromagnetic, weak, and strong forces) in terms of groups (specifically, the unitary group of degree and the special unitary groups and of degree and , respectively).
In recent years, the notion of generalized symmetry has been introduced to further describe graph symmetry [
1,
2]. The generalized symmetries of a graph are a generalization of the notion of the automorphism group of a graph and are derived from the application of Green’s equivalence relations to the endomorphism monoid of the graph (the automorphism group is a subgroup of the graph’s endomorphism monoid). Since these symmetries and invariant properties are strictly associated with a single graph, they do not address properties that remain fixed when the connection topology of the graph changes.
An important problem in network theory is identifying those properties of networks that remain fixed (invariant) as the network’s connection topology changes with time. It was shown in [
3] that the set of all networks (i.e., all connection topologies) on a fixed number of nodes also forms a semigroup. There it was also shown that the application of Green’s equivalence relations to this semigroup partitions the associated set of networks into equivalence classes, each of which contains many fixed node number networks with various connection topologies, such that all networks within each class share some identifiable invariant connectivity property. If the connection topology of a network changes such that its initial and final configurations are in the same equivalence class, then the initial and final configurations share a common invariant property. It follows that, in this context, Green’s equivalence classifications can be useful for identifying invariant properties of networks which evolve within an equivalence class. Such connectivity invariants can be used, for example, to identify important actors in evolving social networks and to select communication network reconfigurations that will retain a desired connectivity between specific node sets.
Transformations between networks within an equivalence class which preserve the associated invariant connectivity properties are called “Green’s symmetries”. Here, in addition to reviewing the Green’s classification of networks [
3], the “Green’s symmetry problem” is introduced and defined. This problem is to determine (by calculation) the ensemble set of all the Green’s symmetries which evolve an initial network configuration into a final configuration within a fixed Green’s
equivalence class or within a fixed Green’s
equivalence class. As discussed below, each such symmetry encodes information about the internal dynamics of the evolution, i.e., how node neighborhoods in the initial network configuration are joined to form node neighborhoods in the final configuration such that the invariant properties are preserved.
Since the cardinality of such ensembles can be large, the statistical notion of propensity is introduced. This quantity provides measures of the overall tendency of node neighborhoods in an initial network configuration to associate and form node neighborhoods in the final network configuration. Propensities are used to define “propensity energies”, which quantify the overall complexity of the internal dynamics of a network evolution, and “energies of evolution”, which quantify the complexity of internal dynamical activity for an evolution produced by a specific ensemble symmetry.
The objective of this paper is to motivate the application of Green’s symmetry principles to network science by demonstrating how Green’s equivalence relations can be applied to: classify networks; identify associated structural invariants; determine symmetries that preserve these invariants; and define associated measures that quantify aspects of the internal dynamics of network evolutions. The remainder of this paper is organized as follows: To make this paper reasonably self-contained, the relevant definitions and terminology from semigroup theory are summarized in the next section (for additional depth and clarification the reader is invited to consult such standard references as [
4,
5]). The semigroup
of all binary relations on a finite set
and the semigroup
of
Boolean matrices are defined and shown to be isomorphic to one another in
Section 3. The semigroup of networks
on a fixed set
of nodes is introduced and is shown to be isomorphic to
in
Section 4. This isomorphism provides for the Green’s equivalence classifications of
given in
Section 5. Green’s evolutions of networks and their associated invariant properties are discussed in
Section 6. The “Green’s symmetry problem” is defined in
Section 7 and its satisfiability and computational complexity are discussed in
Section 8. The information encoded in symmetries as internal dynamics is detailed in
Section 9. Symmetry “ensembles” and their “propensities” and “energies” are introduced in
Section 10. A simple example illustrating aspects of the theory is presented in
Section 11. Concluding remarks comprise the final section of this paper.
2. Semigroups
A semigroup is a set and an associative binary operation called multiplication defined upon the set (contrast this with the above definition of a group and note that a group is a semigroup endowed with the additional special properties given by items (iii) and (iv)). The one-sided right (one-sided left) multiplication of by is the product . An element is an identity if for . An identity can be adjoined to by setting and defining for . Semigroup and the semigroup on set with associative binary operation are isomorphic (denoted ) when there is a bijective map (i.e., an isomorphism) such that for all .
The well-known and Green’s equivalence relations on a semigroup partition into a highly organized “egg box” structure using their relatively simple algebraic properties. In particular, the equivalence relation on is defined by the rule that if and only if for and the equivalence relation is similarly defined so that if and only if and . The relations and commute under the composition of binary relations and is the smallest equivalence relation containing and .
For and denote the class containing by where or when or , respectively. Thus, if and only if . If and , then there exist elements in such that (hereafter the juxtaposition will also be used for the multiplication ).
3. The Semigroups and
The semigroup
of Boolean matrices is the set of all
matrices over
with Boolean composition
defined by
as the semigroup multiplication operation. Here
, where
is a counting number,
denotes Boolean multiplication (i.e.,
), and
denotes Boolean addition (i.e.,
).
The rows (columns) of any are Boolean row (column) , vectors, i.e., row (column) , tuples over , and come from the set of all Boolean row (column) -vectors. These vectors can be added coordinate-wise using Boolean addition. If , then when the ith coordinate implies the ith coordinate , ( is a partial order).
Let
be either the zero (unit) row or zero (unit) column vector (the context in which
is used defines whether it is a row or column vector). The matrix with
in every row, i.e., the zero matrix, is denoted by
and the matrix with
in every row is denoted by
. For
, the row space
of
is the subset of
consisting of
and all possible Boolean sums of (one or more) nonzero rows of
.
is a lattice
under the partial order
. The row (column)
basis (
) of
is the set of all row (column) vectors in
that are not Boolean sums of other row (column) vectors in
. Please note that each vector in
must be a row (column) vector of
. The vector
is never a basis vector and the empty set
is the basis for the
matrix [
6,
7].
The semigroup
of binary relations on a set
of cardinality
(denoted
) is the power set of
with multiplication
being the “composition of binary relations” defined by
It is easy to see that a bijective index map
induces an isomorphism
defined by
, where
if
and is
if
.
is therefore the Boolean matrix representation of
[
8].
4. The Semigroup
A network of order is the pair , where is a nonempty set of nodes with , and the binary relation is the set of directed links connecting the nodes of the network. Thus, is both a digraph and a binary relation. If , then node is an in(out)-neighbor of node . The in-neighborhood of is the set of all in-neighbors of and the out-neighborhood of is the set of all out-neighbors of .
Let
be the set of networks on
and define “multiplication of networks” by
, where
,
, and
Lemma 1. is a semigroup that is isomorphic to.
Proof. The operation “multiplication of networks” is the same as the operation “composition of binary relations”. Since it is clearly an associative binary operation on , then is a semigroup under the operation “multiplication of networks”. Also, the bijective map defined by preserves multiplication. Thus, is a semigroup isomorphism and . ☐
Lemma 2. If, then.
Proof. This follows from the facts that
(Lemma 1) and
[
8]. ☐
Thus, is also a Boolean matrix representation of .
5. Green’s Equivalence Classifications of
Let
be the isomorphism of Lemma 2 and
be an associated index bijection. If
is the
ith Boolean row vector and
is the
jth Boolean column vector in the matrix
corresponding to network
, then
encodes the out-neighbors of node
in
as the set
and
encodes the in-neighbors of node
in
as the set
When and , then is a basis out-neighborhood and is a basis in-neighborhood for network . Thus, a basis neighborhood in is a nonempty neighborhood in which is not the set union of other neighborhoods in .
Let be the set of basis out-neighborhoods and be the set of basis in-neighborhoods in network . Also, define as the set whose elements are and the sets generated by the closure under set union of the out-neighborhoods in and let be the poset ordered by the set inclusion relation . Thus, when , it may be formally stated that:
Lemma 3. is a lattice that is isomorphic to.
Proof. The proof for this Lemma is the same as that given as the proof of Lemma 3.3 in [
3]. ☐
In what follows, will be referred to as the lattice for .
The following major theorem provides complete and equivalence classifications of all fixed-order networks:
Theorem 1. Let. Then
- i.
if and only if
- ii.
if and only if
- iii.
if and only ifand;
- iv.
if and only ifandare lattice isomorphic.
Proof. The proof of this result is the same as the proof of Theorem 3.4 in [
3]. ☐
Thus, the Green’s
and
equivalence classifications of the networks in
depend entirely upon their having (generally distinct) nodes with identical out-neighborhoods, identical in-neighborhoods, and both identical out-neighborhoods and in-neighborhoods, respectively, whereas the
equivalence classification of networks in
depends entirely upon their having isomorphic
lattices which are generated by their out-neighborhoods. As an illustration of this theorem the reader is invited to consult the simple example given in [
3] which corresponds to the complete Green’s equivalence classification of (and the associated “egg box” structure for) all order two networks.
6. Green’s Evolutions of Fixed-Order Networks
For , let denote the evolution of a network during a time interval , where is the initial network at and is the final network at . If , then the evolution is a Green’s evolution. It is important to note that since and , then and evolutions are also evolutions, whereas evolutions are both and evolutions, as well as evolutions.
Theorem 2. The following statements are true for network evolutions in:
- i.
evolutions preserve basis out-neighborhood sets andlattice isomorphism;
- ii.
evolutions preserve basis in-neighborhood sets andlattice isomorphism;
- iii.
evolutions preserve basis out-neighborhood and in-neighborhood sets andlattice isomorphism;
- iv.
evolutions preservelattice isomorphism.
Proof. This is a direct and obvious consequence of the definitions of Green’s evolutions and Theorem 1. ☐
To illustrate this theorem, consider the order two networks
and
in the example in [
3], where
and
. As can be seen from the associated Green’s equivalence classification performed there, since
and
, the evolution
is both a Green’s
evolution and a Green’s
evolution. Theorem 2 (
i) is satisfied, since, from Table 1 and the discussion in [
3], it is also seen that
and that the
lattices are isomorphic undirected paths of length
.
7. The Green’s Symmetry Problem
In general, a symmetry associated with a “situation” is defined as an “immunity to change” for some aspect of the “situation”. For a “situation” to have a symmetry: (a) the aspect of the “situation” remains unchanged when a change is performed; and (b) it must be possible to perform the change, although the change does not actually have to be performed [
9].
Recall from
Section 2 that for an
evolution
in
, there exists at least one
such that
). Although
(
) does not have to be applied to
, it can produce the desired evolution when applied as a right (left) multiplication of
. In so doing, this multiplication not only preserves
, but also
’s
lattice structure. Thus, (a) and (b) above are satisfied and both
and the associated
lattice structure can be considered as the invariant properties associated with the
symmetries which produce the evolution. Symmetries such as
are
Green’s symmetries.
The “Green’s symmetry problem” is defined here as the determination of all symmetries that produce an evolution from an initial to a final network within an or an class such that each symmetry preserves the structural invariants required by Theorem 2. As will be discussed below, such symmetries encode information about which node neighborhoods in the initial network can be joined to form neighborhoods in the final network such that the structural invariants required by the evolution are preserved.
8. Satisfiability and Computational Complexity of the Green’s Symmetry Problem
The Green’s symmetry problem for an evolution is satisfiable if there are symmetries which can produce the evolution.
Theorem 3. The Green’s symmetry problem is at leastsatisfiable for both Green’sandevolutions.
Proof. Semigroup theory guarantees the existence of at least one Green’s symmetry in that can produce a Green’s evolution and at least one Green’s symmetry in that can produce a Green’s evolution. ☐
8.1. Green’s Evolutions
The isomorphism established in Lemma 2 provides for computational solutions to the Green’s symmetry problem. In particular, if
is a Green’s
evolution, then, since
and
are known, the equation
can be solved for
for each
using the disjunctive normal form logical expression
where use is now made of the Boolean matrix representations of
and
. This expression for fixed
and all
defines a system of
equations for node
.
This system of equations is column- satisfied if there exists a column vector for which (1) is a true statement for each . For each , let be the set of all for which the associated system of equations is satisfied and define . Clearly, if , then is column- satisfied for each and the evolution is -satisfiable. Each instantiation of is represented by a Boolean matrix in which has an as its jth column.
Let index the unit valued entries in the row vector .
Lemma 4. Letfor someand. Ifcolumn-satisfies, thenhaswhen.
Proof. Assume for some that column- satisfies . If and for some , then (1) is true and zero valued for and that value, and the following implication chain is valid: However, since column- satisfies , it must also satisfy (1) for all has when . ☐
Corollary 1. If, then
Proof. ☐
The computational complexity of the Green’s symmetry problem for Green’s evolutions is the number of remaining combinations of values which must be checked for satisfiability after the assignments specified by Lemma 4 have been made. Assume that and for each let index the zero valued Boolean equations of form (1).
Theorem 4. Proof. For each , the set (which can possibly be empty) indexes all row locations in for which in every that column- satisfies . The set indexes all for which must be evaluated to determine the column- satisfiability of an associated . Since there are such evaluations for each , then for all there are a total of evaluations required to determine all which column- satisfy . ☐
8.2. Green’s Evolutions
If
is a Green’s
evolution, then, since
, it can be solved for
for each
using the disjunctive normal form logical expression
which, for fixed
and all
, defines a system of
equations for node
. This system is
row- satisfied if there exists a row vector
for which (2) is a true statement for each
. For each
, let
be the set of all
for which the associated system of equations is row-
satisfied and define
. If
, then
is row-
satisfied for each
and
the evolution is satisfiable. Each instantiation of
is represented by a Boolean matrix in
which has a
as its
ith row.
Let index the unit valued entrees in the column vector .
Lemma 5. Letfor someand. Ifrow-i satisfies, thenhaswhen.
Proof. Assume for some that row, satisfies . If for some and , then (2) is true and zero valued for and that value, and the following implication chain is valid: However, since row- satisfies , it must also satisfy (2) for all has when . ☐
Corollary 2. If, then.
Proof. . ☐
The computational complexity of the Green’s symmetry problem for Green’s evolutions is the number of remaining combinations of values which must be checked for satisfiability after the assignments specified by Lemma 5 have been made. Assume that and for each let index the zero valued Boolean equations of form (2).
Theorem 5. Proof. For each , the set (which can possibly be empty) indexes all column locations for which in every that row- satisfies . The set indexes all for which must be evaluated to determine the row- satisfiability of an associated . Since there are such evaluations for each , then for all there are a total of evaluations required in order to determine all which row- satisfy . ☐
9. Symmetries: Instantiations of Internal Dynamics
Since Green’s symmetries are themselves effectively elements of , they correspond to special binary relations between network nodes that encode aspects of the internal dynamics of a Green’s evolution . In particular, they generally identify many-to-one correspondences between neighborhood sets in that are joined by set union to produce a neighborhood in . Each of these correspondences occurs in such a way as to preserve the structural invariants required by Theorem 2. These correspondences are the internal dynamics of the evolution.
Consider a Green’s
evolution
where each symmetry
satisfies
and is one instantiation of a possible set of symmetries which produce the evolution and preserve the required invariants. If
is a column in
with a
in each of the rows in the set
and zeros in every other row location (i.e., there are
1’s and
0’s), then this column encodes an internal dynamic of the evolution where the in-neighborhoods of nodes
in
are joined together as
and associated with the in-neighborhood
in
according to
This expression is called a
internal dynamic of the evolution and the set
is
the associated motion of the dynamic. Clearly, for the special case where
,
If
is a Green’s
evolution, a symmetry
which produces the invariant preserving evolution satisfies
. If
is a row in
with a 1 in each of the column locations in
, then this row encodes an internal dynamic of the evolution where the out-neighborhoods of nodes
in network
are joined by set union and associated with the out-neighborhood
in network
according to
This expression is a
internal
dynamic of the evolution and the set
is the associated motion of the dynamic. When
, then
These notions will be clarified below using a simple example.
10. Symmetry Ensembles, Propensities, and Energies
Since the symmetry which produces a Green’s evolution is not necessarily unique, it can be unclear as to how to assign a specific symmetry to an evolution. However, the collection of symmetries obtained from Green’s symmetry problem, i.e., the symmetry ensembles, can be used to construct propensities. Propensities can be viewed as weighted symmetries which, in some sense, represent their respective ensembles.
Let
index the symmetries which are solutions to the Green’s symmetry problem for some Green’s
(
) evolution
. The sets
and
are the associated symmetry ensembles. The propensities associated with each ensemble are defined as
and
Thus, is a measure of the tendency of the nodes in column in network to form motions that associate in-neighborhoods in with in-neighborhoods in network according to the internal dynamic (3). Similarly, is a measure of the tendency of nodes in row in to form motions that associate out-neighborhoods in with out-neighborhoods in according to the internal dynamic (4).
Propensities can be used to associate energies with both ensembles and specific symmetries. These energies quantify in a directly proportional manner the complexity level of the internal dynamical activity that is associated with an evolution. The
propensity energies provide a representative measure of the “overall” complexity of internal dynamical activity for an evolution based upon ensemble propensity. The propensity energies for ensembles
and
are defined as
and
respectively.
The
energies of evolution for the specific symmetries in an ensemble measure the complexity of internal dynamical activity for an evolution produced by a specific symmetry in an ensemble. In particular, if
and
, then the associated energies of evolution are defined as
and
The following Lemma guarantees that the energy of evolution for a symmetry never exceeds the propensity energy for the associated ensemble.
Lemma 6. For any Green’sorevolution,, whereorwhenor.
Proof. . ☐
Recall that internal
and
dynamics are strictly defined by their motions. These motions also have energies that provide a measure of the level of internal dynamical activity induced by the motion. Since the symmetries
and
encode
and
internal dynamics with motions
and
, respectively, then the associated
energies of motion are the quantities
and
The energies of motion are related to their energies of evolution by the following theorem:
Theorem 6 (Conservation of Energy of Evolution). The energy of evolution of a Green’s symmetry is conserved by the energies of motion of its internal dynamics.
Proof. Let and set index all the internal dynamics encoded by . Then , where use has been made of the fact that is equivalent to because when and . It is similar for the dynamics. ☐
11. Example
Let
be a Green’s
evolution in
(or equivalently in
), where (in
)
with
(note that this evolution corresponds to the
Green’s
evolution in [
3]). Theorem 3 guarantees the existence of at least one
such that
The disjunctive normal form logical expression (1) for this equation yields the following system of equations
which can be used to solve the associated Green’s symmetry problem.
For the two equations in the second row of this system to be satisfied requires the assignment
. By inspection it is seen that the complete system is satisfied when, in addition to these assignments,
and
each assume both values from the set
. Thus,
so that
and the evolution
is 4-satisfiable. The associated symmetry ensemble is the set
To calculate the computational complexity of this Green’s symmetry problem, refer to
Section 8.1 and observe that
and
. Application of Theorem 4 yields
, i.e., four combinations of value assignments must be checked for each
since, according to Lemma 4,
values cannot be assigned when
because
.
The propensity and propensity energy for the ensemble are
and
, respectively, and the energies of evolution are
and
Please note that this validates Lemma 6. These energies also indicate that
produces the least energy of evolution in the sense that the evolution involves simpler internal dynamical activity than evolutions produced by the other symmetries in the ensemble.
To illustrate this further, first observe that and (here the jth column vector is set directly equal to the nodes in the in-neighborhood of node ). It is also easily determined that the motions of the dynamics for: are ; are and ; are and ; and are . By inspection it is found that each of these motions satisfies (3). Using as an example, it is seen that (3) yields the correct set theoretic relationship or or for both and . Also note that the internal dynamics for are simpler than those for the other symmetries in the ensemble, in the sense that both of the motions are singleton sets, whereas at least one of the motions for the other symmetries is a doubleton set. This is consistent with the fact mentioned above that produces the least energy of evolution.
Now consider the energies of motion for each ensemble symmetry. They are easily calculated from the theory and are found to be:
and
Thus, the motions associated with an
evolution are the least energetic since
This is also consistent with the fact that an induced evolution is the least energetic and involves the least complex internal dynamics.
Finally, observe that these results validate Theorem 6. In particular,
and
12. Concluding Remarks
The research documented in [
3] was inspired by earlier research performed by Konieczny [
6] and Plemmons et al. [
7]. This paper has reviewed the results developed in [
3], i.e., that the set of all networks on a fixed number of nodes can be classified using the Green’s equivalence relations of semigroup theory and that all networks within a Green’s equivalence class have a common structural invariant (neighborhoods or poset relationships between node sets generated by neighborhoods). By extension, it was deduced in this paper from these results that if a network evolves from an initial network configuration to a final network configuration such that both the initial and final networks are in the same Green’s equivalence class, then the structural invariants for the class are preserved by the evolution. In addition, the Green’s symmetry problem was also defined in this paper. This problem is to determine by computation all symmetries which produce a network evolution within a Green’s
or a Green’s
equivalence class (i.e., a symmetry ensemble). These symmetries were shown to be solutions to special Boolean equations whose form is dictated by semigroup theory. Each such symmetry encodes information about the internal dynamics of the associated evolution and an ensemble associated with an evolution was used to define propensities and energies which quantify aspects of the internal dynamics of the evolution. However, it should be noted that a practical limitation exists for solving the Green’s symmetry problem. This occurs because the cardinality of symmetry ensembles associated with large real networks can be quite large, thereby requiring the use of considerable computational resources to solve such problems (see future research suggestions below).
In conclusion, it is believed that the results of this paper are new and not in general use (perhaps having the closest resemblance to these results are the applications of Green’s relations to social networks [
10] and automata theory, e.g., [
11]). However, the results of this paper are important and should be of general interest to network science researchers and those working in areas of applied network theory. In addition to applications similar to those mentioned in
Section 1 (actor identification in social networks and communication network reconfiguration), contemporary areas of frontier research, such as identifying emerging scientific disciplines, e.g., [
12], analyzing brain connectivity, e.g., [
13,
14,
15], and finding symmetries in engineering processes [
16], could also benefit from the results of this paper.
Before closing it is worthwhile to mention several directions for related future research. First, because of the computational resources required to solve the Green’s symmetry problem, it would be useful to investigate how sampling and statistics can be used to obtain symmetry sub-ensembles that effectively yield the same information about propensities and energies as the associated full ensemble. A second research area involves understanding symmetries and their computation for network evolutions occurring within Green’s
and
equivalence classes. A third and potentially very interesting research area concerns determining the relationships (if any) between the theory developed in this paper and the relatively new theory of persistence that is used to analyze large data sets, e.g., [
17].