A multipole-expanded effective field theory for vortex ring-sound interactions

The low-energy dynamics of a zero temperature superfluid or of the compressional modes of an ordinary fluid can be described by a simple effective theory for a scalar field — the superfluid ‘phase’. However, when vortex lines are present, to describe all interactions in a local fashion one has to switch to a magnetic-type dual two-form description, which comes with six degrees of freedom (in place of one) and an associated gauge redundancy, and is thus considerably more complicated. Here we show that, in the case of vortex rings and for bulk modes that are much longer than the typical ring size, one can perform a systematic multipole expansion of the effective action and recast it into the simpler scalar field language. In a sense, in the presence of vortex rings the non-single valuedness of the scalar can be hidden inside the rings, and thus out of the reach of the multipole expansion. As an application of our techniques, we compute by standard effective field theory methods the sound emitted by an oscillating vortex ring.


Introduction
Vortex lines are the only allowed vorticose configurations in zero-temperature superfluids. This is because superfluidity requires the fluid flow to be irrotational, and so vorticity can only be localized on one-dimensional defects. These are curves that from the viewpoint of the long-distance hydrodynamical description have zero thickness. In practice, their thickness is microscopic, e.g. of order of atomic size for liquid helium, and the hydrodynamical description breaks down within their core.
For ordinary finite-temperature fluids, such as water, one can have much more general vorticose configurations. Still, one can set up fairly easily vortex line-like solutions. Their thickness in this case is not microscopic and the hydrodynamical description still holds at those distance scales, but as long as the thickness is much smaller than the other length JHEP02(2018)022 scales in the system -say the radius of curvature of the vortex line, or the typical wavelength of sound waves in the surrounding fluid -one can take in first approximation the zero-thickness limit, and parametrize finite-thickness corrections according to the standard philosophy of effective field theories.
In the incompressible limit, the fact that vortex lines are consistent solutions of the hydrodynamical equations for ordinary fluids is guaranteed by Kelvin's theorem. To see this, consider an initial configuration with vorticity localized on a one-dimensional defect and apply Kelvin's theorem to a very small loop that wraps around the defect and follows the fluid flow. Its circulation Γ ≡ v·d stays constant, while ∇· v = 0 implies that the loop does not grow, thus conserving vorticity in the same infinitesimal spatial (comoving) region. Therefore, one finds that the defect stays one-dimensional -i.e., it does not dissolve and spread out vorticity -and that it moves along with the surrounding fluid flow. This of course ignores the dissipative effects associated with viscosity, which are not contemplated in Kelvin's theorem and in fact make all fluid flows die eventually. However, compared to perfect-fluid hydrodynamics, viscosity is associated with higher derivative corrections, which are suppressed for long-distance and low-frequency phenomena, and thus negligible in first approximation. Notice also that there is a relativistic generalization of Kelvin's theorem (see e.g. [1]), and so the same conclusions apply to vortex lines in relativistic fluids as well.
The above argument involving Kelvin's theorem in principle can fail for compressible fluids. Indeed, a non-zero divergence ∇ · v = 0 could a priori make the loop around the defect grow, thus allowing the thickness of the defect to spread out and become large over time. 1 In appendix A we prove that a zero-thickness vortex line evolves in an ordinary (perfect) fluid exactly like it does in a superfluid: it keeps its circulation and its onedimensional nature, it does not generate other forms of vorticity, and it moves along with the surrounding fluid. This is consistent with the fact that, in the absence of more general forms of vorticity, superfluids and ordinary fluids obey the same equations of motion. Indeed, there exists a duality between vorticity-free ordinary fluids and superfluids that holds directly at the level of the action, in the relativistic case as well [1]. It then follows that, for the more realistic case where vorticity is localized on a tube of small (yet finite) thickness a, the rate at which a can grow is suppressed by some positive power of a itself, and can thus be neglected in first approximation. In the following we will therefore restrict ourselves to the superfluid case. The discussion above tells us that, to lowest order in the vortex line thickness, our results will be valid for ordinary fluids as well.
For a review of the general properties of vortex lines and vortex rings, we refer the reader to refs. [2,3]. Recently, considerable progress has been made in studying their dynamics via effective field theory techniques [4][5][6][7][8]: one can couple systematically the bulk degrees of freedom that parametrize generic fluid flows and long-wavelength perturbations such as sound waves to the embedding coordinates X(t, σ) of a zero-thickness stringthe vortex line -with σ being an arbitrary coordinate along the string. This is done at the level of a derivative expansion for a long-distance/low energy effective action, which is valid at distance scales much bigger than the vortex line's thickness.

