Abstract :
We analyse inclusive photoproduction of polarized 3P1 quarkonium in the framework of QCD. To separate non-perturbative and perturbative parts in the density matrix of the produced quarkonium we use a method, which is equivalent to the diagrammatic expansion widely used in analyzing deeply inelastic scattering. A systematic expansion in a small velocity v, with which a heavy quark moves inside the quarkonium in its rest frame, is performed for the non-perturbative parts, and they are expressed as matrix elements in non-relativistic QCD. At the leading order of v there are four matrix elements representing non-perturbative physics. The perturbative parts are calculated at the leading order of coupling constants. Some numerical results, especially numerical results for HERA, are given.