# Bound states in a superfluid vortex: A detailed study along the BCS-BEC crossover

###### Abstract

The bound states that can occur in a superfluid vortex have recently called for attention owing to the capability of detecting them experimentally. However, a detailed theoretical account for the presence of these vortex bound states is still lacking, for all temperatures in the superfluid phase and couplings along the BCS-BEC crossover. Here, we fill this gap and present a systematic theoretical study based on the Bogoliubov-de Gennes equations for the bound states that occur over the two characteristic (inner and outer) spatial ranges in which the extension of a superfluid vortex can be partitioned. It is found that the total number of bound states decreases from the BCS (weak-coupling) side of the crossover toward the intermediate-coupling region where they are still present, whereas the bound states disappear upon entering the BEC (strong-coupling) side. A scaling relation is also obtained that connects the number of bound states in the inner spatial range of the vortex to the depth and width of the vortex itself. A criterion is finally provided in terms of the local density of states, to distinguish where a given fermionic superfluid is located in the coupling-temperature phase diagram of the BCS-BEC crossover.

## I Introduction

Recently, interest has arisen in the internal structure of vortices in type-II superconductors, owing to the presence of bound states in the vortex core region of a Y123 superconductor that were detected experimentally Berthod-2017 (). The relevance of this finding has been highlighted SX-2017 (), as an indication that the superconducting state of high-temperature superconducting materials should be well described by the conventional BCS pairing theory BCS-1957 (), in spite of the fact that the normal state of these materials is highly unconventional.

In Ref. Berthod-2017 () the experimental findings were supported by theoretical calculations, based on the approach to the solution of the Bogoliubov-de Gennes (BdG) equations DeGennes-1966 () developed in Ref. Berthod-2016 (), to come up with the local density-of-states as the relevant quantity to be compared with the experimental data. Although these calculations considered the effects of a realistic band structure as well as the influence of disorder and nearby vortices on a given vortex, they were limited to zero temperature and to a single value of the inter-particle coupling. To understand in a more complete fashion the properties of bound states in superconducting vortices, it would then be desirable to characterize them also as a function of temperature and inter-particle coupling, at least for the simplest case of an isolated vortex. Purpose of this paper is to address this aspect of the problem, thereby complementing to some extent the theoretical information provided in Ref. Berthod-2017 ().

The occurrence of bound states in a superconducting vortex was originally proposed in Ref. CdGM-1964 () by solving the BdG equations. That calculation was restricted to zero temperature and to the case of weak coupling for which , where is the bulk value of the gap parameter away from the center of the vortex and the Fermi energy. In more modern language, this would correspond to the BCS (weak-coupling) limit of the BCS-BEC crossover Physics-Reports-2018 (). Later works, that considered the self-consistent solution of the BdG equations for a vortex line, have either relaxed the restriction to zero temperature but still at weak coupling Schluter-1991 (), or spanned the whole BCS-BEC crossover but still at zero temperature Sensarma-2006 (); Levin-2006 (). To obtain a more complete information on the properties of bound states in a superconducting vortex, however, it would be worthwhile to address both their coupling and temperature dependences at the same time.

From a technical point of view, in Refs. Schluter-1991 (); Sensarma-2006 (); Levin-2006 () the study of a single vortex was performed in terms of the BdG equations by constraining the vortex into a cylinder of radius with infinite walls, such that all single-particle wave functions were to vanish at . In this way, however, all single-particle energies turn out to be discrete, in such a way that no sharp distinction there exists between truly bound states with finite spatial extension and unbound continuum states past a well-defined threshold. This limitation was later removed in Ref. SPS-2013 (), where the BdG equations were solved by introducing suitable Òfree boundary conditionsÓ through which the vortex profile is matched to its asymptotic behavior at a large distance from the center. In this way, a clear distinction can be marked between bound (localized) states and unbound states extending to infinity. By this approach, the temperature dependence of the healing length associated with the size of the vortex could be determined even for temperatures quite close to the critical temperature at which diverges. In Ref. SPS-2013 () this analysis was carried out from weak to strong coupling across the whole BCS-BEC crossover, to determine not only the profile of the gap parameter but also those of the local number density and current. In addition, in Ref. SPS-2013 () an advanced regularization procedure was developed for the self-consistent gap equation, which permits to increase considerably the accuracy of the numerical calculations.

