Bilinear R-parity violation in neutrinoless double beta decay

We discuss some phenomenological issues of the effective quark-lepton operators emerging from the bilinear lepton-Higgs couplings in the superpotential and in the soft supersymmetry breaking sector of the supersymmetric models with explicit R-parity violation. The contribution of these operators to the neutrinoless double beta decay (0νββ) is derived. The corresponding nuclear matrix elements are calculated within the renormalized quasiparticle random phase approximation, which includes the Pauli effect of fermion pairs and does not collapse for the physical values of the nuclear force strength. On this basis we extract from the experimental data new stringent limits on the first generation mass parameter characterized the lepton-Higgs bilinear coupling and on the electron sneutrino vacuum expectation value.