Let's define $P_k \in \{ I, X, Y, Z \}^{\otimes n}$ and called each of these $P_k$ as a Pauli string (or word) then given that $$U = \sum_{k=1}^L c_kP_k $$ with the following conditions:
- $\sum_{k=1}^L |c_k|^2 = 1 $
- ${\rm Im}(c_i^* c_j) = 0$ for all $i\neq j$
- $\{P_i, P_j\} = 0$ for all $i\neq j$ where $\{A,B \}:=AB+BA$ is the anti-commutator. In other words, the Pauli strings $P_i$ in the sum that make up $U$ are anti-commute.
With the above three conditions, $U$ is a unitary operator and thus can be implemented as a quantum circuit.
In the following paper at the end of page 3, it shows that $U$ can be implemented as the product of $2L-1$ exponents of these Pauli strings (or words) in the following way:
$$ U = \prod_{k = 1}^L e^{i \theta_k P_k/2} \prod_{k = L}^1 e^{i \theta_k P_k/2} $$
where $\theta_k = \arcsin(\frac{c_k}{\sum_{j=1}^k c_j^2 } ) $. Is there a neat proof of this fact or is it just full out algebraic manipulation which gets really messy? That is, is there an intuitive way to see why this is true? Since the authors seem to indicate this is a pretty trivial fact.
Bonus: What is the best implementation for $U$ in term of the quantum circuit outside of the proposed method above?