Structures of aqueous [Be(H2O)(4)](2+), its outer-sphere and inner-sphere complexes with F-, Cl-, and SO42-, and dinuclear complexes with a [Be-2(kappa-OH)(kappa-SO4)](+ )core have been studied through Car-Parrinello molecular dynamics (CPMD) simulations with the BLYP functional. According to constrained CPMD/BLYP simulations and pointwise thermodynamic integration, the free energy of deprotonation of [Be(H2O)(4)](2+) and its binding free energy with F- are 9.6 and -6.2 kcal/mol, respectively, in good accord with available experimental data. The computed activation barriers for replacing a water ligand in [Be(H2O)(4)](2+) with F- and SO42-, 10.9 and 13.6 kcal/mol, respectively, are also in good qualitative agreement with available experimental data. These ligand-substitution reactions are indicated to follow associative interchange mechanisms with backside (S(N)2-like) attack of the anion relative to the aquo ligand it is displacing. Outperforming static density functional theory computations of the salient kinetic and thermodynamic quantities involving simple polarizable continuum solvent models, CPMD simulations are validated as a promising tool for studying the structures and speciation of beryllium complexes in aqueous solution.