JHEP02(2018)022
With this paper we would like to go one step further in this program, and write down an effective field theory for small vortex rings 2 interacting with the surrounding fluid, valid at distance scales much bigger than the typical ring size, organized as a multipole expansion.
The advantage of doing so is twofold: on the one hand, if indeed the surrounding fluid only has perturbations of very long wavelengths, our effective field theory replaces the infinitely many string degrees of freedom of a vortex line with a few collective coordinates for the vortex ring (such as its center's position) and a handful of multipole moments. The multipole expansion of course involves infinitely many terms; but, as usual, these are naturally organized in decreasing order of importance, and so for long distance phenomena only the first few multipoles will be relevant.
On the other hand, the multipole expansion for a vortex ring allows one to avoid a technical complication that is needed to describe a more general vortex line: because of the topological defect nature of a vortex line -the superfluid scalar phase is not single valued in the presence of a line -one cannot use a scalar description for the fluid bulk modes, but rather one has to resort to a magnetic-type dual description involving a two-form field, which is technically more tedious [5,6,9,10]. In the case of our multipole expansion for a small vortex ring, such a complication disappears. Essentially, this is because the non-single valuedness of the scalar phase can only be probed at short distances, by going through the ring. Any process that does so is however not contemplated by our multipole expanded effective theory, which can only describe long-distance phenomena. After performing the multipole expansion, we are thus able to revert to the scalar description, which is much easier to use for computations. As a concrete application of our formalism, we will compute the sound emitted by an oscillating vortex ring.
For technical convenience, we will use a manifestly relativistic notation, with metric signature (−, +, +, +), and natural units ( = c = 1). Our results are fully relativistic, but if one wishes to take the non-relativistic limit, one can do so at any stage in our computations. Since a vortex line or ring has to move very slowly anyway [5,6], the resulting modifications are minor. For instance, for the sound emitted by an oscillating vortex ring (section 7), the expression for the emission power (7.27) remains the same in the non-relativistic limit, but with the relativistic enthalpy density w replaced by the mass density mn.

Scalar description
From a QFT viewpoint, a superfluid at equilibrium can be defined as a system with a conserved U(1) charge Q in a homogeneous state |ψ such that (i) Q has a nonzero density, and (ii) Q is spontaneously broken:

JHEP02(2018)022
This is the QFT analog of the statement that the ground state is a Bose-Einstein condensate (see e.g. [11] for a recent review of superfluids), but, unlike that statement, it has the advantage of not relying on a weakly coupled microscopic model for superfluidity. At zero temperature, which is what we will assume here, the conditions above are enough to determine the low-energy dynamics of the system. 3 The condition of having a finite density is usually implemented by introducing a chemical potentialμ and demanding that the equilibrium state |ψ be the ground state of the modified Hamiltonian [13] (the eigenvalue can always be set to zero by a c-number offset of H.) Then, if Q is spontaneously broken, so must be H, i.e., the internal symmetry generated by Q and the time translations generated by H are broken down to the diagonal symmetry generated by H , which can be thought of as the appropriate time translations in the superfluid phase.
Adopting an EFT point of view, the low-energy/long-distance dynamics must be governed by the Goldstone mode associated with the above symmetry breaking pattern [12,14]. In particular, the low-energy effective action must be a functional of a real scalar field φ on which the U(1) symmetry acts as a shift symmetry 3) The effective action can thus depend on φ only through its derivatives so, to lowest order in a derivative expansion, it must read for some function p. The energy momentum tensor is then that of a perfect fluid if we identify p with the rest-frame pressure, its Legendre transform ρ ≡ µp − p with the rest-frame energy density, and with the four-velocity field.
Notice that the Noether current associated with the U(1) symmetry (2.4) is and so p is physically the rest-frame number density n. We thus have that the enthalpy density is w ≡ ρ + p = µn , (2.8)

JHEP02(2018)022
which shows that our variable µ is in fact the local value of the chemical potential -hence the symbol chosen for it. Finally, the speed of sound is In this EFT, the desired symmetry breaking pattern is realized as a nontrivial expectation value for φ at equilibrium. The requirement that φ be invariant under H but not under H and Q separately, implies which is indeed a solution of the equations of motion following from (2.4), with µ(x) =μ. It corresponds to the fluid being at rest, U µ = δ µ 0 , and with a non-zero charge density, n ≡ p (μ) = 0. The Goldstone field π(x), can then be thought of as a perturbation of this background solution φ(x) =μ t + π(x) . (2.11) Expanding the action in powers of π one gets and all 'barred' quantities are evaluated on the background solution (2.10), that is at µ =μ. From the quadratic part of this action we see that π propagates at the speed of sound, and can thus be identified with the phonon field. The interaction terms involve derivatives of π and thus become stronger and stronger at higher and higher energies. The momentum cutoff for our phonon effective field theory is of order (see [15]) corresponding to an energy cutoff of order

Two-form description
To introduce the dual two-form description, it is convenient for what follows to use Legendre transform techniques [16]. Consider thus the following action for two independent fieldsa 1-form P µ and a 2-form A µν = −A νµ : where p(µ) is the same function as above, and ε 0123 ≡ +1.

JHEP02(2018)022
A µν appears at most linearly in the action, and so it is effectively a Lagrange multiplier. Varying the action w.r.t. it yields the equation of motion ∂ [µ P ν] = 0, which is solved by 4 for some scalar field φ (the minus sign is conventional). Plugging this back into (2.16) reproduces the scalar action (2.4) up to a boundary term, so the two actions are completely equivalent descriptions of the same bulk dynamics.
On the other hand, we can proceed in the opposite order and, starting again with eq. (2.16), integrate out P µ first. Its equation of motion reads which we can interpret as an algebraic non-linear equation for P µ , to be solved in terms of F µ . In particular, we get P µ −F µ , and their squares are related by which we had already identified with the number density, i.e. the thermodynamic variable conjugate to µ. Replacing P µ with its solution in the action, we get an equivalent action which is now a functional of A µν only [6]: This describes exactly the same bulk physics as eq. (2.4), but it does so using a twoform rather than a scalar field. Notice that above we had already identified the Legendre transform of p with the energy density ρ, which is naturally a function of n.
The action (2.20) has an abelian gauge symmetry so one also needs to add a gauge fixing term, a convenient choice for our purposes being where ξ > 0 is an arbitrary parameter. Notice that, despite A µν 's self-interactions, the abelian nature of its gauge transformations implies that a gauge-fixing term like the one above is all we need to make sense of an unconstrained path integral over A µν . That is, we don't need any ghost fields, even if we were to consider quantum/loop computations. The gauge symmetry implies that there is only a single propagating degree of freedom, in agreement with the dual scalar formulation. Indeed, A 0i appears in the action with no time-derivatives and therefore does not feature propagating wave solutions (but can JHEP02(2018)022 mediate long-range Coulomb-type interactions), while ε ijk A jk has a pure-gauge transverse part, leaving the longitudinal part to carry the sound degrees of freedom. Now the background solution describing the superfluid at equilibrium, i.e. the analog of the scalar ground state (2.10), is, up to a gauge transformation, Perturbations about this background are parametrized by two three-vector fields A and B, In the language of [4,6], A is the "hydrophoton" and B is the sound field. In the following we will work in the Landau gauge ξ → 0 exclusively, which for perturbative computations is equivalent to imposing ∇ · A = ∇ × B = 0 in interaction terms and using the transverse and longitudinal propagators for A and B, respectively:

