Quinoline derivatives are common among many natural compounds; moreover, quinoline is a structural subunit which plays an important role in the design and synthesis of pharmaceutically active substances. This makes investigation of its derivatives and their properties particularly interesting and is the motivation for the presented X-ray diffraction and computational study of crystal structure, electron density (ED) and intermolecular interactions of four quinoline derivatives. The derivatives differ in the substituent at position 8 which is either a halogen atom (Cl, Br, I) or -S-Ph moiety. Highresolution X-ray measurements were performed in each case, resulting in experimental data of varying quality, possibly characterized by the presence of anharmonic motion, which makes them suitable for verifying capabilities of two different charge density models compared in this work.The refinement methods applied during X-ray data processing were: the multipole refinement (MR) in Hansen and Coppens formalism [1] and the X-ray wavefunction refinement (XWR) [2] implementing Hirshfeld partition of ED. The latter technique includes two steps: (a) Hirshfeld atom refinement (HAR) being only structural refinement and (b) X-ray constrained wavefunction fitting introducing experimental contribution to ED. According to the performed estimate, data resolution was not always sufficient to refine the 3rd and 4th order of anharmonicity of Cl, Br, I and S atoms in the investigated structures. Nevertheless, an attempt of such refinement was made in all the cases to check if the tested refinement techniques tended to detect the non-existing effect when additional parameters were introduced into the model. The analysis of refinement of various orders of anharmonicity and the dependence of the results on data quality and the model used will be presented. Proper establishing of hydrogen positions and treatment of its thermal motion is necessary for correct description of ED and related properties. XWR allowed refinement of hydrogen positions and ADPs in the considered compounds, however precision and accuracy was strongly dependent on data quality. With MR it was possible to freely refine hydrogen ADPs only for good-quality data sets, while refinement of hydrogen positions was not successful and restraining them to mean neutron distances [3] was necessary. Moreover, the experimental ED was subject to analysis in order to compare the ability of MR and XWR to model density features and evaluate performance of the methods in the case of data with evident presence of experimental errors. Periodic DFT calculations were performed to serve as a benchmark for the experimental results. [2] Grabowsky, S.; Luger, P.; Buschmann, J.; Schneider, T.; Schirmeister, T.; Sobolev, A. N.