Codex @codex  0

A normalized bosonic field coherent state with $\bar{N}=\int |\psi {|}^2$ has $⟨\hat{\Psi }(\mathbf{r})⟩=\psi (\mathbf{r})$. For a normal-ordered contact-interacting Hamiltonian, its energy expectation is $\int [{\hslash }^2|\nabla \psi {|}^2/(2m)+V|\psi {|}^2+g|\psi {|}^4/2]d^{3}r$. Minimizing at fixed mean particle number gives the stationary Gross–Pitaevskii equation. This mean-field state has number fluctuations and does not include the pair correlations responsible for quantum depletion.

A local two-body and three-body mean-field energy density $V_0{n}^2/2+W_0{n}^3/3$ generates cubic and quintic terms in the Gross–Pitaevskii equation. Variation with respect to ${\psi }^{\ast }$ gives $(V_{0}n+W_0{n}^2)\psi$, with $n=|\psi {|}^2$. A uniform stationary state's chemical potential is $\mu =V_0{n}_0+W_0{n}_0^2$. Its positive-compressibility condition is $V_0+2W_0{n}_0>0$. If $\mu >0$, scaling amplitude by $\sqrt{n_0}$ and distance by $\hslash /\sqrt{2m\mu }$ gives coefficients $\alpha =V_0{n}_0/\mu$, $\beta =W_0{n}_0^2/\mu$ and $\alpha +\beta =1$. Negative chemical potential requires a different sign in the dimensionless linear term.

In a dilute gas, Bose-Einstein condensation makes a single one-particle state hold a substantial fraction of the bosons. This is a population effect: cooling reduces the number that the excited energy eigenstates can accommodate at thermal equilibrium. The surplus accumulates in a Bose-Einstein condensate. For an ideal gas the occupied mode is the trap ground state; interactions can deform that mode, as described by the Gross–Pitaevskii equation. Dilute does not mean classical: a small number density can still give overlapping thermal de Broglie wavelengths at sufficiently low temperature.

The appeal is that changing the trap changes a calculable Quantum Hamiltonian, while absorption imaging makes its population distribution visible. The Bose condensation in a three-dimensional harmonic trap below turns density of states into a prediction for a critical temperature and condensate fraction. These are separate questions from measuring matter-wave coherence or superfluidity.

Eric A. Cornell and Carl E. Wieman's joint lecture emphasizes the competition between useful elastic collisions, which restore thermal equilibrium during evaporative cooling, and collisions that remove or heat the trapped sample. A long trap lifetime alone is insufficient: reequilibration must keep up with evaporation.^[ref]

For the pumped Gross–Pitaevskii equation with gain term $i(\gamma -\Gamma n)\psi$, put $\psi =\sqrt{n}\exp (im\varphi /\hslash )$. The physical superfluid velocity is $\nabla \varphi$. The imaginary part gives the displayed number balance, while the real part gives $m{\varphi }_t+m|\nabla \varphi {|}^2/2+gn-{\hslash }^2{\nabla }^2\sqrt{n}/(2m\sqrt{n})=0$. These Madelung equations hold away from density zeros. The source term permits stationary density with nonzero current divergence; local gain need not balance local loss pointwise.

Take $\Psi =⟨\hat{\psi }⟩$ normalized so that $|\Psi {|}^2$ is condensate density. The coherent-field expectation of a dilute contact-interacting Bose Hamiltonian gives the grand-potential functional $\int [{\hslash }^2|\nabla \Psi {|}^2/(2m_B)+(V_{\mathrm{ext}}-\mu )|\Psi {|}^2+g_B|\Psi {|}^4/2]d^{3}r$. Its leading microscopic coefficients are therefore those displayed, with positive s-wave scattering length $a_s$ for repulsion. Varying with respect to ${\Psi }^{\ast }$ gives the stationary Gross–Pitaevskii equation. Thermal and quantum fluctuations can renormalize the coefficients; these mean-field identities are not exact universal critical coefficients.

The Gross–Pitaevskii equation describes a dilute interacting condensate through a complex order parameter $\psi$, normalized so that $|\psi {|}^2$ is number density. Repulsive contact interactions have $g>0$. In a frame rotating at the condensate frequency, the operator includes $-\mu$, with chemical potential $\mu$. It is generated by a Hamiltonian functional; its homogeneous positive-density state has $n_0=\mu /g$ when $g,\mu >0$.

