The agreement in workfunction predictions between RISM and AIMD simulations, despite the failure to capture the characteristic double peak is intriguing. The fine features of interfacial charge polarization, heavily influenced by water structure, are generally considered crucial for accurate workfunction predictions. The observed agreement might suggest that workfunctions are more sensitive to the average density distribution, which RISM adequately captures, rather than the fine details of water structure. Alternatively, this agreement could arise from a fortuitous cancellation of errors within the RISM framework.