Created: 10 Mar 2014  Modified: 23 Jun 2017  BibTeX Entry  RIS Citation 
The goal is to understand how nauty
output and symmetry metrics behave when we consider a forest with multiple disconnected trees. This situation would occur if we calculated symmetries across the entire set of trait trees when capturing statistics.
Some statistics, such as orbits, might crosscut trait trees if two trait trees shared adjacency relationships. In other situations, I’d expect them to be a disjoint union of results, and thus the statistics might simply be additive (and thus needing to be normalized by the number of configured trait trees). In other cases, such as \(\textit{Aut}(G)\), the relationship might be multiplicative and combinatorial.
I proceed by constructing a union graph of 2 balanced trait trees, then 3 and 4, and examining how the nauty
output statistics scale. Then I look at a combination of different balanced trees. In the code sections below, I rely upon the axelrodct simulation model available from my Github repository, Brendan McKay’s nauty
library, and the Python networkx
module.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 

Below is the output for a forest with two trees with \(r=4, h=4\). We can see that the orbit number is the same, the multiplicities are simply different because the members of each orbit in each tree are interchangeable, so each multiplicity is twice the value in a single tree. The group size is much larger than twice, however, but we should expect there to be a multiplicative effect given that automorphisms scale as \(n!\) for the complete graphs \(K_n\).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 

Adding a third balanced tree follows the same pattern. The only concern is the scaling of group size. Adding a fourth tree to give me some data points to work with.
1 2 3 4 5 6 7 8 9 10 11 12 

The exponent is definitely going up as the power of the number of trees: \({10^{117}}^2 = 10^{234}\), and so on. But the other values are rising as well, and not by a simple multiplicative factor (which looked like 2). It looks like \(n!\).
So, with a forest of identical trees, with n trees, and \(s = \textit{Aut}(g)\) for each tree, the total automorphism group size for the entire forest is thus \(n! s^n\). This is true for identical tree copies, but obviously more complex for nonidentical components. We would not also see orbit multiplicities which are a simple factor of the base orbit multiplicity, and we’d expect that the orbit number would not be constant, of course.
This turns out to be a wellknown result in graph theory, and applicable to any set of repeated graphs, not just trees. Frucht (1949) gave the theorem but no explicit proof. But it’s easy to see why this is true. Take a graph \(G\) with \(n\) copies of a connected component, \(H\). For each of the \(n\) copies of \(H\), there are \(s = \textit{Aut}(H)\) automorphisms which map vertices within each copy. Thus there are \(s^n\) total automorphisms which perform mappings within a single copy of \(H\). Then, we consider automorphisms between the identical copies, and find that there are \(n!\) ways of permuting each automorphism with vertices from the \(n\) copies. The total is thus \(n! s^n\). Another way to think of this is that we are doing the cartesian product of the symmetric group \(S_n\) which gives all permutations of \(n\) elements, with the direct product of all the copies. I believe this is called an “external wreath product”.
So how do the orbit numbers etc behave with multiple nonidentical components?
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 

So…since the graph components are different, the orbits are disjoint and restricted to their own component. The total orbit number is simply the sum of the orbit numbers for the two components, and the orbit multiplicites and vertex membership are obviously restricted to each component (easy to see given the sequential vertex numbering).
The group size is greater than that of a single tree, but many orders of magnitude smaller than the combined group size of two identical trees. This likely reflects there being very few automorphisms between trees. In fact, the total group size is simply the product of the individual components: \(s = \textit{Aut}(G_1)\textit{Aut}(G_2)\ldots\textit{Aut}(G_n)\). This means that the automorphism group of nonidentical graph components taken together as a single graph with multiple components is the “direct product” of the individual components. What’s missing from the simple direct product is the interaction across components, which I suspect is what makes the above calculation a “wreath product,” although I don’t fully understand the latter yet.
Frucht, Robert. 1949. “On the Groups of Repeated Graphs.” Bulletin of the American Mathematical Society 55 (4): 418–20.