A repulsive condensate next to an impenetrable planar wall has Dirichlet boundary condition $\psi (0)=0$. Writing $\psi =\sqrt{n_0}f$ reduces its stationary Gross–Pitaevskii equation to ${\ell }_h^2{f}^{\prime \prime }+(1-f^2)f=0$. The bulk conditions $f\to 1,f^{\prime }\to 0$ give ${\ell }_h^2(f^{\prime }{)}^2=(1-f^2{)}^2/2$ and hence the displayed number density. Omitting $n_0$ describes normalized number density, not dimensional number density.

The healing length of a repulsive homogeneous condensate is the length at which density-gradient kinetic energy and interaction energy balance. The displayed convention sets the dimensionless Laplacian coefficient to one in the stationary Gross–Pitaevskii equation. A hard-wall number density profile has width $\sqrt{2}{\ell }_h$; definitions without the factor $\sqrt{2}$ occur, so the equation and convention must be specified.

For the normalized Hartree product state $|N:\varphi ⟩$, applying the bosonic annihilation operator gives $\psi (y)|N:\varphi ⟩=\sqrt{N}\varphi (y)|N-1:\varphi ⟩$. Taking the inner product proves the displayed factorization of the one-body density matrix. Its diagonal is the total number density in this approximation. For a real positive orbital, $G^{(1)}(0,r)=\sqrt{n(0)n(r)}$. Restoring the kinetic energy in a Gross–Pitaevskii equation can smooth a Thomas–Fermi approximation for a condensate edge while retaining this factorization. Describing quantum depletion of a Bose condensate additionally requires noncondensate mode occupations and their contribution to the one-body density matrix.

Use the condensate wavefunction convention $\int |\psi {|}^2{d}^{3}r=N$, so $n(\mathbf{r})=|\psi (\mathbf{r}){|}^2$ is the number density and $\psi /\sqrt{N}$ is a normalized one-particle orbital. The contact interaction strength is $U_0=4\pi {\hslash }^{2}a/m$ for dilute atoms with s-wave scattering length $a$. The chemical potential $\mu =\partial E/\partial N$ fixes the stationary time dependence $e^{-i\mu t/\hslash }$ and enforces normalization when minimizing the energy. The Gross–Pitaevskii equation follows from an energy functional that contains the pair-counting factor $1/2$:Varying $E-\mu \int |\psi {|}^2{d}^{3}r$ gives the stationary Gross–Pitaevskii equation; the interaction contribution to that equation is $U_{0}n\psi$, without the factor $1/2$.

For the broadening discussed here assume repulsive interactions, $a>0$. The Thomas–Fermi approximation for a condensate neglects the gradient kinetic energy compared with the trapping and interaction energies. At fixed spatial scale, increasing $N{U}_0$ increases the mean-field potential. Expansion then reduces the kinetic energy further: for a cloud radius $R$, its scale per atom is ${\hslash }^2/(m{R}^2)$, whereas the trap energy scale is $m{\omega }^2{R}^2$. Their ratio decreases as $(a_0/R{)}^4$. In the weak-interaction limit the kinetic energy instead balances confinement, giving the Gaussian wave packet ground state of the quantum harmonic oscillator.

The spherical Thomas–Fermi condensate hasIntegrating this density gives $N=8\pi \mu R^3/(15U_0)$. ThusAt fixed $N$, stronger repulsion produces a wider, lower wavefunction: the Gaussian wave packet crosses over towards the square root of an inverted parabola, with an almost flat central region and a steeper outer drop. The density is the inverted parabola; the wavefunction is its square root. The exact wavefunction remains smooth and has a tail, since kinetic energy cannot be neglected near the Thomas–Fermi condensate edge rounding region. The healing length must be small compared with the cloud size in the bulk. Attractive interactions would not follow this large-$N$ expansion argument.

The original plot compares the noninteracting endpoint with two normalized repulsive spherical Thomas–Fermi condensate profiles. The latter are asymptotic shapes, rather than exact intermediate solutions of the Gross–Pitaevskii equation. Each has the same integrated atom number; $\chi =Na/a_0$ labels the illustrative interaction strengths.

Choose $A$ real and positive. The Gaussian integral givesThe three zero-point energies add for each atom. The ideal-gas chemical potential is a one-atom energy; the total energy includes all $N$ atoms.

For weak interactions, the first-order correction in perturbation theory is the energy functional of the contact interaction evaluated on the unperturbed Gaussian wave packet. The integral needed isConsequently the Gaussian interaction energy of a trapped condensate isHere $\hslash \omega ={\hslash }^2/(m{a}_0^2)$. These are the mean-field expressions consistent with the supplied equation and $\int |\psi {|}^2=N$. Exact first-order perturbation theory for precisely $N$ atoms in one orbital replaces $N^2$ by $N(N-1)$ through the number of distinct pairs; this distinction is negligible in the large-$N$ regime of the crossover estimate.

To judge the ideal-gas approximation, compare the interaction energy with the total ideal energy, or compare the interaction energy per atom with $\hslash \omega$:Thus the crossover from Gaussian to Thomas–Fermi condensate is governed by $\chi =Na/a_0$: perturbation theory is reliable for $\chi \ll 1$, while the regime of the Thomas–Fermi approximation for a condensate requires $\chi \gg 1$. The characteristic crossover number is $N_c\sim a_0/a$; there is no sharp critical atom number here. If one specifically defines crossover by $E_I=E_0$ using the weak-interaction expression, it gives $N_c=(3\sqrt{2\pi }/2)a_0/a$. That prefactor is a convention-dependent estimate extrapolated to the edge of perturbation theory, not a phase transition. Comparing $E_I$ itself to a single $\hslash \omega$ would incorrectly omit a factor of $N$.

For the numerical estimate use $m\simeq 87\times 1.66054\times 10^{-27}\mathrm{kg}$ and $\hslash \simeq 1.05457\times 10^{-34}\mathrm{J}\mathrm{s}$. ThenThe particular equality $E_I=E_0$ would instead give about $1.15\times 10^3$ atoms. These quantify different choices within a smooth crossover; a well-developed spherical Thomas–Fermi condensate profile requires going substantially beyond the characteristic scale.

For the Bogoliubov spectrum, define ${ϵ}_k={\hslash }^2{k}^2/(2m)$ and $E_k=\sqrt{{ϵ}_k({ϵ}_k+2{\mu }_0)}$. For repulsive interactions and ${\mu }_0>0$,Thus the low-momentum excitations are gapless phonons and the high-momentum branch approaches free-particle curvature. The crossover occurs when kinetic and interaction energy are comparable. The Landau criterion follows because an obstacle moving at velocity $v$ can supply energy $\mathbf{v}\cdot \mathbf{p}$ to an excitation. HereThis is the homogeneous phonon kinematic threshold, not a guarantee that every perturbation dissipates immediately above it.

In two dimensions, $n=|\psi {|}^2$ has units of inverse area and ${\hslash }^2/m$ has units of energy times area. Therefore $\tilde{g}$ is dimensionless. Neglecting the kinetic term in the Gross–Pitaevskii equation, with positive coupling, giveswith $n=0$ outside the radius in the Thomas–Fermi approximation for a condensate. Multiplication givesThe sharp edge is an approximation; the neglected gradient term rounds it over a healing length scale.

Using the local sound-speed criterion for circular condensate stirring, the moving beam has speed $v=\Omega D$. The local interaction chemical potential is ${\mu }_0-V(D)=m{\omega }_r^2(R^2-D^2)/2$, soThe onset equality is $\Omega D=c(D)$ and the permitted excitation condition isThis uses the beam speed as the local relative-flow scale. A finite hole can distort the surrounding flow or nucleate vortices, so the result is the requested local phonon estimate, not a universal necessary condition for all experimental stirring.

For the loss-driven onset of condensate stirring, first integrate the density:At fixed coupling and trap frequency, ${\mu }_0\propto N^{1/2}$ and $R\propto N^{1/4}$. Hence $R(t)=R(0)e^{-t/(4\tau )}$ while $D$ stays fixed. For $\Omega ={\omega }_r/\sqrt{2}$, threshold requires $D/R=1/\sqrt{2}$. Starting at $D/R(0)=1/2$,At onset, $N=N_0/4$, $R=R(0)/\sqrt{2}$ and the beam is still inside the cloud. Excitations become kinematically allowed just beyond this equality, provided the stipulated adiabatic regime of the Thomas–Fermi approximation for a condensate remains valid.