Vortex lines
Although the scalar and two-form descriptions are physically equivalent as far the bulk dynamics are concerned, in the presence of vortex lines one should really choose the twoform one. The reason is that the scalar φ is not single-valued in the presence of a vortex line, while the two-form is. For instance, for an infinitely long, straight vortex line in an otherwise unperturbed superfluid one has where ϕ is the angle around the string, and Γ is the string's circulation, Γ = v · d x. Whenever ϕ winds once around the string, ∆ϕ = 2π, φ(x) undergoes a nonzero change, ∆φ = −μΓ. On the other hand, the corresponding solution for the two-form reads (up to gauge transformations) where we are using the standard form notation, r is the distance from the string, z a cartesian coordinate along it, and for simplicity we neglected the r-dependence of µ and JHEP02(2018)022 n, which is appropriate at large distances from the string, or in fact at all distances in the non-relativistic limit. 5 Finally, r 0 is an arbitrary scale needed to compensate dimensions in the log; changing r 0 just amounts to performing a gauge transformation on A. The solution above does not depend explicitly on ϕ, and is thus perfectly single-valued around the string. Related to this, if one tries to describe the dynamics of the vortex line by coupling the embedding coordinates X µ (σ, τ ) (σ and τ being world-sheet coordinates) to the bulk degrees of freedom, one discovers that in the scalar language all interactions with φ must involve at least one derivative acting on each φ (because of shift invariance). On the other hand, the two-form language enables one to write down a non-derivative coupling between the line and A µν , which is thus more relevant at low energies than all the couplings that are allowed in the scalar field case. This is completely analogous to having a magnetic monopole in E&M: the electric gauge field A µ is not single valued, and cannot couple locally to the monopole without derivatives, that is, only higher-multipole couplings are allowed; conversely the magnetic dual fieldÃ µ is single valued and, compatibly with all the symmetries, admits a local world-line coupling of the form In our case, the symmetries are: • Poincaré invariance, under which A µν (x) and X µ (τ, σ) transform in the obvious way; • World-sheet reparametrization invariance.
The leading string-bulk interaction compatible with all these symmetries is given by a Kalb-Ramond term (λ is a generic coupling constant), which is a direct generalization of (3.3). Notice that, like (3.3), this is gauge-invariant only up to total derivatives. One can easily check that the equations of motion deriving from the action where S bulk is the bulk action (2.20) discussed above, admit a straight-string solution with A µν given by (3.1). At the next order in a derivative expansion for A µν , one can write down a Nambu-Goto-like term (" NG ") in which the tension is an arbitrary function of the two invariants JHEP02(2018)022 that can be built out of the superfluid's and string's degrees of freedom: 6 Although (3.6) is formally a higher-derivative correction to (3.4), the two action terms are comparable as far as their contributions to the string's dynamics are concerned. The reason is that if in (3.4) one sets A µν to its background value (2.23) and ignores its fluctuations, the string has degenerate dynamics.
This says that the velocity of the string projected perpendicularly to the string itself is zero. Since the parallel component of the velocity is unphysical to begin with -because of σ-reparametrization invariance -we reach the conclusion that the only source of string dynamics in (3.4) comes from the interaction of the string's degrees of freedom with A µν 's fluctuations, through diagrams like (3.8) where the plane represents the string world-sheet, the dotted line represents A and the wavy line represents B. Such diagrams are log-divergent in the UV, 7 and -in agreement with standard renormalization theory in QFT -the structure of the UV divergences is identical to what one would get from the local couplings in (3.6). As a result, these diagrams induce a running of the coefficients appearing in the function T . For instance, for the tension T ≡ T (n,ū) one gets the RG equation where k is a renormalization (momentum) scale, with solution where k 0 is a UV momentum scale, possibly of order of the inverse string thickness, but logically separate from it. In conclusion, (3.6) is subleading as a source for A µν , but it is JHEP02(2018)022 as important as (3.4) as far as the string dynamics are concerned, and we should thus keep it whenever we are interested in those. Things simplify substantially in the non-relativistic limit. In fact, for the string's motion to be within the regime of validity of our effective theory, the string has to move much slower than the speed of sound, and thus than that of light. So, we can expand our full action for slow string speeds (while allowing for a relativistic equation of state for the surrounding superfluid, i.e., for a relativistic sound speed). The leading-order terms are where we chose the gauge τ = X 0 ≡ t for the world-sheet time coordinate, and the function T implicitly depends on the renormalization scale k, which -to minimize contributions coming from virtual A µν fluctuations -can be conveniently identified with the typical momentum scale of the process under consideration. In fact, we can use (3.10) to reconstruct T (n) and its scale dependence: eq. (3.10) depends on the background value of n explicitly throughn, and implicitly throughw =ρ +p and possibly the UV scale k 0 . 8 This tells us that the generalized tension T (n) at scale k is where we used the zero-temperature thermodynamic identity ρ + p = µn. As a check, notice that this matches precisely the results of [6] for the running of T (n), which were derived directly via Feynman diagrams.

Back to the scalar?
As clear from (2.17) and (2.18) and as customary for Hodge dualities, our scalar/two-form duality is local at the level of shift/gauge invariant field strengths (∂ µ φ and F µ ) but nonlocal at the level of field potentials (φ and A µν ). This suggests that all action terms that are written directly in terms of field strengths and derivatives thereof can be processed through the duality and written in a local fashion in either language, while terms that are invariant only up to total derivatives -like our Kalb-Ramond coupling above -will look non-local if written in the wrong language (the scalar one, in this case.) However, it is instructive to try anyway and rewrite the action above in the scalar language: this will show us explicitly where the obstacle is, and how it will be overcome by the multipole expansion of the next section. Roughly speaking, when we consider a vortex ring and we Taylor expand the A µν in the Kalb-Ramond coupling about a suitably-defined center of the ring, the monopole term -that involving A µν without derivatives -vanishes. All higher multipoles involve gaugeinvariant combinations of derivatives of A µν , that is, the field strength F µ and derivatives thereof, and can thus be rewritten in the dual scalar language as local couplings. However, the situation is complicated by the fact that our duality involves a Legendre transform,

