In Geometric quantization theory the Moyal product is one of main tools. We know Moyal product for the smooth functions $f$ and $g$ on $ℝ^{2n}$ takes the form
$f\star g = fg + \sum_{n=1}^{\infty} \hbar^{n} C_{n}(f,g)$
where each $C_n$ is a certain bidifferential operator of order $n$ with the following properties and also where $\hbar$ is the reduced Planck constant.
1.$\quad f\star g = fg + \mathcal O(\hbar)$
$\quad f\star g-g\star f = \mathrm i\hbar\{f,g\} + \mathcal O(\hbar^2) \equiv \mathrm i\hbar \{\{f,g\}\}$
$\quad f\star 1=1\star f=f$
4.$ \quad \overline{f\star g} = \overline{g}\star \overline{f}$
My question is how can we find $C_n$. I am looking for a method , not explicit formula based on Poisson bivector. Also when can we use of this product?
Let me give you the naïve picture at least, which I hope partially answers your question; from what I understand, in the actual research literature on deformation quantisation, i.e., in noncommutative geometry à la Kontsevich, one takes a somewhat more flexible but technically more delicate and sophisticated approach to the same basic strategy. In any event, the basic idea is to translate the properties of $\star$ into a system of equations in the $$C_k : C^\infty(\mathbb{R}^{2n},\mathbb{C}) \otimes_{\mathbb{C}} C^\infty(\mathbb{R}^{2n},\mathbb{C}) \to C^\infty(\mathbb{R}^{2n},\mathbb{C}),$$ which one could try to solve recursively.
For simplicity, write $$ C_0(f,g) := fg, \quad C_1(f,g) = i\{f,g\}, $$ so that $$ f \star g = \sum_n \hbar^n C_n(f,g). $$ Let us look at the properties of $\star$ in turn:
So, constructing your deformation quantisation, at least in principle, proceeds inductively: