Further improvements to stabilizer simulation theory: classical rewriting of CSS-preserving stabilizer circuits, quadratic form expansions of stabilizer operations, and framed hidden variable models

Simulation of stabilizer circuits is a well-studied problem in quantum information processing, with a number of highly optimized algorithms available. Yet, we argue that further improvements can arise from the theoretical structure of stabilizer operations themselves. We focus on the subclass of stabilizer circuits composed of Calderbank-Shor-Steane (CSS)-preserving stabilizer operations, which naturally appear in fault-tolerant computations over CSS stabilizer codes. Using elementary circuit-transformation techniques, we show that such circuits can be exactly rewritten as classical probabilistic circuits that reproduce measurement statistics. This rewriting introduces no computational overhead, in contrast to the general case of stabilizer circuits. To clarify the origin of this simplification, we introduce the standard quadratic-form representation of general stabilizer operations (Clifford channels). It provides an efficient way to describe compositions of stabilizer operations and thus to simulate stabilizer circuits. CSS-preserving operations correspond to purely linear forms, which under a Walsh-Hadamard-Fourier transform yield a noncontextual hidden variable model, providing an alternative proof of the introduced rewriting. Finally, we develop a theory of reference frames for multiqubit systems, where frames are encoded by quadratic forms. This allows us to express stabilizer operations as probabilistic maps for proper reference frames. Non-CSS-preserving stabilizer circuits require dynamical modifications of reference frames, embodying a contextuality resource that leads to the computational overhead. This framework provides a new perspective on simulating stabilizer and near-stabilizer circuits within dynamically evolving quasiprobability models.