nature is filled with fluid flow , from small scales to large scales .
our daily life on a human scale constantly experiences currents of air , water , and the like .
the earth sustains the ocean flows and atmospheric circulation , leading to an abundance of life . in space ,
stars and galaxies flow , forming many large - scale structures .
most flow in nature is actually turbulent .
as the velocity increases , flow generally changes from laminar to turbulent .
turbulence has been a mystery for a long time .
it has been investigated not only in basic science , including physics and mathematics research , but also in applied sciences , such as fluid engineering and aeronautics . as feynman said , turbulence is a problem `` that is common to many fields , that is very old , and that has not been solved . ''
@xcite this is chiefly because turbulence is a complicated dynamical phenomenon with strong nonlinearity far from an equilibrium state .
`` turbulence '' reminds one of some sketches by leonardo da vinci .
he observed the turbulent flow of water and drew pictures showing that turbulence has a structure comprised of vortices of different sizes ( fig .
vortices may be a key to understanding turbulence .
however , vortices are not well defined for a classical fluid , and the relationship between turbulence and vortices remains unclear .        independently of these studies in classical fluid dynamics , quantum fluids such as superfluid helium and atomic bose - einstein condensates ( becs ) have been investigated .
these systems are subject to quantum restrictions ; the appearance of order parameters makes the rotational motion exist only in the presence of quantized vortices .
a quantized vortex is a stable topological defect with quantized circulation .
such vortices give rise to quantum turbulence ( qt ) .
since quantized vortices are well defined as elements composing a turbulent flow , qt is expected to be an easier system to study than classical turbulence ( ct ) and has a simpler model of turbulence .
this article addresses recent developments in quantum fluid dynamics and qt in superfluid helium and atomic becs . in sec .
2 , basic concepts such as superfluidity , bose - einstein condensation , and quantized vortices are briefly reviewed .
section 3 is devoted to quantized vortices in a rotating bec , addressing typical phenomena of quantum hydrodynamics .
section 4 describes present important topics of qt in superfluid helium . applying these ideas , qt in atomic becs
is discussed in sec .
section 6 is devoted to conclusions .
quantum fluid dynamics has been studied in connection with superfluid helium for over 50 years .
liquid @xmath0he enters a superfluid state below the @xmath1 point ( 2.17 k ) through bose
einstein condensation of the atoms @xcite .
the characteristics of superfluidity were discovered experimentally in the 1930s by kapitza _
the hydrodynamics of superfluid helium is well described by the two - fluid model , in which the system consists of an inviscid superfluid ( of density @xmath2 ) and a viscous normal fluid ( of density @xmath3 ) with two independent velocity fields @xmath4 and @xmath5 .
the mixing ratio of the two fluids depends on temperature .
as the system is cooled down below the @xmath1 point , the ratio of the superfluid component increases , and the entire fluid becomes superfluid below about 1 k. the bose - condensed system exhibits the macroscopic wavefunction @xmath6 as an order parameter .
the superfluid velocity field representing the potential flow is given by @xmath7 with boson mass @xmath8 .
since the macroscopic wavefunction must be single - valued for space coordinate @xmath9 , the circulation @xmath10 for an arbitrary closed loop within the fluid is quantized in terms of the quantum value @xmath11 . a vortex with such quantized circulation
is called a quantized vortex . any rotational motion of a superfluid is only sustained by quantized vortices .
a quantized vortex is a topological defect characteristic of a bose
einstein condensate and is different from a vortex in a classical viscous fluid .
first , its circulation is quantized and conserved , in contrast to a classical vortex whose circulation can have any value and is not conserved .
second , a quantized vortex is a vortex of inviscid superflow .
thus , it can not decay by the viscous diffusion of vorticity that occurs in a classical fluid .
third , the core of a quantized vortex is very thin , on the order of the coherence length , which is only a few angstroms in superfluid @xmath0he and submicron in atomic becs . since the vortex core is thin and does not decay by diffusion , it is possible to identify the position of a quantized vortex in the fluid .
these properties make a quantized vortex more stable and definite than a classical vortex .
early experimental studies of superfluid turbulence focused primarily on thermal counterflow , in which the normal fluid and superfluid flow in opposite directions .
the flow is driven by an injected heat current , and it is found that the superflow becomes dissipative when the relative velocity between the two fluids exceeds a critical value @xcite .
feynman proposed that it is a superfluid turbulent state consisting of a tangle of quantized vortices @xcite .
vinen later confirmed feynman s proposal experimentally by showing that the dissipation arises from mutual friction between vortices and the normal flow @xcite .
the quantization of circulation was also observed by vinen using a vibrating wire technique @xcite .
subsequently , many experimental studies have examined superfluid turbulence ( st ) in thermal counterflow systems , revealing a variety of physical phenomena @xcite .
the dynamics of quantized vortices are nonlinear and nonlocal , so it has not been easy to quantitatively understand the experimental results on the basis of vortex dynamics .
a breakthrough was achieved by schwarz , who clarified the picture of st consisting of tangled vortices by a numerical simulation of the quantized vortex filament model in the thermal counterflow @xcite . however , since the thermal counterflow has no analogy in conventional fluid dynamics , these experimental and numerical studies are not helpful in clarifying the relationship between st and classical turbulence ( ct ) .
superfluid turbulence is often called quantum turbulence ( qt ) , which emphasizes the fact that it is comprised of quantized vortices .    comparing qt and
ct reveals definite differences .
turbulence in a classical viscous fluid appears to be comprised of vortices , as pointed out by da vinci .
however , these vortices are unstable , repeatedly appearing and disappearing .
moreover , the circulation is not conserved and is not identical for each vortex .
quantum turbulence consists of a tangle of quantized vortices that have the same conserved circulation .
looking back at the history of science , _ reductionism _ , which tries to understand the nature of complex things by reducing them to the interactions of their parts , has played an extremely important role .
the success of solid state physics owes much to _ reductionism_. in contrast , conventional fluid physics is not reducible to elements , and thus does not enjoy the benefits of _
reductionism_. however , quantum turbulence is different , being reduced to quantized vortices .
thus _ reductionism _ is applicable to quantum turbulence .
consequently , qt should lead to a simpler model of turbulence than ct .
based on these considerations , research into quantum fluid dynamics has opened up new directions since the mid 1990s .
one new direction has occurred in the field of low temperature physics by studying superfluid helium .
it started with the attempt to understand the relationship between qt and ct @xcite .
the energy spectrum of fully developed classical turbulence is known to obey the kolmogorov law in the inertial range .
recent experimental and numerical studies support a kolmogorov spectrum in qt . following these studies ,
qt research on superfluid helium has moved to important topics such as the dissipation process at very low temperatures , qt created by vibrating structures , and visualization of qt @xcite .
another new direction is the realization of bose - einstein condensation in trapped atomic gases in 1995 , which has stimulated intense experimental and theoretical activity @xcite .
as proof of the existence of superfluidity , quantized vortices have been created and observed in atomic becs , and numerous efforts have been devoted to a number of fascinating problems @xcite .
atomic becs have several advantages over superfluid helium .
the most important is that modern optical techniques enable one to directly control condensates and visualize quantized vortices .
because a series of experiments on becs clearly show the properties of quantum fluid dynamics @xcite , we shall start with this topic .
what happens if we rotate a cylindrical vessel with a classical viscous fluid inside ? even if the fluid is initially at rest , it starts to rotate and eventually reaches a steady rotation with the same rotational speed as the vessel . in that case
, one can say that the system contains a vortex that mimics solid - body rotation .
a rotation of arbitrary angular velocity can be sustained by a single vortex . however
, this does not occur in a quantum fluid .
because of quantization of circulation , superfluids respond to rotation , not with a single vortex , but with a lattice of quantized vortices .
feynman noted that in uniform rotation with angular velocity @xmath12 the curl of the superfluid velocity is the circulation per unit area , and since the curl is @xmath13 , a lattice of quantized vortices with number density @xmath14 ( `` feynman s rule '' ) arranges itself parallel to the rotation axis @xcite .
such experiments were performed for superfluid @xmath0he : packard _ et al .
_ visualized vortex lattices on the rotational axis by trapping electrons along the cores @xcite .
this idea has also been applied to atomic becs .
several groups have observed vortex lattices in rotating becs @xcite . among them , madison _ et al_. directly observed nonlinear processes such as vortex nucleation and lattice formation in a rotating @xmath15rb bec @xcite . by sudden application of a rotation along the trapping potential ,
an initially axisymmetric condensate undergoes a collective quadrupole oscillation to an elliptically deformed condensate .
this oscillation continues for a few hundred milliseconds with gradually decreasing amplitude .
then the axial symmetry of the condensate is recovered and vortices enter the condensate through its surface , eventually settling into a lattice configuration .
this observation has been reproduced by a simulation of the gross - pitaevskii ( gp ) equation for the macroscopic wavefunction @xmath6 in two - dimensional @xcite and three - dimensional @xcite spaces . the corresponding gp equation in a frame rotating with frequency @xmath16
is given by @xmath17 \psi({\bf x},t ) .
\label{gpeqrot}\ ] ] here @xmath18 is a trapping potential , and @xmath19 is the angular momentum along the rotational axis .
the interparticle potential @xmath20 is approximated by a short - range interaction @xmath21 , where @xmath22 is a coupling constant characterized by the s - wave scattering length @xmath23 .
the term @xmath24 indicates phenomenological dissipation . in this simulation @xcite . ]
figure 2 shows the profile of the condensate density @xmath25 and the phase @xmath26 when there is a quantized vortex in a trapped bec .
the density has a hole representing the vortex core .
the phase has a branch cut between 0 and @xmath27 , and the edge of the branch cut corresponds to the vortex core around which the phase rotates by @xmath27 as the superflow circulates .
one can therefore clearly identify the vortex both in the density and the phase .    a typical two - dimensional numerical simulation of eq . ( [ gpeqrot ] ) @xcite for the vortex lattice formation
is shown in fig .
3 , where the condensate density and the phase are displayed together . the trapping potential is @xmath28,\ ] ] where @xmath29 hz , and the parameters @xmath30 and @xmath31 describe small deviations from axisymmetry corresponding to experiments @xcite .
we first prepare an equilibrium condensate trapped in a stationary potential ; the size of the condensate cloud is determined by the thomas - fermi radius @xmath32 .
when we apply a rotation with @xmath33 , the condensate becomes elliptic and performs a quadrupole oscillation [ fig .
3(a ) ] . then
, the boundary surface of the condensate becomes unstable and generates ripples that propagate along the surface [ fig .
as stated previously , it is possible to identify quantized vortices in the phase profile also .
as soon as the rotation starts , many vortices appear in the low - density region outside of the condensate [ fig .
3(a ) ] . since quantized
vortices are excitations , their nucleation increases the energy of the system . because of the low density in the outskirts of the condensate , however , their nucleation contributes little to the energy and angular momentum . since these vortices outside of the condensate are not observed in the density profile , they are called `` ghost vortices ''
. their movement toward the thomas - fermi surface excites ripples [ fig .
it is not easy for these ghost vortices to enter the condensate , because that would increase both the energy and angular momentum .
only some vortices enter the condensate cloud to become `` real vortices '' wearing the usual density profile of quantized vortices [ fig .
3(d ) ] , eventually forming a vortex lattice [ fig .
3 ( e ) and ( f ) ] .
the number of vortices forming a lattice is given by `` feynman s rule '' @xmath34 .
the numerical results agree quantitatively with these observations .
note the essence of the dynamics .
the initial state has no vortices in the absence of rotation .
the final state is a vortex lattice corresponding to rotational frequency @xmath12 . in order to go from the initial to the final state
, the system makes use of as many excitations as possible , such as vortices , quadruple oscillation , and surface waves .
we refer the readers to ref .
26 for details .
these experimental and theoretical results demonstrate typical behavior of quantum fluid dynamics in atomic becs .
63 ms , ( b ) 305 ms , ( c ) 350 ms , ( d ) 410 ms , ( e ) 450 ms , and ( f ) 850 ms after the start of the rotation .
the phase varies continuously from 0 ( red ) to @xmath27 ( blue ) . ]
this section reviews research on qt since the mid 1990s . before considering qt , we begin with classical fluid dynamics and statistical properties of ct @xcite .
classical viscous fluid dynamics is described by the navier  stokes equation , @xmath35 where @xmath36 is the velocity of the fluid , @xmath37 is the pressure , @xmath38 is the density of the fluid , and @xmath39 is the kinematic viscosity .
the flow can be characterized by the ratio of the second term on the left - hand side of eq .
( [ eq - navier - stokes ] ) , referred to as the inertial term , and the second term on the right - hand side , called the viscous term .
this ratio defines the reynolds number @xmath40 , where @xmath41 and @xmath42 are a characteristic velocity and length scale , respectively .
when @xmath41 increases and the reynolds number exceeds a critical value , the system changes from a laminar to a turbulent state , in which the flow is complicated and contains eddies .
such turbulent flow shows characteristic statistical behavior @xcite .
we assume a steady state of fully developed turbulence of an incompressible classical fluid .
energy is injected into the fluid at a rate @xmath43 in an energy - containing range . in an inertial range , this energy is transferred to smaller length scales without dissipation . in this range ,
the system is locally homogeneous and isotropic , leading to energy spectral statistics described by the kolmogorov law , @xmath44 the energy spectrum @xmath45 is defined by @xmath46 , where @xmath47 is the kinetic energy per unit mass and @xmath48 is the wavenumber from the fourier transform of the velocity field .
the spectrum of eq .
( [ eq - kolmogorov ] ) is derived by assuming that @xmath45 is locally determined only by the energy flux @xmath43 and by @xmath48 .
the energy transferred to smaller scales is dissipated at the kolmogorov wavenumber @xmath49 in an energy - dissipative range via the viscosity of the fluid at a dissipation rate @xmath43 in eq .
( [ eq - kolmogorov ] ) , which is equal to the energy flux @xmath50 in the inertial range .
the kolmogorov constant @xmath51 is a dimensionless parameter of order unity .
the kolmogorov spectrum has been confirmed experimentally and numerically in turbulence at high reynolds numbers .
the inertial range is thought to be sustained by a self - similar richardson cascade in which large eddies break up into smaller eddies through vortex reconnections . in ct , however , the richardson cascade is not completely understood because it is impossible to definitely identify each eddy .      as described in sec . 2 ,
most experimental studies on superfluid turbulence have examined thermal counterflow .
however , the nonlinear and nonlocal dynamics of vortices delayed progress in the microscopic understanding of the vortex tangle .
schwarz overcame these difficulties @xcite by developing a direct numerical simulation of vortex dynamics connected with dynamical scaling analysis , enabling the calculation of physical quantities such as the vortex line density , anisotropic parameters , and mutual friction force .
the observable quantities calculated by schwarz agree with typical experimental results .
two formulations are generally available for studying the dynamics of quantized vortices .
one is the vortex filament model and the other is the gross  pitaevskii ( gp ) model .      as discussed in sec .
1 , a quantized vortex has quantized circulation .
the vortex core is extremely thin , usually much smaller than other characteristic length scales of the vortex motion .
these properties allow a quantized vortex to be represented as a vortex filament . in classical fluid dynamics ,
the vortex filament model is an idealization .
however , the vortex filament model is accurate and realistic for a quantized vortex in superfluid helium .
the vortex filament model represents a quantized vortex as a filament passing through the fluid , having a definite direction corresponding to its vorticity . except for the thin core region ,
the superflow velocity field has a classically well - defined meaning and can be described by ideal fluid dynamics @xcite .
the velocity at point @xmath52 due to a filament is given by the biot
savart expression , @xmath53 where @xmath54 is the quantum of circulation .
the vortex moves with the superfluid velocity . at finite temperatures ,
mutual friction occurs between the vortex core and the normal flow .
a numerical simulation method on this model has been described in detail elsewhere @xcite .
a vortex filament is represented by a single string of points separated by distance @xmath55 .
the vortex configuration at a given time determines the velocity field in the fluid , thus causing the vortex filaments to move .
vortex reconnection needs to be included when simulating vortex dynamics .
a numerical study of a classical fluid shows that close interaction of two vortices leads to their reconnection , chiefly because of the viscous diffusion of vorticity .
schwarz assumed that two vortex filaments reconnect when they come within a critical distance of each other and showed that statistical quantities such as the vortex line density are not sensitive to how the reconnections occur @xcite . even after schwarz s study , it remained unclear as to whether quantized vortices can actually reconnect .
however , koplik and levine directly solved the gp equation to show that two closely quantized vortices reconnect even in an inviscid superfluid @xcite .
more recent simulations have shown that reconnections are accompanied by the emission of sound waves having wavelengths on the order of the healing length @xcite .
schwarz numerically studied how remnant vortices develop into a statistical steady vortex tangle under thermal counterflow @xcite .
the tangle is self - sustained by the competition between the excitation due to the applied flow and the dissipation by mutual friction .
the numerical results were quantitatively consistent with experimental results , confirming feynman s original picture that qt is a tangle of quantized vortices .
the following two quantities are characteristic of a vortex tangle .
the line density @xmath56 is defined as the total length of vortex cores per unit volume .
the mean spacing @xmath57 between vortices is given by @xmath58 .      in a weakly interacting bose system
, the macroscopic wavefunction @xmath59 appears as the order parameter of the bose
einstein condensate , obeying the gross  pitaevskii ( gp ) equation @xcite , @xmath60 writing @xmath61 , the squared amplitude @xmath62 is the condensate density and the gradient of the phase @xmath63 gives the superfluid velocity @xmath64 , corresponding to frictionless flow of the condensate .
this relation causes quantized vortices to appear with quantized circulation .
the only characteristic scale of the gp model is the coherence length defined by @xmath65 , which determines the vortex core size .
the gp model can explain not only the vortex dynamics but also phenomena related to vortex cores , such as reconnection and nucleation .
however , the gp equation is not applicable to superfluid @xmath0he , which is not a weakly interacting bose system .
but it is applicable to bose
einstein condensation of a dilute atomic bose gas @xcite , as discussed in sec .
most older experimental studies on qt were devoted to thermal counterflow .
since this flow has no classical analog , these studies did not significantly contribute to the understanding of the relation between ct and qt .
since the mid 1990s , important experimental studies have been published on qt in the absence of thermal counterflow , differing significantly from previous studies @xcite .
the first important contribution was made by maurer and tabeling @xcite , who confirmed experimentally the kolmogorov spectrum in superfluid @xmath0he .
a turbulent flow was produced in a cylinder by driving two counter - rotating disks .
the authors measured the local pressure fluctuations to obtain the energy spectrum .
the experiments were made at three different temperatures : 2.3 k , 2.08 k , and 1.4 k. both above and below the @xmath1 point , the kolmogorov spectrum was confirmed .
next was a series of experiments on grid turbulence performed for superfluid @xmath0he above 1 k by the oregon group @xcite .
flow through a grid is usually used for generating turbulence in classical fluid dynamics @xcite . at sufficient distance behind the grid , the flow displays homogeneous isotropic turbulence .
this method has also been applied to superfluid helium . in the oregon experiments ,
the helium was contained in a channel along which a grid was pulled at constant velocity .
a pair of second - sound transducers was implanted in the walls of the channel to observe vortex tangles . in combining their observations with the decay of the turbulence , the authors had to make several assumptions .
the analysis is too complicated to be described in detail here .
the important point is that the coupling between the superfluid and the normal fluid by mutual friction makes a quasiclassical flow appear at length scales much larger than @xmath57 and causes the fluid to behave like a one - component fluid @xcite .
the line density is found to decay as @xmath66 .
a simple analysis shows that the @xmath66 decay can be derived from a quasiclassical model with a kolmogorov spectrum @xcite .
subsequently many experimental , theoretical , and numerical works appeared .
the present areas of study of qt are @xcite : the energy spectra and the dissipation mechanism at zero temperature ; qt created by vibrating structures ; and visualization of qt .
the pioneering experiments by maurer _
et al . _ and by the oregon group were performed at relatively high temperatures where the normal fluid component was significant , and so the situation is complicated because of the two - fluid behavior .
here we confine ourselves to the case of pure superfluid at zero temperature .
what happens to qt at zero temperature is not trivial @xcite .
the first problem is the nature of the energy spectrum of the turbulence for the pure superfluid component .
the second problem is the dissipation in this system .
since there is no normal fluid component , any dissipative mechanism does not work at large scales .
however , some dissipative mechanism should operate at small length scales .
one possibility is acoustic emission during vortex reconnections . in classical fluid dynamics ,
vortex reconnections cause acoustic emission . in quantum fluids , numerical simulations of the gp model
show acoustic emission at every reconnection event @xcite .
however , this mechanism is thought not to be important in superfluid helium because of its short coherence length .
another possible mechanism is the radiation of sound ( phonons ) by the oscillatory motion of the vortex cores .
the third problem is how energy is transferred from large to small scales at which the dissipative mechanism operates .
if the energy spectrum obeys the kolmogorov law , the energy should be transferred through a richardson cascade of quantized vortices in contrast to the classical case .
vortex reconnections sustaining a richardson cascade should be less effective at size scales shorter than the mean vortex spacing @xmath57 .
no experimental studies have directly addressed the energy spectra of qt at zero temperature , although three numerical studies have been made .
the first study was done by nore _
_ using the gp model @xcite .
they numerically solved the gp equation starting from taylor  green vortices and followed their time development .
the quantized vortices became tangled and the energy spectra of the incompressible kinetic energy appeared to obey the kolmogorov law for a short period , though eventually deviated from it .
the second study was done by the vortex filament model @xcite .
made a vortex tangle arising from taylor  green vortices and obtained an energy spectrum consistent with the kolmogorov law .
the third used the modified gp model by kobayashi and tsubota @xcite .
the kolmogorov spectra were confirmed for both decaying @xcite and steady @xcite qt by the modified gp model .
the normalized gp equation is @xmath67\phi(\mathbf{x},t),\label{eq - gp}\ ] ] which determines the dynamics of the macroscopic wavefunction @xmath68 $ ] .
the hydrodynamics in the gp model are compressible .
the total number of condensate particles is @xmath69 and the total energy is @xmath70\phi(\mathbf{x},t),\ ] ] as represented by the sum of the interaction energy @xmath71 , the quantum energy @xmath72 , and the kinetic energy @xmath73 @xcite , @xmath74 ^ 2 , e_{kin}(t)=\frac{1}{n}\int d\mathbf{x}\:[f(\mathbf{x},t)\nabla\phi(\mathbf{x},t)]^2.\label{eq - energy - definition } \end{aligned}\ ] ] the kinetic energy is further divided into a compressible part @xmath75 due to compressible excitations and an incompressible part @xmath76 due to vortices .
the kolmogorov spectrum is expected for @xmath76 .
the failure to obtain a kolmogorov law in the pure gp model @xcite is attributable to the following .
the simulations showed that @xmath76 decreases and @xmath75 increases while the total energy @xmath77 is conserved because many compressible excitations are created through vortex reconnections @xcite and disturb the richardson cascade of quantized vortices .
_ overcame these difficulties and obtained a kolmogorov spectrum in qt that revealed an energy cascade @xcite .
they made numerical calculations of the fourier - transformed gp equation with dissipation , @xmath78\tilde{\phi}(\mathbf{k},t ) + \frac{g}{v^2}\sum_{\mathbf{k}_1,\mathbf{k}_2}\tilde{\phi}(\mathbf{k}_1,t)\tilde{\phi}^\ast(\mathbf{k}_2,t ) \times\tilde{\phi}(\mathbf{k}-\mathbf{k}_1+\mathbf{k}_2,t ) .
\label{eq - fourier - gp1 } \end{aligned}\ ] ] here @xmath79 is the spatial fourier component of @xmath80 and @xmath81 is the system volume .
the healing length is @xmath82 .
the dissipation has the form @xmath83 with step function @xmath63 , which only dissipates excitations smaller than @xmath84 .
this form of dissipation can be justified by a coupled analysis of the gp equation and the bogoliubov  de gennes equations for thermal excitations @xcite .
kobayashi _ et al .
_ confirmed the kolmogorov spectra for decaying turbulence @xcite . to obtain a turbulent state
, they started the calculation from an initial configuration in which the density was uniform and the phase of the wavefunction had a random spatial distribution .
the initial state was dynamically unstable and soon developed turbulence with many vortex loops .
the spectrum @xmath85 was then found to obey the kolmogorov law .
a more elaborate analysis was made for steady qt by introducing energy injection at large length scales in addition to energy dissipation at small scales @xcite .
energy injection at large scales was effected by moving a random potential @xmath86 that had a characteristic spatial scale @xmath87 .
once steady qt occurs as a balance between energy injection and dissipation , it has an energy - containing range @xmath88 , an inertial range @xmath89 , and an energy - dissipative range @xmath90 .
a typical simulation of steady turbulence was performed for @xmath91 .
the dynamics started from the uniform wavefunction .
the moving random potential nucleates sound waves as well as vortices , but the incompressible kinetic energy @xmath76 due to vortices becomes dominant in the total kinetic energy @xmath73 .
the energies soon become almost constant , and steady qt is obtained .
the energy cascade can be confirmed quantitatively by checking that the energy dissipation rate @xmath43 of @xmath76 is comparable to the energy flux @xmath92 via a richardson cascade in the inertial range .
the energy flux is found to be approximately independent of @xmath48 and comparable to @xmath43 . as shown in fig .
[ fig - steady - kolmogorov ] ( b ) , the energy spectrum is quantitatively consistent with the kolmogorov law in the inertial range @xmath89 , equivalent to @xmath93 .     for qt .
the plotted points are from an ensemble average of 50 randomly selected steady states .
the solid line is the kolmogorov law @xcite.,title="fig : " ] + ( a )     for qt .
the plotted points are from an ensemble average of 50 randomly selected steady states .
the solid line is the kolmogorov law @xcite.,title="fig : " ] + ( b )    the above arguments are chiefly limited to length scales larger than the mean spacing @xmath57 of a vortex tangle , in which the richardson cascade is effective .
what happens at smaller scales where the richardson cascade is less effective ?
the most probable scenario is the kelvin - wave cascade .
a kelvin wave is a deformation of a vortex line into a helix , with the deformation propagating as a wave along the vortex @xcite .
kelvin waves were first observed by making torsional oscillations in uniformly rotating superfluid @xmath0he @xcite .
the approximate dispersion relation for a rectilinear vortex is @xmath94 with a constant @xmath95 and a cutoff parameter @xmath96 equal to the vortex core size . here
@xmath48 is the wavenumber of an excited kelvin wave , which is different than the wavenumber used for the energy spectrum above . at a finite temperature ,
a significant fraction of the normal fluid contains damped kelvin waves because of mutual friction .
at very low temperatures , however , mutual friction does not occur and the only possible mechanism of dissipation is the radiation of phonons @xcite .
phonon radiation becomes effective only when the frequency is very high , typically of the order of ghz ( @xmath97 nm@xmath98 ) .
therefore , a mechanism is required to transfer energy to such high wavenumbers for kelvin - wave damping . an early numerical simulation based on the vortex filament model showed that kelvin waves are unstable to the buildup of sidebands @xcite .
consequently nonlinear interactions between different kelvin wavenumbers can transfer energy from small to large wavenumbers , which constitutes a kelvin - wave cascade .
this idea was first suggested by svistunov @xcite and later developed and confirmed through theoretical and numerical work by kivotides _
@xcite , vinen _ et al . _
@xcite , and kozik and svistunov @xcite .
consider the nature of the transition between the richardson and the kelvin - wave cascades .
turbulence at scales larger than the mean vortex spacing obeys the kolmogorov law with a richardson cascade . in this regime , the nature of individual vortices is not relevant , and the large - scale velocity field created collectively by the vortices mimics classical turbulent behavior .
this regime may therefore be called classical . on the other hand , at length scales smaller than the mean vortex spacing the energy
is transferred along each vortex through a kelvin - wave cascade ; this regime has no classical analog and can be called quantum .
several theoretical considerations of the classical - quantum crossover have been proposed @xcite , but this topic is not yet settled and remains controversial .
recently , vibrating structures including discs , spheres , grids , and wires , have been used for research into qt @xcite . despite differences between the structures ,
the experiments show surprisingly similar behavior . when a classical viscous fluid steadily flows past a structure , the flow changes from laminar to turbulence in the wake as the reynolds number increases .
the net drag force on the structure is often expressed in terms of the drag coefficient @xmath99 as @xmath100 , where @xmath101 is the projected area of the structure onto a plane normal to the flow .
for laminar viscous flow , the drag is approximately proportional to @xmath102 so that @xmath103 , while turbulent flow with high reynolds number yields a drag coefficient @xmath99 of order unity .
this drastic reduction in the drag force is called the drag crisis , and is closely related to the formation of the boundary layer @xcite .
a boundary layer in a high velocity gradient produces eddies . however , such boundary layers do not appear in the superfluid component .
another classical aspect is the effect of oscillations .
the oscillatory case is more complicated than that of steady flow because a second length scale appears in addition to the linear size of the structure , namely the viscous penetration depth @xmath104 where @xmath105 is the oscillation frequency of the structure @xcite .
owing to the importance of the topic in marine engineering , a number of experimental studies have been published on oscillatory classical flow past a cylinder .
as laminar flow becomes unstable , the drag coefficient @xmath99 reduces to a value of order unity . however , the values of @xmath99 oscillate within the range 0.5 to 2.0 , showing a more complicated transition to turbulence than in the case of steady flow .
visual observations show that vortices are formed in the wake of the oscillating cylinder .
a few experimental studies have reported classical flow around an oscillating sphere .
a typical experiment shows that a single vortex ring is generated and shed during each half period of oscillation @xcite .
typical behavior appeared in the pioneering observation of qt for an oscillating microsphere in superfluid @xmath0he by jger _
the sphere they used had a radius of approximately 100 @xmath106 m .
it was made from a strongly ferromagnetic material and was magnetically levitated in superfluid @xmath0he .
its response to an alternating drive was observed . at low drives ,
the velocity response @xmath102 was proportional to the drive @xmath107 , taking the `` laminar '' form @xmath108 with a temperature - dependent coefficient @xmath109 . at high drives ,
the response changed to the `` turbulent '' form @xmath110 above a critical velocity @xmath111 . at low temperatures , the transition from laminar to turbulent response
was accompanied by hysteresis .
subsequently , several groups have experimentally investigated the transition to turbulence in superfluid @xmath0he and @xmath112he - b by using grids , wires , and tuning forks @xcite .
these experimental studies reported some common behaviors independent of the details of the structures , including its type , shape , and surface roughness .
the observed critical velocities were in the range from 1 to 200 mm / s .
since the velocity was lower than the landau critical velocity of approximately 50 m / s , the transition to turbulence should come not from intrinsic nucleation of vortices but from extension or amplification of remnant vortices.he at rest usually contains remnant vortices @xcite .
they are metastably pinned between the walls of the container and are established when the liquid is cooled down through the @xmath1 point to the superfluid phase . ]
such behavior is seen in the numerical simulation using the vortex filament model @xcite .
figure 5 shows how the remnant vortices initially attached to a sphere develop into turbulence under an oscillating flow .
many unresolved problems remain , such as the nature of the critical velocity and the origin of the hysteresis in the transition between the laminar and turbulent response .
m in an oscillating superflow of 150 mm / s at 200  hz .
[ hnninen , tsubota , and vinen , phys . rev .
b * 75 * ( 2007 ) 064502 , reproduced with permission .
copyright 2007 by the american physical society . ] ]      there has been little direct experimental information about flow in superfluid @xmath0he because standard flow visualization techniques are not applicable to it
. however , the situation is changing @xcite .
for qt , one can seed the fluid with tracer particles in order to visualize the flow field .
quantized vortices are observable by appropriate optical techniques .
a significant contribution was made by zhang and van sciver @xcite . using the particle image velocimetry ( piv ) technique with 1.7-@xmath106m - diameter polymer particles , they visualized a large - scale turbulent flow both in front of and behind a cylinder in counterflowing superfluid @xmath0he at finite temperatures . in classical fluids ,
such turbulent structures are seen downstream of objects such as cylinders , and the structures periodically detach to form a vortex sheet . in the case of @xmath0he counterflow , however ,
the locations of the large - scale turbulent structures were relatively stable ; they did not detach and move downstream , although local fluctuations in the turbulence were evident .
do the tracer particles follow the normal flow , the superflow , or a more complex flow ?
_ studied this question theoretically and numerically to show that the situation changes depending on the size and mass of the tracer particles @xcite .
another important contribution was the visualization of quantized vortices by bewley _
_ @xcite . in their experiments ,
liquid helium was seeded with solid hydrogen particles smaller than 2.7 @xmath106 m in diameter at a temperature slightly above @xmath113 , following which the fluid was cooled to below @xmath113 . when the temperature is above @xmath113 , a homogeneous cloud of the particles was seen that disperses throughout the fluid . however , on passing through @xmath113 , the particles coalesced into web - like structures .
_ suggested that these structures represent decorated quantized vortex lines .
they reported that the vortex lines form connections rather than remaining separated , and were homogeneously distributed throughout the fluid .
the observed fork - like structures may indicate that several vortices are attached to the same particle . applying this technique to thermal counter flow , paoletti
_ have directly visualized the two - fluid behavior @xcite .
the long history of research into superfluid helium has uncovered two main cooperative phenomena of quantized vortices : vortex lattices under rotation and vortex tangles in qt .
however , almost all studies of quantized vortices in atomic becs have been limited to vortex lattices .
this section briefly discusses qt in atomic becs .    kobayashi and tsubota proposed an easy , powerful method to make a steady qt in a trapped bec , by using precession @xcite .
the dynamics of the wavefunction are described by the gp equation with dissipation .
first , one traps a bec cloud in a weakly elliptical harmonic potential , @xmath114,\ ] ] where the parameters @xmath115 and @xmath116 exhibit elliptical deformation in the @xmath117- and @xmath118-planes .
second , to transform the bec into a turbulent state , a rotation is applied along the @xmath119-axis followed by a rotation along the @xmath120-axis , as shown in fig .
the rotation vector is @xmath121 where @xmath122 and @xmath123 are the frequencies of the first and second rotations , respectively . consider the case where the spinning and precessing rotational axes are perpendicular to each other .
then the two rotations do not commute and can not be represented by their sum .
this form of precession is also used for turbulence in water @xcite .
starting from a stationary solution without rotation or elliptical deformation , kobayashi _ et al . _
numerically calculated the time development of the gp equation by turning on the rotation @xmath124 and elliptical deformation @xmath125 . by monitoring the kinetic energy and anisotropic parameters ,
the system eventually becomes statistically steady and isotropic . in the steady state ,
the spectrum @xmath85 of the incompressible kinetic energy per unit mass is consistent with the kolmogorov law [ fig .
the inertial range which sustains the kolmogorov law is determined by the thomas - fermi radius @xmath32 and the coherence length @xmath84 .
the application of a combined precession around three axes enables one to obtain more isotropic qt @xcite .
the velocity field in a bec cloud can be observed by bragg spectroscopy @xcite , enabling one to obtain the energy spectrum .
is the thomas - fermi radius amd @xmath126 is the characteristic length scale of the trap . ]
there are several advantages of studying qt in atomic becs compared to in superfluid helium .
first , one can observe the vortex configuration , probably even the richardson cascade process of vortices .
in fact one can directly study the relation between the real - space richardson cascade and the wavenumber - space cascade ( the kolmogorov spectrum ) .
second , one can control the transition to turbulence by changing the rotational frequencies or other parameters .
for example , rotation along one axis forms a vortex lattice .
if one applies another rotation , it may rotate the lattice if the frequency is low . if the frequency is high , however , the second rotation can melt the lattice into a vortex tangle .
that would enable one to investigate in detail the entire transition to turbulence .
third , by changing the shape of the trapping potential , one can study the effect of the anisotropy on turbulence ; a typical question is how the kolmogorov spectrum is changed when the bec becomes anisotropic .
this issue is related to the transition between 2d and 3d turbulence .
fourth , one can create rich qt in multicomponent becs by controlling the interaction parameters .
multicomponent becs have been extensively studied and found to create many kinds of vortex structures depending on the intercomponent interaction @xcite . by using such systems
, we will be able to make multicomponent qt coupled through the interaction .
in this article , we have reviewed recent research on quantum hydrodynamics and qt in superfluid helium and atomic becs .
research on qt is currently a key branch in low - temperature physics , attracting the attention of many scientists .
it is comprised of quantized vortices as definite elements , in contrast to conventional turbulence .
thus , the investigation of qt may lead to a breakthrough in understanding a mystery of nature pondered since the era of da vinci .
there are many related topics not addressed in this article , regarding which we refer the readers to other review articles @xcite .
the author thanks kenichi kasamatsu and michikazu kobayashi for help with the manuscript .
the author acknowledges a grant - in - aid for scientific research from jsps ( grant no .
18340109 ) and a grant - in - aid for scientific research on priority areas from mext ( grant no .
17071008 ) .
r. p. feynman , r. b. leighton , and m. sands , _ the feynman lectures on physics _ , addison - wesley , san francisco , 2006 .
d. r. tilley and j. tilley , _ superfluidity and superconductivity _ , 3rd ed .
, institute of physics publishing , bristol and philadelphia , 1990 . c. j. gorter and j. h. mellink , _ on the irreversible processes in liquid helium ii _ ,
physica 15 ( 1949 ) , pp .
r. p. feynman , _ application of quantum mechanics to liquid helium _ , in _ prog .
c. j. gorter , north - holland , amsterdam , 1955 , vol .
i , pp . 1753 .
w. f. vinen , _ mutual friction in a heat current in liquid helium ii : i. experiments on steady heat currents _ , proc .
soc . a 240 ( 1957 ) ,
. w. f. vinen , _ mutual friction in a heat current in liquid helium ii : ii .
experiments on transient effects _ , proc .
soc . a 240 ( 1957 ) ,
. w. f. vinen , _ mutual friction in a heat current in liquid helium ii : iii .
theory of the mutual friction _ , proc .
soc . a 242 ( 1957 ) ,
. w. f. vinen , _ mutual friction in a heat current in liquid helium ii : iv .
critical heat currents in wide channels _ , proc .
soc . a 243 ( 1958 ) ,
. w. f. vinen , _ the detection of a single quantum of circulation in liquid helium ii _ , proc .
soc . a 260 ( 1958 ) ,
400413 . j. t. tough , _ superfluid turbulence _
, in _ prog
. low temp .
_ , ed . c. j. gorter , north - holland , amsterdam , 1982 , vol
viii , pp .
. k. w. schwarz , _ three - dimensional vortex dynamics in superfluid 4he : line - line and line - boundary interactions _ , phys .
b 31 ( 1985 ) , pp . 57825803 .
k. w. schwarz , _ three - dimensional vortex dynamics in superfluid 4he : homogenous superfluid turbulence _ , phys . rev .
b 38 ( 1988 ) , pp . 23982417 .
w. f. vinen and j. j. niemela , _ quantum turbulence _ , j. low temp .
( 2002 ) , pp .
. w. f. vinen , _ an introduction to quantum turbulence _ , j. low temp .
145 ( 2006 ) , pp .
, eds . w. p. halperin and m. tsubota , elsevier , amsterdam , 2008 , vol .
m. tsubota , _ quantum turbulence _ , j. phys .
77 ( 2008 ) , 111006 . c. j. pethick and h. smith , _ bose - einstein condensation in dilute gases _ , cambridge university press , cambridge , 2002 .
k. kasamatsu and m. tsubota , _ quantized vortices in atomic bose - einstein condensates _ , in _ prog .
w. p. halperin and m. tsubota , elsevier , amsterdam , 2008 , vol .
16 , pp . 351403
. j. r. abo - shaeer , c. raman , j. m. vogels , and w. ketterle , _ observation of vortex lattices in bose - einstein condensates _ ,
science 292 ( 2001 ) , pp . 476479 .
k. w. madison , f. chevy , w. wlhlleben , and j. dalibard , _ vortex formation in a stirred bose - einstein condensate _ ,
lett . 84 ( 2000 ) , pp .
. k. w. madison , f. chevy , w. wlhlleben , and j. dalibard , _ statonary states of a rotating bose - einstein condensate : routes to vortex nucleation _ , phys .
86 ( 2001 ) , pp .
44434446 . g. a. williams and r. e. packard , _ photographs of quantized vortex lines in rotating he ii _
33 ( 1974 ) , pp .
e. j. yarmchuck and r. e. packard , _ photographic studies of quantized vortex lines _ , j. low temp
46 ( 1982 ) , pp . 479515
. m. r. matthews , b. p. anderson , p. c. haljan , d. s. hall , c. e. wieman , and e. a. cornell , _ vortices in a bose - einstein condensate _ ,
83 ( 1999 ) , pp .
m. tsubota , k. kasamatsu , and m. ueda , _ vortex lattice formation in a rotating bose - einstein condensate _ ,
a65 ( 2002 ) , 023603 .
k. kasamatsu , m. tsubota , and m. ueda , _ nonlinear dynamics of vortex lattice formation in a rotating bose - einstein condensate _ , phys .
a67 ( 2003 ) , 033610 .
k. kasamatsu , m. machida , n. sasa and m. tsubota , _ three - dimensional dynamics of vortex - lattice formation in bose - einstein condensate _ phys .
a71 ( 2005 ) , 063616 .
u. frisch , _ turbulence _ , cambridge university press , cambridge , 1995 .
a. n. kolmogorov , _ the local structure of turbulence in incompressible viscous fluid for very large reynolds number _ , dokl .
nauk sssr 30 ( 1941 ) , pp .
299 - 303 [ reprinted in proc .
soc . a 434 ( 1991 ) ,
a. n. kolmogorov , _ on degeneration ( decay ) of isotropic turbulence in an incompressible viscous liquid _ ,
nauk sssr 31 ( 1941 ) , pp .
538 - 540 [ reprinted in proc .
a 434 ( 1991 ) , pp .
15 - 17 ] .
m. tsubota , t. araki , and s. k. nemirovskii , _ dynamics of vortex tangle without mutual friction in superfluid @xmath0he _ , phys . rev .
b 62 ( 2000 ) , pp .
11751 - 11762
. j. koplik and h. levine , _ vortex reconnection in superfluid helium _ , phys .
71 ( 1993 ) , pp .
1375 - 1378 .
m. leadbeater , t. winiecki , d. c. samuels , c. f. barenghi , and c. s. adams , _ sound emission due to superfluid vortex reconnections _ ,
86 ( 2001 ) , pp .
1410 - 1413 .
s. ogawa , m. tsubota , and y. hattori , _ study of reconnection and acoustic emission of quantized vortices in superfluid by the numerical analysis of the gross - pitaevskii equation _
, j. phys .
soc . jpn .
71 ( 2002 ) , pp .
813 - 821 .
j. maurer and p. tabeling , _ local investigation of superfluid turbulence _
, europhys . lett .
43 ( 1998 ) , pp .
s. r. stalp , l. skrbek , and r. j. donnelly , _ decay of grid turbulence in a finite channel _ , phys .
82 ( 1999 ) , pp .
4831 - 4834 .
l. skrbek , j. j. niemela , and r. j. donnelly , _ four regimes of decaying grid turbulence in a finite channel _ , phys .
85 ( 2000 ) , pp .
2973 - 2976 . l. skrbek and s. r. stalp , _ on the decay of homogeneous isotropic turbulence _ , phys .
fluids 12 ( 1997 ) , pp .
1997 - 2019 .
s. r. stalp , j. j. niemela , w. f. vinen , and r. j. donnelly , _ dissipation of grid turbulence in helium ii _ ,
fluids 14 ( 2002 ) , pp .
1377 - 1379
. w. f. vinen , _ classical character of turbulence in a quantum liquid
b 61 ( 2000 ) , pp . 1410 - 1420 . c. nore , m. abid , and m. e. brachet , _ kolmogorov turbulence in low - temperature superflows _ , phys . rev .
78 ( 1997 ) , pp .
3296 - 3299 . c. nore , m. abid , and m. e. brachet , _ decaying kolmogorov turbulence in a model of superflow _ , phys .
fluids 9 ( 1997 ) , pp .
2644 - 2669
. t. araki , m. tsubota and s. k. nemirovskii , _ energy spectrum of superfluid turbulence with no normal - fluid component _ , phys .
( 2002 ) , 145301 .
m. kobayashi and m. tsubota , _ kolmogorov spectrum of superfluid turbulence : numerical analysis of the gross - pitaevskii equation with a small - scale dissipation _
94 ( 2005 ) , 065302 .
m. kobayashi and m. tsubota , _ kolmogorov spectrum of quantum turbulence _
, j. phys .
74 ( 2005 ) , pp .
3248 - 3258 . m. kobayashi and m. tsubota , _ thermal dissipation in quantum turbulence _ ,
97 ( 2006 ) , 145301 .
j. thomson , _ vibrations of a columnar vortex _ , phil mag .
10 ( 1880 ) , pp .
155 - 168 .
h. e. hall , _ an experimental and theoretical study of torsional oscillations in uniformly rotating liquid helium ii _ , proc .
london a 245 ( 1958 ) , pp .
. h. e. hall , _ the rotation of helium ii _ , phil . mag .
suppl . 9 ( 1960 ) , pp .
w. f. vinen , _ decay of superfluid turbulence at a very low temperature : the radiation of sound from a kelvin wave on a quantized vortex _ ,
b 64 ( 2001 ) , 134520 .
d. c. samuels and r. j. donnelly , _ sideband instability and recurrence of kelvin waves on vortex cores _ , phys .
64(1990 ) , pp . 1385 - 1388
. b. v. svistunov , _ superfluid turbulence in the low - temperature limit , _ phys .
b 52 ( 1995 ) , pp .
3647 - 3653 .
d. kivotides , j. c. vassilicos , d. c. samuels , and c. f. barenghi , _ kelvin waves cascade in superfluid turbulence _
86 ( 2001 ) , pp .
3080 - 3083
. w. f. vinen , m. tsubota , and m. mitani , _ kelvin - wave cascade on a vortex in superfluid @xmath0he at a very low temperature _ , phys .
91 ( 2003 ) , 135301 .
e. kozik , and b. v. svistunov , _ kelvin - wave cascade and decay of superfluid turbulence _ ,
92 ( 2004 ) , 035301 .
e. kozik and b. v. svistunov , _ scale - separation scheme for simulating superfluid turbulence : kelvin - wave cascade _ , phys .
94 ( 2005 ) , 025301 .
e. kozik and b. v. svistunov , _ vortex - phonon interaction _ phys .
b 72 ( 2005 ) , 172505 . v.
s. lvov , s. v. nazarenko , and o. rudenko , _ bottleneck crossover between classical and quantum superfluid turbulence _ ,
b 76 ( 2007 ) , 024520 . e. kozik and b. v. svistunov , _ kolmogorov and kelvin - wave cascades of superfluid turbulence at t=0 : what lies between _ , phys .
b 77 ( 2008 ) , 060502 . e. kozik and b. v. svistunov , _ scanning superfluid - turbulence cascade by its low - temperature cutoff _ , phys .
( 2008 ) , 195302 .
l. skrbek and w. f. vinen , _ the use of vibrating structures in the study of quantum turbulence _ , in _ prog
. low temp .
w. p. halperin and m. tsubota , elsevier , amsterdam , 2008 , vol .
195 - 246 .
l. d. landau and e. m. lifshitz , _ fluid mechanics _ , 2nd ed . , pergamon press , oxford , 1987 .
r. j. donnelly , private communication .
j. jger , b. schuderer , and w. schoepe , _ turbulent and laminar drag of superfluid helium on an oscillating microsphere _ ,
74 ( 1995 ) , pp .
566 - 569 .
d. d. awschalom and k. w. schwarz , _ observation of a remanent vorte - line density in superfluid helium _ , phys .
52 ( 1984 ) , pp .
r. hnninen , m. tsubota , and w. f. vinen , _ generation of turbulence by oscillating structures in superfluid helium at very low temperatures _ ,
b 75 ( 2007 ) , 064502 .
s. w. van sciver and c. f. barenghi , _ visualization of quantum turbulence _ , in _ prog
. low temp .
w. p. halperin and m. tsubota , elsevier , amsterdam , 2008 , vol .
247 - 303 .
t. zhang and s. w. van sciver , _ large - scale turbulent flow around a cylinder in counterflow superfluid 4he ( he(ii ) ) _ , nature phys . 1 ( 2005 ) , pp .
d. r. poole , c. f. barenghi , y. a. sergeev , and w. f. vinen , _ motion of tracer particles in he ii _ , phys .
b 71 ( 2005 ) , 064514 .
d. kivotides , c. f. barenghi , and y. a. sergeev , _ collision of a tracer particle and a quantized vortex in superfluid helium : self - consistent calculations _
b 75 ( 2007 ) , 212502 .
d. kivotides , c. f. barenghi , and y. a. sergeev , _ interactions between particles and quantized vortices in superfluid helium _ , phys .
b 77 ( 2008 ) , 014527 .
g. p. bewley , d. p. lathrop , and k. r. sreenivasan , _ visualization of quantized vortices _ ,
nature 441 ( 2006 ) , 588 .
m. s. paelotti , r.b .
firrito , k. r. sreenivasan , and d. p. lathrop .
_ visualization of superfluid helium flow _ , j. phys .
77 ( 2008 ) , 111007 .
m. kobayashi and m. tsubota , _ quantum turbulence in a trapped bose - einstein condensate _
a76 ( 2007 ) , 045603 . s. goto , n. ishii , s. kida , and m. nishioka , _ turbulence generator using a precessing sphere _
fluids 19 ( 2007 ) , 061705 .
m. kobayashi and m. tsubota , _ quantum turbulence in a trapped bose - einstein condensate under combined rotations around three axes _
, j. low temp .
( 2008 ) , pp .
587 - 592 . j. stenger , s. inouye , a. p. chikkatur , d. m. stamper - kurn , d. e. pritchard , and w. ketterle , _ bragg spectroscopy of a bose - einstein condensate _
82 ( 1999 ) , pp .
4569 - 4573 . k. kasamatsu , m. tsubota , and m. ueda , _ vortices in multicomponent bose - einstein condensates _ , int .
b 19 ( 2005 ) , pp . 1835 - 1904 .