Starting from an ab initio-derived two-site dimer Hubbard hamiltonian on a triangular lattice, we calculate the superconducting gap functions and critical temperatures for representative κ-(BEDT-TTF)_2X superconductors by solving the linearized Eliashberg equation using the Two-Particle Self-Consistent approach (TPSC) extended to multi-site problems. Such an extension allows for the inclusion of molecule degrees of freedom in the description of these systems. We present both, benchmarking

more »

... s for the half-filled dimer model as well as detailed investigations for the 3/4-filled molecule model. Remarkably, we find in the latter model that the phase boundary between the two most competing gap symmetries discussed in the context of these materials -d_xy and the recently proposed eight-node gap s+d_x^2-y^2 symmetry- is located within the regime of realistic model parameters and is especially sensitive to the degree of in-plane anisotropy in the materials as well as to the value of the on-site Hubbard repulsion. We show that these results provide a more complete and accurate description of the superconducting properties of κ-(BEDT-TTF)_2X than previous Random Phase Approximation (RPA) calculations and, in particular, we discuss predicted critical temperatures in comparison to experiments. Finally, our findings suggest that it may be even easier to experimentally switch between the two pairing symmetries as previously anticipated by invoking pressure, chemical doping or disorder effects.