In the present paper, we take advantage of the above procedures developed in Ref. SPS-2013 () and concentrate our analysis on the bound-state part of the spectrum for an isolated vortex embedded in an otherwise homogeneous fermionic superfluid of infinite extent. In this context, we will determine the local density of states as a function of energy and of the spatial position from the vortex core, for varying temperature in the superfluid phase from to and for couplings along the BCS-BEC crossover. We will thus also be able to study how the total number of bound states evolves as a function of temperature and coupling, paying additional attention to how the bound states are distributed over the two characteristic spatial ranges in which the vortex (including its tail) can be partitioned.

The main results obtained in this paper are as follows:

(i) By relying on the numerical procedure to solve the BdG equations for a fermionic superfluid vortex developed in Ref. SPS-2013 (), an energy threshold is identified that clearly separates truly bound states from continuum states. In this way, the counting of bound states for given coupling and temperature becomes a meaningful process.

(ii) In addition, by partitioning the spatial extent of the vortex into an inner and an outer region (to be identified below), the number of bound states that hinge on each region is determined as a function of energy below threshold. It is found that the number of bound states in the inner region scales in a universal way on the depth and width of the vortex itself.

(iii) The number of bound states is found to be non-negligible even in the intermediate-coupling region of the crossover, where the Cooper pair size becomes comparable with the inter-particle distance. To the extent that this small size is compatible with the occurrence of high-temperature superconductivity Dagotto-1994 (), there is thus no a priori reason to associate the occurrence of bound states in a vortex with the superconductivity being of the conventional BCS type as recently asserted SX-2017 ().

(iv) Finally, it is proposed that the shape of the local density of states vs energy, taken about the center of the vortex, can serve as a guide to distinguish where a given fermionic superfluid stands in the coupling-temperature phase diagram of the BCS-BEC crossover. In this context, one should recall that the local density of states has accurately been measured with a scanning tunneling microscope also in high-temperature superconductors RMP-2007 (), as well as in more complex SNS structures made with conventional superconductors Esteve-2008 ().

The plan of the paper is as follows. Section II recalls the main features of the solution of the BdG equations drawn from Ref. SPS-2013 () for an isolated vortex embedded in an infinite superfluid. Emphasis is there given to the occurrence of an inner and outer region with different spatial behaviors of the vortex profile where the bound states can reside. Section III presents the numerical results for the number of bound states and the local density of states, from which the scaling relation and the criterion for spotting the position along the BCS-BEC phase diagram mentioned above are derived. Section IV gives our conclusions. Appendix A discusses the asymptotic profile of an isolated vortex in the context of the Gross-Pitaevskii equation, which holds for the composite bosons that form in the BEC limit of the crossover.

## Ii Formal aspects

We begin by briefly recalling the method used in Ref. SPS-2013 () to solve the BdG equations in the presence of an isolated vortex embedded in an otherwise infinite homogeneous superfluid, as a function of both coupling across the BCS-BEC crossover and temperature.

The BdG equations read DeGennes-1966 ():

(1) |

Here, where is the fermion mass and the chemical potential (we set troughout). The local gap parameter in Eq. (1) is determined via the self-consistent condition:

(2) |

where is the Fermi function at temperature ( being Boltzmann constant) and is the bare coupling constant of the contact interaction. A suitable regularization of the inhomogeneous gap equation (2) was implemented in Appendix B of Ref. SPS-2013 (), by drawing elements from the derivation of the Gross-Pitaevskii equation Pitaevskii_Stringari-2003 () for composite bosons that form in the BEC limit, which was obtained in Ref. Pieri-2003 () starting from the BdG equations. In the process, the bare coupling constant that enters Eq. (2) gets replaced by the (dimensionless) coupling parameter , where is the Fermi wave vector associated with the (bulk) number density and is the scattering length of the two-fermion problem. This coupling parameter enables one to span the whole BCS-BEC crossover Physics-Reports-2018 (). In practice, the crossover between the BCS and BEC regimes is exhausted within the range about the unitary limit (UL) where .

