The standard QCD action is improved by the addition of irrelevant operators built with chiral composites. An effective Lagrangian is derived in terms of auxiliary fields, which has the form of the phenomenological chiral Lagrangians. Our improved QCD action appears promising for numerical simulations as the pion physics is explicitely accounted for by the auxiliary fields. QCD in the high energy limit is efficiently described by the perturbative expansion in the gauge coupling constant. For the low energy properties instead, we do not have an equally satisfactory theory.At the phenomenological level a description of the low energy physics related to the symmetry breaking of the chiral invariance is obtained by using the phenomenological chiral Lagrangians [1,2], but a direct derivation from the basic theory of quarks and gluons is still lacking.At a more fundamental level the lattice formulation has allowed us to investigate the quark confinement [3] and the spontaneous breaking of chiral invariance, but due to the complexities related to the definition of chiral fermions [3,4], in numerical simulations the chiral limit is achieved only through a fine-tuning procedure. Even though these difficulties can be overcome, we remain unable to unify the treatment of high energy and low energy properties.For these reasons we have developed an approach based on the use of quark composites as fundamental variables. The idea behind it is that a significant part of the binding of the hadrons can be accounted for in this way, so that the "residual interaction" is sufficiently weak for a perturbative treatment.The quark-composites approach is in principle fairly general, since it allows us to treat all the hadrons composite of quarks, but for technical reasons the composites with the quantum numbers of mesons and barions are treated in a different way. The nucleonic composites, for instance, naturally satisfy the Berezin integration rules and we derived the substitution rules which allow us to replace polynomials of the quark fields by appropriate polynomials of these composites [5] in the partition function. The mesonic composites instead, due to the complexity of the integral over even elements of a Grassmann algebra, are replaced by auxiliary fields [6] by means of the Stratonovich-Hubbard transformation [7]. Even though our effective action is not renormalizable by power counting, due to the renormalizability of QCD only a finite number of free parameters can be generated by the counterterms because of the BRS identities.We assume the modified partition functionwhere S Y M is the Yang-Mills action, S q is the action of the quark field and S C is a four fermions irrelevant operator which provides the kinetic terms for the quark composites with the quantum numbers of the chiral mesons. λ is the quark field while the gluon field is associated to the link variables V . Differentials in square brackets are understood to be the product of the differentials over the lattice sites and the internal indices. All the fields live in a...