Improved upper bounds on the stabilizer rank of magic states

In this work we improve the runtime of recent classical algorithms for strong simulation of quantum circuits composed of Clifford and T gates. The improvement is obtained by establishing a new upper bound on the stabilizer rank ofmcopies of the magic state|T⟩=2−1(|0⟩+eiπ/4|1⟩)in the limit of largem. In particular, we show that|T⟩⊗mcan be exactly expressed as a superposition of at mostO(2αm)stabilizer states, whereα≤0.3963, improving on the best previously known boundα≤0.463. This furnishes, via known techniques, a classical algorithm which approximates output probabilities of ann-qubit Clifford + T circuitUwithmuses of the T gate to within a given inverse polynomial relative error using a runtimepoly(n,m)2αm. We also provide improved upper bounds on the stabilizer rank of symmetric product states|ψ⟩⊗mmore generally; as a consequence we obtain a strong simulation algorithm for circuits consisting of Clifford gates andminstances of any (fixed) single-qubitZ-rotation gate with runtimepoly(n,m)2m/2. We suggest a method to further improve the upper bounds by constructing linear codes with certain properties.