JHEP02(2018)022
which relates the form of the duality relation between the fields to the action: every time we write down a new term in the action, the duality transformation changes. As a warm-up, consider our non-relativistic action above without the Kalb-Ramon term, where we have left the x-and X-dependence of f implicit, since they will not matter in the duality transformation. This action depends on A µν only through the gauge invariant combination n = −F µ F µ . We can thus follow the simple example of section 2, reintroduce the one-form P µ , and write an equivalent action that depends on A µν at most linearly: where µ ≡ −P µ P µ , and g(µ) is related to f (n) by a Legendre transform: The logic is exactly the same as before: on the one hand, integrating out P µ from (4.4), gives us back exactly the action (4.1). On the other hand, solving first the equation of motion coming from varying A µν , and plugging back into the action gives us the scalar action where µ now is This result is formally correct, but it does not make immediately obvious what the string-localized term looks like in the scalar formulation. The reason is that the deltafunction term in f (n) (see eqs. (4.2)-(4.3)) contributes to g(µ) both directly, via g = µn−f , and indirectly, because it affects the relationship between µ and n. To understand its overall JHEP02(2018)022 effect, it is useful to treat the delta-function term as a small perturbation to the action and work perturbatively in it. We formally write where the subscript 'B' refers to bulk relationships, and so in particular the functions µ B (n) and n B (µ) are those associated with the superfluid's equation of state, while δf, δµ, δn are their string-localized perturbations. The Legendre transform of f , g(µ), thus is where we kept up to first order in the perturbations, since as we will see in section 6.1 higher order terms can be set to zero using a suitable regularization scheme. The zeroth order term is just the Legendre transform of the bulk action, using bulk relationships: it is the pressure p(µ). The combination multiplying δn(µ) vanishes, because We are left with −δf (n B (µ)). We thus see that, to this order, the scalar action is just the original bulk one, p(µ), supplemented by the same NG term appearing in the twoform formulation, where now the tension should be expressed as a function of the chemical potential µ -in particular, using the bulk equation of state to relate n, ρ, and p to µ: Note that such a simple result cannot hold if we keep the u dependence in T as well; as a consequence, away from the non-relativistic limit our duality will have to be modified. Consider now adding the KR term (3.4). We can work directly at the P -A level, since the KR term is already linear in A µν : Again, this action is equivalent to (3.11) upon integrating out P µ . Now however switching to the scalar formulation is more complicated. Varying with respect to A µν one gets the eom which has 0i and ij components

JHEP02(2018)022
where we denote ∂ σ by a prime from now on. Equation (4.20) can be solved by taking a second curl and inverting a Laplacian for the transverse part of P , yielding the standard Biot-Savart law The undetermined total gradient ∇Φ is, as we already know, the scalar field of the scalar formulation. Now plugging (4.22) into the divergence of (4.21) and using the fact that the Biot-Savart term is divergence-free, we can again invert the Laplacian to get We thus have that (2.17) is modified as and plugging this expression back into the action we get (up to a boundary term) Now, since ∇ · J = 0, equation (4.22) is the Helmholtz decomposition of P , meaning that J carries the rotational information of the fluid flow, while Φ carries the compressional part. It is thus clear that the physical interpretation of the Φ fluctuations is sound. Indeed, comparing with in the two-form formulation, we see that we have reorganized the degrees of freedom as follows: Unfortunately, due to the non-local, integral structure of J µ (x) the action (4.25) is not very useful -this was precisely the reason to introduce the two-form language in the first place. From this viewpoint, it's not clear how in the case of vortex rings the multipole expansion is going to help: it looks like one will simply get a multipole expansion of J µ , with long-range multipole tails ∼ 1/| x − X| +2 , which are still non-local. One has to keep in mind, however, that the Φ-J split of P µ is not unique, since one can transfer a fourgradient from one part to the other. This is precisely what allows one to recover a local (i.e., non-integral) expression for P µ in the multipole-expanded case. This is easier to see if one imagines performing the multipole expansion before solving eqs. (4.20), (4.21). For example, the first nontrivial term for the first equation is a dipole term,

JHEP02(2018)022
where q is a constant vector -the dipole moment -and x c is an arbitrary point in a neighborhood of the vortex ring. This equation tells us that the curl of P is localized at x c . Correspondingly, there is a solution for P such that the non-gradient part of P is localized at x c : The difference between this and the lowest-order multipole of (4.22) is that the non-gradient term now is not divergence free, that is, it's not a pure curl. So, some gradient part has been moved from Φ to J , with the end result of making P local. The same magic works for all higher multipoles, and for P 0 as well, as is clear from the results of the next section.

The multipole expansion
Consider now a vortex ring of typical size R, interacting with much longer bulk modes. One can perform a multipole expansion of the interaction terms, which, as usual, just amounts to a Taylor expansion of the bulk modes about an arbitrarily defined "center" x c of the ring: To keep the notation light, we are suppressing the (σ, t)-dependence of our variables. Since x c is a collective coordinate for the ring, it only depends on t. On the other hand, Y depends both on t and on σ. Notice that, as usual, the multipole expansion is done at fixed time, that is, it corresponds to a Taylor expansion in space only. Let's focus first on the Kalb-Ramond term in the action. Before proceeding with the multipole expansion, it is useful to separate A µν 's background from its fluctuations A and B, as in (2.24). The Kalb-Ramond term thus splits as When we decompose X as X = x c + Y , we can formally perform the σ integral inS KR exactly and end up with a simple point-particle action. To lowest nontrivial order (O(Y 2 )), we getS where q -in general a function of t -is defined as