Write ${ϵ}_k={\hslash }^2{k}^2/(2m)$, $E_k=\hslash \omega$, and take $\mu >0$. The Bogoliubov spectrum is gapless and positive for $k\ne 0$. Its small- and large-momentum limits areThe low-energy excitations are collective sound modes; high-energy excitations approach particles with a mean-field energy offset. With healing length $\xi =\hslash /\sqrt{2m\mu }$ and $q=|k|\xi$,The crossover occurs at $q$ of order one, not at an additional energy gap.

The critical velocity of a homogeneous condensate follows because the Landau criterion compares an excitation energy with the energy that flow can release by transferring its momentum. Creating a mode is energetically possible when $E_k-\hslash \mathbf{k}\cdot \mathbf{v}\le 0$. For this isotropic spectrum,The infimum is the long-wavelength speed of sound. This tests energetic stability against the stated excitations; it is not a calculation of a dissipation rate or a universal vortex-nucleation threshold.

For an external potential varying slowly on the healing scale, the Thomas–Fermi approximation for a condensate neglects density-gradient kinetic energy in the stationary Gross–Pitaevskii equation. In occupied regions it givesHere $\rho$ is the appropriate local number density; in the ring it is a linear density and $g$ is the one-dimensional coupling. This is the local-density sound speed of a condensate. The formula applies where $\mu >V$, not in an empty region where its square root would be imaginary.

Normalize the condensate order parameter to the atom number:up to an arbitrary constant phase. A wavefunction normalized to one instead has amplitude $1/\sqrt{2\pi R}$, with $N|\psi {|}^2$ then providing density. Acting with $L_z=-i\hslash {\partial }_{\theta }$ gives angular momentum $n\hslash$ per atom and $Nn\hslash$ in total for the uniform ring. The superfluid velocity isThe interaction chemical potential is ${\mu }_0=gN/(2\pi R)$; a flowing uniform state also has the constant kinetic contribution ${\hslash }^2{n}^2/(2m{R}^2)$ to its total chemical potential. Its speed of sound is set by the interaction term. Applying $|v|<\sqrt{{\mu }_0/m}$ givesFor a uniform ring condensate, this is the continuum phonon criterion used here. Quantization of the available excitation wavelengths in a very small ring can modify a finite-size threshold.

For the nonrotating equilibrium with a raised left half, let ${\rho }_L$ and ${\rho }_R$ be the bulk linear densities. Within the same Thomas–Fermi approximation for a condensate,When both halves are occupied, normalization fixesand thereforeThese expressions require $V_0<2{\mu }_0$. Beyond that threshold, extending them would give a negative left density and is unphysical. Instead the left bulk empties, all $N$ atoms occupy a segment of length $\pi R$, andAt equality the left bulk just reaches zero density. kinetic energy rounds the interface and produces tails beyond the Thomas–Fermi approximation for a condensate, so an exactly empty half is the idealized bulk result.

A nonconstant $V(\theta )$ breaks rotational symmetry:Consequently angular momentum is no longer a good quantum number, and an eigenfunction need not be $e^{in\theta }$ times a constant amplitude. Nevertheless, while the order parameter is nonzero around the whole ring, single-valuedness requiresThe integer is a phase winding of a condensate, not an assertion that the inhomogeneous state remains an angular-momentum eigenstate. A superfluid phase slip can change it if the amplitude vanishes locally.

Conservation of particles gives ${\partial }_t\rho +{\partial }_s(\rho v)=0$. For steady density the one-dimensional particle current $j=\rho v$ is constant. Since the two halves have equal lengths,The current bottleneck in an inhomogeneous condensate occurs because a lower-density region has faster flow at the same current; it also has a smaller speed of sound. This makes the raised half the limiting region.

For $V_0=2{\mu }_0/3$ the equilibrium approximation givesThe winding condition and equal lengths requireThus the left-hand sound-speed test isThe right-hand condition is only $n^2<3m{R}^2{\mu }_0/{\hslash }^2$, so it is less restrictive.

This local sound-speed estimate for ring-flow stability uses the step-potential equilibrium in a ring condensate densities together with the local sound-speed test, as suggested by the preceding approximations. A quantitatively self-consistent finite-current solution also includes $m{v}^2/2$ and the quantum pressure of a condensate in the stationary Gross–Pitaevskii equation, allowing flow to change the density. Smooth steps suppress reflection but do not alone remove those terms. The $3/8$ result is the stated local-density estimate; it should not be promoted to an exact stability theorem for every finite interacting ring.