When looking for bound states inside a vortex like in the present context, coupling values on the BCS side of unitarity should primarily be considered, since it will turn out that the number of bound states rapidly vanishes upon entering the BEC side of unitarity. In addition, on the BEC side of the crossover at finite temperature it would be necessary to include pairing fluctuations beyond mean field Physics-Reports-2018 (), a task which is beyond the purposes of the present paper. On the other hand, on the BCS side of the crossover we shall find it necessary to extend the calculations down to the (numerically rather demanding) coupling value , in order to recover what would be expected on the basis of the standard BCS approach.

For an isolated vortex with unit circulation and directed along the axis, the spatially dependent gap parameter is written in cylindrical coordinates as with real. Correspondingly, the eigenfunctions of the BdG equations (1) take the form:

(3) | |||||

(4) |

where is the radial quantum number, is an integer (both positive and negative), and are real functions for the bound states of primary interest here. Within this approach, for given values of and the BdG equations are numerically integrated outwards, starting from with suitable indicial conditions, up to a maximum value (values ranging from for up to for proves sufficient for practical purposes). At that point, “free” boundary conditions are enforced at by using suitable linear combinations of Bessel, Neumann, and Hankel functions, in terms of which the solutions of the BdG equations can be expressed for . This procedure offers a definite advantage, in that it avoids the common practice of constraining the vortex in a cylinder with infinite walls at . For this reason, this procedure allows one to clearly distinguish truly bound states (which decay exponentially for ) from continuum states (with oscillating behavior extending up to infinity). In practice, we have considered values of not smaller than and values of up to .

In this way, for given and irrespective of the value of , a first energy threshold between bound and continuum states occurs at the (temperature-dependent) value recovered by deep in the bulk region when (where ), while a second threshold occurs at when . Accordingly, for only bound states can be found, for bound states embedded in the continuum (with different values of ) can also be found, and for only continuum states occur.

Once the eigenfunctions have been determined together with the corresponding eigenvalues according to the above procedure, one can obtain the local density of states given by:

(5) | |||||

which has dimensions of [volumeenergy]. This quantity, which is of primary experimental interest RMP-2007 (); Esteve-2008 (), contains contributions from both bound and continuum states footnote (). We are now in a position to calculate as a function of coupling along the BCS-BEC crossover and of temperature in the superfluid phase. In addition, by our approach, we can also distinguish the separate contributions to from bound and continuum states.

It is further relevant to point out that, contrary to a common assumption that the radial vortex profile approaches asymptotically the bulk value in an exponential way Berthod-2016 (); Tempere-2017 (), the gap parameter has instead a long tail with a power-law dependence of the type . This property, which went apparently unnoticed in the literature, was found to be valid for all couplings throughout the BCS-BEC crossover by the detailed numerical analysis carried out in Ref. SPS-2013 (). As shown in Appendix A, this property can be also captured analytically in terms of the Gross-Pitaevskii equation Pitaevskii_Stringari-2003 (), to which the BdG equations have been shown to reduce in the BEC (strong-coupling) limit Pieri-2003 (). It is actually in the more limited radial range only that the vortex profile has an exponential dependence, from which a characteristic (temperature and coupling dependent) coherence length can be determined, where marks the position of the maximum value of the radial current SPS-2013 (). As is turns out from the numerical analysis of Ref. SPS-2013 () (and confirmed by from the present one), the two length scales and coincide with each other within numerical error for all couplings and temperatures.

Accordingly, in the following we shall find it convenient to partition the vortex profile into an “inner” region (where the exponential decay applies) and an “outer” region (where the power-law behavior takes over), the boundary between the two regions being (approximately) taken at . The counting of the number of bound states will similarly be partitioned. We are then going to distinguish whether a given bound state with quantum numbers () belongs to the inner or to the outer region, by looking at the (dimensionless) partial normalization condition:

(6) |

