Molecular dynamics simulations have been carried out for aqueous electrolyte solutions between model electrode surfaces. The effect of solvent model flexibility on bulk and double layer properties was observed for electrode surface charge densities of 0, +/-0.1, and +/-0.2 C/m(2) and ion concentrations of 0, 0.5, and 1 M. Two flexible models were used to isolate the effects of flexibility from the effects of a change in the condensed-phase dipole moment. Model flexibility increases the pure water self-diffusion coefficient while a larger liquid dipole moment substantially decreases it. There is an increase in ion contact adsorption and counter ion affinity with the flexible models, suggesting that the ions are less tightly solvated. This conclusion is consistent with observed enhancements of solvated ion densities near uncharged electrodes for the flexible water case. Mobile ions in high concentration quickly damp out the electric field even at high electrode charge densities, but for dilute ion concentrations the field may extend to the center of the cell or beyond. In these cases it is more appropriate to integrate Poisson's equation from the electrode surface outward instead of the common method of assuming zero field at the center of the simulation cell. Using this methodology, we determine the voltage drop across the half-cell for both the rigid and flexible models. The half-cell voltage drop shows some dependence on ion concentration, but solvent flexibility has little effect on that behavior. (C) 2003 American Institute of Physics.