Numerical studies together with asymptotic and spectral analysis establish regimes where soliton pairs in degenerate optical parametric oscillators fuse, repel, or form bound states. A novel bound state stabilized by coupled internal oscillations is predicted.