A new method for constructing empirical valence bond potential energy surfaces for reactions is presented. Building on the generalized Gaussian approach of Chang-Miller, V12(2)(q) is represented by a Gaussian times a polynomial at the transition state and generalized to handle any number of data points on the potential energy surface. The method is applied to two model surfaces and the HCN isom...