## Jeans Instability in Superfluids

Itamar Hason and Yaron Oz Raymond and Beverly Sackler School of Physics and Astronomy, Tel-Aviv University, Tel-Aviv 69978, Israel (Dated: August 4, 2014)

## Abstract

We analyze the effect of a gravitational field on the sound modes of superfluids. We derive an instability condition that generalizes the well known Jeans instability of the sound mode in normal fluids. We discuss potential experimental implications.

PACS numbers: 47.37.+q, 26.60.-c

## I. INTRODUCTION

The Jeans instability of the sound mode in normal fluids that is caused by a gravitational field is a well known phenomenon [1]. It exhibits itself in astrophysical scenarios of aggregation of masses and galaxy formation. The Jeans dispersion relation of a normal fluid sound mode is obtained by linearizing the fluid equations in the presence of a gravitational field. It modifies the sound mode dispersion relation ω 2 = u /vector 2 1 k 2 to

$$\omega ^ { 2 } = u _ { 1 } ^ { 2 } \vec { k } ^ { 2 } - 4 \pi G \rho \.$$

ω is the frequency, /vector k is the momentum vector, u 1 is the speed of sound in a normal fluid u 2 1 = ∂p ∂ρ evaluated at fixed entropy, p is the pressure, ρ the local mass density and G is the gravitational coupling. Instability occurs when the RHS of (1) is negative

$$\vec { k } ^ { 2 } < \frac { 4 \pi G \rho } { u _ { 1 } ^ { 2 } } \.$$

It determines a Jeans length scale λ J , where all scales larger than that being unstable to a gravitational collapse.

In this letter we will analyze the effect of a gravitational field on the sound modes in superfluids. Superfluids are quantum fluids, i.e. fluids at temperatures close to zero where quantum effects are of primary importance [2]. They exhibit remarkable properties such as the ability to flow without viscosity in narrow capillaries. Superfluidity/superconductivity is also expected to be realized in Neutron stars matter (see e.g. [3]), as well as in high density phases of QCD [4].

The hydrodynamic description of superfluids consists of two separate motions, a normal flow and a super flow with densities ρ n and ρ s respectively, ρ n + ρ s = ρ [5, 6]. The super flow moves without viscosity, the normal flow is viscous, and the two flows do not exchange momentum between them. We denote by /vector v n and /vector v s the velocities of the normal and super flows, respectively. The superfluid velocity /vector v s corresponds to the gradient of the condensate phase that breaks spontaneously the global/local symmetry, which results in superfluidity/superconductivity, respectively. Thus, the superfluid motion is irrotational

$$\vec { \nabla } \times \vec { v } _ { s } = 0 \.$$

There are two sound modes in superfluids: the first sound u 1 which is a density wave as in normal fluids, and the second sound u 2 , which is a temperature wave and is unique to

superfluids. Its dispersion relation reads

$$\omega ^ { 2 } = u _ { 2 } ^ { 2 } \vec { k } ^ { 2 }, \quad u _ { 2 } ^ { 2 } = \frac { \rho _ { s } s ^ { 2 } } { \rho _ { n } \frac { \partial s } { \partial T } }$$

where T is the temperature, s is the entropy density per particle and the derivative is taken at fixed pressure. As expected, the second sound vanishes in the limit ρ s → 0. In general the waves can be a superposition of the two.

As we will show, in the absence of density fluctuations the pure second sound (4) remains stable, while density fluctuations imply a stability criterion for both sound modes and their superpositions. It reads

$$\vec { k } ^ { 2 } < \frac { 4 \pi G \rho } { u _ { 1 } ^ { 2 } + \frac { 1 } { ( 1 - J ) } u _ { 2 } ^ { 2 } }$$

where J = ∂s ∂ρ ∂p ∂T ∂s ∂ρ ∂T ∂p . This condition generalizes the ordinary Jeans instability (2) and reduces to it in the limit ρ s → 0. The RHS of (5) depends on the details of the system thermodynamic properties. An interesting limit to take is that of very low temperatures T → 0, where the super flow component is dominant. In this limit one has u 2 1 = 3 u 2 2 [6] and J → 0 [7]. Thus, at very low temperaure we have a larger Jeans length than that of a normal fluid by a factor of 2 √ 3 .

The letter is organized as follows. We will introduce a gravitational potential to the superfluid hydrodynamics equations and study fluctuations at the linear order. We will then analyze the instability conditions and derive (5), which is our main result. Finally, we will discuss potential experimental implications.

## II. GRAVITATIONAL INSTABILITY IN SUPERFLUID HYDRODYNAMICS

We will consider the evolution equations of a superfluid in the presence of a gravitational field. We denote the gravitational potential by φ . It satisfies Gauss's law

$$\nabla ^ { 2 } \phi = 4 \pi G \rho \.$$

## A. Self-gravitating superfluid hydrodynamics

We will consider an ideal superfluid. The superfluid density current reads

$$\vec { j } = \rho _ { n } \vec { v } _ { n } + \rho _ { s } \vec { v } _ { s }$$

and it satisfies a continuity equation, which is not affected by the gravitational field

$$\frac { \partial \rho } { \partial t } + \vec { \nabla } \cdot \vec { j } = 0 \.$$

Similarly, the entropy conservation is not affected by the gravitational field and reads

$$\frac { \partial } { \partial t } \left ( \rho s \right ) + \left ( \rho s \right ) \vec { \nabla } \cdot \vec { v } _ { n } = 0 \.$$

Note, that entropy is carried only by the normal flow and not by the super flow.

The energy and momentum conservation equations, however, are modified in the presence of the gravitational field and take the form

$$\frac { \partial j _ { i } } { \partial t } + \frac { \partial \Pi _ { i k } } { \partial x _ { k } } + \rho \frac { \partial \phi } { \partial x _ { i } } = 0$$

$$\frac { \partial E } { \partial t } + \vec { \nabla } \cdot \vec { Q } + \vec { j } \cdot \vec { \nabla } \phi = 0 \.$$

Π ik is the momentum flux tensor and /vector Q is the energy flux as calculated without a gravitational field in [6].

In order to derive the the modified equations in the presence of the gravitational field, one uses thermodynamics and the Galilean principle [6]. One denotes by K 0 the reference frame, where the super flow velocity is zero. The velocity of the normal flow in this frame is /vector v n -/vector v s . In this frame one has

$$E _ { 0 } = - p + T \rho s + \mu \rho + \phi \rho + ( \vec { v } _ { n } - \vec { v } _ { s } ) \cdot \vec { j } _ { 0 }$$

and

$$d \left ( \mu + \phi \right ) = - s d T + \frac { 1 } { \rho } d p - \frac { \rho _ { n } } { \rho } \left ( \vec { v } _ { n } - \vec { v } _ { s } \right ) d \left ( \vec { v } _ { n } - \vec { v } _ { s } \right ) \.$$

E 0 and j 0 are the energy and density current in the system K 0 , respectively. There are two effects of the gravitational field. First, to replace the chemical potential µ by µ + . Second, φ to introduce a new term in the fluid flow: in addition to the force term -∇ /vector p , we have - ∇ ρ /vector φ .

Finally, the super flow being a potential flow satisfies

$$\frac { \partial \vec { v } _ { s } } { \partial t } + \vec { \nabla } \left ( \frac { 1 } { 2 } \vec { v } _ { s } ^ { 2 } + \mu + \phi \right ) = 0 \. \quad \ \ \ ( 1 4 ) \\. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

It is straightforward to see that the hydrodynamics equations together with the thermodynamic relations constitute a complete set that determines all the charge densities and velocities.

## B. Self-gravitating superfluid sound

Linearizing the above ideal superfluid hydrodynamics equations in a presence of the gravitational field we get

$$\mathfrak { t } \quad & \frac { \partial \rho } { \partial t } + \vec { \nabla } \cdot \vec { j } = 0 \\ & \frac { \partial } { \partial t } \left ( \rho s \right ) + ( \rho s ) \, \vec { \nabla } \cdot \vec { v } _ { n } = 0 \\ & \frac { \partial \vec { j } } { \partial t } + \vec { \nabla } p + \rho \vec { \nabla } \phi = 0 \\ & \frac { \partial \vec { v } _ { s } } { \partial t } + \vec { \nabla } \mu + \vec { \nabla } \phi = 0 \, \\ \Delta ^ { 2 } \alpha$$

from which we derive

$$\frac { \partial ^ { 2 } \rho } { \partial t ^ { 2 } } & = \nabla ^ { 2 } p + \rho \nabla ^ { 2 } \phi, \\ \frac { \partial ^ { 2 } s } { \partial t ^ { 2 } } & = \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \nabla ^ { 2 } T + \frac { \rho _ { s } } { \rho _ { n } } s \nabla ^ { 2 } \phi \. \\ \intertext { a / 1 \in \bigvee \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \delta }$$

Perturbing 1 around (16), using Gauss's law for the gravitational potential, and expressing all quantities in terms of the pressure p and the temperature T we get

$$\begin{array} { c c c c c } \cdot & & \cdot & \cdot & \cdot & \cdot & \cdot & \searrow \\ & \frac { \partial \rho } { \partial p } \frac { \partial ^ { 2 } \delta p } { \partial t ^ { 2 } } + \frac { \partial \rho } { \partial T } \frac { \partial ^ { 2 } \delta T } { \partial t ^ { 2 } } = \nabla ^ { 2 } \delta p + 4 \pi G \rho \left ( \frac { \partial \rho } { \partial p } \delta p + \frac { \partial \rho } { \partial T } \delta T \right ) \\ & \frac { \partial s } { \partial p } \frac { \partial ^ { 2 } \delta p } { \partial t ^ { 2 } } + \frac { \partial s } { \partial T } \frac { \partial ^ { 2 } \delta T } { \partial t ^ { 2 } } = \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \nabla ^ { 2 } \delta T + 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \left ( \frac { \partial \rho } { \partial p } \delta p + \frac { \partial \rho } { \partial T } \delta T \right ) \. \\ \searrow \cdots & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \end{array} \cdots \quad \begin{array} { c } \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \end{array} \cdots \quad \begin{array} { c } \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \\ \end{array}$$

This system of equations for a wave of the form exp ( i /vector k · /vector x -iωt ) reads

$$& \left ( \frac { \partial \rho } { \partial p } \left ( \omega ^ { 2 } + 4 \pi G \rho \right ) - \vec { k } ^ { 2 } \right ) \delta p + \frac { \partial \rho } { \partial T } \left ( \omega ^ { 2 } + 4 \pi G \rho \right ) \delta T = 0 \\ & \left ( \frac { \partial s } { \partial p } \omega ^ { 2 } + \frac { \partial \rho } { \partial p } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \right ) \delta p + \left ( \frac { \partial s } { \partial T } \omega ^ { 2 } - \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \vec { k } ^ { 2 } + \frac { \partial \rho } { \partial T } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \right ) \delta T = 0 \. \\ \text{In the absence of a variational field } ( C \rightarrow \partial ) \text{ the equations reduce to the known analysis}$$

In the absence of a gravitational field ( G → 0) the equations reduce to the known analysis of superfluid sound modes [6].

In the general case, for a solution of these system of equations we require the determinant to vanish

$$H \omega ^ { 4 } + \left ( 4 \pi G H \rho - \left ( \frac { \partial \rho } { \partial p } \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } + \frac { \partial s } { \partial T } \right ) \vec { k } ^ { 2 } \right ) \omega ^ { 2 } + \left ( \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \vec { k } ^ { 2 } - 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \left ( \frac { \partial \rho } { \partial p } \rho s + \frac { \partial \rho } { \partial T } \right ) \right ) \vec { k } ^ { 2 } = 0$$

1 As in the ordinary Jeans instability analysis, we perturb around a stable state and therefore ignore the zeroth order contribution of the Laplacian of the gravitational potential ∇ 2 φ 0 .

where H = ∂ ρ,s ( ) ∂ p,T ( ) is the Jacobian for changing variables from ( ρ, s ) to ( p, T ), and we assume that H &gt; 0. Solving for ω 2 we get

$$2 H \omega ^ { 2 } \, & = \, - \left ( 4 \pi G H \rho - \left ( \frac { \partial \rho } { \partial p } \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } + \frac { \partial s } { \partial T } \right ) \vec { k } ^ { 2 } \right ) \\ & \, \pm \, \sqrt { \left ( 4 \pi G H \rho + \left ( \frac { \partial \rho } { \partial p } \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } - \frac { \partial s } { \partial T } \right ) \vec { k } ^ { 2 } \right ) ^ { 2 } + 4 \frac { \partial \rho } { \partial T } \frac { \rho _ { s } } { \rho _ { n } } s \vec { k } ^ { 2 } \left ( 4 \pi G H + \frac { \partial s } { \partial p } s \vec { k } ^ { 2 } \right ) } \,. \, ( 2 0 ) \\ \text{In the absence of a gravitational field. we have the plus and minus sign solutions reducing}$$

In the absence of a gravitational field, we have the plus and minus sign solutions reducing to the first and second sounds, respectively.

Consider the leading effect of the super flow in the dimensionless parameter ρ s ρ n /lessmuch 1. At first order we get

$$2 H \omega ^ { 2 } & = \, - \left ( 4 \pi G H \rho - \frac { \partial s } { \partial T } \vec { k } ^ { 2 } \right ) \pm \left ( 4 \pi G H \rho - \frac { \partial s } { \partial T } \vec { k } ^ { 2 } \right ) \\ & \quad + \frac { \rho _ { s } } { \rho _ { n } } \left ( \frac { \partial \rho } { \partial p } s ^ { 2 } \pm \frac { \partial \rho } { \partial p } s ^ { 2 } \pm \frac { 2 \frac { \partial \rho } { \partial T } s \vec { k } ^ { 2 } \left ( 4 \pi G H + \frac { \partial s } { \partial p } s \vec { k } ^ { 2 } \right ) } { \left ( 4 \pi G H \rho - \frac { \partial s } { \partial T } \vec { k } ^ { 2 } \right ) } \right ) \,. \\ \intertext { o r e w e h e f i s t s o n d } \text{solution}$$

Therefore, we have the first sound solution

$$H \omega ^ { 2 } = & - \left ( 4 \pi G H \rho - \frac { \partial s } { \partial T } \vec { k } ^ { 2 } \right ) - \frac { \rho _ { s } } { \rho _ { n } } \left ( \frac { \frac { \partial \rho } { \partial T } s \vec { k } ^ { 2 } \left ( 4 \pi G H + \frac { \partial s } { \partial p } s \vec { k } ^ { 2 } \right ) } { \left ( 4 \pi G H \rho - \frac { \partial s } { \partial T } \vec { k } ^ { 2 } \right ) } \right ) \, \quad \ \ ( 2 2 ) \\ \text{and the second sound solution}$$

and the second sound solution

$$\cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \\ H \omega ^ { 2 } & = \frac { \rho _ { s } } { \rho _ { n } } \left ( \frac { \partial \rho } { \partial p } s ^ { 2 } + \frac { \frac { \partial \rho } { \partial T } s \vec { k } ^ { 2 } \left ( 4 \pi G H + \frac { \partial s } { \partial p } s \vec { k } ^ { 2 } \right ) } { \left ( 4 \pi G H \rho - \frac { \partial s } { \partial T } \vec { k } ^ { 2 } \right ) } \right ) \. \\ \\ \sim \quad \cdot \quad \cdot \quad \cdot \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots$$

## C. Comparison to the Non-Gravitational Case

It is instructive to compare special solutions to the superfluid sound equations in the presence of the gravitational interaction (18) to the non-gravitating case.

## 1. The Pure Second Sound Limit

In the limit ∂ρ ∂T = 0, the non-gravitational sound equations (setting G = 0 in (18)) read

$$\begin{array} { c } \psi \end{array} \left ( \frac { \partial \rho } { \partial p } \omega ^ { 2 } - \vec { k } ^ { 2 } \right ) \delta p = 0 \\ \left ( \frac { \partial s } { \partial p } \omega ^ { 2 } \right ) \delta p + \left ( \frac { \partial s } { \partial T } \omega ^ { 2 } - \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \vec { k } ^ { 2 } \right ) \delta T = 0 \. \\ \hat { \ } \end{array} \right ) \hat { \ } \end{array}$$

One has a solution where δp = 0, δT arbitrary, and the dispersion relation (4) of the pure second sound describing the temperature wave. There is another solution, where the dispersion relation is that of a first sound mode, and the pressure and temperature variations are related by

$$\left ( \frac { \partial s } { \partial p } \right ) \delta p = \left ( \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \frac { \partial \rho } { \partial p } - \frac { \partial s } { \partial T } \right ) \delta T \. \\ \intertext { w h e o r a v i t v } \text{ still having } \underline { \partial \rho } - \partial \text{ } w o o t$$

