We consider a nonlinearly coupled electromechanical system, and develop quantitative theory for two-phonon cooling. In the presence of cooling, mechanical Hilbert space is effectively reduced to its ground first excited states, allowing quantum operations at level individual phonons preparing nonclassical states with negative Wigner functions. propose scheme performing arbitrary Bloch sphere ro...