Nonstabilizerness in open XXZ spin chains: Universal scaling and dynamics

Magic, or nonstabilizerness, is a crucial quantum resource, yet its dynamics in open quantum systems remain largely unexplored. We investigate magic in the open XXZ spin chain under either boundary gain and loss or bulk dephasing using the stabilizer Rényi entropy M_{2}. To enable scalable simulations of large systems, we develop a unique, highly efficient algorithm for computing M_{2} within the matrix product state formalism while maintaining constant bond dimension—an advancement over existing methods. For boundary driving, we uncover universal scaling laws, M_{2}(t)∼t^{1/z}, linked to the dynamical exponent z for several distinct universality classes. We also disentangle classical and quantum contributions to magic by introducing a mean-field approximation for magic, thus emphasizing the prominent role of quantum-critical fluctuations in nonstabilizerness. For bulk dephasing, dissipation can transiently enhance magic before suppressing it and drive it to a nontrivial steady-state value. These findings position magic as a powerful diagnostic tool for probing universality and dynamics in open quantum systems.