Ernesto García-Alfonso

Francois Coppens

Manuel Barranco

Martí Pi

Frank Stienkemeier
email: stienkemeier@uni-freiburg.de

Nadine Halberstadt
email: nadine.halberstadt@irsamc.ups-tlse.fr

Alkali atoms attached to vortex-hosting Helium nanodroplets

Light absorption or accordingly fluorescence excitation spectroscopy of alkali atoms attached to 4 He droplets is investigated as a possible way of detecting the presence of vortices. To this end, we have calculated the equilibrium configuration and energetics of alkali atoms attached to a 4 He1000 droplet hosting a vortex line using 4 He density functional theory. We use them to study how the dipole absorption spectrum of the alkali atom is modified when the impurity is attached to a vortex line. Spectra are found to be shifted to the blue (higher frequencies) and broadened compared to vortex-free droplets, because the dimple in which the alkali atom sits at the intersection of the vortex line and the droplet surface is deeper. This effect is smaller for lighter alkali atoms, and all the more so when using a quantum description since in this case, they sit further away from the droplet surface on average due to their zero-point motion. Spectral modifications due to the presence of a vortex line are minor for np ← ns excitation, and therefore insufficient for vortex detection. In the case of higher n p ← ns or n s ← ns (n > n) excitations the shifts are larger, as the excited state orbital is more extended and therefore more sensitive to changes in the surrounding helium density.

I. INTRODUCTION

4 He droplets built in supersonic beam experiments can be characterized by a temperature of about 0.37 K. [1] At that temperature they are superfluid, which makes them ideal ultra-cold matrices for spectroscopy studies of captured dopants. They are also ideal laboratories for addressing superfluidity at the nanoscale. [2][3][4][5][START_REF] Callegari | Handbook of High Resolution Spectroscopy[END_REF] In particular, they may host quantized vortices, a dramatic fingerprint of the superfluid state. [START_REF] Donnelly | Quantized vortices in helium II[END_REF][8][9] Since their existence was proposed by Onsager and Feynman, the study and detection of quantized vortices in superfluid and superconducting systems has not ceased and attracts the interest of physicists and chemists working in low temperature processes, semiconductors, quantum optics, and atomic and molecular physics and chemistry. It is hard to find a more widespread subject of study.

Quantized vortices are ubiquitous in quantum systems under rotation, as e.g. liquid helium samples at temperatures below 2.17 K (the superfluid transition temperature), [START_REF] Donnelly | Quantized vortices in helium II[END_REF] cold-gas Bose-Einstein condensates, [10] and electron droplets in low-dimensional quantum dots in a magnetic field. [11] However, the only small selfbound system where vortices are likely to occur is a 4 He droplet. [12,13] So far, no clear experimental signature of vortex formation in small 4 He droplets has been found, although this has been the subject of several experimental and theoretical proposals and studies. Many of them involve light absorption or laser induced fluorescence excitation spectroscopy, which usually give similar spectra, the former being easier to simulate theoretically and the latter easier to measure experimentally: we will refer to them as light absorption spectroscopy, or simply absorption spectroscopy, throughout this article. Close et al. [14] suggested the light absorption spectrum of Ag atoms in the bulk of superfluid droplets as a possible candidate: if Ag atoms are attached to the vortex line, their environment changes and this could be reflected in the absorption spectrum. To our knowledge, this possibility has never been explored. Along the same line, a theoretical study has been carried out for electron bubbles attached to vortex lines in superfluid liquid helium as a function of pressure. It was shown that the electron absorption spectrum is shifted by the presence of a vortex line, but the shift is too small to be experimentally detectable. [15] Another theoretical proposal, [16] which has not been tested experimentally either, suggested that microwave absorption spectroscopy of a Ca dopant could show a difference when a vortex was present. More recently, absorption spectroscopy was also suggested [17] for the same dopant. The rationale is basically the same as for Ag atoms: the change induced in the dopant environment due to its binding to a vortex line would translate as a change in the absorption spectrum. Ca impurities are known to reside in a deep dimple at the surface of the helium droplet. [17] Their binding to a vortex line would result in their further sinking inside the droplet. As a result, their absorption spectrum would evolve from that of an impurity on the droplet surface to that of an im-This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0008923

