My reference is Infinite series in p-adic fields by Keith Conrad.
Corollary 5.6. If $f(x)=\sum_{n≥0} a_nx^n$ has a positive radius of convergence in the $p$-adic field $\mathbb Q_p$ then $f$ is infinitely differentiable on its disc of convergence in $K$ and for all $n ≥ 0$ and $a_n=\frac{f^{(n)}(0)}{n!}$. In particular, $f$ has only one choice of power series coefficients.
Proof. Let $D$ be the disc of convergence of $f(x)$. By Theorem 5.3 (of the same note), $f'(x)=\sum_{n \geq 0} n a_n x^{n-1}$ for all $x \in D$. This is again a power series, so by Theorem 5.3, it is differentiable on $D$ with derivative $f^{''}(x)=\sum_{n \geq 2}n(n-1)x^{n-2}$. By induction, the $k$th derivative is $$f^{(k)}(x)=n(n − 1) \cdots (n − (k − 1))a_nx^{n−k}$$ for all $x \in D$. In particular, the constant term is $f^{(k)}(0)= k(k−1) \cdots (k−(k−1))a_k= k!a_k$, so $a_k =\frac{f^{(k)}(0)}{k!}$. This shows the coefficients of the power series are determined by $f$ as a function around $0$.
Uniqueness follows from Corollary 3.15 of the same Note, and it is not required to use derivatives of $f$.
My question: Is there a generalization of the above results for vector-tuple power series $f: \mathbb Q_p^n \to \mathbb Q_p^n$?
For example, take the simple case $n=2$ and consider the 2-tuple power series $f(x,y)=(f_1(x,y),f_2(x,y))$ defined in the following way $$f(x,y)=\begin{pmatrix} f_1(x,y) \\ f_2(x,y) \end{pmatrix}=\sum_{m,n \geq 0} \begin{pmatrix} a_{mn}^{11} & a_{mn}^{12} \\ a_{mn}^{21} & a_{mn}^{22} \end{pmatrix} X^{m+n}, $$ where $X^{l+m}=\begin{pmatrix} x^i y^j \\ x^r y^s \end{pmatrix}$ is column vector and $i,j,r,s \geq 0$ such that $i+j=m$ and $r+s=n$. Also note that $a_{mn}^{uv}$ is the $(u,v)$-th entry of the 2-by-2 matrix corresponding to the coefficient of $X^{l+m}$.
Of course, there can be more different expositions of $f(x,y)=(f_1(x,y),f_2(x,y))$, for example: $$f(x,y)=\begin{pmatrix} f_1(x,y) \\ f_2(x,y) \end{pmatrix}=\begin{pmatrix} \sum_{m,n \geq 0} a_{m,n} x^my^n \\ \sum_{m,n \geq 0} a'_{m,n}x^my^n \end{pmatrix};~a_{m,n}, a_{m,n}^{'} \in \mathbb Q_p.$$
Can we generalise Corollary 5.6?
I see the challenge in handling the derivatives of $f(x,y)$. The coefficients $a_{mn}^{11}, a_{mn}^{12}, a_{mn}^{21}, a_{mn}^{22}$ are determined by the mixed partial derivatives of $f(x,y)$ at $(0,0)$.
Can someone shed some light on it? Thanks