“…The MD simulations carried out in this work were performed with the GROMACS 2019 [39][40][41] code considering the OPLS-AA (optimized potentials for liquid simulations -all atoms) [22] force field for all compounds. This force field was selected due to its versatility, existence of a vast number of parameters in the literature, compatibility and easiness of implementation, and also because it was previously found to enable the calculation of diffusivities of other compounds, such as, D MD 12 of amoxicillin in water and ethanol [27], D MD 11 of different solvents like n-decane, n-hexadecane, noctacosane, ethanol, 1-decanol, cyclohexane, benzene, 2methylpentane, 2,2-dimethylbutane, ethyl butanoate, and nhexanoic acid [42], and D MD 12 of C3-C5 ketones in SC-CO 2 [19], which were in good agreement with available experimental results. Preliminary calculations for some selected systems also supported our choices as shown below.…”