JHEP02(2018)022
This result and the following ones are straightforward consequences of some simple identities, which follow from integration by parts: On the other hand, to perform the σ integral in δS KR we first have to Taylor-expand A and B as in (5.1), and this is what yields the multipole expansion. Working again up to O(Y 2 ) only and using the identities above, we get with the understanding that A's and B's derivatives are to be evaluated at x c . As anticipated, at least at this order, the interactions of the ring with the bulk modes organize themselves into gauge-invariant combinations. Indeed, the two-form's field strength F µ is so that, at this order, the effective point-particle action for the vortex ring coupled to bulk modes can be written as In appendix B we generalize this result to all orders in the multipole expansion. We get where are the ring's multipole moments, characterizing its size, shape, and orientation. Eq. (5.10) just corresponds to the m = 0, O(Y 2 ) terms in the sum. Notice that so far we have been deliberately evasive about how to actually define the ring's center position x c : none of the above results depends on the precise definition of x c . As usual for a multipole expansion, a redefinition of the origin about which one expands keeps the multipole expansion intact, yielding only a reshuffling of the multipole moments, with lower order ones contaminating the higher order ones. For us, a possible definition of x c is that of a center of mass-type coordinate, 14)

JHEP02(2018)022
which has the advantage of being purely geometrical, that is, reparametrization invariant. However, this is totally arbitrary, and in certain cases one might find it convenient to use other choices. The multipole expansion of the non-relativistic generalized Nambu-Goto term is at the same time more straightforward and more tedious, since in that term A µν couples to the string directly through the gauge invariant combination n = −F µ F µ , but it does so in a non-linear fashion. The Taylor expansion of the generalized tension T (n(X)) about x c reads where the complexity of the higher order terms quickly escalates, due to the aforementioned non-linear structure of the coupling. In fact, for what follows it is probably better to keep the derivatives in an unexpanded form, and expand them only when needed for perturbative computations. We can thus define new multipole moments, so that the full multipole expansion of (3.11) reads Notice that with the choice (5.14) for x c , the q -dipole, vanishes.
The expression above includes all multipoles. In the limit of very small rings, or, equivalently, very long bulk modes, the dynamics will be dominated by the lowest multipole terms: where F µ and n are evaluated at x c , and we have defined the ring's typical size R by

JHEP02(2018)022
6 Back to the scalar! Now that everything is expressed in terms of F µ , n = −F µ F µ , and derivatives thereof, we can pass to the scalar formulation straightforwardly. We first rewrite the string's couplings as The action (4.18) now reads with µ = √ −P 2 . The equation of motion associated with A is, in exterior calculus notation, d (P + j) = 0, whose solution is for some new arbitrary scalar field φ. Plugging this back into the action we get Let us now compare this result to (4.25): as anticipated, j µ is not the same as J µ . For instance, J is divergence-free but j is not. As a consequence, φ is not the same as Φ. The advantage of the φ-j decomposition of P over the Φ-J one, is that j µ is local, since it is a series of terms each supported at x c only. The same holds for j of course. In fact, even when the full j and j infinite sums are taken into account, one gets mild non-localities that extend only to distances of order R, being associated with having an extended source, i.e., with the fact that the exact location of the string is somewhere around x c , not precisely at it. Put another way, the magic at work here does not really have to do with the multipole expansion, but rather with having a closed string: only in that case can one rewrite its couplings purely in terms of the gauge-invariant field strength F µ .

UV divergences, new and old
An apparently problematic feature of (6.6) is that it is singular: the sources j and j are series of Dirac-deltas and derivatives thereof. Since the action depends non-linearly on j, any perturbative expansion in j will yield products of deltas and derivatives thereof, all localized at the same point, which are notoriously singular, even as distributions.
For definiteness, consider truncating the j µ and j sums to the lowest order (m = 0) terms, which corresponds to the approximation (5.20) for the effective point-particle action:

JHEP02(2018)022
In this limit, the scalar action (6.6) reduces to The first line of the final expression encodes the familiar bulk dynamics of φ (zeroth order in j, j ) as well the world-line lowest-order multipole interactions between the vortex ring and φ (first order in j, j ). The second line contains the divergences we alluded to: at second order in j, j one has a δ 3 (0); at higher and higher orders one gets higher and higher powers of δ 3 (0). Had one included higher multipoles (higher m) in j and j , one would have found derivatives of δ 3 (0) as well.
What are we to do with these divergences? From a field theory viewpoint, UV divergences are nothing new. For instance, in the case of a long string, self-energy diagrams in which the string exchanges bulk modes with itself are logarithmically divergent [6]. Such UV divergences are harmless -they can be renormalized away -but physical: being logarithmic, they signal that certain local couplings such as the string tension will "run" logarithmically with scale. When we take the limit of a small vortex ring and we treat it as a point-like source, we are effectively going from an effective theory with the string thickness a as UV cutoff, to an effective theory with the ring size R as UV cutoff. As a consequence, those logs get saturated at log R/a, and now the same self-energy diagrams have power-law UV divergences, cut off at distances of order R (more on this below). Power-law divergences are both harmless and un-physical: after renormalization, there is no low-energy remnant of power-law divergences. In fact, if one uses dimensional regularization they don't even show up.
Since δ 3 (0) and its derivatives are pure power-law divergences, we can decide once and for all to use dimensional regularization, so that we can just set all these divergent terms to zero. In that case, the scalar Lagrangian -now including all higher multipoles -stops at first order in j and j , and reduces to 9) where J µ is the conserved U(1) current,

