Here's another answer to the question about multiplicities which can be generalized to other partitions. First, for any partition $\lambda$, let $S_\lambda$ denote the irreducible representation (Schur functor) with highest weight $\lambda$. I'll explain how to get the following:
Theorem: Let $\lambda$ be a partition of even size and let $\mu$ be the result of scaling all column lengths of $\lambda$ by $k$. Then $Sym^k(S_\lambda({\bf C}^n))$ contains $S_\mu({\bf C}^n)$ with multiplicity 1 if $\mu$ has at most $n$ rows (otherwise $S_\mu({\bf C}^n)=0$).
For $\lambda = (2,2)$, scaling all column lengths by $k$ just means the partition $(2,2,2,\dots,2)$ ($2k$ instances of 2).
I'm going to assume basic fluency between symmetric functions and Schur functors.
Then $S_\lambda \otimes S_\mu$ always contains $S_{\lambda + \mu}$ with multiplicity 1: the tensor product of highest weight vectors gives the unique vector (up to scalar multiple) of this weight and it's also a highest weight vector (easy check). In particular, $S_{\lambda}^{\otimes k}$ contains $S_{k\lambda}$ with multiplicity 1. From the product description, the highest weight vector is symmetric so it also belongs to $Sym^k(S_\lambda)$.
Next, we use the involution $\omega$ on symmetric functions. One reference Section I.2 of Macdonald, Symmetric Functions and Hall Polynomials. This has the property that $\omega s_\lambda = s_{\lambda^T}$ where $\lambda^T$ is the transpose partition. It also behaves well with respect to plethysm (there are several kinds, here I refer to the one that corresponds to composition of Schur functors):
If $f$ and $g$ are homogeneous, then
$\omega (f \circ g) = \omega f \circ \omega g$ if $\deg g$ is odd, and
$\omega (f \circ g) = f \circ (\omega g)$ if $\deg g$ is even.
Reference: Macdonald, Example I.8.1(a).
In particular, this means that if $|\lambda|$ is even, then
$\omega (s_k \circ s_\lambda) = s_k \circ s_{\lambda^T}$.
Translating this back to representations:
The right side is $Sym^k(S_{\lambda^T})$ and hence contains $S_{k\lambda^T}$ with multiplicity 1. In particular, $\omega$ is an involution so $Sym^k(S_\lambda)$ contains $S_{(k\lambda^T)^T}$.
We can simplify further: $(k\lambda^T)^T$ is obtained from $\lambda$ by scaling all column lengths by $k$.