  In this article, the stabilizer free weak Galerkin (SFWG) finite element
method is applied to the Ciarlet-Raviart mixed form of the Biharmonic equation.
We utilize the SFWG solutions of the second elliptic problems to define
projection operators, build error equations, and further derive the error
estimates. Finally, numerical examples support the results reached by the
theory.
