Afterwards, either 500 iso-octane or 300 n-hexadecane molecules are randomly distributed by using the Amorphous Builder plugin in MAPS, which roughly corresponds to the expected experimental liquid density for the simulation cell size. A schematic of the simulation box is shown in Figure . The different elements of the system are modelled using the all-atom optimised potential for liquid simulations (OPLS-AA) force field. The long-chain parametrisation is needed to overcome the well-documented overestimate of the liquid-gel transition temperature that occurs for long-chain alkanes, e.g. n-hexadecane. Another shortcoming of additive force fields based on point charges, such as OPLS-AA, is that they underestimate the strength of important intermolecular interactions for the current study. Notably, these force fields neglect the formation of ion-induced dipoles, which underpin the cation-π interactions between the zwitterionic surfactant and the deposit. We overcome this limitation through the explicit consideration of π-electrons with virtual particles using parameters from the INTERFACE force field (IFF). Following the IFF approach, the multipoles along the cloud of conjugated π-electrons are more accurately described, allowing one to capture the physical behaviour surrounding aromatic groups. An additional advantage is that no further parametrisation is needed for its integration with the rest of the OPLS-AA parameters.