The one-dimensional dilute Kondo lattice model is investigated by means of bosonization for different dilution patterns of the array of impurity spins. The physical picture is very different if a commensurate or incommensurate doping of the impurity spins is considered. For the commensurate case, the obtained phase diagram is vertified using a non-Abelian density-matrix renormalization-group algorithm. The paramagnetic phase widens at the expense of the ferromagnetic phase as the f spins are diluted. For the incommensurate case, short-range antiferromagnetic correlations are found to dominate at low doping, which distinguishes the dilute Kondo lattice model from the standard Kondo lattice model. DOI: 10.1103/PhysRevB.69.174425 PACS number͑s͒: 75.10.Jm, 75.40.Mg, 05.50.ϩq Heavy fermion systems have been of great theoretical interest since their discovery some 20 years ago. 1 The central problem posed by heavy fermion materials is to understand the interaction between an array of localized moments ͑gen-erally f electrons in lanthanide or actinide ions͒ and conduction electrons ͑generally s or d band͒. This situation is well described by an antiferromagnetically coupled Kondo-type model.The solution of Kondo-type models is well understood in two limiting cases; the single-impurity limit 2 which can be reduced to a one-dimensional problem and solved via Bethe ansatz, and second the Kondo lattice model ͑KLM͒, which was solved via bosonization 3 and numerous numerical approaches 4,5 in one dimension for half filling and partial conduction-band filling. For half filling the results indicate the existence of a finite spin and charge gap. Accordingly in this case the Kondo lattice model is an insulator with welldefined massive solitonic excitations of the spin sector.For partial conduction-band filling, the conduction electrons form a Luttinger liquid, with spin and charge separation. 5 The localized spins, however, exhibit ferromagnetism, due to an effective double-exchange coupling. 3,4 The double exchange is driving the system toward ferromagnetism, while the fluctuations generated by Kondo singlets compete against this tendency. As a consequence, the paramagnetic to ferromagnetic phase transition is of the quantum order-disorder type, typical to models with an effective random field. 3 However, for small Kondo coupling and close to half filling a Ruderman-Kittel-Kasuya-Yosida liquid state and polaronic regime are always present. 4 For additional properties, see earlier reviews of Ref. 6. Beyond these two solvable limits, no rigorous results exist for the intermediate cases, where the number of impurities are neither one nor equal to the number of sites. This is the focus of our study. We concentrate on the one-dimensional case, and start from the Kondo lattice limit introducing impurity spin holes, that is, we will be dealing with a dilute Kondo lattice model ͑DKLM͒:where L is the number of sites and tϾ0 is the conduction electron hopping. We measure the Kondo coupling J in units of the hopping t. We denote by ...