\section{Introduction}

Cayley graphs have been well-studied as a topology for interconnection networks due to their high symmetry, low diameter, high connectivity, and embeddable properties.  Perhaps the most well-studied topology is the hypercube graph, for which there are many equivalent definitions.  One definition is to view the hypercube as a Cayley graph.

Given a group $H$ and a subset $S \subseteq H$, the Cayley graph of $H$ with respect to $S$, denoted $\Cay(H,S)$, is the graph whose vertex set is $H$, and with arc set $\{ (h, sh): h \in H, s \in S \}$.  If the identity group element $e$ is not in $H$, then $\Cay(H,S)$ has no self-loops, and if $H$ is closed under inverses, then $(h,g)$ is an arc iff $(g,h)$ is an arc; in this  case, the Cayley graph can be viewed as an undirected graph.  The Cayley graph $\Cay(H,S)$ is connected iff $S$ generates $H$.

Let $\mathbb{Z}_2^n$ denote the direct product of $n$ copies of the cyclic group of order 2.  Let $e_i$ denote the unit vector in the vector space $\mathbb{Z}_2^n$ that has a 1 in the $i$th coordinate and zero in the remaining coordinates. The hypercube graph is defined to be the Cayley graph of $\mathbb{Z}_2^n$ with respect to the set of generators $S = \{e_1,e_2, \ldots, e_n\}$.

The hypercube graph has been well-studied in the literature as a topology for interconnection networks.  In order to improve its connectivity and fault-tolerance properties, the folded hypercube graph was proposed \cite{ElAmawy:Latifi:1991}.   The folded hypercube graph is obtained by taking the hypercube graph and joining each vertex $x$ to the unique vertex $\overline{x}$ that is farthest away from $x$ in the hypercube.  Thus, the folded hypercube is the Cayley graph of $\mathbb{Z}_2^n$ with respect to the set of generators $S = \{e_1,e_2,\ldots,e_n, 1\cdots11\}$.  The addition of these edges to the hypercube graph yields a new graph whose diameter is smaller.

The augmented cube graph was introduced in Choudum and Sunitha \cite{Choudum:Sunitha:2002} as a topology for consideration in interconnection networks.  The augmented cube graph of dimension $n$, denoted $AQ_n$, is the Cayley graph of $\mathbb{Z}_2^n$ with respect to the set of $2n-1$ generators $\{e_1,e_2,\ldots,e_n\}$ $ \dot{\cup} \{00\cdots0011, $ $00\cdots0111, 00\cdots1111, 11\cdots1111\}$. Recall that the hypercube graph is the Cayley graph of $\mathbb{Z}_2^n$ with respect to the generator set $\{e_1,\ldots,e_n\}$.  The additional $n-1$ generators in the definition of $AQ_n$ augment the edge set of the hypercube graph, so that the augmented graph has smaller diameter and better connectivity and embeddable properties. For equivalent definitions of the augmented cube graph and a study of some of its properties, see  \cite{Choudum:Sunitha:2002} \cite{Choudum:Sunitha:2008} \cite{Choudum:Sunitha:2001}.

Let $X=(V,E)$ be a simple, undirected graph.  The automorphism group of $X$ is defined to be the set of all permutations of the vertex set that preserves adjacency: $\Aut(X):=\{g \in \Sym(V): E^g = E\}$.    While one can often obtain some automorphisms of a graph, it is often difficult to prove that one has obtained the (full) automorphism group. Automorphism groups of graphs that have been proposed as the topology of interconnection networks have been investigated by several authors; see \cite{Deng:Zhang:2011} \cite{Deng:Zhang:2012} \cite{Ganesan:JACO} \cite{Ganesan:DM:2013}  \cite{Feng:2006} \cite{Mirafzal:2011} \cite{Zhang:Huang:2005} \cite{Zhou:2011}.

Let $G:=\Aut(AQ_n)$ denote the automorphism group of the augmented cube graph.  Since $AQ_n = \Cay(\mathbb{Z}_2^n,S)$ is a Cayley graph, its automorphism group can be expressed as a rotary product $G = \mathbb{Z}_2^n \times_{rot} G_e$, where $G_e$ is the stabilizer in $G$ of the identity vertex $e$ (cf. Jajcay \cite{Jajcay:1994} \cite{Jajcay:2000}).  In  \cite{Choudum:Sunitha:2008} \cite{Choudum:Sunitha:2001}, it was shown that the order of the automorphism group of $AQ_n$ is exactly $2^{n+3}$, for all $n \ge 4$.

In the present paper, the structure of the automorphism group $G:=\Aut(AQ_n)$ is determined to be $G \cong \mathbb{Z}_2^n \rtimes D_8$, for all $n \ge 4$.  Hence, in addition to proving that $G_e \cong D_8$, we also strengthen the rotary product to a semidirect product.  The clique structure of $AQ_4$ is investigated further, and we show that the automorphism group of $AQ_4$ is isomorphic to that of $AQ_3$: $\Aut(AQ_4) \cong \Aut(AQ_3) \cong (D_8 \times D_8) \rtimes C_2$.  We also determine all the nontrivial blocks of $AQ_4$.

An open problem in the literature is to determine, given a group $H$ and a subset $S \subseteq H$, whether $\Cay(H,S)$ is normal \cite{Xu:1998}.  Every Cayley graph $\Cay(H,S)$ admits $R(H) \Aut(H,S)$ as a subgroup of automorphisms; here  $R(H)$ denotes the right regular representation of $H$,  and $\Aut(H,S)$ is the set of automorphisms of $H$ that fixes $S$ setwise (cf. \cite{Biggs:1993} \cite{Godsil:Royle:2001}).  A Cayley graph $\Cay(H,S)$ is said to be normal if its full automorphism group is $R(H) \Aut(H,S)$.  In the present paper, we prove that the Cayley graph $AQ_3$ is non-normal, and we prove that $AQ_n$ is normal for all $n \ge 4$.

\bigskip \noindent \textbf{Notation.} We mention the notation used in this paper. Given a graph $X=(V,E)$ and a vertex $v \in V$, $X_i(v)$ denotes the set of vertices of $X$ whose distance to vertex $v$ is exactly $i$. Thus, if $X=\Cay(H,S)$, then $X_1(e)=S$. The subgraph of $X$ induced by a subset of vertices $W \subseteq V$ is denoted $X[W]$, and so $X[X_i(v)]$ denotes the subgraph of $X$ induced by the $i$th layer of the distance partition of $X$ with respect to vertex $v$.   The dihedral group of order $2n$ is denoted $D_{2n}$.

Given a Cayley graph $X = \Cay(H,S)$, $e$ denotes the identity element of the group $H$ and also denotes the corresponding vertex of $X$.  Throughout this paper, $H$ will be the abelian group $\mathbb{Z}_2^n$.  In this case, $R(H) = R(\mathbb{Z}_2^n)$ is the set $\{\rho_z: z \in \mathbb{Z}_2^n\}$ of translations by $z$, i.e. $\rho_z: \mathbb{Z}_2^n \rightarrow \mathbb{Z}_2^n, x \mapsto x+z$, and $\rho_z$ is an automorphism of the Cayley graph.  We also let $\mathbb{Z}_2^n$ denote the vector space $\mathbb{F}_2^n$.

Let $G:=\Aut(X)$.  Then $G_e$ is the set of automorphisms of $X$ that fixes the vertex $e$, and $L_e$ denotes the set of automorphisms of $X$ that fixes the vertex $e$ and each of its neighbors.

\section{$AQ_2$ and $AQ_3$}

\begin{Lemma}
 Let $AQ_2$ be the augmented cube graph $\Cay(\mathbb{Z}_2^2,S)$, where $S = \{e_1,e_2,11\}$. Then $\Aut(AQ_2) \cong S_4$ and $AQ_2$ is normal.
\end{Lemma}

\noindent \emph{Proof}:
Note that $AQ_2$ is the complete graph on 4 vertices. Hence, $\Aut(AQ_2) \cong S_4$ and the stabilizer in $G:=\Aut(AQ_2)$ of the identity vertex $e$ is $G_e \cong S_3$.  Recall that the Cayley graph $\Cay(\mathbb{Z}_2^2,S)$ is normal iff $G_e \subseteq \Aut(\mathbb{Z}_2^n,S)$ (cf. Xu \cite[Proposition 1.5]{Xu:1998}). The set $\Aut(\mathbb{Z}_2^3,S)$ of automorphisms of the group $\mathbb{Z}_2^2$ that fixes $S$ setwise is exactly the set of invertible linear transformations of the vector space $\mathbb{Z}_2^2$ that fixes $S$ setwise. The elements of $\Aut(\mathbb{Z}_2^2,S)$ are uniquely determined by their action on the basis $\{e_1,e_2\}$ of the vector space. Let $\phi \in \Aut(\mathbb{Z}_2^2,S)$. The image of $\phi(e_1)$ can be chosen from $S$ in 3 ways, and the image $\phi(e_2)$ can be chosen from $S - \{\phi(e_1)\}$ in 2 ways. Hence, $|\Aut(\mathbb{Z}_2^2,S)|=6$, which equals $|G_e|$.  Thus, $AQ_2$ is normal, i.e. it has the smallest possible full automorphism group $R(\mathbb{Z}_2^2) \Aut(\mathbb{Z}_2^2,S) \cong \mathbb{Z}_2^2 \rtimes S_3$.
\qed

\begin{Theorem}
 Let $AQ_3$ be the Cayley graph $\Cay(\mathbb{Z}_2^3,S)$, where $S:=\{e_1,e_2,e_3,011,111\}$. Let $G:=\Aut(AQ_3)$. Then $G_e \cong D_8 \times C_2$, $\Aut(\mathbb{Z}_2^3,S) \cong D_8$, and $AQ_3$ is non-normal.
\end{Theorem}

\noindent \emph{Proof}:
Let $X:=AQ_3$ and $G:=\Aut(AQ_3)$.  We show that the stabilizer in $G$ of the vertex $e$ is $G_e \cong D_8 \times C_2$ and that $\Aut(\mathbb{Z}_2^3,S) \cong D_8$, whence $\Aut(\mathbb{Z}_2^3,S) \ne G_e$ and $AQ_3$ is non-normal.

The distance partition of $AQ_3$ with respect to the identity vertex $e$ is shown in Figure~\ref{fig:distance:partition:AQ3}.  The vertex $e$ is adjacent to each vertex in $X_1(e)$ and each $u \in X_1(e) - \{011\}$ is adjacent to both vertices in $X_2(e)$ (these edges between the layers are not shown in the figure). The graph $AQ_3$ is 5-regular.

\begin{figure}
\centering
\begin{tikzpicture}[scale=1.5]
\vertex[fill] (e) at (0,0) [label=below:$e$] {};

\vertex[fill] (v100) at (2,-1) [label=below:$111$] {};
\vertex[fill] (v111) at (2,1) [label=above:$e_1$] {};

\vertex[fill] (v011) at (3,0) [label=right:$011$] {};

\vertex[fill] (v010) at (4,-1) [label=below:$e_2$] {};
\vertex[fill] (v001) at (4,1) [label=above:$e_3$] {};

\vertex[fill] (v110) at (6,-1) [label=below:$101$] {};
\vertex[fill] (v101) at (6,1) [label=above:$110$] {};

\node at (3,-2) {${X_1(e)}$};
\node at (6,-2) {${X_2(e)}$};
\path
    (v100) edge (v011)
    (v111) edge (v011)
    (v010) edge (v011)
    (v001) edge (v011)

    (v100) edge (v111)
    (v001) edge (v010)

    (v110) edge (v101)
;
\end{tikzpicture}
\caption{Distance partition of $AQ_3$ (edges between different layers are not drawn).}
\label{fig:distance:partition:AQ3}
\end{figure}

The subgraph of $X$ induced by the $i$th layer of the distance partition is denoted $X[X_i(e)]$ .  The stabilizer $G_e$ acts on the subgraph $X[X_i(e)]$ as a group of automorphisms. Observe from Figure~\ref{fig:distance:partition:AQ3} that an automorphism $g \in G_e$ can act on $X_1(e)$ by permuting the endpoints of the edge $\{e_1,111\}$, by independently permuting the endpoints of the edge $\{e_2,e_3\}$, and by thereafter possibly also interchanging these two edges.  Hence, the automorphism group of the induced subgraph $X[X_1(e)]$ is that of two independent edges (or its complement graph - the square) and hence is isomorphic to $D_8$.

The automorphism group of the induced subgraph $X[X_2(e)]$ is clearly $C_2$. Since each vertex in $X_2(e) - \{011\}$ is adjacent to all vertices in $X_1(e)$ and since $G_e$ fixes the vertex $011$, $G_e  \cong D_8 \times C_2$.  This also proves that $\Aut(AQ_3) = R(\mathbb{Z}_2^3) \times_{rot} (D_8 \times C_2)$.

