The six- and four-quark systems are studied in emphasized that the color confinement used in multiquark the framework of constituent quark models. It is system should be different from the one used in two- or three-quark system. For six-quark system, we look for △△ and △△ dibaryon resonances by calculating NN scattering phase shifts with explicit coupling to these dibaryon channels in quark delocalization and color screening model. The model gives a good description of low-energy NN properties and predicts IJP = 03＋ and 01＋ △△resonances, which can be promising candidates for the isoscalar ABC structure reported by the CELSIUS-WASA Collaboration. For tetraquark system, a flux-tube quark model with multi-body confinement interaction is employed to study Y（2175） as a tetraquark state. The Y（2175） with diquark-antidiquark structure has energy 2174 MeV which is very consistent with experimental data. The calculation shows that multi-body confinement potential may play a vital role in the multiquark system.
The NΔ and ΔΔ dibaryon resonances are studied by calculating the NN scattering phase shifts with explicitly coupling these dibaryon channels in a multi-channel coupling calculation with two quark models. These quark models,the chiral quark model and quark delocalization color screening model,describe the NN S-,D-wave phase shifts below the π production threshold quantitatively well. Both quark models predict the 1D2 resonance discovered in NN partial wave phase shift analysis and the J=1 or 3 isoscalar resonance recently reported by CELSIUS-WASA Collaboration are NΔ5S2 and ΔΔ7S3 resonance,respectively. Further measurements at COSY,LEPS and Lanzhou Cooling Storage Ring(CSR) to check the ΔΔ resonance are discussed.
