If $f: \mathbb{R}^n \to \mathbb{R}$ is a Lipschitz function and $X$ is a standard $n$-dimensional Gaussian vector with $\mathbb{E} f(X) = 0$, then $f(X)$ is subgaussian (in a way that does not depend on $n$). If $f$ is $\mathcal{C}^1$, this is equivalent to saying that $|\nabla f|$ bounded implies $f(X)$ is subgaussian.
There seem to be two natural generalizations of this. The first is to ask for weaker bounds on $|\nabla f|$. For example, if $|\nabla f|$ is subgaussian, then $f$ should be subexponential. The second generalization concerns functions $f: \mathbb{R}^n \to \mathbb{R}^k$. If I want to control $|f|$ independently of $k$, it is no longer enough to assume that $f$ is Lipschitz, since for the function $f(x) = (x_1, \dots, x_k)$, $|f|$ concentrates around $\sqrt k$. The natural condition seems to be a bound on the Frobenius norm of $D f$ (the matrix of partial derivatives).
The following statement contains both generalizations simultaneously (and is not hard to prove): If $f: \mathbb{R}^n \to \mathbb{R}^k$ is continuously differentiable and $\mathbb{E} f(X) = 0$ then $$ \big(\mathbb{E} |f(X)|^p\big)^{1/p} \le c \sqrt p \big(\mathbb{E} \|Df\|_F^p\big)^{1/p}. $$
My question is whether a statement like this is known and (if so) where I can find a reference.