purity in its bulk. Hence a satellite peak should appear, which would be blue-shifted from the main peak and with an intensity depending on the number of vortex-hosting droplets. Inspection of existing experimental results [18] showed no such signature.

It is worth stressing that these proposals and experiments involve droplets made of several thousand He atoms created by adiabatic expansion of a helium gas into vacuum. [1] All attempts to detect the presence of vortices in droplets of this size have failed so far. Very recently, large 4 He droplets made of 10 8 -10 11 atoms have been created using hydrodynamic instability of a cryogenic liquid helium jet passing through the nozzle of a molecular beam apparatus, as reviewed in Refs. 9 and 19. Helium drops, which are initially produced in the normal, non-superfluid phase, can acquire angular momentum during the passage of the fluid through the nozzle. Such droplets could be analyzed as single objects by x-ray and XUV light using intense high harmonics sources. [START_REF] Gomez | [END_REF]21] This has made it possible to determine their shapes and, doping them with Xe atoms, [START_REF] Gomez | [END_REF]22] the presence of vortices, confirming their superfluid nature.

The following question therefore arises: why have vortices not been detected in smaller droplets (10 3 -10 4 atoms)? Is it because no clear signature of their presence has been found yet? Or is it because vortices are too scarce due to the droplet formation and doping mechanisms? Let us recall that droplets containing a few thousand He atoms are produced in subcritical conditions by atom condensation and evaporative cooling, [1] which does not favor deposition of large amounts of angular momentum in the droplet. Some angular momentum is transferred to the droplet during the impurity pickup, which does nucleate vortex rings and loops. [23][24][25][26] However, they are dynamically unstable. At variance, large droplets such as those detected as single objects using x-ray or XUV light are created by hydrodynamical instabilities of a very cold but non-superfluid He jet. They can acquire a large amount of angular momentum during their passage through the nozzle of the experimental apparatus. When the droplets become superfluid, this angular momentum helps nucleate vortex arrays stable enough to be experimentally detected. [START_REF] Gomez | [END_REF]22] Since light absorption spectroscopy is very sensitive to the impurity environment, it is a good candidate as a tool for detecting impurities attached to vortices. We have therefore calculated the absorption spectrum of alkali (Ak) atoms from Li to Cs when attached to a vortex line in a superfluid droplet, and compared to the same spectrum with no vortex. The goal is to provide a quantitative answer to the size of the shift that can be expected, and therefore to the possibility for photon absorption experiments to detect the presence of vortices. The absorption spectrum of alkali atoms attached to He droplets has been thoroughly studied experimentally in the past, see e.g. Refs. 27-32, and earlier works have simulated photon absorption spectra of alkali atoms attached to vortex-free helium nanodroplets using different functionals and/or He-alkali pair potentials. [27,28,[33][34][35][36][37] We have repeated them in order to ensure comparison with vortex-hosting droplets with exactly the same parameters. To the best of our knowledge, the spectrum of alkali atoms attached to vortex-hosting droplets has not been previously addressed theoretically.

This article is organized as follows. In Sec. II we outline the 4 He Density Functional Theory ( 4 He-DFT) approach used to simulate doped 4 He droplets and their photon absorption spectrum. The results are presented and discussed in Sec. III, and a Summary and Outlook are given in Sec. IV.

II. THE 4 HE-DFT APPROACH

We have used the 4 He-DFT method to describe superfluid helium droplets. This approach has been recently reviewed and comprehensively described in Ref. 35, where the interested reader may find all the details for its application to pure and doped 4 He drops hosting or not vortex lines. We only give here the main features and some details to help following the discussion of the results.

Within 4 He-DFT, the total energy of a 4 He N droplet at zero temperature is written as a functional of the 4 He atomic density ρ(r)

E[ρ] = T [ρ] + E c [ρ] = h2 2m dr |∇Ψ(r)| 2 + dr E c [ρ]
(1) where T [ρ] is the kinetic energy, Ψ(r) an effective wave function defined by ρ(r) = |Ψ(r)| 2 , and the functional E c contains the He-He interaction term (within Hartree approximation) as well as additional terms describing nonlocal correlation effects. Strictu sensu, for vortex-free droplets Ψ(r) is a real function and hence Ψ(r) = ρ(r)

The results presented in this work have been obtained using the most accurate He functional to date, namely the Orsay-Trento functional (OT). [38] It includes a nonlocal correction to the kinetic energy and a backflow contribution. The backflow term only plays a role in the presence of density currents, which is the case in vortex-hosting droplets. Both contributions (especially the backflow term) make the OT functional computationally quite expensive. In addition to the computational cost, numerical problems arise in the case of helium samples doped with very attractive impurities. [35] Several variants of the OT functional have been designed and used to circumvent this problem. The so-called "solid functional" [39] has become very popular in static and time-dependent applications because it contains a "penalty term" to prevent the He density to become unphysically large when the helium sample is doped with very attractive impurities. The results obtained with this simplified functional are presented in the Supplementary Material. Because the backflow term should be the main difference between both functionals (the alkali-helium in-This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0008923

teraction is only weakly attractive), comparing the results in the Supplementary Material with those presented here will give an idea of the importance of this term for absorption spectra in the presence of a vortex.

For alkali-doped droplets Eq. ( 1) has to be modified to include the interaction of the droplet with the dopant. For the heavier alkalis, Ak ≡ K, Rb or Cs, which can be described as point-like, classical particles, this is simply done by adding the impurity-liquid interaction energy to the energy functional E [ρ]. This term, which acts as an external field for the helium density, is constructed within the pairwise sum approximation by integrating the impurity-He atom interaction, V Ak , over the liquid density,

E[ρ] → E[ρ] + dr ρ(r) V Ak (|r -r Ak |) (2)
where r Ak is the location of the impurity.

The droplet equilibrium configuration is obtained by solving the Euler-Lagrange equation deduced from functional minimization of Eq. ( 2) with respect to Ψ(r)

- h2 2m ∇ 2 + δE c δρ + V Ak (|r -r Ak |) Ψ(r) = µΨ(r) (3)
where µ is the The results presented in this work have been obtained using the 4He-DFT BCN-TLS computing package. [START_REF] Pi | 4He-DFT BCN-TLS: A Computer Package for[END_REF] Details on how Eq. ( 2) is solved can be found in Refs. 35, 42 and references therein. In short, Ψ(r) is written on a Cartesian grid and extensive use is made of the Fast-Fourier Transform to compute the convolutions entering the definition of E c . The droplet equilibrium configuration is obtained by imaginary-time relaxation for Ψ(r).

Including vortices does not require further developments in the formalism. A vortex line along the z axis for instance is obtained by the imprinting procedure, in which imaginary-time relaxation starts from the effective wave function [35]

Ψ(r) = ρ 1/2 0 (r) x 2 + y 2 (x + iy) (4)
where ρ 0 (r) is the density corresponding to a vortex-free droplet (note that in the presence of a vortex Ψ(r) must be complex). Imaginary-time relaxation eventually converges to a configuration which is an eigenstate of the angular momentum operator Lz with eigenvalue N h. We have defined the following energies to characterize the alkali-droplet equilibrium configurations: [35] Ak Ak@ 4 He1000 S Ak (Ak@V )@ 4 He1000 I: Influence of the presence of a vortex line V on the characteristics of an alkali doped 4 He1000 droplet: total energy E, and distance d of the impurity to the center of mass of the droplet. Ak@ 4 He1000 stands for the vortex-free droplet and (Ak@V )@ 4 He1000 for the droplet hosting a linear vortex. S Ak is the (vortex-free) impurity solvation energy defined in Eq. ( 5), and B Ak is the binding energy of the impurity to the vortex line defined in Eq. ( 7). The results presented in this table are for the OT functional. The energy of the pure 4 He1000 droplet is -5440.79 K when it is vortex-free, and -5318.08 K when it hosts a linear vortex, giving EV = 122.71 K for the energy for creating a vortex, Eq. ( 6). The asterisk in the first 3 lines of the table denotes a quantum treatment of the dopant.

B Ak E (K) d ( Å) (K) E (K) d ( Å) (K) Li * -
Solvation energy of the alkali:

S Ak = E(Ak@ 4 He N ) -E( 4 He N ) (5)
Vortex (V ) energy:

E V = E(V @ 4 He N ) -E( 4 He N ) (6)
Binding energy of the alkali to the vortex line:

B Ak = E[(Ak + V )@ 4 He N ] -E[(Ak@V )@ 4 He N ] (7) = S Ak + E V -{E[(Ak@V )@ 4 He N ] -E( 4 He N )}
A negative solvation energy implies that the impurity is attached to the droplet, and a positive binding energy of the Ak to the vortex implies that it is attached to the vortex line. This binding energy is the result of a delicate balance between the contributing terms, and the resulting values are rather small. Intuitively, it corresponds to the kinetic energy of the superfluid flow in the volume excluded by the impurity. [START_REF] Donnelly | Quantized vortices in helium II[END_REF] Figure 1 shows as an example a two-dimensional (2D) cut of the He density in a symmetry plane for a Cs-doped 4 He 1000 droplet, without or with a vortex. The appearance is qualitatively the same for all the alkali atoms. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. The vortex is visible as an empty line along the vertical axis. Its presence makes the dopant distance d to the droplet center of mass shorter, which is the results of two effects: the presence of a vortex flattens the droplet, and it makes deeper the dimple in which the alkali atom sits. This is general for all alkalis, as can be checked in Table I.

PLEASE CITE THIS ARTICLE AS

Table I collects the most relevant characteristics of the calculated equilibrium configurations. As can be seen there, the binding energy to the vortex line is positive for all the alkalis, indicating that all alkali-doped vortexhosting droplets are more stable if the alkali sits at the intersection of the vortex core with the surface. The binding energy increases with atomic number, except for the heaviest one studied here, Cs, which is slightly less bound to the vortex than Rb.

Lighter alkalis require a quantum treatment, as performed earlier in vortex-free droplets. [33,34,36] This was performed here for Li, Na, and K. Because of zero-point delocalization, they sit further away from the droplet surface when described quantum mechanically. This is quite a sizeable effect, especially for Li as can be seen in Table I. As a consequence, the Ak solvation energy in the vortex-free droplet (S Ak ) is reduced (in absolute value), as already observed. [36] The effect becomes smaller with increasing atomic number, as expected. The binding energy to the vortex line (B Ak ) increases in the quantum treatment, likely because the excluded volume already mentioned [START_REF] Donnelly | Quantized vortices in helium II[END_REF] increases due to zero point delocalization of the impurity.

As indicated before, we have also used the solid functional for this study. The resulting impurity distances to the droplet center of mass, solvation and vortex-binding energies differ only very slightly from the OT ones, as can be checked in Table I in the Supplementary Material, and the qualitative conclusions are not affected.

B. Absorption spectra

We now turn to the influence of a vortex on the absorption spectrum of an alkali-doped helium droplet. Once the equilibrium configuration of the system has been determined, the absorption spectrum is calculated as explained in Ref. 17. The required excited Π and Σ Ak-He pair potentials were taken from Pascale, [START_REF] Pascale | [END_REF] except in the case of Li, for which we have used those from Ref. 34. The method, semiclassical in nature, essentially yields the peak energies of the absorption spectrum and their relative intensity.

The width of the spectrum comes from fluctuations in the helium density and/or, to a lesser extent, from the probability distribution given by the square of the Ak wave function when it is treated quantum mechanically. To incorporate the helium density fluctuations, we use the DF sampling method described in detail in Refs. 42 and 44. The additional spectral width originating from the dopant zero-point delocalization was negligible in comparison and was not calculated.

np ← ns absorption spectra

Figure 2 shows the effect of a vortex on the np ← ns absorption spectrum of the alkali atoms attached to a 4 He 1000 droplet. It can be seen that the absorption peaks are further blue-shifted and broadened in the presence of a vortex line. In addition, treating the dopant quantum mechanically decreases the shift and the width as expected, since the impurity is on average further away from the droplet due to its zero point motion. This is illustrated in the figure for the case of Na and also happens This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS

PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0008923

for Li and K. The quantitative analysis of the corresponding spectral characteristics is presented in Table II.

Let us mention that the semiclassical approach is less justified in the case of Li, for which the absorption spectrum is dominated by bound-bound transitions. [27,34] Including them properly would require e.g. to carry out a Fourier analysis of the time-correlation function of the Li atom moving in the full 3D 2 Π 1/2 , 2 Π 3/2 , and 2 Σ 1/2 potentials. [34] Hence the corresponding spectral characteristics are omitted in Table II. The spectrum obtained semiclassically is still presented in Fig. 2 for the sake of completeness.

It is tempting to make a direct comparison with experimental spectra which have been measured for all the calculated transitions. However, quantitative agreement of peak positions and line shapes from the DFT as well as the quality of available interaction potentials is very limited and deviations from measured absorption lines are beyond the vortex-induced differences calculated in this work. In particular, shifts due to the presence of a vortex are significantly smaller than the widths of the peaks. This is all the more clear when examining Table III, where the difference ∆ between the shift with and without vortex is compared to the average FWHM. For all cases studied here, the average width is of the order of twice the difference in shifts. Having said this, comparing line positions and shapes is not suited to experimentally asses the presence of vortices for these transitions. Moreover, for very large droplets, coherent diffraction imaging results suggest that the proportion of droplets hosting a vortex is only on the order of a few percent, [START_REF] Gomez | [END_REF][21][22] implying that changes in the overall peak shape are at most on the same level. Although experimental spectra can readily be obtained on that level of signal-to-noise ratios, only methods being able to detect differences between vortexfree and vortex-hosting droplets would contain valuable information. As long as there are no methods at hand to specifically control the content of vortices of produced droplets in spectroscopic measurements, spectral changes will remain concealed.

(n + 1)s ← ns or (n + 1)p ← ns absorption spectra

The spectral characteristics of alkali-doped helium droplets are mostly governed by the usually repulsive interaction of their excited state orbital with the surrounding helium at the equilibrium configuration of the ground state. Since this orbital is more diffuse and more extended in higher (n > n) electronic states, its sensitivity to changes in the surrounding helium density is increased. Hence, absorption spectra to n > n excited states could reveal more differences between vortexhosting and vortex-free droplets than np ← ns ones.

We have tested this idea on three cases which have already been studied before in vortex-free droplets, namely K(5s ← 4s), Rb(6p ← 5s), and Cs(7s ← 6s) absorption spectra. [32,36,45,46] Fig. 3 shows the corresponding This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0008923 This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

(∆E) vortex-f (∆E) vortex-h FWHM (cm -1 ) Ak (cm -1 ) (cm -1 ) 2 Π 1/2 2 Π 3/2 2 Σ 1/2 2 Π 1/2 2 Π 3/2 2 Σ 1/2 2 Π 1/2 2 Π 3/2 2 Σ 1/2 vortex-f vortex-h vortex-f vortex-h vortex-f vortex-h Na * (3p ←
PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0008923

absorption spectra, and the shifts and FWHM are collected in the lower part of Table II. As expected, the difference in shifts between spectra originating from vortexhosting and vortex-free droplets are much larger than in the np ← ns spectra. They are now of the same order of magnitude as the spectral widths, as can be seen more clearly in Table III. This could make them distinguishable in the most favorable cases. Yet, several groups have examined (n + 1) ← n and higher transitions and did not report any unassigned shoulder or separated peak. [29][30][31][32] This is not surprising, providing, as discussed above, an expected contribution of vortices only on the percent level. Since the peak shapes and consequently also the far-out wings cannot be reproduced by the simulations on that level of accuracy, without the ability of measuring difference spectra it will be very difficult to confirm vortices in spectroscopic measurements also at these transitions.

Checking which case would be the most favorable one for vortex-line detection would require a systematic search for all the alkalis and their (n + 1) ← n or even (n + 2) ← n transitions, which in turn would call for systematic improvement of the commonly used Pascale's helium-alkali interactions in the excited electronic states; for instance, the He-Rb * (6p) interactions are notoriously too repulsive, resulting in too large spectral shifts. This systematic study is beyond the scope of the present work.

IV. SUMMARY AND OUTLOOK

We have investigated the possibility of using alkali absorption spectroscopy as a means of detecting the presence of vortices in a superfluid helium droplet. To this end, we have conducted 4 He-DFT simulations using the most accurate functional to date, namely the OT functional. For the sake of completeness, the absorption spectra of Figs. 2 and 3 obtained using the solid functional of Ref. 39, which lacks the backflow and non-local ki- The results show a blue shift and broadening of the absorption peaks in a vortex-hosting compared to a vortexfree droplet. However, these spectral modifications are very modest for np ← ns excitation: the shifts are significantly smaller than the peak widths. Since the proportion of droplets hosting a vortex is expected to be small, these modifications will be difficult to extract from experimental observations. In addition, the experimental droplet size distribution also affects the spectrum, which would make the identification of vortex presence even more difficult.

∆ Γ ∆ Γ (cm -1 ) (cm -1 ) (cm -1 ) (cm -1 ) (cm -1 ) (cm -1 ) Na * (3p ←
The situation is more favorable when turning to higher n p ← ns or n s ← ns (n > n) excitations. In this case the excited state orbital of the alkali atom is more diffuse and more extended, which makes it more sensitive to changes in its local environment. As a consequence, the corresponding n p or n s ← ns absorption spectra are more clearly affected, with larger blue shifts and broadening of the peaks, the shifts being then of the same order of magnitude as the widths. This could make experimental detection of vortex-hosting droplets feasible, at least if one finds ways to increase their proportion.

Note that in previous works simulating the capture of dopants by vortex-hosting droplets, dopants were found to spin around vortex lines rather than being attached to them. [26,47,48] Spectral modifications in the alkali absorption spectrum would be smaller in that case. However, these simulations were describing the first several hundred picoseconds of the capture dynamics. The alkali location used in the current work is that corresponding to its equilibrium location, i.e., at the vortex line-droplet surface intersection. The flight time of the droplet between the pickup cell and detection is of the order of millisecond, which gives plenty of time for the system to relax after the capture process and to reach its equilibrium configuration.

In droplets formed by supersonic expansion of helium gas into vacuum, vortex nucleation is expected to be scarce. When occurring it is due to collisions with impurities, as observed in studies of impurity capture by helium nanodroplets. [23,24] The faster the impurity, the more probable vortex creation is. On the other hand, if the impurity is too fast it can escape the droplet, even if the He-impurity interaction is attractive. Combining these arguments, one could imagine a two-step experiment where the amount of vortices could be changed: collision with fast heliophilic impurities to create vortices, then usual pickup of alkalis for detection through n s ← ns or n p ← ns (n ≥ n) absorption spectroscopy. Additional doping from the first step could be discriminated by spectroscopy. Certainly, the first step also induces changes in the droplet size distribution. However, corresponding effects are reasonably understood and can be modelled. [49] Moreover, cluster size dependencies of electronic spectra of alkalis appear to be weak beyond a This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset. PLEASE CITE THIS ARTICLE AS DOI:10.1063/5.0008923 few thousand helium atoms per droplet and can be accurately characterized. [50] Direct diffractive imaging of small droplets being impossible at present, the ultimate detection of vortices in helium droplets made of a few thousand atoms will be likely based on finding an observable sensitive to the changes caused by vortices in the impurity environment, that might be simultaneously studied experimentally and theoretically. Let us point out that femtosecond pump-probe photoionization of alkali atoms attached to helium droplets provides such observables, as e.g. desorption/fall-back times and final velocities of desorbed atoms and ions. [45] However, it does not seem that these observables will be more sensitive to the presence of vortices than the alkali absorption spectrum.

SUPPLEMENTARY MATERIAL

See supplementary material for the results obtained with the solid density functional of Ref. 39. This is the author's peer reviewed, accepted manuscript. However, the online version of record will be different from this version once it has been copyedited and typeset.

PLEASE CITE THIS ARTICLE AS

  FIG.1: 2D density in a plane of symmetry of a 4 He1000 droplet doped with a Cs atom (green dot). The color bar shows the He atom density in units of Å-3 . Top, vortex-free droplet; bottom, vortex-hosting droplet.

  FIG. 2: np ← ns absorption spectrum for the alkali-doped 4 He1000 droplet obtained with the OT functional. Li, Na, and K atoms are treated quantum-mechanically ([Q]), whereas Rb, and Cs are treated classically ([Cl]). For the sake of comparison, the results for Na obtained with the classical description are also shown. Black dashed line: vortex-free (vortex-f) droplet; red solid line: vortex-hosting (vortex-h) droplet. Gas-phase atomic transitions are represented by thin vertical lines. Units are arbitrary for absorption intensity, and cm -1 for excitation energy.

FIG. 3 :

 3 FIG. 3: Dipole absorption spectrum for K(5s ← 4s) (top), Rb(6p ← 5s) (middle), and Cs(7s ← 6s) (bottom) transitions in a vortex-free (black dashed line) or vortex-hosting (red solid line) 4 He1000 droplet obtained with the OT functional. The K atom is treated quantum mechanically and the Rb and Cs atoms classically. The thin vertical lines represent gasphase atomic transitions. Units are arbitrary for absorption intensity, and cm -1 for excitation energy.

  DOI:10.1063/5.0008923

  5443.32 27.70 -2.53 -5325.24 23.66 4.63 Na * -5446.07 26.55 -5.28 -5327.30 23.28 3.94 K * -5447.95 26.83 -7.16 -5328.93 23.51 3.69

	Li -5450.62 25.06 -9.83 -5330.66 21.55 2.75
	Na -5450.47 25.37 -9.68 -5330.57 21.99 2.81
	K -5451.08 26.07 -10.29 -5331.11 22.96 2.74
	Rb -5451.99 26.06 -11.20 -5332.19 22.97 2.91
	Cs -5451.29 26.60 -10.50 -5331.30 23.55 2.72
	TABLE

TABLE II :

 II Comparison of the np ← ns spectral characteristics of an alkali atom sitting on a vortex-free (vortex-f) and vortexhosting (vortex-h) 4 He1000 nanodroplet. 2 nd to 4 th columns: energy shifts (∆E) vortex-f for the vortex-free droplet with respect to the gas-phase atomic lines, calculated from the absorption maxima for the 2 Π 1/2 , 2 Π 3/2 , and 2 Σ 1/2 final states; 5 th to 7 th columns: energy shifts (∆E) vortex-h for the vortex-hosting droplet; last 6 columns: full width at half maximum (FWHM) for each peak for the vortex-free or vortex-hosting droplet. The top part of the table is for np ← ns absorption spectra and the bottom part, separated by a horizontal line, is for (n + 1)s ← ns or (n + 1)p ← ns absorption spectra. The asterisk denotes a quantum treatment of the dopant.

	3s) 6.4	3.1	58.3 12.9 18.5	76.5	10	21	18	43	111	132
	K * (4p ← 4s) 7.6	-4.5	66.4 14.8	7.0	97.8	13	22	14	34	120	135
	Na (3p ← 3s) 17.4 29.9 140.0 33.0 65.9 171.1	28	51	63	99	204	212
	K (4p ← 4s) 13.4	3.9 116.8 25.1 26.7 151.2	21	36	34	62	178	185
	Rb (5p ← 5s) 51.2 11.0 138.6 74.0 37.1 169.7	43	59	37	63	176	185
	Cs (6p ← 6s) 62.6	9.5 131.3 88.9 34.0 168.5d	57	71	31	53	158	167
	K * (5s ← 4s) -	-	618.9	-	-	880.6	-	-	-	-	242	270
	K (5s ← 4s)	-	-	839.4	-	-1079.3	-	-	-	-	282	307
	Rb(6p ← 5s) 942.6 1101.3 3205.7 1518.0 1766.8 3535.6 314	421	326	428	743	761
	Cs(7s ← 6s)	-	-1070.3 -	-1402.1	-	-	-	-	384	434

TABLE III :

 III Comparison of the energy shift difference ∆ between the vortex-hosting and vortex-free droplet and Γ, the half sum of the FWHM of the spectra with or without vortex, calculated from data in TableII.netic energy terms, are presented in the Supplementary Material, together with the corresponding table of their spectral characteristics.

Π 1/2 2 Π

3/2 2 Σ 1/2 Ak ∆ Γ

Acknowledgments

This work has been performed under Grant No FIS2017-87801-P (AEI/FEDER, UE). M.B. and E.G.A. thank the Université Fédérale Toulouse Midi-Pyrénées for financial support throughout the "Chaires d'Attractivité 2014" Programme IMDYNHE. A computer grant from CALMIP high performance computer center (grant P1039) is gratefully acknowledged.

Data availability statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.