We then attribute the location of the bound state to the inner (outer) region when exceeds (falls behind), say, as compared to the value of the complete normalization when the upper limit of the integral in Eq. (6) extends to infinity. Owing to the () power-law tail of the gap parameter, it turns out that the vortices residing in the outer region by far exceed in number those residing in the inner region, although most of them lie far away from the vortex center.

In practice, the bound states in the inner region (and therefore their number) are numerically much more under control than those in the outer region, especially when the latter ones hinge mostly close to the boundary at where the connection with decaying exponentials is enforced through the boundary conditions. In addition, when several vortices would be packed up together to form a vortex lattice, the nature of the bound states in the remote external part of the outer region of what would have been a vortex in isolation is expected to be strongly modified by the presence of the surrounding vortices. For this reason, it may anyway not be worth to insist in extending the size of the outer region beyond those values of that we have considered in practice for the numerical calculations.

## Iii Numerical results

We pass to characterize the presence of bound states in an isolated vortex, for all temperatures in the superfluid phase and couplings along the BCS-BEC crossover, by taking advantage of our accurate numerical calculations in terms of the BdG equations following the methodologies developed in Ref. SPS-2013 (). We shall find that the number of bound states rapidly decreases on the BEC side of unitarity, so that, in practice, only the coupling region from the BCS regime up to unitarity will be relevant to our purposes. Several features will be highlighted by looking at the details of the local density of states, or simply by counting the number of bound states that can be identified for given coupling and temperature. Our main findings are summarized by the following features.

A. Contribution from bound and scattering states to the local density of states

One of the main questions to be asked at the outset is where the various contributions, in which the local density of states of Eq. (5) can be partitioned, are spatially located over the extension of an isolated vortex. These contributions can be identified as those corresponding to the bound states residing essentially in the inner and outer regions of the vortex, and as that corresponding to the scattering states. Here, the partial normalization criterion (6) is used to assign a given bound state either to the inner or outer regions of the vortex (see also subsection III-B below).

Accordingly, we express the local density of states of Eq. (5) as the sum of three contributions, namely, , , and , where “ib” stands for inner-bound, “ob” for outer-bound, and “sc” for scattering. These three contributions are reported in Figs. 1, 2, and 3, respectively, as a function of both and , in each case for three characteristic couplings and two significant temperatures . [One should, however, be aware that, due to the progressive importance of pairing fluctuations beyond mean field at finite temperature for increasing coupling, the data at unitarity and finite temperature are expected to be quantitatively less reliable than those at weaker couplings.]

0.0 | 0.6864 | 0.3445 | 0.4965 |

-1.5 | 0.1001 | 0.0523 | 0.0579 |

-3.0 | 0.0097 | 0.0051 | 0.0055 |

From Figs. 1-3 one notices that the contribution is mostly spatially concentrated over the inner portion of the vortex, while the contribution is spread over a much wider spatial region. On the other hand, at low energies above (the temperature dependent) threshold is depressed in the central region of the vortex owing to the orthogonality requirement of quantum mechanical states with different energies, while at high energies above threshold is almost uniformly spread over all space owing to the occurrence of spatial oscillations with small wavelength. The situation remains qualitatively the same also for increasing temperature from to . More marked differences appear instead for increasing coupling, from deep in the weak-coupling (BCS) region to unitarity where , along which one notices a progressively enhanced asymmetry between negative and positive energies. We shall return below to a more detailed discussion of this asymmetry.

From Figs. 1-3 one also notices that, for each coupling and temperature, a sharp threshold for the continuum states occurs at (cf. Fig. 3), and that bound states “embedded” in the continuum are also found when (cf. Figs. 1-2), although they give a small overall contribution to . This can be seen, in particular, in the left-bottom panel of Fig. 1, where the two arrows mark the values of the continuum threshold beyond which bound states embedded in the continuum appear. For convenience, the numerical values of relevant to Figs. 1-3 are reported in Table 1.

We have further verified that, past unitarity upon approaching the BEC regime, there occurs a progressive depletion of in the central region of the vortex, implying that bound states tend to disappear in this coupling regime. In addition, for given coupling this depletion is essentially complete for , while it becomes only partial upon increasing the temperature. We shall return below to a more detailed discussion of this feature.

B. Counting the number of bound states in the inner and outer regions

