## Thermodynamics of random walking vortex loops in counterflowing superfluids

Sergey K. Nemirovskii Institute of Thermophysics, Lavrentyev ave, 1, 630090, Novosibirsk, Russia and Novosibirsk State University, Novosibirsk (Dated: August 2, 2024)

Based on the theory of the thermodynamic equilibrium in a system of quantum vortices in superfluids in the presence of a counterflow, the influence of a vortex tangle on various thermodynamic phenomena in quantum liquids is studied. Using the early calculated partition function we study some of the properties of He II related to counterflow, such as the distribution of vortex loops in their length, the suppression of the superfluid density ρ s and the shift T λ . Good agreement with the early obtained results is a fairly strong argument in favor of the point of view that the gas of string-like topological excitations can indeed be considered as a additional kind of quasi-particles having the inner structure at high temperatures, especially near the phase transition. The application of the developed formalism to the theory of quantum turbulence is briefly discussed.

## I. INTRODUCTION AND SCIENTIFIC BACKGROUND

## II. GIBBS DISTRIBUTION AND PARTITION FUNCTION

Long ago, in pioneering papers by Onsager [1] and by Feynman [2] the idea that the quantum vortices might be responsible for the superfluid phase transition had been suggested. Since then, many authors developed this idea by studying the gas of string like topological excitations in various media. Of course, the main achievements were in the two-dimensional case, where the BerezinskiiKosterlitz-Thouless transition was predicted and thoroughly investigated [3] As for 3D case, the situation here is far from a more less developed theory. Nevertheless, there are quite a few works that use various approximations [4],[5],[6], [7], [8],[9],[10]). In these papers, the effect of an ensemble of vortex filaments on the phase transition, on the behavior of the superfluid component, on the behavior of the heat capacity, and other issues were discussed. It can be stated [8],[7],[6] that the theory of string-like structures can provide a physical understanding of phenomena in superfluid liquids, superconductors, models of cosmic strings, crystalline bodies (melting), pulsars (glitches).

In the letter we apply this ideas to the counterflowing superfluids. There are many works in which various properties of superfluidity in the presence of counterflow are investigated. It is of interest to consider these phenomena from the standpoint of an ensemble of quantum vortices. In the paper we will focus on the effects related to counterflow. In particular we will study the distribution of vortex loops in the space of their lengths, the vortex line density (VLD) L total length per unit volume, the Lamb impulse and the superfluid flow induced by vortices. We also investigate suppression of superfluid density ρ s and shift T λ (which we associate with Hagedorn temperature T H ) due to the counterflow velocity. Additionally, we discuss possible applications of thermodynamically equilibrium vortices to the theory of quantum turbulence (QT).

In paper by the author [11] it was shown that a vortex filaments in superfluids under the action of a random Langevin force in the presence of a counterflow with a relative velocity v ns come into a state of thermodynamic equilibrium with the Gibbs distribution

$$\mathcal { P } ( \{ \mathbf s ( \xi ) \}, t ) = \mathcal { N } \exp ( - \frac { H \{ \mathbf s \} } { k _ { B } T } ). \quad \quad ( 1 )$$

Here s ( ξ,t ) are the radius vectors of the vortex line elements, ξ is a label parameter, in this case it coincides with the arc length, and, accordingly, runs from 0 to the length of the loop l ; Quantity N is a normalization factor. The Hamiltonian in the presence of relative velocity has the form

$$H \{ \mathbf s \} = & E \{ \mathbf s \} - \mathbf P \cdot \mathbf v _ { n s }. \quad \quad ( 2 )$$

The energy E { } s and the Lamb impulse P are defined as (see, for example, [12])

$$E ( \mathbf s ) = \frac { \rho _ { s } \kappa ^ { 2 } } { 8 \pi } \int _ { 0 } ^ { l } \int _ { 0 } ^ { l } \frac { \mathbf s ^ { \prime } ( \xi ) \mathbf s ^ { \prime } ( \xi ^ { \prime } ) } { | \mathbf s ( \xi ) - \mathbf s ( \xi ^ { \prime } ) | } d \xi d \xi ^ { \prime }, \quad ( 3 )$$

