The chiral phase transition and color superconductivity in an extended NJL model with eight-quark interactions are studied. The scalar-type nonlinear term hastens the chiral phase transition, the scalar-vector mixing term suppresses effects of the vector-type linear term and the scalar-diquark mixing term makes the coexisting phase wider.Quantum Chromodynamics (QCD) has non-perturbative properties. First-principle lattice QCD simulations are useful to study thermal systems at zero or small density [1,2]. At high density, however, lattice QCD is still not feasible due to the sign problem. Therefore, effective models are used in finite density region. One of the models is the Nambu-Jona-This model has the mechanism of spontaneous chiral symmetry breaking, but it has not the confinement mechanism. However, this model has been widely used [4,5] with the mean field approximation (MFA), for example, for analyses of the critical endpoint of chiral phase transition [6,7,8,9,10,11].As for the NJL model, only a few studies were done so far on roles of higher-order multiquark interactions [12,13], except for the case of the six-quark interaction coming from the