The paper analyzes the Morley element method for the Cahn-Hilliard equation. The objective is to derive the optimal error estimates and to prove the zero-level sets of the Cahn-Hilliard equation approximate the Hele-Shaw flow. If the piecewise L ∞ (H 2 ) error bound is derived by choosing test function directly, we cannot obtain the optimal error order, and we cannot establish the error bound which depends on 1 polynomially either. To overcome this difficulty, this paper proves them by the following steps, and the result in each next step cannot be established without using the result in its previous one. First, it proves some a priori estimates of the exact solution u, and these regularity results are minimal to get the main results; Second, it establishes L ∞ (L 2 ) and piecewise L 2 (H 2 ) error bounds which depend on 1 polynomially based on the piecewise L ∞ (H −1 ) and L 2 (H 1 ) error bounds; Third, it establishes piecewise L ∞ (H 2 ) optimal error bound which depends on 1 polynomially based on the piecewise L ∞ (L 2 ) and L 2 (H 2 ) error bounds; Finally, it proves the L ∞ (L ∞ ) error bound and the approximation to the Hele-Shaw flow based on the piecewise L ∞ (H 2 ) error bound. The nonstandard techniques are used in these steps such as the generalized coercivity result, integration by part in space, summation by part in time, and special properties of the Morley elements. If one of these techniques is lacked, either we can only obtain the sub-optimal piecewise L ∞ (H 2 ) error order, or we can merely obtain the error bounds which are exponentially dependent on 1 . The approach used in this paper provides a way to bound the errors in higher norm from the errors in lower norm step by step, which has a profound meaning in methodology. Numerical results are presented to validate the optimal L ∞ (H 2 ) error order and the asymptotic behavior of the solutions of the Cahn-Hilliard equation. m 0 := 1 |Ω| Ω u 0 (x) dx.(2) There exists a nonnegative constant σ 1 such that