Mixed Biharmonic Problem with zero normal derivative on internal facets

Thanks to both of you!

I think the approach with RT and DG spaces could lead to a solution of my problem. However just changing the spaces in the code does not work because an error is thrown. The error is related to dolfinx_mpc.

I made a MWE in this new question where I am going a step back and solve the mixed poisson problem.