In a first step, the standard non-perturbative numerical renormalization group method is employed to demonstrate the hierarchical structure of correlations unveiled by natural orbitals. This simplification brought new insights for simulating quantum impurity problems, and a new algorithm is developed to generate an optimized subset of natural orbitals independently of existing methods, going beyond their usual limitations. This algorithm is presented in detail in the book, and a careful benchmark on known results is carried out to guarantee the validity of the method. It is then used to study spatial entanglement structures under various conditions that were not accessible with previous methods, such as representing the electron bath by a realistic 2D square lattice or taking account of static disorder in the metallic host. In the last chapter, the non-interacting problem in the presence of disorder is studied through random matrix theory, reproducing some of the results presented in the previous chapters. The main original result of this chapter lies in the analytical calculation of the joint distribution of one-particle orbitals energies and amplitudes of the impurity, which makes it possible to calculate any disordered averaged local correlation functions. Starting from this result, calculations in the large-N limit are compared with numerical simulations.