I have two states, $|\psi\rangle$ and $|\phi\rangle$. I have in mind that they live on a length $L$ spin chain with finite local Hilbert space dimension.
I know that for every Schmidt decomposition bipartitioning the system into a region $A$ and a region $B$, the following is true:
$$|\psi\rangle = \sum_{i=1}^n \lambda_i|a_i\rangle|b_i\rangle $$ $$|\phi\rangle = \sum_{i=1}^n \lambda_i|a'_i\rangle|b'_i\rangle $$
That is, while the states are not necessarily equal, their Schmidt coefficients are equal across every cut in physical space. The Schmidt coefficients may depend on the choice of cut. I imagine that all of my cuts are in physical space, but I do allow regions $A$ and $B$ to contain sites that are not contiguous; for example, $A$ could contain all even sites and $B$ can contain all odd sites.
Given this, am I guaranteed that there exists a unitary that is a tensor product of single-site unitaries, $U = \otimes_{i=1}^L U_i$, such that $ U|\psi\rangle = |\phi\rangle$?
Here are my thoughts. If I had the weaker statement that the Schmidt coefficients were equal $$|\psi\rangle = \sum_{i=1}^n \lambda_i|a_i\rangle|b_i\rangle $$ $$|\phi\rangle = \sum_{i=1}^n \lambda_i|a'_i\rangle|b'_i\rangle $$ for some specific regions $A$ and $B$, then I know I can make a unitary $U = U_A \otimes U_B$ with $U_A = \sum_i |a'_i\rangle \langle a_i|$ and $U_B = \sum_i |b'_i\rangle \langle b_i|$ that takes $|\psi\rangle$ to $|\phi\rangle$.
However, it's less clear to me how to use the information from all of the bipartite cuts together. I was thinking I could consider all $L$ cuts where $A$ contains a single site, and then attempt to argue that $U$ can be written as a product of single-site operators, but I'm not sure that will work. In particular, I'm getting suspicious that perhaps all of the bipartite cuts aren't enough, and that I'll need to know things about multipartite decompositions of the state.