$$\mathbf P = & \frac { \rho _ { s } \kappa } { 2 } \int _ { 0 } ^ { l } \mathbf s ( \xi ) \times \mathbf s ^ { \prime } ( \xi ) \ d \xi, \quad \quad ( 4 )$$

where κ is the quantum of circulation, and ρ s is superfluid density. Relations (1) - (3) should be used to calculate the partition function (PF), and, accordingly, to determine the various properties of the vortex tangle (VT). In the ordinary statistical mechanics when P is the true momentum of a particle (or quasiparticle), Eqs. (1-2) are obvious and follows from the Galilean transformation. Unlike this, the situation for vortex loops is more complicated, since the Lamb impulse is not the 'real' momentum, and Eqs. (1) -(4) need the verification, which had been accomplished in [11].In the next work

by the author devoted to this topic [13] the one loop PF Z 1 , corresponding the stated problem (1) - (4) was calculated. In this calculation the vortex lines were supposed to be closed loops having a random walking structure on the cubic lattice with edge equal to a and obeying the Wiener distribution. The starting expression for Z 1 reads:

$$Z _ { 1 } = \int d \mathbf r \sum _ { l } \left ( \frac { a } { l } \right ) \int _ { s ( 0 ) = \mathbf r } ^ { s ( l ) = \mathbf r } D s ( \xi \} \exp \left [ - \frac { 3 } { 2 a } \oint _ { ( } s ^ { \prime } ( \xi ) ) ^ { 2 } d \xi + \ln 5 \frac { l } { a } \right ] \\ \times \exp \left [ - \beta \frac { \rho _ { s } \kappa ^ { 2 } } { 8 \pi } \int _ { 0 } ^ { l } \int _ { 0 } ^ { l } \frac { s ^ { \prime } ( \xi _ { 1 } ) \cdot s ^ { \prime } ( \xi _ { 2 } ) } { | s ( \xi _ { 1 } ) - s ( \xi _ { 2 } ) | } d \xi _ { 1 } d \xi _ { 2 } + \beta \nu _ { n s } \frac { \rho _ { s } \kappa } { 2 } \int _ { 0 } ^ { l } s ( \xi ) \times s ^ { \prime } ( \xi ) \, d \xi \right ].$$

The first line is the contribution into PF due to configurations of vortex loops. The path integral take into account all curves starting and ending in the same point r . Term ( l/a ) ln 5 appears due to fact that l/a links, placed on the cubic lattice, have (2 D -1) options for which the string can move at each step, so the number of line states N l, a ( ) increases with its length as N l, a ( ) = (2 D -1) l/a ( D is the dimension of space). The quantity ln 5 can be considered as the entropy (normalized by the Boltzman constant k B ) per link. It is interesting that for polymer chains the entropy per leg, obtained in numerical simulations (see [10]) is about 1 6, which is very close to ln 5. . The quantity a is the length parameter of model. For the cubic lattice model, it is the cube edge; for polymer chains it is an elementary step. In the case of quantum vortices, quantity a is the coherence length (see for details, e.g., [4],[5], [7]). Further the integration over r will be omitted, since the spacial homogeneity is assumed and, accordingly, Z 1 is considered as it is taken per unit volume.

In work [13] the calculation of the PF was implemented on the base of the local induction approximation (LIA) for the energy E ( s ). In LIA, the energy of the vortex loop is E loc = ε V l , where ε V is the energy of a vortex line per unit length (tension, see, e.g. [14]). The main achievement of the LIA is that the path integral in 5 becomes Gaussian and can be calculated exactly. Realizing that a complete description should be carried out on the basis of a complete expression for the energy E (3), we would like note, however, that the LIA is capable of describing many important features of the vortex ensemble. It can be motivated by that the main part of the dynamics is described precisely by the LIA (See, e.g. review article [12]). And, more importantly, since the main goal of our investigations is to study the effects related to the counterflow, the LIA is convenient option for this purpose.

