A molecular model for an electron transfer reaction at a solution-electrode interface is developed. The solvent diabatic free energy curves for the reaction Fe+3+e(-)-->Fe+2 are calculated using an umbrella sampling procedure. These are used to calculate the rate of electron transfer as a function of the electrode-solution potential difference in the electronically nonadiabatic region. A model Hamiltonian for the adiabatic case is also developed and used to calculate the adiabatic free energy surface. Reactive flux correlation functions are used to determine the solvent dynamical corrections to the rate. A comparison between the molecular dynamics and the Kramers and Grote-Hynes theories is made.