Another relevant piece of information can be obtained by a simple count of the number of bound states, separately in the inner and outer regions as defined above. To this end, we recall that the bound-state wave functions are localized in the - plane orthogonal to the axis of the vortex, so that a whole branch spanned by the wave vector is associated with each bound state with given angular momentum and radial quantum number . Accordingly, in the following the numbers of bound states will be given per unit length along the direction of the vortex axis, since they will be summed (integrated) over the wave vector (as it is also done in Eq. (5)). With this procedure, the line density of bound states and , that reside, respectively, in the inner and outer regions of an isolated vortex, can be obtained as a function of coupling and temperature.

Figure 4 shows both (lower panel) and (upper panel) for three characteristic temperatures , with the coupling ranging from deep in the BCS regime up to unitarity. In addition, for the calculation has been extended to on the BEC side of unitarity.

A few general features can be evidenced from these plots: (i) Both and rapidly decrease from the BCS regime to the UL and quickly vanish once the BEC side past unitarity is reached, thus implying that the bound states have no relevance for the composite bosons that form out of fermion pairs in the BEC regime; (ii) For given coupling, increases upon increasing the temperature, whereby the depth of the vortex decreases but at the same time its width increases; (iii) For given coupling and temperature, by far exceeds , although appears to be essentially independent of temperature.

It should be remarked that the values of we have obtained are numerically less reliable that those of , because they depend on the value of the radius where the boundary conditions are enforced on the solutions of the BdG equations (cf. Section II). In this context, we found it numerically too demanding to make this radius exceed for the smaller coupling we could reach deep in the BCS regime. For these reasons, in the following we shall limit ourselves to further examine the dependence of on and , leaving aside similar considerations for .

C. Symmetry vs asymmetry of the local density of states about the vortex center

A characteristic feature, that distinguishes whether a homogeneous attractive Fermi gas lies in the weak-coupling (BCS) or the unitary regime, is the degree of asymmetry of the single-particle spectral function, with the two peaks of this spectral function being much more symmetric about zero frequency in the BCS regime than at the UL PPPS-2012 (). The occurrence of particle-hole asymmetry is a characteristic feature of the BCS-BEC crossover, whereby the underlying Fermi surface gets progressively modified when passing from the BCS to the BEC limits across UL, as the ratio increases. Our interest here is to verify to what extent this characteristic feature is maintained also for the present inhomogeneous problem, such that a measurement of the local density of states about the vortex center would be able to detect whether or not the superfluid properties of the Fermi gas deviate from those of an ordinary BCS superconductor.

As we have already noted while commenting Figs. 1-3, all three contributions to the local density of states that we have considered become progressively more asymmetric when passing from the BCS to the unitary regime. A more quantitative characterization of this asymmetry can be obtained as follows. Let us consider a spatial region in the inner part of the vortex, say, of linear extension around the vortex center at (in such a way to concentrate our attention close to the center of the vortex). We then look for the behavior of the integrated quantity

(7) |

where is the local density of states (5). The quantity (7) has dimensions of [lengthenergy] and contains, in principle, contributions from both bound and continuum states, depending on the value of . In the following, we shall be interested in the energy interval where is the bulk value of the gap parameter at temperature , such that only bound states contribute to in this energy interval.

Plots of vs are shown in Fig. 5 for three characteristic couplings and three characteristic temperatures for each coupling (with the exception of the coupling for which is lacking). Irrespective of temperature, drastic changes are seen to occur in the shape of , which becomes progressively more asymmetric when passing from (upper panel) to (lower panel).

An even sharper characterization, about the development of this asymmetry when approaching the UL from the BCS side, can be obtained by integrating over either the positive or negative portions of the gap region, thereby defining

(8) | |||||

(9) |

which have dimensions of [length]. In this way, a suitable asymmetry parameter can then be introduced, which can be analyzed as a function of coupling and temperature. The results are reported in Fig. 6. One sees that, essentially for all temperatures, there occurs a steep increase of this asymmetry parameter, from the value zero in the BCS regime to almost unity at UL, which at has reached half the way at about the coupling . This quantity thus provides one with a quick clue for the closeness of the superfluid Fermi gas to the UL, by looking at the distribution of the bound states about the vortex center.

