  Ultracold polar molecules in multilayered systems have been experimentally
realized very recently. While experiments study these systems almost
exclusively through their chemical reactivity, the outlook for creating and
manipulating exotic few- and many-body physics in dipolar systems is
fascinating. Here we concentrate on few-body states in a multilayered setup. We
exploit the geometry of the interlayer potential to calculate the two- and
three-body chains with one molecule in each layer. The focus is on dipoles that
are aligned at some angle with respect to the layer planes by means of an
external eletric field. The binding energy and the spatial structure of the
bound states are studied in several different ways using analytical approaches.
The results are compared to stochastic variational calculations and very good
agreement is found. We conclude that approximations based on harmonic
oscillator potentials are accurate even for tilted dipoles when the geometry of
the potential landscape is taken into account.