The further calculations are close to the problem of the motion of a charged particle in a constant magnetic field [4],[5], although there are some distinctions and a bit different methods have been used. Result of according calculations (when v ns is directed along z ) is that the following expression for the one loop PF (see [13])

$$\int _ { \substack { t o a c - \\ \text{ind} } } \int _ { \substack { Z _ { 1 } = \int \frac { d l } { l } \left ( \frac { 3 } { 2 \pi l a } \right ) ^ { 3 / 2 } \left [ \frac { \sin ( l \kappa \beta a \rho _ { s } v _ { n s } / 6 ) } { l \kappa \beta a \rho _ { s } v _ { n s } / 6 } \right ] ^ { - 1 } \exp \left [ - \beta \tilde { \varepsilon } l \right ] \\ \mathfrak { h } \text{ the } \quad \text{in} \, \quad \text{in} \, \quad \text{and} \, \quad \text{comment} \, \quad \tilde { \ } s = \tilde { \ } w \cdots \quad \ln 5 \, \underbrace { T h }$$

Here we have introduced quantity ˜ = ( ε ε V -ln 5 βa ). The striking difference from the Kleinert's result is that we have an ordinary sine instead of a hyperbolic one. Formally this difference arises from the imaginary unit in the (Euclidean) path integral in the expression for the motion of a particle in a magnetic field. Physically the sine is responsible for instability of the flow with respect to formation of large loops. Indeed, the combination E s { } -P · v ns is just an analog of the Landau criterion (but not for quasiparticles, but for large vortex loops).

Developing the factor in square brackets into a power series in the variable v ns (up to second term, the first one is obviously equal to unity), we split the one loop PF Z 1 (6) in two parts - the velocity-free PF Z 1 free and the velocity dependent part Z 1 V ns .

The quantity Z 1 ts the PF for a single closed string of arbitrary length. Consistently assuming that in LIA vortex loops do not interact, we calculate the grand canonical PF Z for an ensemble of many loops using the following simple formula

$$Z = \sum _ { N = 0 } ^ { \infty } \frac { Z _ { 1 } ^ { N } } { N! } = \exp \left [ Z _ { 1 } \right ]. \text{ \quad \ \ } \text{(7)}$$

The quantities Z Z , 1 , Z 1 free and Z 1 V ns are the key ones for calculating various properties of superfluids due to the presence of a VT. Let's discuss a series of issues.

## III. VARIOUS APPLICATIONS

Two populations Loosely speaking, there are two ingredients that affect PF. The first one is associated with thermal activation. The dynamics of such vortices, which we will call thermodynamic vortices, is determined mainly by the exponential (Gibbs) factor in (6), and

their sizes are concentrated in the vicinity of the correlation length a. The role of the velocity enclosed in square brackets in (6), is reduced to small corrections to the structure of the VT and to the thermodynamic properties of the system.

There exists one more population of vortex loops, which we call hydrodynamic vortices. Their appearance is quite similar to nucleation of free vortices from bound pairs in the 2D case [15]. From the form of the Hamiltonian H { } s and from the saddle-point approximation, it follows that this population is concentrated near the scales where the energy E { s } competes with the Lamb impulse P multiplied by v ns (see Eq. (1)-(3)). It is seen that the most preferred scale is of the order κ/v ns . This result is close to what takes place in the theory of quantum turbulence. A full consideration of such a situation seems very complicated and is beyond the scope of this article. In this article, we will focus only on thermodynamic vortices and their influence on the properties of superfluids. Questions related to hydrodynamic vortices, as well as other questions relating to the relationship of equilibrium quantum vortices with quantum turbulence (QT), are of great interest and will be investigated in the future.

The Hagedorn temperature It can be seen that for ˜ ε &lt; 0. the PF Z 1 (6) diverges. The corresponding limiting temperature T H , when ˜( ε T H ) = 0 is called the Hagedorn temperature [4],[5], [7], [16] and is defined as

$$T _ { H } = \frac { \varepsilon _ { V } a } { k _ { B } \ln 5 }. \quad \quad \quad ( 8 )$$

In paper [16] it was shown that the T H is practically equal to the phase transition temperature ( T λ in the case of He II). A small difference arises due to the finiteness of the system and disappears at large volumes. The physical meaning of T H is that the configurational entropy (multiplied by temperature) is comparable to the energy of a vortex filament link. Therefore, the energy entering the system is spent on the increasing the length of the lines, and the phase transition is accompanied by the appearance of very long filaments.

Vortex line density. A velocity-free state v ns = 0 . A velocity-free state describes thermodynamically equilibrium vortices in helium at rest. Previously, this problem was studied under various assumptions. The most popular were studies in which vortex filaments were considered as ideal rings (see, for example, [8],[9],[10]). Another approach used in our article, when vortex filaments were either flexible polymers with an effective link, or lattice models were used [4],[5], [7], [17]. One more direction, was the study a statistics of topological defects of the underlying fields (see, e.g., [16]). In these papers the authors predicted that VLD L at the Hagedorn temperature T H has a value of the order of the inverse squared coherence length, that is, L ˜ (1 /a ) 2 . Then L is exponentially suppressed below T H . All these facts are immediately follows from Eq. (6) for PF (if we put v ns = 0).

FIG. 1: (Color online) The loop length distribution for the Hagedorn temperature T H and for T &lt; T H , see text, In this figure, we do not take into account the divergence of the correlation length near the lambda point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00680v1-figure-b28e2af8-f677-495b-a004-e1afeba4c813.jpg)

Indeed, the integrand is the production of the power and exponential ingredients. Therefore, at the Hagedorn temperature T H , when ˜ = 0 the loops distribution becomes ε a pure power function n l ( ) = ( 1 l ) 5 / 2 ( 3 2 πa ) 3 / 2 , which indicates the random walking nature of thevortex line structure. For the temperatures T &lt; T H , when ˜ ε &gt; 0, the Gibbs factor prevails and long loops are exponentially suppressed.

Thus, we can state, that our results in the velocityfree state are in good agreement with those previously obtained.

̸

Vortex line density due to counterflow L ( v ns ) Let's come to effects related to presence of the counterflow, v ns = 0. We begin with the correction for the VLD L ( v ns ). It can be evaluated from the velocity dependent part Z 1 V ns of PF Z 1 (6). The most important result is that the VLD is proportional v ns squared, i.e. L therm ( v ns )= γ 2 therm v 2 ns . It is noteworthy that this result is also valid while using the sine (in Eq. 6) without a power series developing. Quantity γ therm , calculated from Z 1 V ns , as a function of the temperature, is shown in the inset of Fig. 2.

It is remarkable fact that γ therm is close to the γ , which takes place in the theory of quantum turbulence, L turb = γ 2 turb v 2 ns (see, e.g., [18]). However, in other aspects, the set of thermodynamic vortex loops related to Z 1 V ns has properties different from the loops observed in the theory of quantum turbulence. In particular, the average loop length ⟨ ⟩ l is of the order of the coherence length a , in contrast to the QT state, where ⟨ ⟩ ∼ l 1 / √ L .This is understandable, since the main source of length is the thermal effect, and k B T/ε V ∽ a . Furthermore, the interline space 1 /L therm ( v ns ) far exceeds the radius of curvature and sizes of loops (which are of the order of a ), in contrast to the case of quantum turbulence. The same can be said about the distribution of the loops in space of there lengths. In the theory of quantum turbulence, the loop sizes are gathered on scales close to the intervortex distance 1 / L .

FIG. 2: (Color online) Quantities ⟨ Pv ns ⟩ /v 2 ns and energy ⟨ E ( { } ⟩ s ) /v 2 ns . The inset shows the quantity γ therm as a function of temperature (see text)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00680v1-figure-cc9d39af-5842-4783-857a-29e51563db86.jpg)

