Extended operator expansion method for neutrinoless double beta decay