You can't.
If you know the energy as a function of bond length E(r), then you can determine the harmonic-oscillator eigenfrequencies from the second-order derivatives of E with respect to r, taken at the equilibrium point, where the first-order derivatives (forces) are zero. Getting anharmonic corrections requires even higher-order derivatives.
Or, to put it another way, the vibrational frequencies depend on (and are in fact a kind of description of) the potential-energy curve E(r) of the bond distance. You need to know the whole curve to know the frequencies. If you only know the energy and bond length, all you have is a single point on the curve. The simplest way to get the information about it is to determine the derivatives as above, in what's essentially a Taylor expansion around the equilibrium point.
You can't get E(r) with anything less than a full quantum-mechanical calculation.