The Lamb impulse Let's study now the Lamb impulse P and the effects associated with it. The quantity ⟨ Pv ns ⟩ can be evaluated from (6), applying it to the velocity-dependent part of PF Z 1 V ns . In Fig.2 there are shown the quantity ⟨ Pv ns ⟩ /v 2 ns and energy ⟨ E ( { } ⟩ s ) /v 2 ns due to extra vortices induced by counterflow. It is seen that they are of the same order as it should be. Formally, the Lamb impulse P and the flow induced by it appear due to the nonzero polarization of the extra vortex loops arising from counterflow.

Suppression of the superfluid density ∆ ρ s ( v ns ) An analysis of the average ⟨ Pv ns ⟩ shows that the VT induces the superfluid flow directed against the external superfluid v s . From macroscopic point of view, in particular in hydrodynamic experiments, this additional superfluid mass current which cancels the part of applied external velocity should display itself as a suppression of the superfluid density ∆ ρ s ( v ns ). Depletion of the superfluid density ∆ ρ s ( v ns ) in counterflow was predicted in Landau's pioneering work [19] It was calculated by Khalatnikov [20] on the base of the phonon and roton model. The results are that ∆ ρ s ( v ns ) /ρ s (0) are estimated as some numerical factor multiplied by v 2 ns /c 2 1 for phonons ( c 1 is the first sound velocity), and v 2 ns / k ( B T/p 0 ) 2 for rotons ( p 0 is the momentum of the roton minimum). In paper [21] the velocity dependence of the superfluid density near the λ -point was investigated. The result is that ∆ ρ s ( v ns ) / ∆ ρ s (0) scales as v 2 ns / κ/a ( λ ) 2 ( a λ is the correlation length, singular neat T λ ).