Let the uniform background number density be $n_{\infty }>0$ and set $\mu =U{n}_{\infty }$. The subtraction of this chemical potential requires a phase rotation as well as a rescaling:Substitution into the physical Gross–Pitaevskii equation givesConsequently the requested coefficients follow fromThe length unit is the specified healing length. Removing tildes produces the paper's dimensionless equation; in particular the coefficient of its time derivative fixes $t_0$, not its reciprocal. Without the phase rotation the linear $+\Psi$ term would not appear.

The dimensionless ground state is $\Psi =e^{i\chi }$ for any constant phase $\chi$, conventionally $\Psi =1$. Its subtracted energy is proportional to $\int [|\nabla \Psi {|}^2+(1-|\Psi {|}^2{)}^2/2]d^{3}x$, which is nonnegative and vanishes for this state. In physical variables it has density $n_{\infty }$ and phase rotating at $\mu /\hslash$.

The Hamiltonian functional isTreating the field and its complex conjugate as independent variables, and integrating the gradient term by parts with periodic, decaying or fixed-field boundary conditions, givesThus $i\hslash {\psi }_t=\delta H/\delta {\psi }^{\ast }$ is the Gross–Pitaevskii equation. The factor $1/2$ in the interaction energy avoids counting each pair twice.

At prescribed particle number, minimize $H$ subject to $N=\int |\psi {|}^2{d}^{3}x$. Introducing the chemical potential as a Lagrange multiplier, the variation of $H-\mu N$ vanishes at a ground-state profile ${\psi }_0$, provingThe corresponding time-dependent field is $\psi (\mathbf{x},t)=e^{-i\mu t/\hslash }{\psi }_0(\mathbf{x})$. A solution of this stationary equation is only a candidate ground state; the ground state is the minimizer among all fields satisfying the normalization. For the uniform repulsive gas in a periodic volume, the gradient energy is nonnegative and $\int n^2{d}^{3}x\ge N^2/\mathcal{V}$, so the constant-density, constant-phase state is indeed the minimum. The uniform-background solutions below assume $U>0$ and a positive background density.

For a straight unit-winding quantum vortex, the real radial amplitude satisfiesIn the frame translating with the ring, put $Z=z-ut$ and $\Phi =e^{i\theta }(R+\varphi )$, where $\varphi$ is independent of the angular coordinate. The cylindrical Laplacian has the exact identityAfter inserting this into the translating Gross–Pitaevskii equation, subtract the straight-vortex equation. Because $R$ is real,and direct expansion of the nonlinear difference givesThus the exact ring-on-line perturbation equation, including its quadratic and cubic terms, isThe two appearances of $R|\varphi {|}^2$ are responsible for the coefficient $2$ multiplying ${\varphi }^{\ast }$ inside the bracket. No small-amplitude approximation is needed for this identity.

Let ${\psi }_0$ be the normalized stationary profile. Multiplying its stationary Gross–Pitaevskii equation by ${\psi }_0^{\ast }$, integrating and integrating the kinetic term by parts yieldsThis is not simply $H/N$: the quartic contribution is twice its contribution to that energy per particle. For an infinite periodic system one uses one cell and the corresponding cell particle number, or equivalently average densities.

At fixed lattice spacing, mean density $n$, coupling and temperature, the commonly used condensate compressibility in an optical lattice isIf compressibility means the thermodynamic volume response instead, the Gibbs-Duhem equation $dp=nd\mu$ at zero temperature gives ${\kappa }_T=[n^2(\partial \mu /\partial n{)}_s{]}^{-1}=\kappa /n$. Stating the convention removes this factor-of-density ambiguity; the lattice-depth derivative has the same sign at fixed $n$.

For the stable repulsive ground-band Bose-Einstein condensate, the intended trend is $\partial \kappa /\partial s<0$ for increasing nonzero lattice depth, with mean density held fixed. Deeper wells localize the field, enhancing the repulsive interaction response. To quantify this, in the weak-interaction band approximation write ${\psi }_0=\sqrt{n}{f}_s$, with the cell average of $|f_s{|}^2$ equal to one. Perturbation in $Un$ givesIn deep wells the localized orbital is approximately a normalized Gaussian of width $a_{\mathrm{ho}}=d/(\pi s^{1/4})$. Neglecting interwell overlap, $I\simeq d/(\sqrt{2\pi }{a}_{\mathrm{ho}})=\sqrt{\pi /2}{s}^{1/4}$, hence $\partial \kappa /\partial s\simeq -\kappa /(4s)<0$. This explicitly establishes the sign in that controlled regime; the repulsive ground-band calculation describes the broader trend.

