On the basis of linear programming, new sets of entanglement witnesses (EWs) for 3 ⊗ 3 and 4 ⊗ 4 systems are constructed. In both cases, the constructed EWs correspond to the hyperplanes contacting, without intersecting, the related feasible regions at line segments and restricted planes respectively. Due to the special property of the contacting area between the hyper-planes and the feasible regions, the corresponding hyper-planes can be turned around the contacting area throughout a bounded interval and hence create an infinite number of EWs. As these EWs are able to detect entanglement of some PPT states, they are non-decomposable (nd-EWs).

One of the problems concerning entanglement witnesses (EWs) is the construction of them by a given set of operators. Here several multi-qubit EWs called stabilizer EWs are constructed by using the stabilizer operators of some given multi-qubit states such as GHZ, cluster and exceptional states. The general approach to manipulate the multi-qubit stabilizer EWs by exact(approximate) linear programming (LP) method is described and it is shown that the Clifford group play a crucial role in finding the hyper-planes encircling the feasible region. The optimality, decomposability and non-decomposability of constructed stabilizer EWs are discussed.

The Gell-Mann λ matrices for Lie algebra su(3) are the natural basis for the Hilbert space of Hermitian operators acting on the states of a three-level system(qutrit). So the construction of EWs for two-qutrit states by using these matrices may be an interesting problem. In this paper, several two-qutrit EWs are constructed based on the Gell-Mann matrices by using the linear programming (LP) method exactly or approximately. The decomposability and non-decomposability of constructed EWs are also discussed and it is shown that the λ-diagonal EWs presented in this paper are all decomposable but by adding λ-non-diagonal terms, one can obtain various non-decomposable EWs.