A new statistical mechanics approach to calculating probabilities of chemical bonds being open is developed. The method uses effective phonon theory and becomes a mean field theory which shows the observed second order phase transition on helix melting. The temperature and width of the transition is accurately calculated without fitting parameters to melting. The method was applied to the melting of a DNA with a hydration spine. The effect on the DNA and the melting of the spine are explored.