Phase-space representations as given by Wigner functions are a powerful tool for representing the quantum state and characterizing its time evolution in the case of infinitedimensional quantum systems and have been widely used in quantum optics and beyond. Continuous phase spaces have also been studied for finite-dimensional quantum systems such as spin systems. However, much less is known for finite-dimensional, coupled systems, and we present a complete theory of Wigner functions for this case. In particular, we provide a self-contained Wigner formalism for describing and predicting the time evolution of coupled spins which lends itself to visualizing the high-dimensional structure of multi-partite quantum states. We completely treat the case of an arbitrary number of coupled spins 1 2, thereby establishing the equation of motion using Wigner functions. The explicit form of the time evolution is then calculated for up to three spins 1 2. The underlying physical principles of our Wigner representations for coupled spin systems are illustrated with multiple examples which are easily translatable to other experimental scenarios.
Prior work on Wigner functions of finite-dimensional quantum systemsFundamental postulates for phase-space representations of finite-dimensional quantum systems were proposed by Stratonovich [21] (see Appendix C.1), and these build the foundations for Wigner functions of finite-dimensional quantum systems. Reflecting the translational covariance of Wigner functions of infinite-dimensional quantum systems, one of these postulates states that the Wigner function has to transform naturally under rotations. The rotational covariance constrains continuous Wigner representations of spins into a spherical coordinate systems. The resulting Wigner functions can then be given by linear combinations of spherical harmonics, which offers a convenient tool for visualizing spins (see Sec. 1.4).The Wigner transformation of single-spin operators was developed by Várilly and Gracia-Bondía [22] and was then further extended by Brif and Mann [23,24]. In particular, [22] provides an explicit formula for the Wigner transformation for a single spin J, which satisfies the Stratonovich postulates. This formula uses a rank-j-dependent kernel which is based on mapping transition operators Jm⟩⟨Jm ′ onto their corresponding Wigner functions W Jm⟩⟨Jm ′ ; the connection between tensor operators J T jm and spherical harmonics Y jm was also mentioned. A more general kernel for the continuous phase-space representation of a single spin was stated in [23,24]. It defines Wigner functions of tensor operators J T jm of single spins as the corresponding spherical harmonics Y jm . We build on these results in Sec. 3.2 while also unifying normalization factors.Parallel to our work, a general approach for phase-space representations was proposed in [27] which is based on the so-called displaced parity operator [28]. The explicit form of the transformation kernel is computed for the special cases of a single spin J (see Eq. (8) in ...