When we introduce gravity, still having ∂ρ ∂T = 0, we get

$$\begin{array} { c } \Lambda _ { \nu } & \Lambda _ { \nu } & \Lambda _ { \nu } & \Lambda _ { \nu } \\ \left ( \frac { \partial \rho } { \partial p } \left ( \omega ^ { 2 } + 4 \pi G \rho \right ) - \vec { k } ^ { 2 } \right ) \delta p = 0 \\ \left ( \frac { \partial s } { \partial p } \omega ^ { 2 } + \frac { \partial \rho } { \partial p } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \right ) \delta p + \left ( \frac { \partial s } { \partial T } \omega ^ { 2 } - \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \vec { k } ^ { 2 } \right ) \delta T = 0 \. \\ \text{see that the pure second sound solution (4) is not affected by gravity. This is expected} \end{array}$$

We see that the pure second sound solution (4) is not affected by gravity. This is expected since the pure second sound is characterized by no variation of the density, thus gravity does not affect it.

However, the first sound solution is changed. The dispersion relation is the same as the Jeans dispersion relation (1), but the pressure and temperature variations are related by a different formula

$$\left ( \frac { \partial s } { \partial p } \omega ^ { 2 } + \frac { \partial \rho } { \partial p } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \right ) \delta p & = \left ( \left ( \frac { \partial \rho } { \partial p } \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } - \frac { \partial s } { \partial T } \right ) \omega ^ { 2 } + \frac { \partial \rho } { \partial p } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } \rho s ^ { 2 } \right ) \delta T \. \quad \left ( 2 7 \right ) \\ \text{$c$-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text$-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\tt{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\tt {-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text	-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text.-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text--\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{--\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text		-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\tt-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text
-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text}-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text -\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text1-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text'-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\textc-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{-\text{--$$

Comparing (25) and (27), we see a frequency dependence of the linear relation between pressure and temperature (27), which vanishes as ρ s → 0.

## 2. The Pure First Sound Limit

Consider next the classical pure first sound limit, ∂s ∂p = 0. The non-gravitational sound equations in this case are (setting G = 0 and ∂s ∂p = 0 in (18)) :

$$\lambda \lambda _ { \lambda } \varphi _ { \varphi } \lambda _ { \lambda } \bar { \ } \left ( \frac { \partial \rho } { \partial p } \omega ^ { 2 } - \vec { k } ^ { 2 } \right ) \delta p + \left ( \frac { \partial \rho } { \partial T } \omega ^ { 2 } \right ) \delta T = 0 \\ \left ( \frac { \partial s } { \partial T } \omega ^ { 2 } - \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \vec { k } ^ { 2 } \right ) \delta T = 0 \. \\ \text{on where } \delta T = 0, \, \delta p \text{ arbitrary, and the dispersion relation of the pure first}$$

We have a solution where δT = 0, δp arbitrary, and the dispersion relation of the pure first sound. There is another solution, where the dispersion relation is (4), and the pressure and temperature variations are related by

$$\left ( \frac { \partial s } { \partial T } \frac { \rho _ { n } } { \rho _ { s } } \frac { 1 } { s ^ { 2 } } - \frac { \partial \rho } { \partial p } \right ) \delta p = \left ( \frac { \partial \rho } { \partial T } \right ) \delta T \.$$

When we introduce gravity, still with ∂s ∂p = 0, we get

$$\left ( \frac { \partial \rho } { \partial p } \left ( \omega ^ { 2 } + 4 \pi G \rho \right ) - \vec { k } ^ { 2 } \right ) \delta p + \frac { \partial \rho } { \partial T } \left ( \omega ^ { 2 } + 4 \pi G \rho \right ) \delta T = 0 \\ \left ( \frac { \partial \rho } { \partial p } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \right ) \delta p + \left ( \frac { \partial s } { \partial T } \omega ^ { 2 } - \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } \vec { k } ^ { 2 } + \frac { \partial \rho } { \partial T } \cdot 4 \pi G \frac { \rho _ { s } } { \rho _ { n } } s \right ) \delta T = 0 \. \\ \text{We see that the pure first sound solution is not pure anvmore.  This is in contrast to the}$$

We see that the pure first sound solution is not pure anymore. This is in contrast to the pure second sound solution that was not modified by gravity.

## D. General Instability Conditions

In the following we will analyze the instability conditions arising from the roots (20).

## 1. Identical roots

Consider the case where the two solutions (20) are identical. In this case we have the same dispersion relation for the first and second sound and we obtain the instability condition

$$\vec { k } ^ { 2 } & < \frac { 4 \pi G H \rho } { \left ( \frac { \partial \rho } { \partial p } \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } + \frac { \partial s } { \partial T } \right ) } \. \\ \text{react } ( 3 1 ) & \text{ as } ( 5 ). \text{ The constraint } ( 5 ) \text{ is a generalization of } ( 2 ) \text{ to }$$

Dividing by H , we can recast (31) as (5). The constraint (5) is a generalization of (2) to superflows, and is valid for both the first and second sounds. It reduces to (2) in the limit ρ s → 0. In limit T → 0, u 2 1 = 3 u 2 2 [6], and one can see from the data presented in [7] that J → 0. Thus we have that the Jeans length of superfluid is larger by a factor of 2 √ 3 than that of a normal fluid.

## 2. Different roots

The conditions for instability of the two sound modes are the upper bound on the momentum (5) as well as a lower bound

$$\vec { k } ^ { 2 } & > 4 \pi G \rho \left ( \frac { \partial \rho } { \partial p } + \frac { \partial \rho } { \partial T } \frac { 1 } { \rho s } \right ) \,, \\. & \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon$$

which together imply a constraint on the system independent of the gravitational field

$$\left ( \frac { \partial \rho } { \partial p } + \frac { \partial \rho } { \partial T } \frac { 1 } { \rho s } \right ) \left ( \frac { \partial \rho } { \partial p } \frac { \rho _ { s } } { \rho _ { n } } s ^ { 2 } + \frac { \partial s } { \partial T } \right ) < H \.$$

## E. Discussion

As the ordinary Jeans instability has been observed experimentally, it is natural to inquire how can we observe the generalized instability conditions for superfluids. It is probably unlikely that current condensed matter experiments can detect the gravitational instability of superfluids. A potential experimental laboratory could be Neutron stars. Superfluidity is expected to occur both in the inner crust and the core of Neutron stars as a consequence of the formation of Cooper pairs of Neutrons that break spontaneously the U (1) Baryon symmetry [3]. Superfluidity has been invoked in order to explain, for instance, pulsar glitches [8]. Therefore the gravitational instability seems relevant to the dynamics of Neutron stars, such as its rotational and vibrational oscillations. It would be interesting to work out the details of the gravitational instability effects in this framework.

## Acknowledgements

This work is supported in part by the I-CORE program of Planning and Budgeting Committee (grant number 1937/12), and by the US-Israel Binational Science Foundation.

- [1] J. H. Jeans, 'The Stability of a Spherical Nebula'. Philosophical Transactions of the Royal Society A 199 , 1 (1902).
- [2] P. Kapitza, 'Viscosity of Liquid Helium Below the λ -Point', Nature 141 , 74 (1938).
- [3] U. Lombardo and H. J. Schulze, Lect. Notes Phys. 578 , 30 (2001) [astro-ph/0012209].
- [4] M. G. Alford, K. Rajagopal and F. Wilczek, 'QCD at finite baryon density: Nucleon droplets and color superconductivity,' Phys. Lett. B 422 , 247 (1998) [hep-ph/9711395].
- [5] L. D. Landau, 'The Theory of Superfluidity of Helium II', J. Phys. USSR, 5 , 71(1941).
- [6] E. M. Lifshitz and L. D. Landau, 'Course of Theoretical Physics, volume 6, Fluid Mechanics', Pergamon Press (1959).
- [7] J. Brooks and R. Donnelly, 'The Calculated Thermodynamic Properties of Superfluid Helium4', J. Phys. Chem. Data, 6 , 1 (1977).
- [8] P.W. Anderson and N. Itoh, 'Pulsar glitches and restlessness as a hard superfluidity phenomenon', Nature 256 , 25 (1975);.
Itamar Hason Yaron Oz
Jeans Instability in Superfluids
We analyze the effect of a gravitational field on the sound modes of
superfluids. We derive an instability condition that generalizes the well known
Jeans instability of the sound mode in normal fluids. We discuss potential
experimental implications.