Let's study what the vortex model gives us. The quantity ⟨ Pv ns ⟩ can be interpreted as ∆ ρ s ( v ns ) v 2 ns . Developing ⟨ Pv ns ⟩ , extracted from (6) in powers of velocity v ns we get the first nonvanishing term is ∝ v 2 ns , therefore there is a velocity independent contribution into ∆ ρ s . The numerical value of ∆ ρ s can be extracted directly from Fig.2. This contribution appears since excessive (due to counterflow) vortices may be considered as quasiparticles (additional to rotons and phonons), contributing into the superfluid density (and other thermodynamic properties, such as heat capacity).

The next term of the expansion (fourth order v 2 ns ) cor- responds to the v 2 ns dependence discussed above. One can be obtain that the relative shift ∆ ρ s ( v ns ) /ρ s scales as v 2 ns / κ/l ( av ) 2 , where l av is average size of loop, close to correlation length a λ . Literally, this result is close to the one proposed in [21], although with a small numerical factor.

Shift T λ There is a number of works, studying an impact of the counterflow on the shift T λ (see, e.g. review [22]). The physics of this process is quite obvious. The kinetic energy of the counterflow changes the free energy of the entire system, which naturally leads to a shift in the transition point. This scheme had been used in various theoretical models, such as the so called model F [23],[24], the phenomenological theory of superfluidity of helium near the λ point [25], a dynamic new look at the lambda transition [22]. On the basis of these considerations as well as on experimental data, cumulative formula was proposed for the relative shift T λ -T cr ,

$$( T _ { \lambda } - T _ { c r } ) / T _ { \lambda } = ( q / q _ { 0 } ) ^ { 1 / 2 \nu }. \text{ \quad \ \ } ( 9 )$$

Here, the index ν = 2 3 is the crucial quantity in the / theory of the lambda transition , and q 0 spans from 1000 W/cm 2 up to 10000 W/cm . 2

