Abstract
A system consisting of two protons in two double minimum potentials coupled by a bilinear interaction was studied. Explicit calculations were performed for local potentials described by a combination of harmonic oscillator and Gaussian terms. To determine eigenfunctions and -values variational procedures were employed which made use of numerically determined 1-particle solutions. Eigenvalues and (graphical representations of) 2d eigenfunctions are presented, classified according to their symmetry, and used to determine the time evolution of the system for given initial conditions.