We perform a precision computation of hybrid static potentials with quantum numbers ΛThe resulting potentials are used to estimate masses of heavycc andbb hybrid mesons in the Born-Oppenheimer approximation. Part of the lattice gauge theory computation, which we discuss in detail, is an extensive optimization of hybrid static potential creation operators. The resulting optimized operators are expected to be essential for future projects concerning the computation of 3-point functions as e.g. needed to study spin corrections, decays or the gluon distribution of heavy hybrid mesons.