A molecular based understanding of beryllium chemistry in the context of biomolecules is necessary for gaining progress in prevention and treatment of chronic beryllium disease. One aspect that has hindered the theoretical progress has been the lack of a simple classical two-body potential for the aqueous beryllium ion (Be2+) to be used with biomolecular simulations. We provide new parameters for Be2+ that capture the structural and reactive properties of this small dication. Using classical molecular dynamics simulations, we show that these parameters reproduce the correct radial distribution function and coordination numbers for this cation in explicit aqueous solution when compared to published diffraction and NMR measurements. The geometrical parameters obtained using classical simulations are also in agreement with ab initio calculations. We successfully predict the vibrational modes of the tetra aqua Be2+ dication from ab initio calculations on solvated structures obtained from the simulations. The calculated vibrational modes show better agreement with experiments compared to any published work. This new potential also produces a well-established hydrogen bonding between the first and second solvation shells. More importantly, when the molecular dynamics (MD) and ab initio results are interpreted in concert, the dynamics and nature of interactions between the first and second shells capture the pivotal role they play on the reactivity of aqua-Be complexes.