To understand the operation of complex electrolytes we must determine which species are present and active at the interface with the electrode. We combine molecular dynamics sampling with continuum modeling to determine the speciation of Ca(BH4)2 in THF next to a model electrode (graphite). Mulitvalent ions in organic solvents with poor dielectric screening may exist in multiple coordinated states. We explore the free energy landscape associated with the interconversion of species in the bulk electrolyte using metadynamics simulations. We find that the chemical equilibrium between neutral and charged species is dominated by disproportionation following dimer formation - consistent with experimental findings. We also make use of clustering analysis to reveal distinct molecular conformations within constrained coordination states. To assess how the chemical equilibrium may be modified by the potential difference at the electrode, we develop a continuum model based on the modified Poisson-Boltzmann equation, with volume exclusion. This model also incorporates molecular-scale details through sampled free energy profiles of each relevant species within the solvent as a function of distance from the electrode. We find that a narrow (approximately 1 nm) region near the electrode, with highly variable speciation driven by potential difference, dominates the electrochemical activity of this electrolyte.