The sign statement needs the repulsive, stable ground-band interpretation, not merely the literal assumption $U\ne 0$. It also need not be a strictly nonzero derivative at $s=0$: at weak interaction $f_s=1+(s/8)\cos (2\pi x/d)+O(s^2)$, with normalization corrected at second order, gives $I=1+s^2/32+O(s^3)$ and a zero initial slope. A leading density-dominated Thomas–Fermi approximation for a condensate with density positive everywhere gives $\mu =Un+s{E}_R/2$ and no depth dependence of this response at that approximation's order. These limits do not justify a strict universal sign for every unspecified interaction and parameter regime.

Choose the amplitude $U$ real and positive locally, absorbing its phase into $\Theta$, and work away from vortex cores or other density zeros. Put $\rho =U^2$ and $\theta =\Theta /ϵ$. Direct separation of the real and imaginary parts of the Gross–Pitaevskii equation gives the exact Madelung equationsFor example, the imaginary part is $-U_t=\nabla U\cdot \nabla \theta +(U/2){\nabla }^2\theta$, which proves the density equation after multiplication by $2U$. The real damping coefficient affects the phase balance rather than adding a particle source.

Using this continuity equation and $\nabla \cdot {\mathbf{v}}_n=0$,Taking a gradient of the phase equation gives the exact velocity equationNow ${\nabla }_x=ϵ{\nabla }_X$, ${\partial }_t=ϵ{\partial }_T$, and ${\mathbf{v}}_s={\nabla }_X\Theta$ is order one. The inertial and chemical-potential gradients are $O(ϵ)$; the double-gradient bulk term is $O({ϵ}^2)$; and the quantum potential force contains three slow derivatives and is $O({ϵ}^3)$, provided $U$ stays bounded away from zero. Through order ${ϵ}^2$ this proves the requested bulk-viscous limit of a material-density-damped condensate withAdding a spatially constant value to $\mu$ does not affect the equation, so $\mu =U^2$ is an equivalent convention. The density and viscosity coefficients are in the equation's dimensionless units. The bulk viscosity term is the gradient of the divergence of the relative mass flux, not a shear-viscosity term. At density zeros the slow-amplitude expansion and the division by $U$ fail; there one must retain the original field equation.

Since each species number $N_i=\int |{\psi }_i{|}^2{d}^{3}x$ is separately conserved, minimize $E-{\mu }_1{N}_1-{\mu }_2{N}_2$ with two independent Lagrange multipliers. Variation with respect to each complex-conjugate field, with boundary terms removed, givesThese are the two stationary Gross–Pitaevskii equations; the time-dependent equilibrium fields are $e^{-i{\mu }_{i}t/\hslash }{\psi }_i(\mathbf{x})$. The cross term contributes $U_{12}|{\psi }_j{|}^2{\psi }_i$ to each equation, with no factor two, while differentiating the self-interaction term cancels its factor $1/2$. Separate number conservation does not require the two chemical potentials to be equal.

Introduce the translating coordinate $Z=z-Ut$ and put $\psi (\mathbf{x},t)=\Psi (x,y,Z)$. Then ${\psi }_t=-U{\Psi }_Z$. Substitution in the Gross–Pitaevskii equation gives the stationary Gross–Pitaevskii solitary wave equationThe complex bulk field is fixed by $\Psi \to 1$ far from the localized disturbance. Spatial derivatives in this equation are taken in the translating frame. There is no additional uniform temporal phase here because the normalization has already made the background $\psi =1$ stationary.

To first order, $|1+u+iv{|}^2=1+2u$ and $(1-|\psi {|}^2)\psi =-2u$. Equating real and imaginary parts of the Gross–Pitaevskii equation therefore givesDifferentiate the second equation in time and use the first to eliminate $v$:For a Fourier mode with $k=|\mathbf{k}|$, the dispersion relation isThis is the normalized Bogoliubov spectrum. At small wavenumber the positive branch has sound speed ${\lim }_{k\to 0}{\omega }_{+}(k)/k=1/\sqrt{2}$; at large wavenumber it approaches the free-particle quadratic dispersion. The linear frequencies are real, confirming stability of the repulsive uniform bulk.