JHEP02(2018)022
This is certainly the most economic approach, and in fact in the two-form case dimreg has the additional advantage of respecting gauge-invariance, unlike for instance a hard momentum cutoff. However, it is also interesting to wonder whether in going from the two-form language to the scalar one, the UV divergences just get shuffled around: some of the UV divergences that one would get from the self-energy diagrams in the two-form language, appear directly as divergent terms in the scalar Lagrangian. We checked that, to lowest order in the multipole expansion and to order qq and q 2 , this is indeed the case.
As alluded to above, another aspect associated with UV divergences is the RG running of certain localized couplings between the ring and bulk modes. For a long string, one gets log divergences in self-energy diagrams that induce a running of the generalized tension as in (3.12). Now, when we consider a closed string of typical size R interacting with much longer bulk modes, in going to lower and lower momenta the running stops at k ∼ 1/R. Related to this, in the point-particle limit for the ring the self-energy diagrams that used to be log-divergent for a long string are now power-law divergent, and thus have no running associated with them. (There are certainly higher-order diagrams that are log-divergent, and thus induce RG running for higher-order coefficients -see e.g. [17,18].) In the pointparticle limit we can thus take expression (3.12) for T and set k ∼ 1/R in the log. As for the UV scale k 0 , strictly speaking one should determine that and its µ dependence from experiment. However, up to order-one factors one can generically expect k 0 ∼ 1/a -the inverse string thickness. In the limit R a one is dominated by the large log R/a, and any corrections due to the mismatch between k 0 and 1/a can be ignored in first approximation, and so can their dependence on µ. One thus gets n(µ) µ log R/a , (6.11) Finally, it is worth pointing out that if the vortex ring has perturbations (Kelvin waves) of typical wavelength much shorter than R, then is the scale that should enter the log, not R.
To summarize: eq. (6.9) is the effective action one should use to describe the interactions between a small vortex ring and much longer bulk modes. The ring's multipole moments are defined in (5.12), (5.13), and (5.17). For perturbative computations, one should expand φ =μt + π and µ = −(∂φ) 2 in powers of π. If the vortex ring's typical size R is much bigger that the string's core size a, one can safely replace T (µ) by eq. (6.11), in which case its µ-(and thus π-) dependence only comes from the prefactor n(µ)/µ. For smaller vortex rings the full expression (3.12) should be used.

Sound emission by an oscillating vortex ring
As a simple application of our formalism, let us now compute the rate at which an oscillating vortex ring emits sound. To lowest order in λ (and thus v/c s , with v the typical speed of the ring -more on this below) and to zeroth order in bulk self-interactions, the corresponding JHEP02(2018)022 diagram is simply (7.1) where the wavy line is sound (π(x)), and the double bold line is the world-line of the vortex ring, to be treated as an external source in first approximation. The relevant part of the action thus comes from expanding (6.9) in powers of π, to second order in the bulk and to first order on the vortex ring's world-line: where the subscript (1) means that we take the part that is linear in the phonon field π. We have kept the first multipoles only, which, as usual, dominate emission at low frequencies.
Recall that where µ = −(∂φ) 2 , φ =μ(t+π), is the typical wavelength of the ring's oscillations, and we used the approximate expression (6.11) for the tension. From these we can compute the relevant terms appearing in (7.2): where we used that λ is related to the vortex line's circulation Γ by λ = (w/n)Γ [6].

Power counting
Our multipole expanded action contains several scales controlling the dynamics of the bulk modes and of the string and their mutual interactions. To understand which term in the action dominates phonon emission, we have to estimate the orders of magnitude of the various multipoles. First, notice that for the string to be captured by the EFT we clearly need its typical size to obey L 1/k * , where k * is the strong coupling momentum scale of the bulk interactions (2.14). Another relevant scale is the core radius of the vortex a L and, for the string approximation to be valid, we need the momenta to be k X , k π 1/a. However, because of the multipole expansion ∼ ∞ n=0 (k π L) n , (7.5) we actually have a stronger bound on the bulk momentum, namely k π L −1 . Next, the parameter λ, which controls the bulk-string interactions, has dimensions of action and is

JHEP02(2018)022
typically of the order of in the superfluid case, so O(1) in natural units. This means than rather then doing perturbation theory in λ, we should expand in some other quantityin our case this will be related to the multipole expansion. Let's then go back to (6.9), and consider the linear order in π, but all orders in the multipole expansion. For our estimates we will write ∇π ∼ k π π,π ∼ ω π π =c s k π π, and so on. The parameters related to the string dynamics are (up to logs) where v is the string's typical speed, and ω X and k X are the frequency and wavenumber of the string's oscillations, i.e. the Kelvin waves (see appendix C). For definiteness let's specialize to the lowest Kelvin wave modes, with k X ∼ 1/L. Notice that, as usual for waves emitted by an oscillating source, the frequency of the wave is the same as the frequency of the source, ω π = ω X ≡ ω. This gives a relationship between v and k π : v ∼c s (k π L) , (7.7) which means that, as usual for a multipole expansion, v/c s and k π L are not independent expansion parameters -they are in fact of the same order. We will choose to parametrize the multipoles in powers of v/c s , which is more readily measurable. We next estimate the multipole moments as follows: where A is the Kelvin waves' amplitude (in units of length). Note that we only keep terms that are linear in A, which is enough, since the terms independent of A don't depend on time and thus cannot emit, while higher powers of A are necessarily suppressed because of the bound A < L. Collecting everything we obtain 9 The explicit factors of v/c s 1 upfront tell us the relative importance of the various multipoles: the m-th q 0 and q multipoles are equally important, and as important as the (m + 1)-st q multipole.

