Absorption and emission lines of the iron-peak species Fe II are prominent in the infrared, optical and ultraviolet spectra of a myriad of astrophysical sources, requiring extensive and highly reliable sets of atomic structure and collisional data for an accurate quantitative analysis. However, comparisons among existing calculations reveal large discrepancies in the effective collision strengths, often up to factors of three, highlighting the need for further steps towards new converged calculations. Here we report a new 20 configuration, 6069 level atomic structure model, calculated using the multi-configurational Dirac-Fock method. Collision strengths and effective collision strengths are presented, for a wide range of temperatures of astrophysical relevance, from substantial 262 level and 716 level Dirac R-matrix calculations, plus a 716 level Breit-Pauli R-matrix calculation. Convergence of the scattering calculations is discussed, and results are critically compared with existing data in the literature, providing us with error estimates for our data. As a consequence, we assign an uncertainty of ±15% to relevant forbidden and allowed transitions encompassed within a 50 level subset of the 716 level Dirac R-matrix dataset. To illustrate the implications of our new datasets for the analysis of astronomical observations of Fe II, they are incorporated into the Cloudy modelling code, sample Fe II spectra are generated and compared.