Radicals in aqueous solution: assessment of density-corrected SCAN functional.

We study self-interaction effects in solvated and strongly-correlated cationic molecular clusters, with a focus on the solvated hydroxyl radical. To address the self-interaction issue, we apply the DC-rSCAN method, with the auxiliary density matrix approach. Validating our method through simulations of bulk liquid water, we demonstrate that DC-rSCAN maintains the structural accuracy of rSCAN while effectively addressing spin density localization issues. Extending our analysis to solvated cationic molecular clusters, we find that the hemibonded motif in the [CHS∴CHSH] cluster is disrupted in the DC-rSCAN simulation, in contrast to rSCAN that preserves the (three-electron-two-center)-bonded motif. Similarly, for the [SH∴SH] cluster, rSCAN restores the hemibonded motif through spin leakage, while DC-rSCAN predicts a weaker hemibond formation influenced by solvent-solute interactions. Our findings demonstrate the potential of DC-rSCAN combined with the auxiliary density matrix method to improve electronic structure calculations, providing insights into the properties of solvated cationic molecular clusters. This work contributes to the advancement of self-interaction corrected electronic structure theory and offers a computational framework for modeling condensed phase systems with intricate correlation effects.