Abstract
We set up and parametrize a Hubbard model for interacting quantum dots in bilayer graphene and study double dots as the smallest multi-dot system. We demonstrate the tunability of the spin and valley multiplets, Hubbard parameters, and effective exchange interaction by electrostatic gate potentials and the magnetic field. Considering both the long- and short-range Coulomb interaction, we map out the various spin and valley multiplets and calculate their energy gaps for different dot sizes and inter-dot separations. For half-filling and large valley splittings, we derive and parametrize an effective Heisenberg model for the quantum dot spins.