
$Q$systems and $T$systems are systems of integrable difference equations
that have recently attracted much attention, and have wide applications in
representation theory and statistical mechanics. We show that certain
$\tau$functions, given as matrix elements of the action of the loop group of
${\rm GL}_{2}$ on twocomponent fermionic Fock space, give solutions of a
$Q$system. An obvious generalization using the loop group of ${\rm GL}_3$
acting on threecomponent fermionic Fock space leads to a new system of 4
difference equations.

We introduce hierarchies of difference equations (referred to as
$nT$systems) associated to the action of a (centrally extended, completed)
infinite matrix group $GL_{\infty}^{(n)}$ on $n$component fermionic Fock
space. The solutions are given by matrix elements ($\tau$functions) for this
action. We show that the $\tau$functions of type $nT$ satisfy bilinear
equations of length $3,4,\dots,n+1$. The $2T$system is, after a change of
variables, the usual $3$ term $T$system of type $A$.
Restriction from $GL_{\infty}^{(n)}$ to a subgroup isomorphic to the loop
group $LGL_{n}$, defines $nQ$systems, studied earlier by the present authors
for $n=2,3$.

We briefly describe what taufunctions in integrable systems are. We then
define a collection of taufunctions given as matrix elements for the action of
$\widehat{GL_2}$ on twocomponent Fermionic Fock space. These taufunctions are
solutions to a discrete integrable system called a $Q$system.
We can prove that our taufunctions satisfy $Q$system relations by applying
the famous "DesnanotJacobi identity" or by using "connection matrices", the
latter of which gives rise to orthogonal polynomials. In this paper, we will
provide the background information required for computing these taufunctions
and obtaining the connection matrices and will then use the connection matrices
to derive our difference relations and to find orthogonal polynomials.
We generalize the above by considering taufunctions that are matrix elements
for the action of $\widehat{GL_3}$ on threecomponent Fermionic Fock space, and
discuss the new system of discrete equations that they satisfy. We will show
how to use the connection matrices in this case to obtain "multiple orthogonal
polynomials of type II".