Since the geometry optimization may converge to various local energy minima based on the initial ion orientation, we repeat the optimization process for each system while varying the initial configuration of the OTF -anion(s) relative to both the sodium cation(s) and the electrode surface. Subsequently, we identify the geometry with the lowest overall energy as the optimal configuration and calculate the adsorption energy based on this configuration as