$0\nu\beta\beta$ to the first $2^+$ state with two-nucleon mechanism for L-R symmetric model

We develop the formalism for calculating the decay rate of neutrinoless double beta decay to the $2^+$ excited states within L-R symmetric model. We consider the effects from induced hadronic currents up to NLO. The QRPA method in a spherical basis is adopted for the nuclear many-body calculation and the corresponding nuclear matrix elements are given. Also, the phase space factors are obtained with numerical electron wave functions. Our results suggest that the nuclear matrix elements are nucleus dependent and they are generally smaller than that of the decay to the ground states. And finally, we give a naive analysis of how current experiment data constrains the L-R symmetric model.