Increasing tensor powers of the k × k matrices M k (ℂ) are known to give rise to a continuous bundle of C ∗ -algebras over I = { 0 } ∪ 1 / ℕ ⊂ [ 0 , 1 ] with fibers A 1 / N = M k (ℂ) ⊗ N and A 0 = C (X k) , where X k = S (M k (ℂ)) , the state space of M k (ℂ) , which is canonically a compact Poisson manifold (with stratified boundary). Our first result is the existence of a strict deformation quantization of X k à la Rieffel, defined by perfectly natural quantization maps Q 1 / N : à 0 → A 1 / N (where à 0 is an equally natural dense Poisson subalgebra of A 0 ). We apply this quantization formalism to the Curie–Weiss model (an exemplary quantum spin with long-range forces) in the parameter domain where its ℤ 2 symmetry is spontaneously broken in the thermodynamic limit N → ∞. If this limit is taken with respect to the macroscopic observables of the model (as opposed to the quasi-local observables), it yields a classical theory with phase space X 2 ≅ B 3 (i.e. the unit three-ball in ℝ 3 ). Our quantization map then enables us to take the classical limit of the sequence of (unique) algebraic vector states induced by the ground state eigenvectors Ψ N (0) of this model as N → ∞ , in which the sequence converges to a probability measure μ on the associated classical phase space X 2 . This measure is a symmetric convex sum of two Dirac measures related by the underlying ℤ 2 -symmetry of the model, and as such the classical limit exhibits spontaneous symmetry breaking, too. Our proof of convergence is heavily based on Perelomov-style coherent spin states and at some stage it relies on (quite strong) numerical evidence. Hence the proof is not completely analytic, but somewhat hybrid. [ABSTRACT FROM AUTHOR]