From the experimental results available from Ref. Berthod-2017 (), no clear evidence for the occurrence of particle-hole asymmetry can apparently be extracted. One has to consider, however, that in Ref. Berthod-2017 () the presence of bound states in vortices was evidenced in a high- cuprate superconductor for the first time, through a delicate subtraction procedure of a significant background. Additional detailed measurements of density-of-states spectra would then be required to evidence a particle-hole asymmetry over and above this subtraction.

D. A scaling relation

In Fig. 4, was found to depend separately on coupling and temperature. However, both the depth (or bulk value of the gap parameter) and the width of the vortex also depend on coupling and temperature. It will then be interesting to eliminate this double dependence and express directly as a function of a single variable obtained by combining and , irrespective of the corresponding values of coupling and temperature. Through this attempt, we have found that the product is the appropriate scaling variable on which the line density of bound states effectively depends.

Accordingly, we have organized the numerical values from Fig. 4 into the single plot shown in Fig. 7, where is reported vs over a double-log scale so as to put on equal footing quite different sets of data (spanning four decades, from to ). Here, is is units of the inverse of the Fermi wave vector and in units of the Fermi energy , such that the variable of Fig. 7 is dimensionless.

The resulting linear scaling dependence of vs (which is evidenced in the figure by the straight line ) appears quite remarkable, to the extent that essentially all data with individually quite different values of and fall on this straight line. We have further verified that this scaling relation between and holds also when the counting of bound states residing in the inner region is enlarged, by shifting , e.g., to in Eq. (6), provided that the same replacement is also made in the variable on which depends.

## Iv Concluding remarks

In this paper, we have given a detailed account about the number, energy location, spatial location, and shape of the large number of bound states that are present within an isolated vortex embedded in an otherwise homogeneous superfluid. We have done this by an accurate solution of the BdG equations based on the methods developed in Ref. SPS-2013 (), for a dilute Fermi gas spanning the BCS-BEC crossover with the temperature ranging from zero to the superfluid critical temperature. In this way, we have been able to provide a criterion for locating where a given Fermi gas lies along this crossover, in terms of the asymmetry between positive and negative energies of the local density of states. From our numerical calculations we were also able to extract a universal scaling relation, that relates the number of bound states in the inner region of the vortex with the depth and width of the vortex itself.

We have further verified that bound states occur in the interior of a vortex for couplings on the BCS side of unitary, and rapidly disappear when entering the BEC regime past unitarity. Consistently with this theoretical finding, a clear experimental finding for the occurrence of bound states in a vortex in a superconducting material Berthod-2017 () can indeed be greeted as a signature that the superfluid phase of that material may be described by the conventional BCS theory SX-2017 (). Yet, it is also known from the analysis of Ref. SPS-2013 () that it is rather the continuum part of the spectrum of a vortex (as obtained by solving the BdG equations) that exhausts, in practice, most part of the contribution to physical quantities, such as the profile of the gap parameter and the local number density.

In this paper, the analysis has been limited to considering an -wave order parameter, and the question may arise about the role of the symmetry of the order parameter on our results. In this context, we may refer to the results of Ref. Kato-2001 (), where a BdG calculation for an isolated vortex with a -wave order parameter was reported, with the conclusion that the spatial dependence of the quasi-particle density of states is similar to the one with an -wave order parameter, except possibly for the contribution from the extended states.

A further comment is in order on the existence of bound states in the continuum (usually referred to as BICs) that we have found in an isolated vortex for symmetry reasons, owing to the different values of the quantum between these bound states and the continuum states in which they are embedded. As a consequence, a slight geometric perturbation of the vortex line along the -direction could turn BIC states into decaying resonances FW-BIC-1985 (). In this context, it is interesting to mention that BIC states have recently gained considerably experimental and technological interest as a generic wave phenomenon which can occur in supercavity lasing Kodigala-2017 ().

