  We show that spherical Whittaker functions on an $n$-fold cover of the
general linear group arise naturally from the quantum Fock space representation
of $U_q(\widehat{\mathfrak{sl}}(n))$ introduced by Kashiwara, Miwa and Stern
(KMS). We arrive at this connection by reconsidering solvable lattice models
known as `metaplectic ice' whose partition functions are metaplectic Whittaker
functions. First, we show that a certain Hecke action on metaplectic Whittaker
coinvariants agrees (up to twisting) with a Hecke action of Ginzburg,
Reshetikhin, and Vasserot. This allows us to expand the framework of KMS by
Drinfeld twisting to introduce Gauss sums into the quantum wedge, which are
necessary for connections to metaplectic forms. Our main theorem interprets the
row transfer matrices of this ice model as `half' vertex operators on quantum
Fock space that intertwine with the action of
$U_q(\widehat{\mathfrak{sl}}(n))$.
  In the process, we introduce new symmetric functions termed
\textit{metaplectic symmetric functions} and explain how they relate to
Whittaker functions on an $n$-fold metaplectic cover of GL$_r$. These resemble
\textit{LLT polynomials} introduced by Lascoux, Leclerc and Thibon; in fact the
metaplectic symmetric functions are (up to twisting) specializations of
\textit{supersymmetric LLT polynomials} defined by Lam. Indeed Lam constructed
families of symmetric functions from Heisenberg algebra actions on the Fock
space commuting with the $U_q(\widehat{\mathfrak{sl}}(n))$-action. We explain
that half vertex operators agree with Lam's construction and this
interpretation allows for many new identities for metaplectic symmetric and
Whittaker functions, including Cauchy identities. While both metaplectic
symmetric functions and LLT polynomials can be related to vertex operators on
the $q$-Fock space, only metaplectic symmetric functions are connected to
solvable lattice models.
