Room-temperature ionic liquids (ILs) are salts, that are characteristically liquid at room temperatures and own remarkable properties, for example negligible vapour pressure, non-flammability, good thermal stability, and wide electrochemical window. Due to the large number and versatility of cation-anion pair combinations, the potential applicability of ILs is limitless. To improve the industrial efficiency of ILs, where it is used as an electrolyte, catalytic environment or as a thin layer of lubricant, ILs are often combined with various solid supports. Such nanoconfined behaviour of the ionic liquids was extensively studied through numerous experiments; however, a lack of a theoretical, fundamental picture is still noticeable. The reason for that can be originated in the complexity of the investigated model systems, with all conformational fluctuations and spatial reorientation. An ideal and consistent method that gives a correct interpretation of the experimental data and accurately predicts the atomic-level behaviour of nanoconfined ILs is molecular dynamic (MD) simulations. Never the less, employing of MD techniques on nanoconfined systems can be challenging. The first step, during the work on PhD, was to develop a complete and transferable methodology for the simulations of the nanoconfined ILs systems. The geometry of the simulated model system was designed to characterize a prototypical setup employed in the novel, experimentally designed Supported Ionic Liquid Phase catalysis (SILP). While there are currently numerous attempts to characterize thick films of molecular liquids on solid support, using molecular modelling techniques, this study stands out by its technical rigor – the approach used in this thesis is unique in that it is systematically tested against available experiments on liquids in the bulk, solid-liquid, and liquid-vacuum interfaces. In that sense, I believe this project has significantly raised the standards of molecular modelling of ILs and established a prototypical protocol for the determination of the correct parametrization of ILs. Using this model, I was able to relate the conformational molecular features with the structural properties of the liquid throughout the film from the solid to the vacuum interface. Most significantly, these features are further associated with the non-homogenous and anisotropic dynamic properties of the film, which is the first analysis of this type in the literature of ILs. The somewhat surprising emergent result is the identification of multiple time scales governing molecular transport. Consequently, the optimal parametrization was then used to predict novel physical properties on the family of thick imidazolium based IL films ([CnMim][NTf2] (n=2,4,6,8,10)) including the spatial organization of the ions at the solid-liquid and/or liquid-gas interfaces as well as correct mobility of ions between and on those interfaces. The densities together with mobilities in-plane and normal to the interfaces were calculated based on statistical physics techniques to account for the loss of symmetry in the film. Interestingly, strong correlations between the IL’s stratification and diffusive transport were found several nanometres deep into IL films. Also, compared to changes in distribution of anions, which are small, the associated remodelling of density profiles becomes more and more pronounced with the size of cations. Moreover, similar interfacial restructuring is also seen with increasing temperature, as evidenced by the associated drop in surface tension. Furthermore, a line of study focused on thin layers has been comprehensively investigated. What is of importance to interfacial science both from a fundamental as well as an application point of view is to understand the structural ordering of mono- and bi-layers of the IL depending on how they are prepared from a drop in wetting simulations. To understand the relation between the emergent film structure, adsorption and establishment of the film, I performed extensive simulations in which a cylinder of [C2Mim][NTf2] is allowed to spread on a hydroxylated alumina surface, as a function of temperature and cylinder size. In all cases a monolayer is obtained on the surface, however, its structure is found to be very much sensitive to temperature. Interestingly, a transition from an ordered to a disordered structure is found between 300K and 350K, as evidenced by the change in the structure factor of both, cations and anions, which form, in principle intercalated lattices. The stability of the pattern is not only sensitive to temperature, but also to the surface density. Namely, disorder is introduced in the monolayer prior to the nucleation of a bilayer, the latter first appearing as islands on the monolayer at the critical density. This behaviour of IL also points to potentially interesting phase behaviour of fluids, which has been predicted theoretically to change upon confinement. Besides changes in thermodynamic conditions, considerable reorganization can be investigated by creating mixtures with three different combinations of altered cations (PFBMIm+) and anions (PF6-), which provided an important tool to regulate SILP catalysis. Namely, strong anti-correlations between the cation stratification and positioning of a particular catalyst within the film were found in the simulations. Specifically, the calculation of the potential of mean force (PMF) between the catalyst and the film shows that the minima in the potential coincide with the minima in cation density. Hence, the evaluation of PMFs, as well as stratification of density and transport coefficients is crucial for future upscaling.