We now determine $\Aut(\mathbb{Z}_2^3,S)$. Recall that the automorphisms of the group $\mathbb{Z}_2^n$ are precisely the invertible linear transformations of the vector space $\mathbb{Z}_2^n$ (cf. \cite[p. 136]{Dummit:Foote:2004}).  Hence $\Aut(\mathbb{Z}_2^3,S)$ is precisely the set of elements of $G_e$ that can be realized as linear transformations of the vector space $\mathbb{Z}_2^3$.  We first show that $\Aut(X[X_1(e)]) \cong \Aut(\mathbb{Z}_2^3,S)$.
Each automorphism $\phi$ of the induced subgraph $X[X_1(e)]$ also specifies the image $\phi(e_i)$ of the vertex $e_i~~ (i=1,2,3)$. The resulting image of the basis $\{e_1,e_2,e_3\}$ extends uniquely to a linear transformation $\phi'$ of the vector space $\mathbb{Z}_2^3$. Let $f: \Aut(X[X_1(e)]) \rightarrow \Aut(\mathbb{Z}_2^3,X_1(e)), \phi \mapsto \phi'$ be the map that takes an automorphism of the induced subgraph $X[X_1(e)]$ to its linear extension $\phi'$.  We show that the map $f$ is well-defined and a bijective homomorphism.

