P.S. I asked the question on MSE more than a week ago, but didn't get any desired answer, so asking here.
Let $m < n \in \mathbb{N}$. Let us equip $\mathbb{R}^m, \mathbb{R}^n $ with their canonical Euclidean (Riemannian) metrics. How can we characterize the extrinsic distance preserving embeddings of $\mathbb{R}^m$ into $\mathbb{R}^n $? To be more precise, I'm looking for sufficiently regular, injective, distance preserving transformations from $\Phi: \mathbb{R}^m \to \mathbb{R}^n$, so that: $d_m(x, y)= d_n(\Phi(x), \Phi(y))$, where $d_m, d_n$ represents the distances in $m, n$ dimensional Euclidean spaces respectively, so $d_m(x,y):=||x-y||_{\mathbb{R}^m}, d_n(x,y):=||x-y||_{\mathbb{R}^n}$. Note that this means: we're looking to preserve the extrinsic distances between the two Euclidean spaces, we don't care about if the first one embeds isometricaly as a submanifold in the second.
I think the answer is:
$$x \mapsto A(Bx, O(Bx))$$ where $A$ is an Euclidean isometry of $\mathbb{R}^n$ (i.e. a rigid motion), $B$ is an Euclidean isometry of $\mathbb{R}^m$ (i.e. a rigid motion), and $O:\mathbb{R}^m \to \mathbb{R}^n$ is "zero padding $n-m$" times, namely: $O(x)= (x, 0,\dots 0)$.
If the above correct/incorrect, how do I go about proving it or characterizing the Euclidean embeddings?
Approaches that immediately come to mind are as follows:
(1) Assume $\phi: \mathbb{R}^m \to \mathbb{R}^n$ is such an isometric embedding. Then we can try to consider the canonical projection $\pi: \mathbb{R}^n \to \mathbb{R}^m$ and try to consider the map: $\phi \circ \pi: \mathbb{R}^n \to \mathbb{R}^n $, but this will not be an isometry (not one to one), which is a problem.
(2) Alternately, we can consider: $\pi \circ \phi: \mathbb{R}^m \to \mathbb{R}^m$. We can't guarantee that $\pi \circ \phi$ is an isometry of $\mathbb{R}^m$, which is again a problem.
(3) As someone suggested in the comments of the MSE question, there exists an isometry, say $A$, of $\mathbb{R}^n$ so that $A \circ \phi (\mathbb{R}^m)= \mathbb{R}^m \times \{0\}^{n-m}$. This intuitively seems correct, as the isometry group of $\mathbb{R}^n$ has dimension $\frac{n(n+1)}{2}$, which means it's "high dimensional enough", i.e. "there're enough elements in it" to render the last $n-m$ components to $0$. This could be made rigorous if $\phi$ was a linear or affine map, but if $\phi$ is nonlinear, I can't see how this argument goes through? In this case we need to show first that $\phi$ is affine, i.e. if $\phi: \mathbb{R}^m \to \mathbb{R}^n$ fixes $0$, then it has to be a linear map from $\mathbb{R}^m$ to $\mathbb{R}^n$.(?)