A QCD multiquark cluster system is studied in the relativistic harmonic oscillator potential model （RHOPM）, and the electromagnetic form factors of the pion, proton and deuteron in the RHOPM are predicted. The calculated theoretical results are then compared with existing experimental data, finding very good agreement between the theoretical predictions and experimental data for these three target particles. We claim that this model can be applied to study QCD hadronic properties, particularly neutron properties, and to find six-quark cluster and/or nine-quark cluster probabilities in light nuclei such as helium 3He and tritium 3H. This is a problem of particular importance and interest in quark nuclear physics.
Based on the Dyson-Schwinger Equations （DSEs）, the two-quark vacuum condensate, the four-quark vacuum condensate, and the quark gluon mixed vacuum condensate in the non-perturbative QCD vacuum state are investigated by solving the DSEs with rainbow truncation at zero- and finite- temperature, respectively. These condensates are important input parameters in QCD sum rule with zero and finite temperature, and in studying hadron physics, as well as predicting the quark mean squared momentum rn02- also called quark virtuality in the QCD vacuum state. The present calculated results show that these physical quantities are almost independent of the temperature below the critical point temperature Tc=131 MeV, and above Tc the chiral symmetry is restored. For comparison we calculate the temperature dependence of the ＂in-hadron condensate＂ for pion. At the same time, we also calculate the ratio of the quark gluon mixed vacuum condensate to the two-quark vacuum condensate by using these condensates, and the unknown quark mean squared momentum in the QCD vacuum state has been obtained. The results show that the ratio m2/0（T） is almost fiat in the temperature region from 0 to To, although there are drastic changes of the quark vacuum condensate and the quark gluon mixed vacuum condensate at the region. Our predicted ratio comes out to be m2/0（T）=2.41 GeV2 at the Chiral limit, which is consistent with other theory model predictions, and strongly indicates the significance that the quark gluon mixed vacuum condensate has played in the virtuality calculations.
