The aggregation in a solution of charged dyes such as Rhodamine B (RB) is significantly affected by the type of counterion, which can determine the self-assembled structure that in turn modulates the optical properties. RB aggregation can be boosted by hydrophobic and bulky fluorinated tetraphenylborate counterions, such as F5TPB, with the formation of nanoparticles whose fluorescence quantum yield (FQY) is affected by the degree of fluorination. Here, we developed a classical force field (FF) based on the standard generalized Amber parameters that allows modeling the self-assembling process of RB/F5TPB systems in water, consistent with experimental evidence. Namely, the classical MD simulations employing the re-parametrized FF reproduce the formation of nanoparticles in the RB/F5TPB system, while in the presence of iodide counterions, only RB dimeric species can be formed. Within the large, self-assembled RB/F5TPB aggregates, the occurrence of an H-type RB-RB dimer can be observed, a species that is expected to quench RB fluorescence, in agreement with the experimental data of FQY. The outcome provides atomistic details on the role of the bulky F5TPB counterion as a spacer, with the developed classical FF representing a step towards reliable modeling of dye aggregation in RB-based materials.