Let's consider the impact of the studying here vortex model on the shift of the λ point. The free energy F of vortices is related to the one loop PF as F = -k B T ln Z = -k B TZ 1 (see (7)). The quantity Z 1 consists of velocity-free PF Z 1 free and the velocity dependent part Z 1 V ns . The former leads to free energy, which is positive and the latter, associated with counterflow is negative. Accordingly, a phase transition occurs when they balance each other. To move on we have to take all the parameters entering Eq. (6) in vicinity of critical temperature T cr . We adopt the following relations for the superfluid density and the correlation length near T λ (see, e.g., [25], [22])

$$\text{model}. \quad a _ { \lambda } = 4. 7 5 \ast 1 0 ^ { - 8 } \, ( T _ { c r } - T ) ^ { - 2 / 3 }, \quad \rho _ { s \lambda } = 0. 2 \, ( T _ { c r } - T ) ^ { 2 / 3 }. \\ \text{onons} \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdelta \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdclots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \ddots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \vdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \dots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \ \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \ \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots$$

Using these data we get Z 1 free ∝ ( T cr -T ) 2 which agrees with the behaviour of free energy near T λ in the resting He II [25]. Thus, at T cr the free energy associated with the vortex system vanishes as it should be. The velocity dependent part of PF Z 1 V ns changes the whole free energy. To move further we rid of the relative velocity v ns , expressing it via the heat flux q with the use v ns = ρq/ρ ST s . Substituting then the values a , λ ρ sλ (10) into Z 1 V ns , we get Z 1 V ns ∝ q 2 ( T cr -T ) -2 / 3 . Equating the obtained values Z 1 free and Z 1 V ns , and taking into account that at zero counterflow q = 0, T cr coincides with T λ we finally have (see Fig. (3))

$$q = 2. 8 3 * 1 0 ^ { 3 } ( T _ { \lambda } - T _ { c r } ) ^ { 4 / 3 }. \text{ \quad \ \ } ( 1 1 )$$

This relation agrees with the cumulative formula (9). As both theoretical estimates and experimental data, although they give close results with respect to the exponent (4 / 3), lead to a fairly large scatter in the heat flux

FIG. 3: (Color online) Shift T λ as a function of applied heat flux (see Eq. 11). Superfluid phase is inside hatched area.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00680v1-figure-5300df3b-45cf-4f74-91fa-90a8d33a9dd1.jpg)

q 0 . One of the possible reasons for this discrepancy is related to the difficulty of interpreting experimental data, in which both heat transfer and gravity have a significant effect.

## IV. DISCUSSION AND CONCLUSION

The problem of thermodynamic quantized vortices in the counterflowing superfluids, has been studied. Based on the calculated PF, some properties of superfluids due to counterflow were properly investigated. The results are in good agreement with the data obtained on the basis of other formalisms and from experiment.

The first important conclusion that follows from this is that this can be considered a fairly strong argument in favor of the point of view that the gas of topological

- [1] L. Onsager, Nuovo cimento 6 , 249 (1949).
- [2] R. P. Feynman, 'Progress in low temperature physics, vol. 1,' (North-Holland, Amsterdam, 1955) p. 17.
- [3] J. M. Kosterlitz, Rev. Mod. Phys. 89 , 040501 (2017).
- [4] H. Kleinert, Gauge Fields in Condensed Matter Physics (World Scientific, Singapore, 1991).
- [5] H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets (World Scientific, Singapore, 1990).
- [6] A. Vilenkin and E. P. S. Shellard, Cosmic strings and other topological defects , Cambridge monographs on mathematical physics (Cambridge Univ. Press, Cambridge, 1994).
- [7] E. Copeland, D. Haws, S. Holbraad, and R. Rivers, Physica A: Statistical Mechanics and its Applications 179 , 507 (1991).
- [8] G. A. Williams, Phys. Rev. Lett. 82 , 1201 (1999).
- [9] F. Lund, A. Reisenegger, and C. Utreras, Phys. Rev. B 41 , 155 (1990).
- [10] A. Chorin, Vorticity and turbulence , Applied mathematical sciences (Springer-Verlag, 1994).
- [11] S. K. Nemirovskii, J. Low. Temp. Phys. 185 , 365 (2016).
- [12] S. K. Nemirovskii, Physics Reports 524 , 85 (2013).
- [13] S. Nemirovskii, J. Low. Temp. Phys. 208 , 386-393 (2022)

