Highly concentrated aqueous lithium chloride solutions were investigated by classical molecular dynamics (MD) and reverse Monte Carlo (RMC) simulations. At first MD calculations were carried out applying twenty-nine combinations of ion-water interaction models at four salt concentrations. The structural predictions of the different models were compared, the contributions of different structural motifs to the partial pair correlation functions (PPCF) were determined. Particle configurations obtained from MD simulations were further refined using the RMC method to get better agreement with experimental X-ray and neutron diffraction data. The PPCFs calculated from MD simulations were fitted together with the experimental structure factors to construct structural models that are as consistent as possible with both the experimental results and the results of the MD simulations. The MD models were validated according to the quality of the fits. Although none of the tested MD models can describe the structure perfectly at the highest investigated concentration, their comparison made it possible to determine the main structural properties of that solution as well. It was found that four nearest neighbors (oxygen atoms and chloride ions together) are around a lithium ion at each concentration, while in the surroundings of the chloride ion hydrogen atom pairs are replaced by one lithium ion as the concentration increases. While in pure liquid water four water molecules can be found around a central water molecule, near the solubility limit nearly all water molecules are connected to two chloride ions (via their hydrogen atoms) and one lithium ion (by their oxygen atoms).Ion-water potential model; Structure nowadays. Unfortunately, for highly concentrated LiCl solutions they are almost absent, and the comparison with experimental data sets (e.g. XRD or ND data) is completely missing.Several interaction models (force fields, FF) have been developed and used for MD simulations in the context of electrolyte solutions. However, for concentrated aqueous LiCl solutions their applicability is very limited, as was shown recently [65]. In that work 29 different models were investigated; it turns out that their predictions about density, dielectric constant and self-diffusion constants scatter widely. Neutron and X-ray weighted structure factors were also calculated and compared with experimental data: it was found that the agreement between the experimental and simulation curves was poor, especially for the highest investigated concentration (19.55 mol/kg, near the solubility limit).The reverse Monte Carlo method [66] is a useful tool to find three dimensional configurations consistent with experimental data sets. However the number of experimental data sets (mostly one or two structure factors from neutron and/or X-ray diffraction measurements) is low compared to the number of the independent partial pair correlation functions (10). This lack of information results in a highly underdetermined system (see the 'pure' RMC results in R...