On physical grounds, the Caroli-de Gennes-Matricon bound states that are present in the interior of a vortex have a similar origin as the Andreev-Saint-James states that show up as subgap states in the context of the Josephson effect Andreev-1964 (); Saint-James-1964 (); Andreev-1965 () (for a more recent review dealing also with cuprate superconductors, see Ref. Deutscher-2005 ()). Also for the Andreev-Saint-James states, a systematic theoretical study of their occurrence has been performed throughout the BCS-BEC crossover by solving the BdG equations at zero temperature SPS-2010 (). It would then be interesting to address experimentally the occurrence of these subgap states also along the way of the BCS-BEC crossover, in terms of suitable spectroscopic measurements Esteve-2013 ().

ACKNOWLEDGMENTS

This work was partially supported by the Italian MIUR under Contract PRIN-2015 No. 2015C5SEJJ001.

## Appendix A Internal Structure of a Vortex

As mentioned in Section II, the accurate numerical solution of the BdG equations performed in Ref. SPS-2013 () has evidenced the feature that the radial profile of the gap parameter has a () power-law tail for all couplings throughout the BCS-BEC crossover (as well as for all temperatures in the superfluid phase). In this Appendix we show that, in the BEC (strong-coupling) limit of the BCS-BEC crossover, whereby the fermionic BdG equations reduce to the bosonic Gross-Pitaevskii (GP) equation for the composite bosons that form in this limit Pieri-2003 (), the () long-range behavior of the condensate wave function can be determined by simple analytic considerations. Although this result has already been reported for a vortex filament in an almost ideal Bose gas described at low temperature by the GP equation LP-1961-2006 (), the reason to briefly discuss it here is that its relevance for a vortex in a fermionic superfluid described by the BdG equations has passed essentially unnoticed in the literature Berthod-2017 (); Berthod-2016 (); Tempere-2017 ().

The (time-independent) GP equation reads Pitaevskii_Stringari-2003 ():

(10) |

where , , and are the bosonic mass, scattering length, and chemical potential, respectively. [For composite bosons built up in terms of superfluid fermions described by the BdG equations, , , and where is the binding energy of the composite bosons Pieri-2003 ().]

For an isolated vortex filament directed along the axis, one sets and writes generically the vortex solution in cylindrical coordinates in the form , where is an integer referred to as the topological charge of the flow (we will set at the end of the calculation). Far away from the vortex axis, one expects the local bosonic density to reach its bulk value , such that where . Introducing at this point the bosonic healing length and the rescaled radial variable , the GP equation (10) acquires the form:

(11) |

To determine the asymptotic behavior of for , one sets

(12) |

with , and obtains the constants and by entering the approximate form (12) into Eq. (11). At the order here considered, one gets from Eq. (11) the algebraic condition:

(13) |

which to the leading order yields and . For it results that tends to unity as , as reported in Ref. LP-1961-2006 (). Note that, physically, this long-range behavior is dominated by the angular kinetic energy of the vortex.

## References