Let $\phi \in \Aut(X[X_1(e)])$. Observe from Figure~\ref{fig:distance:partition:AQ3} that $\phi$ fixes the vertex $011$ and permutes the vertices $X_1(e)-\{011\} = \{e_1,e_2,e_3,111\}$ among themselves.  In particular, if $\phi$ takes the basis element $e_i$ to $\phi(e_i)$ $(i=1,2,3)$, then $\phi$ must take the fourth vertex $111=e_1+e_2+e_3$ to the remaining vertex (i.e. to the unique vertex in the singleton set $\{e_1,e_2,e_3,111\} - \{\phi(e_1), \phi(e_2), \phi(e_3) \}$. This remaining vertex is the sum $\phi(e_1)+\phi(e_2)+\phi(e_3)$ because every 3 vectors in the set $\{e_1,e_2,e_3,111\}$ sum to the fourth vector in the set.  Also, the subsets $\{e_2,e_3\}$ and $\{e_1,111\}$ are blocks of $\phi$, and the sum of the two vertices in each subset is $011$.  In other words, if $\phi \in \Aut(X[X_1(e)])$, then $\phi$ acts linearly on the vertices in $X_1(e)$, and so the extension of $\phi$ to a linear transformation $\phi'$ is well-defined.

An automorphism of the induced subgraph $X([X_1(e)])$ must take vertices $e_1,e_2$ and $e_3$ to 3 vertices in $\{e_1,e_2,e_3,111\}$. Any 3 of these 4 vectors are linearly independent, whence the linear extension of $\phi$ is an invertible linear transformation. Thus, the map $f$ is well-defined.

Different automorphisms of $X([X_1(e)])$ induce different images of the basis $\{e_1,e_2,e_3\}$ and hence different linear extensions.  Thus, $f$ is injective. Since $\Aut(\mathbb{Z}_2^3,S) \subseteq G_e$, every element $\phi' \in \Aut(\mathbb{Z}_2^3,S)$, when restricted to $S$, is an automorphism $\phi$ of $X([X_1(e)])$.  The linear extension of $\phi$ is unique and hence equals $\phi'$. Thus, $f$ is surjective. The linear extension map is also a homomorphism.  Hence $f$ is an isomorphism.

Since $\Aut(X[X_1(e)]) \cong D_8$, we have that $\Aut(\mathbb{Z}_2^3,S) \cong D_8$. Thus, $\Aut(\mathbb{Z}_2^3,S)$ is a proper subgroup of $G_e = D_8 \times C_2$ and $AQ_3$ is non-normal.
\qed

\begin{Lemma}
 Let $AQ_3$ denote the augmented cube graph $\Cay(\mathbb{Z}_2^3,S)$ of dimension 3. Then $\Aut(AQ_3) \cong (D_8 \times D_8) \rtimes C_2$.
\end{Lemma}

\noindent \emph{Proof}:
Since $|S|=5$, $AQ_3$ is a 5-regular graph on 8 vertices. Its complement graph $\overline{AQ_3}$ is $\Cay(\mathbb{Z}_2^3,S')$, where $S':=\mathbb{Z}_2^3-S-\{0\} = \{110,101\}$. It is easy to see that $\overline{AQ_3}$ is the union of two disjoint squares. Hence $\Aut(\overline{AQ_3}) \cong (D_8 \times D_8) \rtimes C_2$.
\qed

\section{Automorphism group and normality of $AQ_n$ $(n \ge 4)$}

Let $AQ_n$ be the augmented cube graph $\Cay(\mathbb{Z}_2^n,S)$ $(n \ge 4)$. Let $G:=\Aut(AQ_n)$. Since $AQ_n$ is a Cayley graph, its automorphism group is the rotary product $R(\mathbb{Z}_2^n) \times_{rot} G_e$. It was shown in \cite{Choudum:Sunitha:2001}~\cite{Choudum:Sunitha:2008} that the number of automorphisms of the graph $AQ_n$ that fix the identity vertex $e$ is exactly 8, i.e. $|G_e|=8$. Thus, it was shown that the order of the automorphism group of $AQ_n$ is $2^n \times 8 = 2^{n+3}$.  In the present section, we determine the structure of the automorphism group of $AQ_n$: we show that $G_e \cong D_8$ and that every one of these 8 automorphisms can be realized as a linear transformation of the vector space $\mathbb{Z}_2^n$. It then follows that $G_e \subseteq \Aut(\mathbb{Z}_2^n,S)$ and that the rotary product expression can be strengthened to the semidirect product $G = R(\mathbb{Z}_2^n) \rtimes D_8$.

We first recall the following result from the literature:

\begin{Lemma} \cite{Choudum:Sunitha:2001} \cite{Choudum:Sunitha:2008}
 Let $X:=AQ_n = \Cay(\mathbb{Z}_2^n,S)$ be the augmented cube graph $(n \ge 4)$ and $G:=\Aut(X)$. Then $|G_e|=8$.
\end{Lemma}

We recall briefly the proof given in \cite{Choudum:Sunitha:2001} for the result $|G_e|=8$.  First, it can be shown that the induced subgraph $X[X_1(e)]$ has exactly 8 automorphisms (cf. \cite[Lemma 3.2]{Choudum:Sunitha:2001}).  It can also be shown that every automorphism of the graph $AQ_n$ that fixes the identity vertex $e$ and each of its neighbors is trivial, i.e. $L_e=1$ (cf. \cite[Theorem 3.1]{Choudum:Sunitha:2001}).  It follows that $|G_e| \le 8$.  Finally, Table 1 in \cite{Choudum:Sunitha:2001} defines a set of 8 permutations of the vertex set of $AQ_n$ that fix the identity vertex $e$ and that are automorphisms of $AQ_n$.  Hence $|G_e|=8$.  The paper \cite{Choudum:Sunitha:2001} doesn't show the verification that the eight permutations given in \cite[Table 1]{Choudum:Sunitha:2001} preserve adjacency of $AQ_n$.  In the present section, we express these 8 permutations as linear extensions of automorphisms of the induced subgraph $X[X_1(e)]$.

\begin{Proposition} \label{prop:normal:AQn}
Let $X = AQ_n (n \ge 4)$ and $G:=\Aut(X)$. Then each of the 8 elements in $G_e$ can be realized as a linear transformation  of the vector space $\mathbb{Z}_2^n$ that fixes $S$ setwise, and so $G_e \subseteq \Aut(\mathbb{Z}_2^n,S)$.
\end{Proposition}

\noindent \emph{Proof}:
The induced subgraph $X[X_1(e)]$ is shown in Figure~\ref{fig:induced:subgraph:X1e:AQn}, where the subgraph is drawn horizontally. It can be seen from Figure~\ref{fig:induced:subgraph:X1e:AQn} that automorphisms of the induced subgraph $X[X_1(e)]$ are realized by permuting the endpoints of the edge $\{e_1,11\cdots1\}$, by independently permuting the endpoints of the edge $\{e_{n-1},e_n\}$, and by possibly also interchanging these two edges (in which case the automorphism flips the graph left to right).  Hence, $\Aut(X[X_1(e)]) \cong D_8$.

\begin{figure}
\centering
\begin{tikzpicture}[scale=2]

\vertex[fill] (ve1) at (0,0) [label=below:$e_1$] {};
\vertex[fill] (ve2) at (1,0) [label=below:$e_2$] {};
\vertex[fill] (ve3) at (2,0) [label=below:$e_3$] {};

\vertex[fill] (ven_2) at (5,0) [label=below:$e_{n-2}$] {};
\vertex[fill] (ven_1) at (6,0) [label=below:$e_{n-1}$] {};

\vertex[fill] (v11dot1) at (-0.5,1) [label=above:$11\cdot\cdot1$] {};
\vertex[fill] (v01dot1) at (0.5,1) [label=above:$01\cdot\cdot1$] {};
\vertex[fill] (v001dot1) at (1.5,1) [label=above:$001\cdot\cdot1$] {};
\vertex[fill] (v0001dot1) at (2.5,1) [label=above:$0001\cdot\cdot1$] {};

\vertex[fill] (v0dot111) at (4.5,1) [label=above:$0\cdot\cdot111$] {};
\vertex[fill] (v0dot011) at (5.5,1) [label=above:$00\cdot\cdot011$] {};
\vertex[fill] (ve_n) at (6.5,1) [label=above:$e_n$] {};

\node at (3.5,0.5) {$\cdots$};

\path
    (ve1) edge (v11dot1)
    (ve1) edge (v01dot1)

    (ve2) edge (v01dot1)
    (ve2) edge (v001dot1)

    (ve3) edge (v001dot1)
    (ve3) edge (v0001dot1)

    (ven_2) edge (v0dot111)
    (ven_2) edge (v0dot011)

    (ven_1) edge (v0dot011)
    (ven_1) edge (ve_n)

    (v11dot1) edge (v01dot1)
    (v01dot1) edge (v001dot1)
    (v001dot1) edge (v0001dot1)
    (v0dot111) edge (v0dot011)
    (v0dot011) edge (ve_n)

;
\end{tikzpicture}
\caption{Induced subgraph of the first layer of the distance partition of $AQ_n$ $(n \ge 4)$.}
\label{fig:induced:subgraph:X1e:AQn}
\end{figure}

Define $f_1$ to be the map that takes vertex $e_1$ to $11\cdots1$, and $e_i$ to $e_i (i=2,\ldots,n)$. Let $f_1'$ be the linear extension of $f_1$.  Then $f_1'$ takes $11\cdots1=e_1+\cdots+e_n$ to $f_1'(e_1)+\cdots+f_1'(e_n)=11\cdots1+e_2+\cdots+e_n=e_1$.  Similarly, $f_1'$ takes $e_k+e_{k+1}+\cdots+e_n$ $(k \ge 2)$ to itself.  Thus, $f_1'$ just interchanges the endpoints of the edge $\{e_1,11\cdots1\}$ in the induced subgraph $X[X_1(e)]$.  Since $f_1'$ is an invertible linear transformation and fixes $X_1(e)$ setwise, $f_1' \in \Aut(\mathbb{Z}_2^n,S) \subseteq G_e$; hence, $f_1'$ is an automorphism of the graph $X$.  We have thus realized an element $f_1'$ of $G_e$ as a linear transformation of the vector space $\mathbb{Z}_2^n$.

Similarly, define $f_2$ to be the map that takes $e_{n-1}$ to $e_n$, $e_n$ to $e_{n-1}$, and $e_i$ to $e_i$ $(i=1,2,\ldots,n-2)$. Define $f_3$ to be the map that interchanges $e_1$ and $e_{n-1}$, interchanges $e_2$ and $e_{n-2}$, and so on, and takes $e_n$ to $11\cdots1$. Let $f_2'$ and $f_3'$ be the linear extensions of $f_2$ and $f_3$, respectively. It can be seen that the group of linear transformations $\langle f_1',f_2',f_3' \rangle$ is contained in $G_e$ and is isomorphic to $\langle f_1,f_2,f_3 \rangle \cong D_8$.  Thus, $G_e \cong D_8$ and $G_e \subseteq \Aut(\mathbb{Z}_2^n,S)$.
\qed

We have proved above that

\begin{Theorem} \label{thm:autgroup:sdp:AQn}
 Let $X=AQ_n = \Cay(\mathbb{Z}_2^n,S) (n \ge 4)$ and $G:=\Aut(X)$. Then $G_e \cong D_8$, $AQ_n$ is a normal Cayley graph, and $\Aut(AQ_n) \cong \mathbb{Z}_2^n \rtimes D_8$.
\end{Theorem}

\section{Automorphism group of $AQ_4$}

In this section, we study the clique structure of $X:=AQ_4 = \Cay(\mathbb{Z}_2^4,S)$, where $S:=\{e_1,e_2,e_3,e_4,0011,0111,1111\}$. In particular, we describe the action of the automorphism group of $AQ_4$ on the set of all its maximum cliques.  We show that the graph $AQ_4$ has 12 distinct maximum cliques and that the action of $G:=\Aut(AQ_4)$ on a set of 8 maximum cliques is isomorphic to the action of the automorphism group of two disjoint squares acts on its vertices.  It is then deduced that $\Aut(AQ_4) \cong (D_8 \times D_8) \rtimes C_2$.

Because $X$ is vertex-transitive, there exists a maximum clique in the graph which contains the identity vertex $e$.  Observe from Figure~\ref{fig:AQ4:neighborhood:of:e} that the size of a maximum clique in the induced subgraph $X[X_1(e)]$ is 3.  Hence the clique number $\omega(X) = 4$.

\begin{figure}
\centering
\begin{tikzpicture}[scale=2]

\vertex[fill] (ve) at (0,0) [label=below:$e$] {};

\vertex[fill] (ve1) at (2,-1) [label=left:$e_1$] {};
\vertex[fill] (ve2) at (2,0) [label=left:$e_2$] {};
\vertex[fill] (ve3) at (2,1) [label=left:$e_3$] {};

\vertex[fill] (v1111) at (3,-1.5) [label=right:$1111$] {};
\vertex[fill] (v0111) at (3,-0.5) [label=right:$0111$] {};
\vertex[fill] (v0011) at (3,0.5) [label=right:$0011$] {};
\vertex[fill] (ve4) at (3,1.5) [label=right:$e_4$] {};

\node at (3.5,0.5) {$\cdots$};

\path
    (ve1) edge (v1111)
    (ve1) edge (v0111)

    (ve2) edge (v0111)
    (ve2) edge (v0011)

    (ve3) edge (v0011)
    (ve3) edge (ve4)

    (v1111) edge (v0111)
    (v0111) edge (v0011)
    (v0011) edge (ve4)
;
\end{tikzpicture}
\caption{The identity vertex $e$ and its neighbors in $AQ_4$.}
\label{fig:AQ4:neighborhood:of:e}
\end{figure}

There are exactly 3 4-cliques in $AQ_4$ containing the identity vertex $e$ (see Figure~\ref{fig:AQ4:neighborhood:of:e}).  Hence, there are a total of $3 \times 16 = 48$ 4-cliques in $AQ_4$; however, in this count, each clique was counted 4 times, and so there are exactly 12 distinct maximum cliques in $AQ_4$, each of size 4.

Let $\mathcal{C}_1 := \{e,e_3,e_4,0011\}$ be the set of 4 vertices of $AQ_4$ that form the maximum clique corresponding to the upper triangle in Figure~\ref{fig:AQ4:neighborhood:of:e}.  We call $\mathcal{C}_1$ an \emph{upper} \emph{clique} of $AQ_4$.  Note that $\mathcal{C}_1$ is also a subspace of $\mathbb{Z}_2^4$ of dimension 2. Let $\mathcal{C}_2 := \mathcal{C}_1 + 1000 = \{e_1,1010,1001,1011\}$ be the coset obtained by translating $\mathcal{C}_1$ by $1000$.  Let $\mathcal{C}_3 := \mathcal{C}_1 + 0100 = \{e_2, 0110, 0101, 0111 \}$ and $\mathcal{C}_4 := \mathcal{C}_1 + 1100 = \{1100, 1110, 1101, 1111\}$.  The set of upper cliques of $AQ_4$ is defined to equal $\{\mathcal{C}_1, \mathcal{C}_2, \mathcal{C}_3, \mathcal{C}_4 \}$.  The set of translations used to obtain these 4 4-cliques is $\{e, 1000, 0100, 1100\}$, a subspace of dimension 2 orthogonal to $\mathcal{C}_1$.  Hence, the set of four upper cliques is a partition of the vertex set $\mathbb{Z}_2^4$ of $AQ_4$.

Similarly, a lower clique of $AQ_4$ is a 4-clique obtained by a translation of the clique $\mathcal{C}_5 := \{e, e_1, 0111, 1111\}$ corresponding to the lower triangle in Figure~\ref{fig:AQ4:neighborhood:of:e}.  The set of lower cliques is $\{ \mathcal{C}_5, \mathcal{C}_6, \mathcal{C}_7, \mathcal{C}_8 \}$, where $\mathcal{C}_6 := \mathcal{C}_5 + e_4 = \{e_4, 1001, 0110,  1110\}$, $\mathcal{C}_7 := \mathcal{C}_5 + 0011 = \{0011, 1011, e_2, 1100\}$, and $\mathcal{C}_8 := \mathcal{C}_5 + e_3 = \{ e_3, 1010, 0101,  1101\}$.

The set of middle cliques of $AQ_4$ of course consists of the clique $\mathcal{C}_9 := \{e,0011, e_2, 0111\}$ corresponding to the middle triangle in Figure~\ref{fig:AQ4:neighborhood:of:e} and its translations $\mathcal{C}_{10} := \mathcal{C}_9 + e_4 = \{e_4, e_3, 0101, 0110\}$, $\mathcal{C}_{11} := \mathcal{C}_9 + e_1 = \{e_1, 1011, 1100, 1111 \}$, and $\mathcal{C}_{12} := \mathcal{C}_9 + 1001 = \{1001, 1010, 1101, 1110 \}$.

We have defined above the set $\{ \mathcal{C}_1, \mathcal{C}_2, \ldots,  \mathcal{C}_{12} \}$ of all maximum cliques of $AQ_4$ in terms of translations (cosets) of 2-dimensional subspaces. This yields three different partitions of the vertex set of $AQ_4$, namely $\{ \mathcal{C}_1, \ldots, \mathcal{C}_4\}$, $\{\mathcal{C}_5, \ldots, \mathcal{C}_8 \}$, and $\{ \mathcal{C}_9, \dots, \mathcal{C}_{12} \}$.  Only two of these three partitions are equivalent under the action of the automorphism group of $AQ_4$:

                                                                                \begin{Lemma} \label{lemma:AQ4:two:orbits}                                                                                Let $G:=\Aut(AQ_4)$, and let $\mathcal{C}_1, \ldots, \mathcal{C}_{12}$ be the 12 distinct maximum cliques of $AQ_4$ defined above.  Then, the action of $G$ on $\{\mathcal{C}_1, \ldots, \mathcal{C}_{12} \}$ has exactly 2 orbits: $ \{ \mathcal{C}_1, \ldots, \mathcal{C}_8 \}$ and $\{ \mathcal{C}_9, \ldots, \mathcal{C}_{12} \}$.                                                                                \end{Lemma}

\noindent \emph{Proof}:
Each of the upper cliques $\mathcal{C}_1, \ldots, \mathcal{C}_4$ is a translation of the upper clique $\mathcal{C}_1$.  Because translations are automorphisms of the graph, all 4 upper cliques $\mathcal{C}_1, \ldots, \mathcal{C}_4$ lie in the same orbit of the action of $G$ on maximum cliques.  Similarly, the lower cliques $\mathcal{C}_5, \ldots, \mathcal{C}_8$ lie in the same $G$-orbit, and the middle cliques $\mathcal{C}_9, \ldots, \mathcal{C}_{12} $ lie in the same $G$-orbit.

The proof of Proposition~\ref{prop:normal:AQn} shows that there exists an element in $G_e$ that takes the upper triangle (of Figure~\ref{fig:AQ4:neighborhood:of:e}) to the lower triangle.  Hence, the upper clique $\mathcal{C}_1$ and the lower clique $\mathcal{C}_5$ lie in the same $G$-orbit, whence $\{C_1, \ldots, C_8\}$ lies in a single $G$-orbit.

It now suffices to show that the middle clique $\mathcal{C}_9$ and upper clique $\mathcal{C}_1$ do not lie in the same $G$-orbit. The middle clique $\mathcal{C}_9 = \{e, 0011, e_2, 0111\}$ is a 2-dimensional subspace and is closed under addition. Hence, given any vertex $x \in \mathcal{C}_9$, there exists a translation in $\{ \rho_z: z \in \mathcal{C}_9 \}$ that takes $x$ to the identity vertex $e$ and that maps $\mathcal{C}_9$ to itself. Thus, if $\exists g \in G$ such that $g: \mathcal{C}_9 \mapsto \mathcal{C}_1$, then $\exists g \in G_e$ such that $g: \mathcal{C}_9 \mapsto \mathcal{C}_1$. But this is impossible since none of the 8 elements of $G_e \cong D_8$ take the middle triangle to the lower triangle (cf. Figure~\ref{fig:AQ4:neighborhood:of:e}).  Hence, the middle clique $\mathcal{C}_9$ and the upper clique $\mathcal{C}_1$ lie in different $G$-orbits.
\qed

The argument above for $AQ_4$ can be generalized to show that the graph $AQ_n$ $(n \ge 4$) has exactly $(n-1) \times \frac{2^n}{4} = (n-1)2^{n-2}$ distinct maximum cliques, and that the action of $\Aut(AQ_n)$ on this set of maximum cliques has exactly $\lfloor n/2 \rfloor$ orbits.  In particular, it is easy to see that the 4-cliques corresponding to the first triangle and the last triangle of Figure~\ref{fig:induced:subgraph:X1e:AQn} are in the same orbit, the 4-cliques corresponding to the second triangle and the last-but-one triangle are in the same orbit, and so on.

We showed in Lemma~\ref{lemma:AQ4:two:orbits} that the set $\{\mathcal{C}_1, \ldots, \mathcal{C}_8 \}$ of 8 cliques of $AQ_4$ forms a single orbit in the action of $G$ on all maximum cliques.  We study this restricted action of $G$ (i.e. the transitive constituent) further:

\begin{Proposition}  \label{prop:G:action:on:cliques:AQ4:faithful}
 Let $\mathcal{C}_1, \ldots, \mathcal{C}_8$ be the maximum cliques of $AQ_4$ defined above.  Then, the action of $G:=\Aut(AQ_4)$ on $\{\mathcal{C}_1, \ldots, \mathcal{C}_8 \}$ is faithful.
\end{Proposition}

\noindent \emph{Proof}:
Suppose $g \in G$ and $g$ maps each $\mathcal{C}_i$ to itself $(i=1,2,\ldots,8)$.  It suffices to show $g=1$.  Because $g$ maps the lower clique $\mathcal{C}_5$ to itself and the upper clique $\mathcal{C}_1$ to itself, $g$ maps their intersection $\{e\}$ to itself. Hence $g$ fixes the vertex $e$.  Similarly, $g$ fixes $\mathcal{C}_5 \cap \mathcal{C}_2 = \{ e_1 \}$ and $g$ fixes $\mathcal{C}_1 \cap \mathcal{C}_8 = \{e_3\}$.  But it is clear from Figure~\ref{fig:AQ4:neighborhood:of:e} that the only element from $G_e$ that fixes $e_1$ and $e_3$ is the trivial automorphism, i.e. $g=1$.
\qed

If a group $G$ acts on a set $\Omega$, then a block of $G$ is defined to be a subset $\Delta \subseteq \Omega$ such that $\Delta \cap \Delta^g = \Delta$ or $\Delta \cap \Delta^g = \phi$, for all $g \in G$.

\begin{Proposition} \label{prop:AQ4:K:fixed:block}
 Let $\mathcal{K}:=\{\mathcal{C}_1, \mathcal{C}_2, \mathcal{C}_3, \mathcal{C}_4 \}$ and $\mathcal{K}' := \{ \mathcal{C}_5, \mathcal{C}_6, \mathcal{C}_7, \mathcal{C}_8 \}$ denote the set of 4 upper cliques and the set of 4 lower cliques of $AQ_4$, respectively.  Then $\mathcal{K}$ (and hence $\mathcal{K}'$) is a block of the action of $G$ on $\{\mathcal{C}_1, \ldots, \mathcal{C}_8\}$.
\end{Proposition}

\noindent \emph{Proof}:
Consider the partial drawing of $AQ_4$ shown in Figure~\ref{fig:AQ4:4:upper:cliques}.  In this drawing, the 4 vertices in any upper clique are grouped together and the edges between vertices in the same upper clique are not drawn.  Also, only all edges from the upper clique $\mathcal{C}_1 = \{e,e_3,e_4,0011\}$ to vertices in the other three upper cliques are shown. The edges incident to consecutive vertices of $\mathcal{C}_1$ alternate between solid and dotted format, in order to distinguish between these edges more clearly. Note that the edges between $\mathcal{C}_1$ and $\mathcal{C}_3$ are twice as many as the edges between $\mathcal{C}_1$ and any other upper clique.   Thus, $\{\mathcal{C}_1, \mathcal{C}_3 \}$ is a block of $G$ acting on $\mathcal{K}$.  Studying the automorphism groups of graphs also helps to construct new drawings of graphs.

\begin{figure}
\centering
\begin{tikzpicture}[scale=1.1]

\vertex[fill] (v0000) at (2,0) [label=below:$e$] {};
\vertex[fill] (v0001) at (3,0) [label=below:$e_4$] {};
\vertex[fill] (v0010) at (4,0) [label=below:$e_3$] {};
\vertex[fill] (v0011) at (5,0) [label=below:$0011$] {};
\node at (3.5,-1) {$\mathcal{C}_1$};

\vertex[fill] (v1100) at (0,2) [label=left:$1100$] {};
\vertex[fill] (v1101) at (0,3) [label=left:$1101$] {};
\vertex[fill] (v1110) at (0,4) [label=left:$1110$] {};
\vertex[fill] (v1111) at (0,5) [label=left:$1111$] {};
\node at (-2,3.5) {$\mathcal{C}_4$};

\vertex[fill] (v0100) at (2,7) [label=above:$e_2$] {};
\vertex[fill] (v0101) at (3,7) [label=above:$0101$] {};
\vertex[fill] (v0110) at (4,7) [label=above:$0110$] {};
\vertex[fill] (v0111) at (5,7) [label=above:$0111$] {};
\node at (3.5,8) {$\mathcal{C}_3$};

\vertex[fill] (v1000) at (7,2) [label=right:$e_1$] {};
\vertex[fill] (v1001) at (7,3) [label=right:$1001$] {};
\vertex[fill] (v1010) at (7,4) [label=right:$1010$] {};
\vertex[fill] (v1011) at (7,5) [label=right:$1011$] {};
\node at (9,3.5) {$\mathcal{C}_2$};

\path
    (v0000) edge (v1111)
    (v0000) edge (v0100)
    (v0000) edge (v0111)
    (v0000) edge (v1000)
;

\path[dotted]
    (v0001) edge (v1110)
    (v0001) edge (v0101)
    (v0001) edge (v0110)
    (v0001) edge (v1001)
;

\path
    (v0010) edge (v1101)
    (v0010) edge (v0101)
    (v0010) edge (v0110)
    (v0010) edge (v1010)
;

\path[dotted]
    (v0011) edge (v1100)
    (v0011) edge (v0100)
    (v0011) edge (v0111)
    (v0011) edge (v1011)
;

\end{tikzpicture}
\caption{A partition of the vertex set of $AQ_4$ into cosets of clique $\mathcal{C}_1$.}
\label{fig:AQ4:4:upper:cliques}
\end{figure}

Observe from Figure~\ref{fig:AQ4:4:upper:cliques} that the neighbors of $e \in \mathcal{C}_1$ in the three other upper cliques are $1111, e_2, 0111$ and $e_1$.  Three of these neighbors, namely $e_1, 0111$ and $1111$, along with $e$, form the lower clique $\mathcal{C}_5 = \{e, e_1, 0111, 1111\}$, whose 4 vertices lie in distinct upper cliques. In other words, there exists a lower clique $\mathcal{C}_5$ that intersects all 4 upper cliques $\mathcal{C}_1, \mathcal{C}_2, \mathcal{C}_3, \mathcal{C}_4$.

Recall that $G$ acts transitively on $\mathcal{K}$.  Also, an automorphism of the graph must take pairwise disjoint cliques $\mathcal{C}_1, \ldots, \mathcal{C}_4$ to pairwise disjoint cliques.  Hence, if $g \in G$ takes an upper clique, $\mathcal{C}_1$ say, to a lower clique $\mathcal{C}_5$, then $g$ must take each upper clique to a lower clique, i.e. $g$ must interchange $\mathcal{K}$ and $\mathcal{K}'$. Thus, each $g \in G$ either fixes  $\mathcal{K} = \{\mathcal{C}_1, \mathcal{C}_2, \mathcal{C}_3, \mathcal{C}_4 \}$ setwise, or completely moves $\mathcal{K}$ to a set $\mathcal{K}'$ disjoint from $\mathcal{K}$, i.e. $\mathcal{K}$ is a block of the action of $G$ on $\{\mathcal{C}_1, \ldots, \mathcal{C}_8 \}$.
\qed

In the next result, we essentially prove that the action of $G:=\Aut(AQ_4)$ on the cliques $\mathcal{C}_1, \ldots, \mathcal{C}_8$ is isomorphic to the action of the automorphism group of two disjoint squares on the 8 vertices of the two squares.

\begin{Theorem}
 $\Aut(AQ_4) \cong (D_8 \times D_8) \rtimes C_2$.
\end{Theorem}

\noindent \emph{Proof}:
By Lemma~\ref{lemma:AQ4:two:orbits}, $\{\mathcal{C}_1, \ldots, \mathcal{C}_8 \}$ is an orbit of $G$ acting on maximum cliques.  Consider the action of $G$ on $\{ \mathcal{C}_1, \ldots, \mathcal{C}_8 \}$.  By Proposition~\ref{prop:G:action:on:cliques:AQ4:faithful}, this induced action is faithful.  Hence $G$ is determined uniquely by the permutation group image of this induced action. By Proposition~\ref{prop:AQ4:K:fixed:block}, $\mathcal{K} = \{\mathcal{C}_1, \ldots, \mathcal{C}_4 \}$ and $\mathcal{K}' = \{ \mathcal{C}_5, \ldots, \mathcal{C}_8 \}$ are blocks of $G$.

Since $e_3 \in \mathcal{C}_1$, the translation $\rho_{e_3} \in G$ maps each upper clique $\mathcal{C}_i \in \mathcal{K}$ (which is a coset of the subspace $\mathcal{C}_1$) to itself.  Hence $\rho_{e_3}$ fixes $\mathcal{K}$ pointwise.  On the other hand, $\rho_{e_3}$ interchanges $\mathcal{C}_5$ and $\mathcal{C}_8$ because $e \in \mathcal{C}_5$ and $e_3 \in \mathcal{C}_8$.  Also, $\rho_{e_3}$ interchanges $\mathcal{C}_6$ and $\mathcal{C}_7$ because $e_4 \in \mathcal{C}_6$ and $e_4+e_3 \in \mathcal{C}_7$. Thus, $\rho_{e_3}$ effects the permutation $(\mathcal{C}_5, \mathcal{C}_8)(\mathcal{C}_6, \mathcal{C}_7) \in \Sym( \{ \mathcal{C}_1, \ldots, \mathcal{C}_8 \} )$.  Similarly, $\rho_{e_4}$ fixes $\mathcal{K}$ pointwise and effects the permutation $(\mathcal{C}_5, \mathcal{C}_6)(\mathcal{C}_7, \mathcal{C}_8)$ on $\mathcal{K}'$.

Let $A$ be the linear extension of the map $(e_1,e_2,e_3,e_4) \mapsto (e_1,e_2,e_4,e_3)$, i.e. $A$ interchanges the last 2 coordinates of a vector. The coset $\mathcal{C}_1$ is closed under action by $A$ and each vector in the span of $\{e_1,e_2\}$ is also fixed by $A$.  So $A$ maps the coset $\mathcal{C}_1 +a$ $(a \in \mbox{Span}\{e_1,e_2\})$ to itself because $(\mathcal{C}_1+a)^A = \mathcal{C}_1^A + a^A = \mathcal{C}_1 + a$. Thus, $A$ fixes $\mathcal{K}$ pointwise.

The linear transformation $A$ just defined is invertible.  Note also that $A$ fixes $S$ setwise.  Hence $A \in \Aut(\mathbb{Z}_2^4,S) \subseteq G_e$.  Thus, $A$ permutes the cliques $\mathcal{C}_1, \ldots, \mathcal{C}_8$ among themselves.  Since $A$ fixes $\mathcal{K}$ pointwise, $A$ fixes $\mathcal{K}'$ setwise. It can be verified that $A$ fixes $\mathcal{C}_5$ since $e \in \mathcal{C}_5$.  Also, $A$ interchanges $\mathcal{C}_6$ and $\mathcal{C}_8$ because $e_4 \in \mathcal{C}_6$ and $e_4^A = e_3 \in \mathcal{C}_8$.  Thus, $A$ effects the permutation $(\mathcal{C}_6, \mathcal{C}_8) \in \Sym(\{\mathcal{C}_1, \ldots, \mathcal{C}_8 \})$.

Thus $N_1 := \langle \rho_{e_3}, \rho_{e_4}, A \rangle$ $=\langle (\mathcal{C}_5, \mathcal{C}_8)(\mathcal{C}_6, \mathcal{C}_7), (\mathcal{C}_5, \mathcal{C}_6)(\mathcal{C}_7, \mathcal{C}_8), (\mathcal{C}_6, \mathcal{C}_8) \rangle \cong D_8$ (see Figure~\ref{fig:AQ4:cliques:square}).

\begin{figure}
\centering
\begin{tikzpicture}[scale=1.5]

\vertex[fill] (v1) at (-1,1) [label=left:$\mathcal{C}_1$] {};
\vertex[fill] (v2) at (1,1) [label=right:$\mathcal{C}_2$] {};
\vertex[fill] (v3) at (1,-1) [label=right:$\mathcal{C}_3$] {};
\vertex[fill] (v4) at (-1,-1) [label=left:$\mathcal{C}_4$] {};
\draw[thick] (v1) -- (v2) -- (v3) -- (v4) -- (v1);
\draw[dashed] (-1.5,1.5) -- (1.5, -1.5);
\node at (1.6, -1.6) {$A'$};
\draw[dashed] (0,1.5) -- (0, -1.5);
\node at (0, -1.6) {$\rho_{e_1}$};
\draw[dashed] (-1.5,0) -- (1.5, 0);
\node at (1.9, 0) {$\rho_{1111}$};
\node at (0, -2.2) {$\mathcal{K}$};

\vertex[fill] (v5) at (4,1) [label=left:$\mathcal{C}_5$] {};
\vertex[fill] (v6) at (6,1) [label=right:$\mathcal{C}_6$] {};
\vertex[fill] (v7) at (6,-1) [label=right:$\mathcal{C}_7$] {};
\vertex[fill] (v8) at (4,-1) [label=left:$\mathcal{C}_8$] {};
\draw[thick] (v5) -- (v6) -- (v7) -- (v8) -- (v5);
\draw[dashed] (3.5,1.5) -- (6.5, -1.5);
\node at (6.6, -1.6) {$A$};
\draw[dashed] (5,1.5) -- (5, -1.5);
\node at (5, -1.6) {$\rho_{e_4}$};
\draw[dashed] (3.5,0) -- (6.5, 0);
\node at (6.9, 0) {$\rho_{e_3}$};
\node at (5, -2.2) {$\mathcal{K}'$};

\end{tikzpicture}
\caption{$G:=\Aut(AQ_4)$ acting on maximum cliques.}
\label{fig:AQ4:cliques:square}
\end{figure}

Similarly, consider $N_2:=\langle \rho_{e_1}, \rho_{1111}, A' \rangle$, where $\rho_{e_1}$ maps each lower clique to itself and effects the permutation $(\mathcal{C}_1, \mathcal{C}_2)(\mathcal{C}_3,\mathcal{C}_4)$ on the upper cliques; $\rho_{1111}$ effects the permutation $(\mathcal{C}_1, \mathcal{C}_4)(\mathcal{C}_2, \mathcal{C}_3)$.  Let $A'$
be the linear transformation that takes $e_1$ to $1111$ and $e_i$ to itself $(i=2,3,4)$.  Then, $A'$ fixes $S$ setwise, whence $A' \in \Aut(\mathbb{Z}_2^4,S) \subseteq G_e$. It can be verified that $A'$ effects the permutation $(\mathcal{C}_2, \mathcal{C}_4) \in \Sym(\{\mathcal{C}_1, \ldots, \mathcal{C}_8 \})$.  Thus, $N_2 \cong D_8$.

The permutations in $N_1$ fix each element of $\mathcal{K}$ and move only elements in $\mathcal{K}'$.  The permutations in $N_2$ fix each element of $\mathcal{K}'$ and move only elements in $\mathcal{K}$.  Hence, the elements of $N_1$ and $N_2$ commute with each other and $N_1 \cap N_2=1$. Thus, $N_1 N_2 = N_1 \times N_2 \cong D_8 \times D_8$.

A permutation of $\mathcal{C}_1, \ldots, \mathcal{C}_8$ that interchanges $\mathcal{K}$ and $\mathcal{K}'$ is one induced by the linear map $B: (e_1,e_2,e_3,e_4) \mapsto (e_3, e_2, e_1, 1111)$ (cf. Figure~\ref{fig:AQ4:neighborhood:of:e}).
Then $B$ interchanges the upper and lower cliques.  The subgroup $\langle B \rangle \cong C_2$ normalizes $N_1 \times N_2$, which has index 2 in $(N_1 \times N_2) \langle B \rangle $.  Thus, $(D_8 \times D_8) \rtimes C_2$ is isomorphic to a subgroup of $G$. By order considerations (cf. Theorem~\ref{thm:autgroup:sdp:AQn}), $(D_8 \times D_8) \rtimes C_2$ is the full automorphism group $G$.
\qed

\section{Nontrivial blocks of $AQ_4$}

In this section we obtain an equivalent condition for a subset of vertices of $AQ_n$ to be a nontrivial block of the graph, and we determine all the nontrivial blocks of $AQ_4$.  We first recall some definitions on blocks (cf. Wielandt \cite{Wielandt:1964}, Biggs \cite{Biggs:1993}).

Let $G \le \Sym(V)$ be a transitive permutation group.  A block of $G$ is a subset $\Delta \subseteq V$ that satisfies the condition $\Delta \cap \Delta^g = \Delta$ or $\Delta \cap \Delta^g = \phi$, for all $g \in G$.
A block $\Delta$ of $G$ is said to be nontrivial if $1 < |\Delta| < |V|$.  The set $\{\Delta^g: g \in G\}$ of all translates of a block is called a complete block system (or a system of imprimitivity) and forms a partition of $V$.   A permutation group $G$ is said to be imprimitive if it has a nontrivial block.

A graph $X=(V,E)$ is imprimitive if it is transitive and its automorphism group is imprimitive.  The blocks of $X$ are defined to be the blocks of its automorphism group. Thus, if $W \subseteq V$ is a block of $X$ and $g \in \Aut(X)$, then $W^g$ and $W$ are either equal or disjoint (i.e. there cannot be overlap that is only partial).  For example, for the hypercube graph, an antipodal pair $\{x, \overline{x} \}$ of vertices is a nontrivial block. The hypercube is bipartite, and the set of all even weight vertices is another nontrivial block. It can be shown (cf. Biggs \cite[Proposition 22.3]{Biggs:1993}) that the hypercube has only two nontrivial block systems.

We will use the following result:

\begin{Theorem} \cite[Theorem 7.4, 7.5]{Wielandt:1964} \label{thm:lattices:blocks:subgroups}
 Let $G \le \Sym(V)$ be a transitive permutation group and let $e \in V$.  Then, the lattice of blocks of $G$ which contain $e$ is isomorphic to the lattice of subgroups of $G$ that contain $G_e$. More specifically, the correspondence between blocks $\Delta$ and subgroups $H$ is as follows: $\Delta := e^H$ (the orbit of $e$ under the action of $H$), and $H:=G_{ \{ \Delta \}}$ (the set of elements of $G$ that fixes $\Delta$ setwise).
\end{Theorem}

To obtain the blocks of the augmented cube graph $AQ_n$, we determine the subgroups of $G:=\Aut(AQ_n)$ that contain $G_e$ and then use Theorem~\ref{thm:lattices:blocks:subgroups}.   For $z \in \mathbb{Z}_2^n$, let $(\rho_z: \mathbb{Z}_2^n \rightarrow \mathbb{Z}_2^n, x \mapsto x+z)$ denote translation by $z$. For a subset $K \subseteq \mathbb{Z}_2^n$, let $R(K):= \{ \rho_z: z \in K\}$.  We say that a subset $K \subseteq \mathbb{Z}_2^n$ is closed under action by $G_e$ if $k^g \in K$, for all $k \in K, g \in G_e$.

\begin{Lemma} \label{lemma:condition:on:subset:K}
 Let $X:=\Cay(\mathbb{Z}_2^n,S)$ be a normal Cayley graph (here, $X$ need not be $AQ_n$), and let $G:=\Aut(X) = R(\mathbb{Z}_2^n) \rtimes G_e$.  Let $K$ be a subset of vertices of the graph which forms a subspace of $\mathbb{Z}_2^n$ and which is closed under action by $G_e$.  Then $R(K)G_e$ is a subgroup of $G$.
\end{Lemma}

\noindent \emph{Proof}:
The group $G$ is finite, so it suffices to show that $R(K)G_e$ satisfies closure.  We need to show that for all $z_1,z_2 \in K$ and $g_1,g_2 \in G_e$, there exist $z_3 \in K$ and $g_3 \in G_e$ such that $(\rho_{z_1} g_1)(\rho_{z_2} g_2) = \rho_{z_3} g_3$.  This latter equation is equivalent to the condition $x^{(\rho_{z_1} g_1)(\rho_{z_2} g_2)} = x^{\rho_{z_3} g_3}$, for all $x \in V(X)$.  The left hand side is $x^{{\rho_{z_1} g_1}{\rho_{z_2} g_2}} = (x+z_1)^{g_1 \rho_{z_2} g_2}$ $=(x^{g_1} + z_1^{g_1} + z_2)^{g_2}$ $=x^{g_1 g_2} + z_1^{g_1 g_2} + z_2^{g_2}$.
Here, we used the fact that $g_i$ is a linear transformation, which is the case because for a normal Cayley graph $G_e \subseteq \Aut(\mathbb{Z}_2^n)$. The right hand side is $x^{\rho_{z_3} g_3} = x^{g_3} + z_3^{g_3}$.  So we get that
\[
 x^{g_1 g_2} + z_1^{g_1 g_2} + z_2^{g_2} = x^{g_3}+z_3^{g_3}, \forall x \in V(X).
\]
Take $g_3=g_1 g_2$. Then $z_3$ must satisfy the equation $z_1^{g_1 g_2} + z_2^{g_2} = z_3^{g_1 g_2}$, or equivalently, the equation $z_1 + z_2^{g_1^{-1}} = z_3$.  By hypothesis, $K$ is closed under addition and under action by $g_1 \in G_e$.  Hence, there exists a $z_3 \in K$ satisfying the equation. It follows that $R(K)G_e$ is a subgroup of $G$.
\qed

We now show that the subgroups $R(K)G_e$ obtained from Lemma~\ref{lemma:condition:on:subset:K} are all the subgroups of $G$ that contain $G_e$:

\begin{Proposition}
 Let $X:=\Cay(\mathbb{Z}_2^n,S)$ be a Cayley graph and let $G:=\Aut(X)$.  If $G_e \le H \le G$, then $H=R(K) G_e$ for some subspace $K$ of $\mathbb{Z}_2^n$ that is closed under action by $G_e$.
\end{Proposition}

\noindent \emph{Proof}:
Suppose $H$ satisfies $G_e \le H \le G$. Then $H$ can be expressed as a disjoint union of left cosets of $G_e$ in $H$, i.e.  $H = h_1 G_e \cup h_2 G_e \cup \ldots \cup h_k G_e$, where $k$ is the index of $G_e$ in $H$.  Note that $h_i \in H \le G = R(\mathbb{Z}_2^n) G_e$, and so each $h_i$ can be expressed as $h_i = \rho_{z_i} g_i$ for some $z_i \in \mathbb{Z}_2^n, g_i \in G_e$. Then, $h_i G_e = \rho_{z_i} g_i G_e = \rho_{z_i} G_e$. Thus, $H=R(K)G_e$, where $K:=\{z_1,\ldots,z_k\}$.  This proves that if $H$ is a subgroup of $G$ which contains $G_e$, then $H=R(K) G_e$ for some subset $K \subseteq \mathbb{Z}_2^n$.

To show $K$ is a subspace, suppose $z_1,z_2 \in K$ (it is clear $K$ is nonempty). We show $z_1+z_2 \in K$.  By closure in the subgroup $R(K)G_e$,  $(\rho_{z_1} g_1)(\rho_{z_2} g_2) \in R(K) G_e$, for all $g_1, g_2 \in G_e$. In particular, taking $g_1 = g_2 = 1$, we get that $\rho_{z_1} \rho_{z_2} \in R(K) G_e$. This implies that there exist $z_3 \in K$ and $g_3 \in G_e$ such that $\rho_{z_1 + z_2} = \rho_{z_3} g_3$. The right hand side takes vertex $z_3$ to $z_3^{\rho_{z_3} g} = e^g =e$, whence $z_3^{\rho_{z_1+z_2}} = e$. But this implies that $z_3 = z_1 + z_2$.  Thus, for closure to hold in $R(K)G_e$, $R(K)$ must also contain $\rho_{z_1+z_2}$, i.e. $K$ is closed under addition and hence is a subspace.

It remains to show that $K$ is closed under action by $G_e$. Let $z_1 \in K, g_1 \in G_e$.  We show $z_1^{g_1} \in K$.  Since $R(K)G_e$ is a subgroup, there exist $z_2 \in K$ and $g_2 \in G_e$ such that $\rho_{z_1} g_1 = g_2 \rho_{z_2}$, i.e. $x^{\rho_{z_1} g_1 \rho_{z_2}} = x^{g_2}$, for all $x \in V(X)$.  Taking $x=e$, we get $z_1^{g_1} = z_2$, i.e. $z_1^{g_1}$ is also in $K$.
\qed

 Thus, the subgroups $H$ of $G$ that contain $G_e$ are in bijective correspondence with the subspaces $K$ of $\mathbb{Z}_2^n$ that are closed under action by $G_e$.  The conditions of Lemma~\ref{lemma:condition:on:subset:K} require $K$ to be closed under action by $G_e$. This implies $K$ is a union of orbits of $G_e$.

Suppose we have a subset $K \subseteq \mathbb{Z}_2^n$ satisfying the conditions of Lemma~\ref{lemma:condition:on:subset:K}. By Theorem~\ref{thm:lattices:blocks:subgroups}, the subgroup $H:=R(K)G_e$ gives rise to the block $\Delta = e^H = e^{R(K)G_e} = K^{G_e} = K$.  Given such a block $\Delta$, a complete block system $\{\Delta^g: g \in G\}$ is the set  of all translations of $\Delta$: if $\rho_z g \in G$ $(z \in \mathbb{Z}_2^n, g \in G_e)$, then $\Delta^{\rho_z g} = K^{\rho_z g} = (K+z)^g = K^g + z^g = K+z^g = \Delta+z^g$. Thus, to determine all nontrivial block systems of a normal Cayley graph $X=\Cay(\mathbb{Z}_2^n,S)$, it suffices to determine all subspaces $K$ of $\mathbb{Z}_2^n$ that are closed under action by $G_e$.

By the arguments above, we get the following characterization for the nontrivial blocks of a normal Cayley graph:

\begin{Theorem} \label{thm:characterization:blocks}
 Let $X:=\Cay(\mathbb{Z}_2^n,S)$ be a normal Cayley graph and let $G:=\Aut(X)$.  Then, the nontrivial blocks of $X$ that contain the identity vertex $e$ are precisely the subspaces of $\mathbb{Z}_2^n$ that are closed under action by $G_e$.
\end{Theorem}

We now apply the above results to the special case where the graph is $AQ_4$.

\begin{Proposition} \label{prop:nontrivial:blocks:AQ4}
Each of the following subsets is a nontrivial block of $AQ_4$:
$\Delta = \{e,e_2\}$; $\Delta' := \{e, e_2, 0011, 0111\}$; $\Delta'' := \{ e, e_2, 0011, 0111,$ $1001, 1010, 1101, 1110 \}$.
\end{Proposition}

\noindent \emph{Proof}:
The vertex $e_2$ is fixed by $G_e$ (cf.  Figure~\ref{fig:AQ4:neighborhood:of:e}).  Hence, the subspace spanned by $\{e_2\}$ is fixed by $G_e$. By Theorem~\ref{thm:characterization:blocks}, this subspace is a block of $AQ_4$.  The four vectors in $\Delta'$ form a 2-dimension subspace $\mbox{Span} \{e_2,0011\}$ of $\mathbb{Z}_2^4$ and this subset is closed under action by $G_e$ (see Figure~\ref{fig:AQ4:neighborhood:of:e}).  By Theorem~\ref{thm:characterization:blocks}, $\Delta'$ is a block of $AQ_4$.

It suffices to show that $\Delta''$ is a union of orbits of $G_e$ and that $\Delta''$ forms a subspace.  It is clear from the drawing of $AQ_4$ (see Figure~\ref{fig:AQ4:neighborhood:of:e}) that $\{e\}$, $\{e_2\}$ and $\{0011, 0111\}$ are fixed blocks (hence union of orbits) of $G_e$.
Since $G_e$ fixes $e_2$, $G_e$ fixes the set $\{1001, 1010, 1101, 1110\}$ of non-neighbors of $e_2$ in $X_2(e)$ setwise. Thus, this set is a fixed block of $G_e$ and $\Delta''$ is a union of orbits of $G_e$. It can be verified that $\Delta''$ is a subspace of dimension 3 spanned by $\{0011, 0100, 1001\}$.  Hence $\Delta'' \subseteq \mathbb{Z}_2^4$ is such that $R(\Delta'') G_e$ is a subgroup, and hence, the corresponding set $\Delta''$ is a block of $G$.
\qed

\begin{Theorem}
 The 3 nontrivial blocks given in Proposition~\ref{prop:nontrivial:blocks:AQ4} are the only nontrivial blocks of $AQ_4$ which contain the vertex $e$.
\end{Theorem}

\noindent \emph{Proof}:
By Theorem~\ref{thm:characterization:blocks}, it suffices to find all the subspaces $K$ of $\mathbb{Z}_2^4$ that are closed under action by $G_e$.  We can find the orbits of $G_e$, and determine which fixed blocks (union of orbits) of $G_e$ are subspaces.  Since $K$ must be nontrivial, $K \ne \{e\}$ and $K \ne \mathbb{Z}_2^4$.

It is clear from Figure~\ref{fig:AQ4:neighborhood:of:e} that $\Delta_1 := \{e\}$, $\Delta_2 := \{e_2\}$, $\Delta_3 := \{0011, 0111\}$ and $\Delta_4 := \{e_1, 1111, e_3, e_4 \}$ are orbits of $G_e$. Note that $\Delta_4$ spans $\mathbb{Z}_2^4$, and so $K$ does not contain $\Delta_4$.  Define the linear transformations $A_1: (e_1,e_2,e_3,e_4): \mapsto (e_1,e_2,e_4,e_3)$, $A_2: (e_1,e_2,e_3,e_4) \mapsto (1111, e_2, e_3, e_4)$, and $A_3: (e_1,e_2,e_3,e_4) \mapsto (e_3,e_2,e_1,1111)$.  Then $\langle A_1,A_2,A_3 \rangle = G_e \cong D_8$.

Let $x:=0101$. Then $x^{A_1} = 0110$, $x^{A_3} = 1011$ and $x^{A_1 A_3} = 1100$.  Hence, $\Delta_5 := \{0101, 0110, 1011, 1100 \}$ lies in a single $G_e$-orbit.  Let $y:=1110$. Then $y^{A_1} = 1101$, $y^{A_2} = 1001$ and $y^{A_1 A_2} = 1010$.  Hence, $\Delta_6 := \{1110, 1101, 1001, 1010 \}$ lies in a single $G_e$-orbit. If $\Delta_5 \dot{\cup} \Delta_6 \subseteq K$, then $|K| \ge |\Delta_5| + |\Delta_6| + |\{e\}| = 9$, but since $K$ is a  subspace, $|K|$ is a power of 2, and so $K = \mathbb{Z}_2^4$, a contradiction. Thus, at most one of $\Delta_5$ or $\Delta_6$ is contained in $K$.  If $K \supseteq \Delta_5$, then $|K| \ge 5$, and so $|K|=8$. But it can be verified that $\Delta_1 \dot{\cup} \Delta_2 \dot{\cup} \Delta_3 \dot{\cup} \Delta_5$ is not a subspace.  The subset $\Delta_1 \dot{\cup} \Delta_2 \dot{\cup} \Delta_3 \dot{\cup} \Delta_6$ is the subspace $\Delta''$ of Proposition~\ref{prop:nontrivial:blocks:AQ4}.  If $K$ contains neither $\Delta_5$ nor $\Delta_6$, then the only choices for $K$ are $\Delta_1 \dot{\cup} \Delta_2$ and $\Delta_1 \dot{\cup} \Delta_2 \dot{\cup} \Delta_3$ (denoted $\Delta$ and $\Delta'$, respectively, in Proposition~\ref{prop:nontrivial:blocks:AQ4}).
\qed

\section{Acknowledgements}

Thanks are due to the authors of \cite{Choudum:Sunitha:2002} \cite{Choudum:Sunitha:2001} for sending copies of their papers.
{
\bibliographystyle{plain}


\section{Introduction}

The understanding of the morphology of galaxies and its relation with their physical properties has evolved during the last decades as new insights are provided by observations. It is well-known that the angular momentum content correlates with the morphology of galaxies along the Hubble sequence.
\citet{fall1983} has found evidences suggesting that the specific angular momentum content of both disc and elliptical galaxies holds a clear correlation with the stellar mass, but with the elliptical systems it shifts to lower values. Recently \citet{romanowsky2012}  have analysed a larger database of galaxies with different morphologies, revisiting and confirming \citet{fall1983}'s findings. These authors estimated the specific angular momentum of the stellar component $j_{\rm s}$ using kinetic and photometric data and the stellar mass obtained from near infrared luminosities. Using this larger database, they found that not only disc-dominated galaxies follow this relation but also elliptical ones. However, the latter has lower $j_{\rm s}$ at a given $M_{\rm star}$. The authors suggested that observed galaxies of all morphological types lie along almost parallel sequences with exponents $\alpha \sim 0.6$ in the $j_{\rm s}-M_{\rm star}$ diagram, but with an off-set between late and early type.   More recently, \
citet{Fall2013} showed that galaxies position in the $j_{\rm s}-M_{\rm star}$ diagram is correlated with their position in the disc-to-bulge-ratio-$M_{star}$ plane. They concluded that the description of galaxies in terms of $j_{\rm s}-M_{\rm star}$ may be regarded as a physically motivated alternative to the traditional morphological classification. The detected  $j_{\rm s}-M_{\rm star}$ relation for spheroidal galaxies suggests an angular momentum efficiency of $10-20\%$ regardless of $M_{\rm star}$  \citep{Fall2013}.

In a hierarchical clustering universe, protogalaxies acquire their angular momenta via tidal torques due to small density fluctuations \citep{peebles1969}. As baryons are initially drawn by the same tidal torques as the dark matter haloes, they are expected to acquire similar specific angular momenta.
The specific angular momentum of the dark matter haloes is supposed to grow as $M_{\rm vir}^{2/3}$. Hence, if baryons collapse conserving their angular momentum, they should also determine a similar relation.
 Based on this hypothesis, \citet{fall1980} postulated the so-called standard  disc model (SDM). This model was able to explain the observed correlation between the angular momentum  of discs and the stellar masses \citep{MMW98}.

Several numerical works have shown that if a disc component is able to assemble then it does so conserving its specific angular momentum. Although within the hierarchical scenario for galaxy formation, this might imply a combination of processes which could contribute to retain or balance the angular momentum content so that discs are formed as expected from the SDM, \citep[e.g.][]{Agertz2011,brook2011,deRossi2012,Padilla2014,Nelson2015}.  In particular, SN feedback plays a critical role in determining the fraction of disc components, as well as mergers and interactions \citep[e.g][]{dominguez1998,Scan2009,guedes2011,Sales2012,Pedrosa2014,Lagos2015}.
Elliptical galaxies  are expected to be affected by major mergers or repetitive minor ones. This tends to decrease the angular momentum  content of the material which ends up dominated by velocity dispersion\citep[e.g.][]{meza2003,parry2009,hopkins2010,bournaud2011,sales2010,Pedrosa2014}. Internal disc instabilities could also contribute to  shape bulges and elliptical galaxies \citep[e.g.][]{efstathiou1982}.

As the final morphology of a galaxy is directly linked to its angular momentum content, there must also be  a connection with its characteristic size: when angular momentum is conserved, more extended structures are able to form. Recently, \citet{kravtsov2013} has found that characteristic sizes of stellar and gas distributions in galaxies scale approximately linearly with the virial radius derived from the abundance matching ansatz. This author found that the relation is in agreement with expectations from  the model of \citet{MMW98}.  Interestingly, he also found that this relation is present not only for late-type galaxies but also for early-type systems. From these findings, he concluded that the angular momentum content plays a crucial role in determining the sizes of galaxies of all morphological types, not only of those dominated by rotation.

In a recent paper, we have made a thorough study of the relation between the morphology of galaxies and their formation history, \citet{Pedrosa2014}. The key role played by the angular momentum content of galaxies was studied. We found that when a disc component is able to form it does so following the expected relation from the standard disc formation model in agreement with previous works \citep[see also][]{governato2010,Sales2012}. More interestingly, in this simulation,  spheroidal components are shown to follow a similar correlation but  off-set with respect to that of  rotational supported systems. It became clear from this work that SN feedback played a major role.

In this paper, we  analyse the angular momentum content of the disc and bulge components of  the  best resolved galaxies in a $\Lambda$-CDM scenario.  \citet{Pedrosa2014} analysed a sister simulation (S230A) from the Fenix Project with a stronger feedback parameter. We explore
a run (S230D) with a slightly weaker SN feedback and a different chemical distribution of the synthetized metals, which better reproduces observational results of the metallicity gradients (Tissera et al. 2015, submitted).
This is a very important point which helps to set constraints on the strength of the SN feedback \citep[e.g.][]{gibson2013}. The simulated galaxies have also larger $D/T$ stellar ratios compared to those measured in S230A.
 We also follow them back in time to
establish the possible existence of a redshift dependence.

This work is organized as follows: In Section \ref{sec:simus}, we describe the initial conditions and free parameters used to run the analysed simulation and  the methods used to identify and decompose galaxies. In Section \ref{sec:angmom}, we analyse the angular momentum content of our galaxies and their dynamical components. In Section \ref{sec:size}, we study  the characteristic sizes of the galaxies, and the connection with their angular momentum content. Finally, in Section \ref{sec:conclu}, we summarize our main findings.

\section{Numerical Experiments}
\label{sec:simus}

We used a version of the  code {\small GADGET-3}, an update of {\small GADGET-2 } \citep{springel2003, springel2005}, optimized for massive parallel simulations of highly inhomogeneus systems. This version of {\small GADGET-3} includes treatments for metal-dependent radiative cooling, stochastic star formation (SF), chemical enrichment, and the multiphase model for the interstellar medium (ISM) and the Supernova (SN) feedback scheme of \citet{scan2005,scan2006}.
The SN feedback model is able to successfully trigger galactic mass-loaded winds without introducing mass-scale parameters.

Our code considers energy feedback by Type II (SNII) and Type Ia (SNIa) Supernovae. These processes  are grafted into a multiphase model for the ISM, both developed by \citet{scan2006}.
  The energy released by each SN event   is distributed between the cold and the hot phases of the  simulated ISM. The adopted multiphase model  allows the coexistence of diffuse and dense gas phases. In this model, each gas particle defines the cold and hot phases by using a local entropy criteria. This definition allows particles to decouple hydrodynamically from particular low-entropy ones if they are not part of a shock front. Each cold gas particle has a reservoir where injected SN energy is stored until they fullfil the conditions to join their local hot phase. It is  worth mentioning that this SN feedback scheme does not include parameters that  depend on the global properties of the given galaxy (e.g. the total mass, size), thus making it suitable for cosmological simulations where systems with different masses have formed in a complex way. Our code also includes the chemical evolution model developed by \citet{mosco2001}. SNII and SNIa contribute with chemical elements estimated by adopting the
chemical yields of \citet{WW95} and \citet{iwamoto1999}, respectively.

As mentioned before, for this study, we used one of the Fenix Project's simulations (S230D). These share the same initial conditions but assume different parameters for the SF and SN feedback sub-grid physics.
The initial conditions are consistent with the concordance model with $\Omega_{\Lambda}=0.7$, $\Omega_{\rm m}=0.3$, $\Omega_{b}=0.04$, a normalization of the power spectrum of $\sigma_{8}=0.9$ and $H_{0}= 100 h \ {\rm km} \ {\rm s}^{-1}\ {\rm Mpc}^{-1}$, with $h=0.7$.
The simulated volume represents a box of $10$ Mpc h$^{-1}$ comovil with $2 \times 230^3$ initial particles, achieving a mass resolution of $5.9\times 10^{6}{\rm h^{-1}\ M_{\odot}}$ and $9.1\times 10^{5}{\rm  h^{-1}\  M_{\odot}}$ for the dark matter and initial gas particles, respectively. The gravitational softening is $0.5$~kpc~h$^{-1}$.   We acknowledge the fact that the initial condition represents a small volume of the Universe. It  was chosen to represent a typical region of the Universe with no large structure nearby. Using the S230A run, \citet{deRossi2013}  verified that haloes' growth is well-reproduced by confronting the results with those obtained from the Millennium simulation by \citet{fakhouri2010}. The S230A was also analysed by  \citet{Pedrosa2014} to study the morphology and angular momentum content of the simulated galaxies. Recently, \citet{artale2015} introduced the effects of high-mass X-ray binary feedback in the very early Universe in
a counterpart run of S230D.

Simulation S230D used a higher gas density threshold for star formation and a  lower energy per SN event ($0.7 \times 10^{51 }$~erg) than in S230A \citep[e.g.][]{deRossi2013,Pedrosa2014}. Simulation 230D injects a larger fraction of new synthetized chemical elements
 into the cold phase, while the fraction of SN energy is equally distributed between the cold and hot phases (S230A injected equal fraction
of energy and chemical elements into the cold and hot phases close to the SF region). This different combination
of SF, SN feedback parameters  and metal distribution produced a good representation of chemical abundances on the disc components (Tissera et al., submitted).

\subsection{The simulated galaxy sample}
\label{sec:sample}

Our galaxy sample was built up   by using a Friends-of-Friends algorithm to identify the virialized structures. Then, the substructures were discerned by using the SUBFIND scheme \citep{springel2001}, which iteratively determines the self-bound substructures within the virial radius, $r_{200}$. {\bf For illustration purposes, an example of a late type galaxy of our simulated catalog is shown in Fig.~\ref{fig:Map_dens}, (upper panels). }

To diminish resolution issues, we restricted our study to systems with a total number of particles greater than $2000$ in  the main  galaxy (which corresponds to $M_{\rm vir} \sim 10^{10.3} M_{\odot}$).  We also tested that the estimated relations are conserved
when this particle number limit is increased to $\sim 10000$.

In order to classify our simulated galaxies into early-type and late-type systems,  we estimated the parameter $\epsilon$ of the  star (gas) particles defined as $\epsilon = J_{\rm z} /J_{\rm z,max}(E)$, where $J_{\rm z}$ is the angular momentum component in the direction of the total angular momentum, and $J_{\rm z,max}(E)$ is the maximum $J_{\rm z}$ over all particles at a given binding energy, $E$. A particle on a perfect prograde circular orbit in the disc plane should have  $\epsilon=1$.
We considered particles with $\epsilon$ higher than $0.5$ to be part of a disc. The particles  that did not satisfy this requirement were considered to belong to the spheroidal component. By inspecting the $\epsilon - E$ plane, we checked whether the adopted limits were suitable to individualize the rotationally supported components.  This threshold is less restrictive than those usually adopted
\citep[e.g.][]{tissera2012}, so if there were stars in a thick disc, they would contribute by decreasing the signals for angular momentum conservation, for example.
 We also distinguished between the central spheroid and the stellar  halo components depending on their bounding energy following \citet{tissera2012}. In Fig.~\ref{fig:Map_dens} (lower panels)  we display the disk and central spheroid components for the galaxy previously shown, after the descomposition.

Therefore, for each simulated galaxy, we determine the  disc (i.e.rotational-dominated) and the central spheroids components (i.e.dispersion-dominated) based on these dynamical criteria. This procedure  was also applied to galaxies identified at $z = 1$ and $z = 2$ in order to study the evolution of the angular momentum content in the galaxies,
their central spheroids and disc components and the dark matter haloes.

\begin{figure*}
\centering
\resizebox{5cm}{!}{\includegraphics{mapden_faceon_full.eps}}
\hspace*{-0.2cm}
\resizebox{5cm}{!}{\includegraphics{mapden_edgeon_full.eps}}\\
\hspace*{-0.2cm}
\resizebox{5cm}{!}{\includegraphics{faceon_bulge.eps}}
\hspace*{-0.2cm}
\resizebox{5cm}{!}{\includegraphics{faceon_disc.eps}}
\hspace*{-0.2cm}
\caption{Upper panels: Density map for a face-on and an edge-on view of a late type galaxy from the 230D simulation catalog. Lower panels: Disc and central spheroid components.}
\label{fig:Map_dens}
\end{figure*}

\section{Angular Momentum }
\label{sec:angmom}

To calculate the properties of the simulated galaxies, we adopted  the half-mass stellar ratio ($r_{\rm HM}$), defined as the radius that enclosed $50\%$ of the total stellar mass. This characteristic radius was chosen in order to make an adequate comparison with observational results. The half-mass radius is usually adopted to estimate galaxy properties in observations. In the case of \citet{kravtsov2013} for example, he uses the stellar half mass radius to estimate his $R_{hm} - R_{vir}$ relation. We estimated the ratio between  the  angular momentum of the stars (gas) in the discs ($J_{\rm D}$) within $r_{\rm HM}$  and the corresponding of the dark matter haloes within the virial radius ($J_{\rm H} $), $j_{\rm d} = J_{\rm D} /J_{\rm H} $, as a function of the ratio between the corresponding masses,  $m_{\rm d} = M_{\rm D} /M_{\rm H} $\footnote{Note that capital letters will be used for the total angular momentum and total mass and to denote the
dynamical components: disc (D), central spheroid (S) or dark matter halo (H).}. We followed the classical notation of \citet{MMW98} to define $j_{\rm d}$ and $m_{\rm d}$. Similarly, we defined $j_{\rm s}$ and $m_{\rm s}$ ratios by using the star particles identified to belong to the central spheroidal components. Those corresponding to the whole stellar mass and the total baryonic mass were also estimated and analysed. These relations can be seen  in Fig.~\ref{fig:jdmd_D} for $z=0$ (upper panel), $z=1$ (middle panel) and $z=2$ (lower panel).  As shown, there is a clear correlation for the disc and central spheroidal components  in the three redshifts \citep[see also][]{dominguez1998,Pedrosa2014}.

Let us first analyse the relations at $z=0$.  As can be seen in Fig.~\ref{fig:jdmd_D}, there is a one-to-one relation between $j_{\rm i}$ and $m_{\rm i}$, where $\rm i$ represents the gas or stellar discs components as expected in the case of angular momentum conservation.   As discussed before, the one-to-one correlation is the consequence of the fact that all the material in the system experienced the same external torques before it collapsed to form a bounded system, and that the collapse occurred with angular momentum conservation.
 Interestingly, for the stellar central spheroidal components a clear correlation is also detected but off-set from that
determined by the disc components. This indicates that
the material which ended up dominated by dispersion lost angular momentum as expected,  but it did so by a similar fraction, regardless of the  final stellar fraction, $m_{\rm i}$. This trend was previously shown in \citet[][ figure 3]{Pedrosa2014} and in \citet[][see also \citet{dominguez1998} for earlier results with lower numerical resolution and a different numerical code and subgrid physics]{Pedrosa2015}.

In Fig.~\ref{fig:jdmd_D} the same relations are estimated at $z=1$ and $z=2$. As can be seen, both the stellar and gaseous discs are always consistent with the one-to-one relation. However, the most interesting result is again provided by the spheroidal components which determine a similar correlation as a function of redshift. Regardless of their formation history, the disc and spheroidal components
 determine   relations with the stellar and dark matter masses which are  statistically unchanged up to $z \sim 2$.  The gap between the angular momentum content of the galaxy respect to that of the dark matter  does not show a systematic change with redshift.

Hence, our simulated galaxies have disc and central spheroidal components consistent with the relation $ j_{\rm i} = A \times  m_{\rm i}^{\alpha} $, where $A  \sim 1$ for the disc components and
$ A \sim 0.10-0.20$ for the central spheroidal components.
We fit a linear regression to the log-log relations finding slopes of $\alpha \sim 0.9 $ and $\alpha \sim 1.1 $ for the stellar discs and central spheroids, respectively at $z=0$.

These relations are found not to evolve significantly with redshift at least up to $z \sim 2$ ($\alpha \sim 0.90,0.87$ and $\alpha \sim 0.96, 0.86$ for the stellar discs and spheroids at $z=1$ and $z=2$, respectively)
Therefore, the spheroidal components, albeit dominated by dispersion, retain a similar fraction of the angular momentum imprinted by tidal torques onto the  small perturbations
during the first stages of the evolution of the structure \citep{peebles1969}, regardless of their assembly histories or  $m_{\rm i}$.

The slope of the relation $j_{\rm i}$-$m_{\rm i}$ changes when the total stellar mass (i.e. total stellar discs and bulges) or the  baryonic mass are considered for the estimation of these ratios  (blue and black lines, respectively), reaching values close to  $\alpha \sim 2 $ due to the combination of  angular momentum content of both disc and bulge components. This is in agreement with the results obtained by previous works where the angular momentum of baryons were studied without discerning between disc and spheroidal components \citep[e.g.][]{Sales2012}.  In this case, the gap is filled with the combination of contributions from both components and the slope of the mean relation gets steeper. We found indications of a flattening in the relation for the total baryonic mass with redshift ($\alpha \sim 1.54, 1.40 $ for $z=1$ and $z=2$, respectively ). However, a larger sample is needed to confirm the level of evolution since the dispersion in our sample increases with increasing
redshift.

\begin{figure*}
\centering
\resizebox{7cm}{!}{\includegraphics{mdjd_bsg_rhm_z0_D.eps}}
\resizebox{7cm}{!}{\includegraphics{mdjd_scat_plot.eps}}\\
\resizebox{7cm}{!}{\includegraphics{mdjd_bsg_rhm_z1_D.eps}}
\resizebox{7cm}{!}{\includegraphics{mdjd_bsg_rhm_z2_D.eps}}
\hspace*{-0.2cm}
\caption{Relation $j_{\rm i}-m_{\rm i}$ medians for the stellar (red dashed lines) and gaseous (green dashed lines) disc components, the stellar spheroidal components (pink dashed lines), the total stellar component (blue line) and the total baryonic component (black lines). The solid black line represents the relation in the case of angular momentum conservation. The stellar and gas disc components are consistent with angular momentum conservation while the spheroidal components  determine a relation with a similar slope but off-set by a constant factor suggesting that spheroids have lost a similar fraction of angular momentum with respect to the total angular momentum of their virial haloes. The relation has been estimated for $z=0$ (upper panel, in the right one we include the scatter data for the stellar disc component), and for $z=1$ and $z=2$ (bottom panel, left and right respectively), finding similar trends. }
\label{fig:jdmd_D}
\end{figure*}

We have also analysed the angular momentum  efficiency $\eta_{\rm i}$ for both components, defined as the ratio
 between  $j_{\rm D}$  ($j_{\rm S}$) of the disc component (central spheroids) and  $j_{\rm H}$ of the dark matter halo, following  \citet{sales2010}.
In Fig.~\ref{fig:Eta_D} we  plot $\eta_{\rm i}$ versus $m_{\rm i}$ for the three
analysed redshifts. As expected, the disc components are found to have $\eta \sim 1$, on average, reflecting the conservation of angular momentum during their formation. Conversely, the  central spheroid has $\eta \sim 0.10-0.20$, on average.
In both cases, there is no clear evolution with redshift (although the dispersion increases slightly with increasing $z$).
The different locations of $\eta $ for the two dynamical components reflect their different evolutionary histories which  introduce dispersion but, nevertheless, take place in an organized fashion.
The narrow range of $m_{\rm i}$ values detected for the spheroids indicate that there are not tiny bulges or massive spheroids in our simulations. Hence, our results should be verified using larger and higher-resolution simulations.

\begin{figure}
\resizebox{8cm}{!}{\includegraphics{Eta_bd_rhm_md_z_D.eps}}
\hspace*{-0.2cm}
\caption{Angular momentum  efficiency  $\eta_{\rm i}$  as a function of  $m_{\rm i}$ ratio  for  disc (squares) and central spheroid (circles) components
 at $z=0$ (red), $z=1$ (green) and, $z=2$ (blue). }
\label{fig:Eta_D}
\end{figure}

\subsection{The specific angular momenta of galaxies and their dark matter haloes}
\label{sec:angspeci}

In this section, we analyse in detail the $j$ of  the  disc and  central spheroids in relation to that of their dark matter haloes.
In Fig.~\ref{fig:jMvir_z0_D}, we show $j$   for the stellar disc and spheroidal components as a function of virial mass of their host galaxies. They define parallel relations consistent with  $ j \propto M^{2/3}$. This is in agreement with  the expected theoretical relation
derived  under the hypothesis that the dark matter halo spin parameter is constant.
For comparison purposes, we have also included the estimations for the gaseous disc components and the simulated dark matter haloes which
agree well with the expected theoretical relation under the hypothesis of angular momentum conservation (solid line).

Both, the spheroidal and disc components determine relations with the same slopes, which indicates that the stellar spheroids have formed from material of lower content of angular momentum (by a similar fraction), regardless of the virial mass of their host galaxy (or the relative fraction of stellar mass accumulated in the central region as shown in Fig.~\ref{fig:jdmd_D}). This trend is similar to that reported by \citet{genel2015}. This agreement is greatly encouraging considering the differences in the codes, sub-grid physics and criteria adopted to classify the simulated galaxies.

\begin{figure}
\resizebox{8cm}{!}{\includegraphics{Jspecific_rhm_Mvir_z0_D.eps}}
\hspace*{-0.2cm}
\caption{Medians of the specific stellar angular momentum of the gaseous (green line) and stellar  disc (red, dashed line) components and the central stellar spheroidal components (pink, dotted line) as a
function of the virial mass of their host haloes. For comparison purposes, the  relation estimated  for the dark matter haloes (black dotted line)  and the expected theoretical  relation (black solid line) have been also included. }
\label{fig:jMvir_z0_D}
\end{figure}

\subsection{Evolution with redshift}
\label{sec:redshift}

A question naturally arises about the possible evolution of $j$ in the disc and central spheroidal components with redshift. This might imply
a change in the slope or fraction of lost angular momentum,  or both, as  galaxies evolve.
To explore this, we estimated  the $j$ values for the stellar discs and central spheroids as a function of $M_{\rm vir}$  for $z=1$ and $z=2$. Recall that at each  redshift,  we  apply the dynamical decomposition explained in Section 2.1.

 As shown in Fig.~\ref{fig:jBD_z_D}, we found no significant
 evolution of these relations since $z \sim 2$. The gap between the $j$ of rotational-dominated and dispersion-dominated components remains approximately constant, implying a lost  of $80-90\%$ of the initial angular momentum content by the material which ends up in the spheroidal component, regardless of  the virial mass of the halo or of  the redshift.

This is a remarkable trend which suggests again that, independently of the assembly history of galaxies,
their disc and bulge components will contain similar amounts of angular momentum as they evolve.

Note that in Fig. \ref{fig:jBD_z_D} we show the $j$ for the stellar disc and bulge components, independently of the global morphology of the galaxies. Hence, some of the simulated galaxies are
disc-dominated while others are spheroid-dominated systems. So, their total $j$ will be distributed filling the gap and determining a steeper mean relation as shown in
in Fig. ~\ref{fig:jdmd_D}. Therefore, these trends of conservation of the $j$ by rotational-dominated components and equally lost of $j$ by dispersion-dominated ones will be erased as the gap is filled with the combination of systems. The separation
of discs and central spheroids  is crucial
to get a deeper insight into the physical problem. Nevertheless, we can not explain the origin of the similar angular momentum
efficiency for central spheroids yet.

\begin{figure}
\resizebox{8cm}{!}{\includegraphics{Jspeci_bd_rhm_md_z_D.eps}}
\hspace*{-0.2cm}
\caption{Specific stellar angular momentum for the stellar disc (empty triangles) and central spheroid (empty circles) components as a function of
the virial mass for $z=0$ (red), $z=1$ (green) and, $z=2$ (blue). The gap between the two dynamical components does not evolve with redshift.
 For comparison purposes $j$ for the  the dark matter haloes have been also included (filled squares). Note that the
angular momenta of the discs and central spheroids are given at $r_{\rm HM}$.}
\label{fig:jBD_z_D}
\end{figure}

\subsection{Comparison with observations}
\label{sec:obser}

Recently, \citet{Fall2013} reported results which support the existence of a relation between the $j$ and stellar mass with the same slope for both ellipticals and disc dominated-galaxies, as found by \citet{fall1980}. The gap between them indicates a lost of angular momentum by the material which formed the ellipticals by a similar fraction, regardless of stellar mass.
Our simulations seem to reproduce naturally this trend as can be deducted from Fig. ~\ref{fig:jdmd_D} and more clearly seen from
Fig.~\ref{fig:jMstar_DT_D_Fall2} where the observational estimations from  \citet{Fall2013} are included.
As shown in this figure, the stellar disc component determine a relation which agrees well with the one estimated for observed discs.
In the case of central spheroids, they are also consistent with observations provided that stellar bulges behave
similarly to elliptical galaxies \citep[as also suggested by][]{Fall2013}. The good agreement between the
observations and simulations suggests that the
latter are able to grasp  well enough relevant aspects of the physics at play.

In order to improve the comparison of our results with observations, we estimated the total $j$ defined by the whole stellar
populations within the galaxy radius\footnote{The galaxy radius is defined as the one which encloses $\sim 80 \% $ of the baryonic mass.}. Previous figures show results estimated at $r_{\rm HM}$.
The morphology of a simulated galaxy is defined  by the ratio  $D/T$ between the stellar disc mass, $D$,  and the total stellar mass of the galaxy, $T=D+B$.
Simulated galaxies are classified according to the $D/T$ in four subsamples from disc-dominated ($D/T > 0.6$) to spherodial-dominated  ($D/T < 0.2$) ones.

  We should keep in mind that the comparison should be taken as indicative since  differences between the observed and simulated parameters are always expected. For example, while observations use the gas information to trace the angular momentum content of the stellar discs under the hypothesis that they should be dynamical comparable, the simulated calculations make use of the information store directly in the stellar populations.  As we show below, the gaseous discs conserve their angular momentum and they compare well to the relation determined by Sc galaxies, where the dispersion-dominated component is small. This indicates that this type of galaxies is missing in the simulations.

In Fig. \ref{fig:jMstar_DT_D}, we show $j$ as a function of galaxy stellar mass with different $ D/T$ ratios. A correlation can be seen in the sense that higher $D/T$ ratios are related to higher contents of specific angular momentum, as expected. It also presents a good agreement with the observational trend found by \citet{Fall2013}. The gap between disc-dominated systems ($D/T > 0.6$)
 and spheroidal-dominated ones ($D/T < 0.2$) is filled with the corresponding intermediate morphological types.
Our results support the claim that this relation opens the possibility to an alternative way of mapping galaxies' morphologies.
However, and more importantly, the agreement between the observational and numerical results suggests
that the hierarchical clustering signature is not erased in spheroidal galaxies as previously assumed.

It is remarkable that the gap between both morphological types does not depend on the  mass of the galaxy or the virial mass of the haloes. In \citet{Pedrosa2014} a detailed study of the $D/T$ fraction was done. The evolution of the morphology was quantified by analysing the D/T ratios as a function of the redshift for  a group of 15 selected galaxies and the correlation of that parameter with their assembly history. As galaxies evolve, the $D/T$ ratios of the progenitor systems can drastically change, reflecting their different formation histories (major and minor mergers, accretions, etc). The $D/T$ fraction is a highly dynamic parameter, which indicates that the disc can be destroyed and re-built several times along its evolution, as shown in Figure 4 of \citet{Pedrosa2014}.

During galaxy assembly, angular momentum is exchanged through different physical processes. Previous works such as \citet{brook2011} have found that preferential remotion of low angular momentum material at early redshift also contributes to the redistribution of the angular momentum. They found that low angular momentum material is accreted and rapidly expelled when star formation peaks at early times, with shallow potential wells that favour early gas ejection. Later accreted material has higher angular momentum and is accreted preferentially onto the disk, contributing to re-build it. Also, gas might flow into galaxies through the dark matter filaments, feeding the galaxy with high angular momentum material, as found by \citet{dekel2009}. \citet{genel2015} studied, using a more observational approach, the angular momentum content of galaxies in the Illustris simulations. They found, in agreement with our results, that early-type galaxies retention factor is approximately 30\%. They also analised the evolution track of a group of selected galaxies, finding that different processes, such as mergers, could drastically redistribute the angular momentum content.

Therefore, bulges are formed by important contributions of low  angular momentum material which might have lost it during the galaxy formation process \citep{Pedrosa2014}.

\begin{figure}
\resizebox{8cm}{!}{\includegraphics{Jspecific_Mstar_z0_Fall_fig2b.eps}}
\hspace*{-0.2cm}
\caption{Relation between $j$  and the total stellar mass for the simulated stellar disc components (blue, filled squares) and the
stellar spheroids (red, open circles).  The observational estimations reported by \citet{Fall2013} are depicted for comparison:  disc components (blue, open squares) and Elliptical galaxies (red filled circles). The theoretical relation $j \propto M^{2/3}$ is also included for comparison (black, continuous line). }
\label{fig:jMstar_DT_D_Fall2}
\end{figure}

\begin{figure*}
\centering
\resizebox{10cm}{!}{\includegraphics{Jspecific_Mstar_Fall_fig3.eps}}
\hspace*{-0.2cm}
\caption{Relation between $j$  and the total stellar mass for simulated galaxies with different disc-to-bulge ratios: $D/T > 0.6$ (red, filled triangles), $0.4 < D/T < 0.6$ (green, filled triangles), $0.4 < D/T <0.2$ (blue, filled triangles) and  $D/T < 0.2$ (pink, filled triangles) calculated at the galaxy radius. The relation for the gaseous discs are also included (cyan crosses).The observational estimations reported by \citet{Fall2013} are depicted for comparison:  Sa-Sab (black triangles), Sb-Sbc (black squares), Sc-Sm (black diamonds), E (black filled circles), S0 (black filled squares, sE (black circles). The theoretical relation $j \propto M^{2/3}$ is also included (black, continuous line). }
\label{fig:jMstar_DT_D}
\end{figure*}

\section{Size relations}
\label{sec:size}

As shown in previous sections, the disc components in our simulations were formed in agreement with the hypothesis of  $j$ conservation, while the central spheroids were formed from material with lower content of $j$, by a similar fraction. From these results, we would expect a clear scaling relation between the galaxy and halo size.

The proportionality between the specific angular momentum of the baryons and that of  the halo assumed by the disc galaxy formation models would also imply  the existence of a linear relation between the galactic and virial radius of the hosting halo \citep{kravtsov2013}. The galaxy size is expected to scale as $r_{\rm gal} \propto \lambda r_{200}$ .

\citet{kravtsov2013} analysed the relation between sizes of stellar systems of galaxies characterized by their half-mass radius and the virial radius of their dark matter haloes. The latter are derived by  using the abundance matching ansatz. They found that the characteristic size of stellar and gas distributions in galaxies spanning several orders of magnitude in stellar mass, scales approximately linearly with a fraction of the the virial radius with a proportionality factor of  $\sim 0.015$. However, this author did not consider the effects of baryons on the dark matter potential well, which are naturally provided by our hydrodynamical simulations \citep[e.g.][]{Pedrosa2010, tissera2010}.   In fact, \citet{diemer2013}  claimed that due to the most common definition of the dark matter halo that considers a spherical overdensity of matter with respect to a reference density, such as the critical one or the mean matter density of the Universe, could lead to a spurious pseudo-evolution of halo mass due to
redshift evolution of the reference density. They found that for all halo masses a significant fraction of the halo mass growth since $z \sim 1$ is due to this so-called pseudo-evolution rather than the actual physical accretion of matter. They obtained a corrected value for the size relation of $0.03 - 0.04$.
 by taking the contraction effect into account.

In Fig.~\ref{fig:Rhm_rvir_fracD} (upper panel) we show the stellar $r_{\rm HM}$ as a function of a fraction of $r_{200}$ for our simulated systems. Interestingly, we determined that the best factor relating both characteristic radii is $\sim 0.04$.
In  this figure, we show the relation for galaxies with different $D/T$ ratios. It is clear that higher $D/T$ values correspond to more extended structures. And this is consistent with Fig.~\ref{fig:jMstar_DT_D} where it can be seen that also higher $D/T$ ratios correspond to higher values of  angular momentum content.
The size correlation is at place since $z \sim 2$ as expected since the angular momentum relation is also present and unchanged globally
since this time Fig.~\ref{fig:Rhm_rvir_fracD}(lower panel).

\begin{figure}
\resizebox{8cm}{!}{\includegraphics{Rhm_rvir_fracD_D_II.eps}}
\resizebox{8cm}{!}{\includegraphics{Rhm_rvir_D_z.eps}}
\hspace*{-0.2cm}
\caption{Size relations for the simulated galaxies and their host dark matter haloes: $r_{\rm HM} =0.04 \times \lambda r_{200}$. Upper panel: Simulated galaxies are colour-coded according to  the $D/T$ ratio:  $D/T > 0.6$ (red, filled squares), $0.4 < D/T < 0.6$ (green, filled diamonds), $0.4 < D/T <0.2$ (blue, crosses) and  $D/T < 0.2$ (pink, filled circles). Lower panel: The colours correspond to the different redshifts, $z = 0$ (red circles), $z = 1$ (green squares), $z = 2$ (blue triangles)}
\label{fig:Rhm_rvir_fracD}
\end{figure}

\section{Conclusions}
\label{sec:conclu}

We used hydrodynamical cosmological simulation to perform a comparative analysis of the angular momentum content in the disc and bulge components of the simulated galaxies.
We focused our analysis on the the best resolved galaxies and used a dynamical criterion to identify the disc and bulge components at $z=0$ and for their progenitor galaxies at $z=1$ and $z=2$.

Our key findings are:

\begin{itemize}
\item{} We found  that both discs and bulges determine a relation in  the $j_{\rm i}-m_{\rm i}$  plane. Both the rotational-supported  and the dispersion dominated components follow a relation with the same slope ($j_{\rm i} \propto m_{\rm i}$). For the disc components, the correlation is the expected one-to-one relation in the case of angular momentum conservation \citep{MMW98}. Remarkably, the dispersion-dominated systems determine a parallel relation but off-set by an approximately constant factor.  Our simulated bulges  conserved $\eta \sim 10-20\% $ of their original angular momentum, regardless of the $D/T$ ratio of the  galaxy \citep[see also ][]{dominguez1998,Pedrosa2014}.  This is in agreement with observational findings \citep[][and referencies there in]{Fall2013}.

\item A similar analysis was carried out for the progenitor systems up to $z \sim  2$, finding that the central spheroids and discs determine similar relations. Hence,   regardless of the formation history of these galaxies, their disc and spheroidal components have angular momentum contents which determine relations with the  stellar and dark matter masses with similar slope up to $z \sim 2$. In the case of the  spheroidal components a similar off-set is detected showing an angular momentum efficiency of $\sim 0.10 - 0.20$. These values
vary slightly with redshift. A larger sample is required to assess robustly the level of evolution.

\item Regarding  the specific angular momentum, we found that the spheroidal and disc components determine relations with the same
slopes, $ j \propto M^{2/3}$. This indicates that the stellar spheroids have been formed by material of low angular momentum content or that have lost it by similar fraction during the assembly process, regardless of the virial mass of their host galaxy. We have analysed this relation at different redshifts, finding no evolution. Remarkably, this trend is suggesting that independently of virial masses, the disc and bulge components will conserve similar amounts of angular momentum as they evolve.

\item In good agreement with the trend found by \citet{Fall2013}, we found a clear correlation between the morphological type of a galaxy and its total specific angular momentum content: higher D/T ratios are related with higher contents of specific angular momentum, as expected. The gap between disc-dominated systems ($D/T > 0.6$) and spheroidal-dominated ones ($D/T < 0.2$) is filled with the corresponding intermediate morphological types. This supports the claim that this relation opens the possibility to an alternative way of mapping the morphology of galaxies. And this also implies that the hierarchical clustering signature is still present in spheroidal galaxies. It is remarkable that the gap between both morphological types does not depend on the mass of the galaxy or the virial mass of the haloes.

\item In agreement with \citet{kravtsov2013}, we found that the characteristic size of stellar distributions in the simulated galaxies scales approximately linearly with the virial radius in the log-log plane. The proportionality factor for our simulated systems is $\sim 0.04$ \citep{diemer2013}. This relation is found not to evolve with time: galaxies move along the relation as one moves to higher redshift.

\end{itemize}

\begin{acknowledgements}
The authors thanks Michael Fall for kindly providing the data electronically.
This work was partially supported by PICT 2011-0959 and PIP 2012-0396 (Mincyt, Argentina). PBT acknowledges partial support from  the Regular Grant UNAB 2014, Nucleo UNAB 2015 of Universidad Andres Bello and Fondecyt 1150334.
\end{acknowledgements}

\bibliographystyle{aa}