Résumé : In multiconfiguration Dirac–Hartree–Fock (MCDHF) calculations, there is a strong coupling between the localization of the orbital set and the configuration state function (CSF) expansion used to determine it. Furthermore, it is well known that an orbital set resulting from calculations, including CSFs describing core–core correlation and other effects, which aims to lower the weighted energies of a number of targeted states as much as possible, may be inadequate for building CSFs that account for correlation effects that are energetically unimportant but decisive for computed properties, e.g., hyperfine structures or transition rates. This inadequacy can be traced in irregular or oscillating convergence patterns of the computed properties as functions of the increasing orbital set. In order to alleviate the above problems, we propose a procedure in which the orbital set is obtained by merging several separately optimized, and mutually non-orthogonal, orbital sets. This computational strategy preserves the advantages of capturing electron correlation on the total energy through the variational MCDHF method and allows to target efficiently the correlation effects on the considered property. The orbital sets that are merged are successively orthogonalized against each other to retain orthonormality. The merged orbital set is used to build CSFs that efficiently lower the energy and also adequately account for the correlation effects that are important for the property. We apply the procedure to compute the hyperfine structure constants for the (Formula presented.) and (Formula presented.) states in (Formula presented.) Li and show that it leads to considerably improved convergence patterns with respect to the increasing orbital set compared to standard calculations based on a single orbital set, energy-optimized in the variational procedure. The perspectives of the new procedure are discussed in a broader context in the summary.