- (1) C. Berthod, I. Maggio-Aprile, J. Bruér, A. Erb, and C. Renner, Observation of Caroli-deGennes-Matricon vortex states in , Phys. Rev. Lett. 119, 237001 (2017).
- (2) C.-L. Song and Q.-K. Xuey, Cuprate superconductors may be conventional after all, Physics 10, 129 (2017).
- (3) J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of superconductivity, Phys. Rev. 108, 1175 (1957).
- (4) P. G. de Gennes, Superconductivity of Metals and Alloys (Benjamin, New York, 1966), Chap. 5.
- (5) C. Berthod, Vortex spectroscopy in the vortex glass: A real-space numerical approach, Phys. Rev. B 94, 184510 (2016).
- (6) C. Caroli, P. G. De Gennes, and J. Matricon, Bound fermion states on a vortex line in a type II superconductor, Phys. Lett. 9, 307 (1964).
- (7) G. C. Strinati, P. Pieri, G. Röpke, P. Schuck, and M. Urban, The BCS-BEC crossover: From ultra-cold Fermi gases to nuclear systems, Phys. Rep. 738, 1 (2018).
- (8) F. Gygi and M. Schlüter, Self-consistent electronic structure of a vortex line in a type-II superconductor, Phys. Rev. B 43, 7609 (1991).
- (9) R. Sensarma, M. Randeria, and T-L. Ho, Vortices in superfluid Fermi gases through the BEC to BCS crossover, Phys. Rev. Lett. 96, 090403 (2006).
- (10) C.-C. Chien, Y. He, Q. Chen, and K. Levin, Ground-state description of a single vortex in an atomic Fermi gas: From BCS to Bose-Einstein condensation, Phys. Rev. A 73, 041603(R) (2006).
- (11) S. Simonucci, P. Pieri, and G. C. Strinati, Temperature dependence of a vortex in a superfluid Fermi gas, Phys. Rev. B 87, 214507 (2013).
- (12) See, e.g., E. Dagotto, Correlated electrons in high-temperature superconductors, Rev. Mod. Phys. 66, 763 (1994).
- (13) Ø. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod, and C. Renner, Scanning tunneling spectroscopy of high-temperature superconductors, Rev. Mod. Phys. 79, 353 (2007).
- (14) H. le Sueur, P. Joyez, H. Pothier, C. Urbina, and D. Esteve, Phase controlled superconducting proximity effect probed by tunneling spectroscopy, Phys. Rev. Lett. 100, 197002 (2008).
- (15) L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Clarendon Press, Oxford, 2003).
- (16) P. Pieri and G. C. Strinati, Derivation of the Gross-Pitaevskii equation for condensed bosons from the Bogoliubov-de Gennes equations for superfluid fermions, Phys. Rev. Lett. 91, 030401 (2003).
- (17) For continuum states above threshold, the quantum number is replaced by the energy eigenvalue and the in Eq. (5) is replaced by an integral over .
- (18) N. Verhelst, S. Klimin, and J. Tempere, Verification of an analytic fit for the vortex core profile in superfluid Fermi gases, Physica C: Superconductivity and its Applications 533, 96 (2017).
- (19) F. Palestini, A. Perali, P. Pieri, and G. C. Strinati, Dispersions, weights, and widths of the single-particle spectral function in the normal phase of a Fermi gas, Phys. Rev. B 85, 024517 (2012).
- (20) M. Kato and K. Maki, Bound states and extended states around vortices in -wave superconductors, J. Magn. Magn. Mater. 226-230, 280 (2001).
- (21) H. Friedrich and D. Wintgen, Interfering resonances and bound states in the continuum, Phys. Rev. A 32, 3231 (1985).
- (22) A. Kodigala, T. Lepetit, Q. Gu, B. Bahari, Y. Fainman, and B. Kanté, Lasing action from photonic bound states in continuum, Nature 541, 196 (2017).
- (23) A. Andreev, The thermal conductivity of the intermediate state in superconductors, Sov. Phys. JETP 19, 1228 (1964) [Zh. Eksp. Teor. Fiz. 46, 1823 (1964)].
- (24) D. Saint-James, Elementary excitations in the vicinity of the surface separating a normal metal and a superconducting metal, J. Phys. France 25, 899 (1964).
- (25) A. Andreev, Electron spectrum of the intermediate state of superconductors, Sov. Phys. JETP 22, 455 (1966) [Zh. Eksp. Teor. Fiz. 49, 655 (1965)].
- (26) G. Deutscher, Andreev-Saint-James reflections: A probe of cuprate superconductors, Rev. Mod. Phys. 77, 109 (2005).
- (27) A. Spuntarelli, P. Pieri, and G. C. Strinati, Solution of the Bogoliubov-de Gennes equations at zero temperature throughout the BCS-BEC crossover: Josephson and related effects, Phys. Rep. 488, 111 (2010).
- (28) L. Bretheau, Ç. Ö. Girit, C. Urbina, D. Esteve, and H. Pothier, Supercurrent spectroscopy of Andreev states, Phys. Rev. X 3, 041034 (2013).
- (29) L. P. Pitaevskii, Vortex lines in an imperfect Bose gas, Sov. Phys. JETP 13, 451 (1961) [Zh. Eksp. Teor. Fiz. 40, 646 (1961)]; see also, E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2, Theory of the Condensed State (Butterworth-Heinemann, Oxford, 2006), Sec. 30.