We investigate a hybrid system composed of ultracold Rydberg atoms immersed in an atomic Bose-Einstein condensate (BEC). The coupling between the Rydberg electrons and BEC atoms leads to the excitation of phonons, the exchange of which induces Yukawa interaction between Rydberg atoms. Due to the small electron mass, the effective charge associated with this quasiparticle-mediated interaction can be large, while its range is equal to the healing length of the BEC, which can be tuned by adjusting the scattering length of the BEC atoms. We find that for small healing lengths, the distortion of the BEC can "image" the wave function density of the Rydberg electron, while large healing lengths induce an attractive Yukawa potential between the two Rydberg atoms that can form a new type of ultra-long-range molecule. We discuss both cases for a realistic system.Impurities in a Bose-Einstein condensate (BEC) have attracted much attention and motivated the investigation of a wide range of phenomena. For example, the motion of a single impurity in a BEC can probe the superfluid dynamics [1][2][3], while an ionic impurity in a BEC can form a mesoscopic molecular ion [4]. Due to the selfenergy induced by phonons (excitations of the BEC), a neutral impurity can self-localize in both a homogeneous and a harmonically trapped BEC [5][6][7], which sheds light on polaron physics [8,9]. Exchanging phonons between multiple impurities induces an attractive Yukawa potential between each pair of impurities [10,11], which leads to the so called "co-self-localization" [12] and is related to forming bipolarons and multipolarons [13]. Recent experiments, where an atom of a BEC is excited into a Rydberg state [14] to study phonon excitations and collective oscillations, open the door to exploration of the electron-phonon coupling in ultracold degenerate gases, a phenomenon responsible for the formation of Cooper pairs of two repelling electrons in BCS superconductivity [15].In this Letter, we study Rydberg atoms immersed in a homogeneous BEC, as sketched in Fig. 1(a). Rydberg atoms consist of an ion core and a highly excited electron with its oscillatory wave function Ψ e extending to large distances of the order of ∼n 2 a 0 (n: principle quantum number, a 0 : Bohr radius). As pointed out by Fermi [16], the interaction between the quasi-free electron at x and a ground state atom at r can be approximated at low scattering energies by a contact interaction parametrized by an energy-dependent s-wave scattering lengthWhile the s-wave approximation is valid for qualitative analysis, we include higher-partial wave contributions for quantitative results [17]. A s (k) depends on the scattering energy via the local wave number k(r) given by where R y is the Rydberg constant, 0 the vacuum permittivity, e and m e the charge and mass, respectively, of the electron with angular momentum e and quantum defect δ e . For low-e state, Eq. (1) gives an effective interaction between Rydberg and ground state atoms aswhich leads to an attraction and formation of ...