Emitted power
Now that we have understood the hierarchy between the various multipoles, we can compute the emitted power at low frequencies.
According to our estimates above, the most important multipole is q, of order v/c s . However, to zeroth order in external perturbations (and in the emitted sound), q must be constant in time. The reason is that q is, up to constants, the conserved momentum for a free ring, that is, the Noether charge associated with translations of x c . This is clear from the results of section 5, where in the absence of external fields the only term in the action for a ring that depends on˙ x c is (5.5). This is true to all orders in the perturbations of the ring (Kelvin waves), but to lowest order in gradients of the bulk modes (phonons). Thus, to the order we are working q i does not oscillate and cannot emit sound. We thus have to look at the q 0 , q , q i,j terms, which are all of order (v/c s ) 2 . Notice however that, to first order in the amplitude A of oscillations, q 0 and q cannot depend on time either. This is clear if we expand Y as Y (σ, t) = Y 0 (σ) + ξ(σ, t), where Y 0 describes an unperturbed circular ring, and ξ its perturbations; we have where ϕ is the angle around the ring andn the normal to the ring's plane. The integrals over ϕ vanish for Kelvin waves on a ring, simply because they are integrals of sines and cosines (see appendix C). Furthermore, the˙ x c · q piece in q 0 is constant, because q is constant to all orders, while x c (t) is not affected by Kelvin waves to linear order in their amplitude (appendix C). Thus, if we restrict to oscillations with small amplitudes for our ring, we are left with one interaction term that dominates emission at low frequencies: (7.14) where we used eq. (7.4) for J i (1) . We can now apply standard QFT techniques to compute classical observables, as discussed for instance in the appendices of [4,6]. The amplitude for diagram (7.1) reads (recall that π is not canonically normalized -see eq. (7.2)) where p is the phonon's momentum, ω p =c s | p | its energy, and Q i,j 's is the "world-line Fourier transform" of q i,j :

JHEP02(2018)022
However, notice that the p · x c (t) piece in the exponent is generically of order pvt, and so to lowest order in v/c s we may just use the usual Fourier transform: For instance, for a vortex ring moving at constant velocity v, we would have Following standard QFT techniques, we relate the amplitude above to the emission probability: This counts the number of phonons emitted in a given infinitesimal d 3 p. To get something with a smooth classical limit, we multiply by the energy ω p of each phonon, thus getting the total energy emitted per unit d 3 p: Integrating over the radial component p ∝ ω p first and using standard properties of Fourier transforms, we get the energy per unit solid angle as an integral over time with instantaneous emission power per unit solid angle To integrate over the solid angle we use The total emitted power thus is where we used that, by construction, q i,j is traceless.
We can now consider specifically the case of small-amplitude Kelvin waves on our ring, whose general properties are worked out in appendix C. Using those results and the definition (5.12) we can compute q i,j to first order in the waves' amplitudes. A somewhat tedious computation yields 10 24) 10 The fact that only the m = ±2 modes contribute can be understood from angular momentum considerations: a circular ring preserves rotations about its axis, and its deformations can be classified as eigenfunctions of the corresponding angular momentum; the label m measures precisely their eigenvalues.
A two-index tensor like q i,j can contain up to spin 2; on the other hand, Kelvin waves start at |m| = 2. So, at linear order, q i,j can only depend on the m = ±2 modes.

JHEP02(2018)022
whereε i is the helicity-one polarization vector, which for a ring orthogonal to the z-direction reduces toε = 1 √ 2 (1, i, 0) , (7.25) and the ψ's are oscillating functions defined in (C.10), both with frequency Plugging all this in (7.23), and neglecting terms that average to zero over a period, we finally arrive at This is the power emitted in sound waves. We can compare it to the oscillation energy of the ring times the frequency. Their ratio will provide a measure of the decay rate of the m = ±2 modes, in units of their frequency. The oscillation energy is given by (C.18), so that we have Apart from logs and numerical factors, in the absence of other sources of dissipation the typical lifetime of the m = ±2 modes is thus of order (c s /v) 5 1 periods.

Discussion and outlook
We have developed an effective field theory for small vortex rings interacting with long wavelength fluid flows and sound waves, which we organized as a multipole expansion. Compared to more general configurations involving vortex lines, one of the advantages that arise in this case is that the fluid bulk dynamics can be described in terms of a single scalar field rather than a two-form, with obvious simplifications for concrete computations.
As an application of the formalism, we computed the sound emitted by an oscillating vortex ring via standard QFT techniques (i.e., Feynman diagrams). The same effective theory can be used to study other phenomena as well, such as the long-distance interactions between vortex rings mediated by hydrodynamical modes, or the dragging of vortex rings by the surrounding fluid flow beyond the point-particle limit. We must also stress that the dynamical variables describing all degrees of freedom of the vortex ring in the final action (6.9) are still the world-sheet fields X(t, σ), and these enter in a somewhat cumbersome way through each multipole moment q i 1 ...im µ and q i 1 ...im . Thus, our effective action is mostly useful in cases of given vorticose sources, where the time-dependence of the multipoles is known from the start, as we have illustrated with the oscillating vortex ring for instance. To turn the action (6.9) into an EFT for both sound and the vortex ring's internal degrees of freedom, which can be used to conveniently JHEP02(2018)022 compute back-reaction effects for instance, one must therefore replace the X(t, σ) description with infinitely many degrees of freedom ψ n (t) that "live" on our point-particle and parametrize its internal configuration. This is precisely what we did in appendix C, where we expressed X(t, σ) in terms of the vortex ring's position (x c (t)), radius (R(t)), and orientation (n(t)), and infinitely many oscillators (ψ m (t)) associated to the Kelvin waves (see eqs. (C.2) and (C.19)). We have worked out the free part of this action, finding in particular that the Kelvin wave spectrum on a ring exhibits some peculiarities compared to the straight vortex line case. This therefore amounts to a self-contained effective field theory that determines the dynamics of both the bulk modes and the effective degrees of freedom of the vortex. We plan to explore this theory more systematically in a forthcoming publication. In particular, in the language of the coset construction, our oscillators ψ m (t) correspond to "matter fields", whose interactions, among themselves and with the "Goldstones" x c (t) andn(t), are determined by symmetry considerations. A first step in this direction has been taken in [19]. In fact, along the lines of that paper, we can also include the gravitational field to study how vortex rings and their excitations respond to gravity.
Finally, with suitable modifications, the same formalism that we have been using here can also be applied to relativistic closed strings in empty space. It would be interesting to understand whether our effective-theory parametrization of string dynamics offers a useful formalism for perturbative (closed) string theory.

Acknowledgments
We are especially indebted to Kate Eckerle and Bart Horn for collaboration in the early stages of this project. We also wish to thank R. Penco and D. Roest for helpful and interesting discussions, and B. Horn for insightful comments on the manuscript. This