excitations near the phase transition can indeed play the same role of quasiparticles as photons and rotons at low temperatures. We refrain here from a tempting discussion of the question of the possible connection between rotons and brownian vortex loops.

One more conclusion, having important physical consequence, is that there exists another population of vortex loops, which we call hydrodynamic vortices. Their evolution is mainly due to the relative velocity v ns , entering into Gibbs distribution (1) -(4) as well as into partition function (6). From the form of the Hamiltonian H { } s it follows that this population is concentrated near sizes where the energy E { } s competes with the Lamb impulse P multiplied by v ns . Using the saddle point approximation one can conclude that most favoured scale both size loop size and for the mean radius of curvature is of the order κ/v ns . This result is close to what takes place in the theory of superfluid turbulence - a very challenging topic, actively discussed at present. Our result gives hope to understand the principal mechanisms responsible for the formation the structure of the vortex tangle in the stochastic regime. This question, as well as other questions concerning the relationship of thermodynamic equilibrium vortices with the quantum turbulence, is of great interest and will be investigated in the future.

Study of Gibbs distribution and partition function had been conducted under state contract with IT SB RAS (No. 17-117022850027-5), study of physical applications was financially supported by the Russian Science Foundation (Grant No. 23-22-00128).

- [14] R. Donnelly, Quantized Vortices in Helium II (Cambridge University Press, Cambridge, UK, 1991).
- [15] V. Ambegaokar, B. I. Halperin, D. R. Nelson, and E. D. Siggia, Phys. Rev. B 21 , 1806 (1980).
- [16] N. D. Antunes, L. M. A. Bettencourt, and M. Hindmarsh, Physical Review Letters 80 , 908 (1998).
- [17] L. Burakovsky, D. L. Preston, and R. R. Silbar, Phys. Rev. B 61 , 15011 (2000).
- [18] L. Kondaurova, V. L'vov, A. Pomyalov, and I. Procaccia, Phys. Rev. B 89 , 014502 (2014).
- [19] L. Landau, J. Phys 5 , 71 (1941).
- [20] I. Khalatnikov, An Introduction to the Theory of Superfluidity (Perseus Publishing, Cambridge, MA, 2000) p. 206pp.
- [21] V. Ginzburg and A. Sobyanin, Physics Letters A 69 , 417 (1979).
- [22] D. Goodstein and A. R. Chatto, American Journal of Physics 71 , 850 (2003)
- [23] P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49 , 435 (1977).
- [24] A. Onuki, J. Low Temp. Phys. 51 , 601 (1983).
- [25] V. L. Ginzburg and A. A. Sobyanin, Usp. Fiz. Nauk 120 , 153 (1976).
Sergey K. Nemirovskii
Thermodynamics of random walking vortex loops in counterflowing superfluids
Based on the theory of the thermodynamic equilibrium in a system of quantum
vortices in superfluids in the presence of a counterflow, the influence of a
vortex tangle on various thermodynamic phenomena in quantum liquids is studied.
Using the early calculated partition function we study some of the properties
of He II related to counterflow, such as the distribution of vortex loops in
their length, the suppression of the superfluid density $\rho _{s}$ and the
shift $T_{\lambda}$. Good agreement with the early obtained results is a fairly
strong argument in favor of the point of view that the gas of string-like
topological excitations can indeed be considered as a additional kind of
quasi-particles having the inner structure at high temperatures, especially
near the phase transition. The application of the developed formalism to the
theory of quantum turbulence is briefly discussed.