Department of Physics2024-11-0920222469-992610.1103/PhysRevA.106.0433152-s2.0-85140247273http://dx.doi.org/10.1103/PhysRevA.106.043315https://hdl.handle.net/20.500.14288/7309We study the bound states of N identical bosons that are described by a multiband Bose-Hubbard model with generic hoppings and an attractive on-site interaction. Using a variational approach, we first derive exact integral equations for the dimers, trimers, tetramers, and other multimers, and then apply them to a one-dimensional sawtooth model that features two bands. In particular we reveal the presence of not only the off-site dimer states which consist of two monÖmers on different sites even in the strong-coupling limit but also the off-site trimer states which consist of either a dimer on one site and a monÖmer on another site or three monÖmers on three different sites. Our variational calculations for the ground states of on-site dimers, on-site trimers and off-site trimers benchmark perfectly well with the density matrix renormalization group (DMRG) simulations. We also present DMRG results for the ground states of on-site tetramers, off-site tetramers, on-site pentamers, off-site pentamers, and other multimers.OpticsPhysicsJournal Article2469-9934876436700006Q23798