A Vortex lines in perfect fluids
Consider a non-relativistic perfect fluid. Its field velocity obeys the Euler equation, Taking the curl of this equation we get the standard evolution equation for the vorticity where we have assumed a barotropic equation of state p = p(ρ), so that the pressuredependent term drops. We now consider an ansatz for the vorticity field localized on some JHEP02(2018)022 vortex line X(σ, t) with constant circulation Γ, It is straightforward to check that, at fixed t, such a vorticity field is divergence-free -as befits the curl of the velocity field -and that the associated velocity field has circulation Γ around the vortex line. We now want to prove that such an ansatz is consistent with the evolution equation (A.2). Plugging our ansatz into the various terms of (A.2), and dropping an overall common factor of Γ from now on, we geṫ where all X's and their derivatives are evaluated at (σ, t).
Our claim now is that the vortex line is comoving with the fluid flow. If true, this means that there is a parameterization of the string (i.e., a gauge choice for the variable σ), such that for all σ and t˙ X = v( X, t) . (A.8) From this, we have˙ Plugging both of these expressions into (A.4), and using the distributional identity We thus reach the conclusion that a zero-thickness, constant circulation vortex line, comoving with the fluid, is a solution of the perfect fluid equations of motion.

B All-orders multipole expansion
In this appendix we provide the detailed derivation of the multipole expansion of the KR action (5.11). So let us consider a localized string and work in the physical gauge X 0 = τ . Defining

JHEP02(2018)022
we find that the Kalb-Ramond term (3.4) reads and it will be convenient to also define We begin by Fourier transforming the space-dependence of the bulk fields, expressing everything in terms of Y and expanding around Y = 0: where in the last step we have simply redefined the indexation of the sum for the A terms since the n = 0 term is zero. Using we have that in the contraction one can subtract from ik i 1 . . . ik i n+1 A i n+2 the totally symmetric part and be left with

JHEP02(2018)022
so that this only depends on the curl ∇ × A. Then, thanks to the total symmetry of Y i 1 . . . Y i n+1 , the contraction simplifies considerably For the B term we first need where in the last term we recognize the expression we started with. We therefore isolate it to get Again, this simplifies a lot once contracted by virtue of the symmetry of k i 1 . . . k in . . ik in ik i 1 B ·˙ q i 1 ...in + n (i k · B) q i 2 ...in , (B.12) where q i 1 ...in is the first term of (5.13). Going now back to the action (B.5), we integrate by parts the time-derivative acting on q and redefine the summation of the ∼ q term to obtain

C Kelvin waves on a ring
Consider an approximately circular vortex ring of radius R, moving in an unperturbed fluid. Its action is where k is the inverse of a typical length scale. For simplicity, we consider a situation in which the ring has deformations of some typical wavelength , rather than a spectrum of deformations that spans many orders of magnitude in wavelengths. In the former case, we can just identify k with 1/ . In the latter, we wouldn't be able to use directly the running tension as a shortcut to take into account the effect of virtual A fields; we should instead compute anew the self-energy diagrams where the ring exchanges A with itself. With these qualifications in mind, from now on we will remove the argument of T , and denote T (1/ ) simply by T . For any given time t, we can decompose the degrees of freedom X(σ, t) as X(σ, t) = x 0 + v 0 tn + Rρ(σ) + ξ(t, σ) . (C. 2) The first term is an initial condition for the center of the ring. The second term describes the uniform motion of the center of a perfectly circular vortex ring, moving along somen with velocity v 0 = T λnR . (C. 3) The third term describes the circular ring itself, withρ(σ) being the radial unit vector associated with any given value of σ. Finally, the fourth term describes deviations from circularity: these will evolve in time in a fashion similar to Kelvin waves on a straight vortex line, but with corrections due to the curvature of the ring. Plugging this decomposition into the action, and focusing on the part of the action that depends on ξ, to quadratic order we get whereφ(σ) is the unit angular vector (n,ρ, andφ thus form the standard cylindricalcoordinate basis of unit vectors.) Looking at eq. (C.2), we see that σ-reparametrization invariance-X → X + f (σ, t)∂ σ X -acts on ξ as ξ → ξ + Rφ f (σ, t) , (C.5)

JHEP02(2018)022
Notice also the peculiarities of the dynamics of the low-lying modes, m = 0 and m = ±1: their eigenfreqencies vanish. To understand what these modes correspond to, it is better to go back to the α, β variables, since the ψ m linear combinations are singular for m = 0, ±1. For m = 0, α 0 and β 0 are real and their action reduces to S 0 = (2π) dt 1 2 nλR (β 0α0 − α 0β0 ) + T R β 2 0 , (C.14) with equations of motionβ The most general solution is β 0 = const and α 0 = const − T nλR 2 β 0 · t. Looking at the geometric interpretation of α and β in (C.6), it's clear that this zero-frequency solution corresponds to keeping the ring circular while changing its radius, its normal position, and its speed in a way consistent with the change in radius.
For the m = ±1 modes, the action is with equations of motionα The most general solution now is α 1 = const and β 1 = const − T nλR 2 α 1 · t. Going back to (C.6) again, we see that this corresponds to keeping the ring circular while shifting it parallel to itself and tilting it by an infinitesimal angle δθ = α 1 /R.
To summarize: for an approximately circular ring the m = 0, ±1 modes should be thought of as small changes in the collective degrees of freedom in (C.2), that is x 0 , R, andn. The oscillating degrees of freedom (Kelvin waves) are instead associated with the |m| ≥ 2 modes, and are described by the action (C.11). Their oscillation energy is just given by the corresponding Hamiltonian, Notice that only in the m → ∞ limit do we have ζ · ζ = ζ * · ζ * → 0 and ζ · ζ * → 1, as appropriate for circularly polarized waves. For the low-lying modes instead we have ζ · ζ = ζ * · ζ * ∼ ζ · ζ * ∼ 1.