We consider the quantum entanglement of the electronic and vibrational degrees of freedom in molecules with tendencies towards double welled potentials. In these bipartite systems, the von Neumann entropy of the reduced density matrix is used to quantify the electron-vibration entanglement for the lowest two vibronic wavefunctions obtained from a model Hamiltonian based on coupled harmonic diabatic potential-energy surfaces. Significant entanglement is found only in the region in which the ground vibronic state contains a density profile that is bimodal (i.e., contains two separate local maxima). However, in this region two distinct types of density and entanglement profiles are found: one type arises purely from the degeneracy of energy levels in the two potential wells and is destroyed by slight asymmetry, while the other arises through strong interactions between the diabatic levels of each well and is relatively insensitive to asymmetry. These two distinct types are termed fragile degeneracy-induced entanglement and persistent entanglement, respectively. Six classic molecular systems describable by two diabatic states are considered: ammonia, benzene, BNB, pyridine excited triplet states, the Creutz-Taube ion, and the radical cation of the "special pair" of chlorophylls involved in photosynthesis. These chemically diverse systems are all treated using the same general formalism and the nature of the entanglement that they embody is elucidated.