Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width

The evolution of a pulse of noninteracting quasiparticles, caused by their different velocities and angular distribution of momenta, is studied theoretically. Equations are found that describe the shape of the pulse surface at any time. The time of the beginning, end and duration of the density o...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2004
Автори: Adamenko, I.N., Nemchenko, K.E., Slipko, V.A., Wyatt, A.F.G.
Формат: Стаття
Мова:English
Опубліковано: Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України 2004
Назва видання:Физика низких температур
Теми:
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/119730
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width / I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt // Физика низких температур. — 2004. — Т. 30, № 6. — С. 579-589. — Бібліогр.: 20 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-119730
record_format dspace
spelling irk-123456789-1197302017-06-09T03:04:28Z Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width Adamenko, I.N. Nemchenko, K.E. Slipko, V.A. Wyatt, A.F.G. Квантовые жидкости и квантовые кpисталлы The evolution of a pulse of noninteracting quasiparticles, caused by their different velocities and angular distribution of momenta, is studied theoretically. Equations are found that describe the shape of the pulse surface at any time. The time of the beginning, end and duration of the density of the quasiparticle energy flux is determined at a general spatial point. The quasiparticle energy density is considered at all times and positions, and it is shown that the region of high energy density, in the middle of the pulse, is equal to the initial energy density under certain conditions. These theoretical results are discussed in relation to experimental data on the evolution of a pulse of noninteracting phonons in superfluid helium. 2004 Article Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width / I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt // Физика низких температур. — 2004. — Т. 30, № 6. — С. 579-589. — Бібліогр.: 20 назв. — англ. 0132-6414 PACS: 67.40.Fd http://dspace.nbuv.gov.ua/handle/123456789/119730 en Физика низких температур Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
topic Квантовые жидкости и квантовые кpисталлы
Квантовые жидкости и квантовые кpисталлы
spellingShingle Квантовые жидкости и квантовые кpисталлы
Квантовые жидкости и квантовые кpисталлы
Adamenko, I.N.
Nemchenko, K.E.
Slipko, V.A.
Wyatt, A.F.G.
Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
Физика низких температур
description The evolution of a pulse of noninteracting quasiparticles, caused by their different velocities and angular distribution of momenta, is studied theoretically. Equations are found that describe the shape of the pulse surface at any time. The time of the beginning, end and duration of the density of the quasiparticle energy flux is determined at a general spatial point. The quasiparticle energy density is considered at all times and positions, and it is shown that the region of high energy density, in the middle of the pulse, is equal to the initial energy density under certain conditions. These theoretical results are discussed in relation to experimental data on the evolution of a pulse of noninteracting phonons in superfluid helium.
format Article
author Adamenko, I.N.
Nemchenko, K.E.
Slipko, V.A.
Wyatt, A.F.G.
author_facet Adamenko, I.N.
Nemchenko, K.E.
Slipko, V.A.
Wyatt, A.F.G.
author_sort Adamenko, I.N.
title Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
title_short Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
title_full Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
title_fullStr Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
title_full_unstemmed Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
title_sort evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width
publisher Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
publishDate 2004
topic_facet Квантовые жидкости и квантовые кpисталлы
url http://dspace.nbuv.gov.ua/handle/123456789/119730
citation_txt Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width / I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt // Физика низких температур. — 2004. — Т. 30, № 6. — С. 579-589. — Бібліогр.: 20 назв. — англ.
series Физика низких температур
work_keys_str_mv AT adamenkoin evolutionofapulseofnoninteractingquasiparticleswithdispersionandinitialangularwidth
AT nemchenkoke evolutionofapulseofnoninteractingquasiparticleswithdispersionandinitialangularwidth
AT slipkova evolutionofapulseofnoninteractingquasiparticleswithdispersionandinitialangularwidth
AT wyattafg evolutionofapulseofnoninteractingquasiparticleswithdispersionandinitialangularwidth
first_indexed 2025-07-08T16:29:57Z
last_indexed 2025-07-08T16:29:57Z
_version_ 1837096971573133312
fulltext Fizika Nizkikh Temperatur, 2004, v. 30, No. 6, p. 579–589 Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width I.N. Adamenko1,2, K.E. Nemchenko1, V.A. Slipko1, and A.F.G. Wyatt3 1Karazin Kharkov National University, 4 Svobody Sq., Kharkov 61077, Ukraine E-mail: nemchenko@pht.univer.kharkov.ua 2Scientific and Technological Center of Electrophysics, NAN of Ukraine 28 Chernyshevskyi Str., P.O. Box 8812, Kharkov 61002, Ukraine 3School of Physics, University of Exeter, Exeter EX4 4QL, UK Received December 18, 2003 The evolution of a pulse of noninteracting quasiparticles, caused by their different velocities and angular distribution of momenta, is studied theoretically. Equations are found that describe the shape of the pulse surface at any time. The time of the beginning, end and duration of the den- sity of the quasiparticle energy flux is determined at a general spatial point. The quasiparticle en- ergy density is considered at all times and positions, and it is shown that the region of high energy density, in the middle of the pulse, is equal to the initial energy density under certain conditions. These theoretical results are discussed in relation to experimental data on the evolution of a pulse of noninteracting phonons in superfluid helium. PACS: 67.40.Fd 1. Introduction Quasiparticle pulses are studied theoretically and experimentally in various fields of modern physics. The variety of phenomena that appear are due to a number of factors, the most important of which is the role played by following: 1. The dependence of the quasiparticle energy on momentum (the dispersion or energy—momentum rela- tion). Classical particles can be considered as quasi- particles with a quadratic energy—momentum relation. 2. The angular width of the pulse, which is deter- mined by the angular range of the quasiparticles’ momenta. 3. Interactions between quasiparticles. 4. The interaction of quasiparticles with the me- dium in which the pulse moves. The above-mentioned factors determine, for exam- ple, phenomena observed in photon pulses propagat- ing in an optically transparent medium with variable refraction index, e.g., optical waveguides, or in a medium where the refractive index depends on photon energy. Another interesting property of anisotropic crystals is phonon focusing, which is caused by pho- nons having different velocities in different directions. Also a large variety of phenomena are observed in pulses of common particles. Superfluid helium (He II) is a unique medium for considering quasiparticle propagation. This is due to a number of factors. 1. In superfluid helium there exist quasiparticles with different energy—momentum relations, phonons with an almost linear relation, R � -rotons with a rela- tion which is close to quadratic but at a finite energy and momentum, and R �-rotons, whose velocity and momentum have opposite directions. The results of ex- periments on phonon pulses can be found in [1,2], R � -rotons in [3,4], and R �-rotons in [5]. 2. By changing the temperature of superfluid he- lium, one can control the interactions of quasiparticles with the ambient medium. So, at relatively high tem- peratures, a pulse of quasiparticles moving in super- fluid helium interacts with the gas of thermal excita- tions of He II. At relatively low temperatures, when the contribution of the thermal excitation can be ne- glected, the pulse propagates in a «superfluid vacu- um», and its evolution is just determined by the pro- perties of the quasiparticles in the pulse. © I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt, 2004 3. Interactions between quasiparticles in superfluid helium can be altered by pressure. The use of pressure allows us to reach two limiting cases: the creation of a pulse of strongly interacting phonons and a pulse of almost noninteracting phonons. At the saturated vapor pressure, low-energy phon- ons (l-phonons) with energy � � 10K have anomalous dispersion [6]. In this case, fast three-phonon pro- cesses are allowed and cause an almost instantaneous equilibrium in the system of l-phonons that form the pulse. This fast relaxation causes a unique phenome- non; a pulse of low energy-phonons creates high-en- ergy phonons (h-phonons). Their energy � � 10K, is an order of magnitude greater than the temperature of the l-phonons. This phenomena was discovered and de- scribed in [7,8], and theory of this process was pre- sented in [9–11]. The results of an experimental study of the evolution of strongly interacting phonons was presented in [12,13] and the theory for these pulses was given in [14]. At a pressure of 19 atm, the dispersion of phonons becomes normal [15]. In this case fast three-phonon processes are prohibited by energy and momentum conservation laws, and four-phonon processes are weak. Thus, it is possible to create a pulse of almost noninteracting phonons, that move in a «superfluid vacuum», formed by He II. Results of experiments on the evolution of such pulses formed by noninteracting phonons are presented in [16,17]. The main aim of this paper is the theoretical study of the evolution in time and space of a pulse of noninteracting quasiparticles propagating in vacuum. This is the first step to solving the inverse problem which would enable us to derive the initial character- istics of a phonon pulse, created by a heater, from de- tected signals at different points in space. The results of this theory are compared to the experimental data of [16,17]. It complements the theory of the evolution of a pulse of strongly interacting phonons (see [14]) and can stimulate new experiments. Moreover, the methods presented in this paper can be used to de- scribe the evolution of pulses formed by weakly inter- acting particles and quasiparticles in other fields of physics. 2. Description of a pulse of noninteracting quasiparticles A pulse of quasiparticles can be described by the quasiparticle distribution function n t( , , )p r that deter- mines the number of quasiparticles in a quantum state, which includes the point of phase space p r, . An important property of a pulse is its energy den- sity E t n p dpd p p p ( , ) ( ) r � � � 0 1 2 32 � � � , (1) where � �� ( )p is the energy of the quasiparticle, p is the solid angle in which the momenta of the quasiparticles lie, and the limits p0 and p1 determine the momentum interval for the quasiparticles under consideration. However, the signal amplitude mea- sured by a bolometer is determined by the energy flux density I through the area da, which is equal to the energy passing per unit time through a unit surface area da situated at point r: I t n p dpd p p p ( , ) ( ) ,r v Na� � � 0 1 2 32 � � gr � (2) where v pgr � � is the group velocity of the quasi- particle, and Na is the unit vector of the area da. The distribution function n, included in (1), (2) can be found from the solution of the kinetic equation � � n t n Jv rgr , (3) where J is the collision integral, which determines the change in the number of quasiparticles per unit time in the quantum state under consideration, caused by their interaction. The kinetic equation (3) should be supplemented by the initial condition n t n( , , ) ( ,p r p r)� �0 0 , (4) which determines the pulse in its initial state at t � 0. In experiments this is the time of the end of the heater pulse. Equation (3) together with the initial condition (4) is the complete formulation of the mathematical problem for deriving the function n t( , , )p r . In the general case such a problem has no analytical solution, since the collision integral J in (3) is a nonlinear integral operator, including distri- bution functions with different momenta. So various model solutions are used to treat definite physical problems. For example, in [14], a model solution for a phonon system with fast relaxation allowed a com- plete description of the evolution of a phonon pulse in superfluid helium and gave results which compared with experimental data [12,13]. Here we consider the opposite limiting case, when the interactions between quasiparticles are so weak that their contribution to the evolution of the pulse can be neglected. In this case, the collision integral J, in the kinetic equation (3), is equal to zero, and the solution of the problem can be easily found and looks like 580 Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt n t n t( , , ) ( , )p r p r v� �0 gr . (5) Using this solution we can rewrite (1) and (2) as fol- lows: E t n t p dpd p p p ( , ) ( ) ( ) r pr v� �� � 0 1 0 2 32 � � gr � , (6) I t n t p dpd p p p ( , ) ( ) ( ) r v N p, r va� �� � 0 1 0 2 32 � � gr gr � . (7) According to Eqs. (6) and (7) the temporal deforma- tion of the pulse from its initial state is determined by two factors: 1. The finite angular range of the momenta in the pulse, which is defined by the solid angle p p� 2�� , (8) where � �p p� �1 cos (9) is the angular-width parameter. If � p 0, quasi- particles in the pulse have different directions, and this results in a deformation of the initial pulse. 2. The dispersion, when quasiparticles with differ- ent momentum have different group velocities, results in a temporal widening of initial pulse. If dispersion is neglected, the group velocity can be written in the form v pgr � c p, and it follows from Eqs. (6) and (7) when n p0( , )r that E t E ct p d p p ( , ) ( )r r p � �� 0 , (10) I t c E ct p p d p p ( , ) ( )r r p p Na� �� 0 , (11) where E E t p p dp p p p 0 0 2 3 0 2 0 1 ( ) ( , ) ( , ) ( ) r r n r� � � � � � � (12) is the energy density in the initial pulse. Equations (10) and (11) have a simple physical meaning which allows us to describe, the temporal and spatial evolution of the phonon pulse without do- ing the integration. In order to find the energy density and the energy flux density at the point r and at time t according to Eqs. (10) and (11), one should consider a sphere of radius ct and center at r. The energy density at the point r is equal to the average value of the ini- tial density on the spherical surface of a segment of this sphere with solid angle p . By analogy, one can interpret the integrands in the general expressions (6) and (7) as the partial contri- bution from quasiparticles with momentum p to the energy E t( , )r and to the flux I t( , )r . 3. The shape of the quasiparticle pulse at different times Let us consider a pulse, moving along the z axis, with the initial distribution function n r n p g p p z p0 0( , ) ~ ( , ( ) ( cos )p r) r� �� � , (13) where g L r L z z( ) ( ) ( ) ( )| |r � � ��� � � (14) where � is the Heaviside step function, which is equal to unity or zero when the argument is positive or ne- gative, respectively. n r0( , )p is nonzero inside the cy- linder (14) for all values of momentum, determined by the limits of integration in Eqs. (1) and (2). This distribution function describes the pulse considered below. It consists of quasiparticles whose momenta lie inside a cone with solid angle p , moving along an axis parallel to the z axis of the cylinder (14). Let us suppose that the moduli of the group veloci- ties of the quasiparticles in this pulse, vary between the maximum value v c vmax � � 1 (15) to the minimum value v c vmin � � 2. (16) Let us obtain the equation that describes the shape of this pulse at different times. At time t, the length of the quasiparticle’s path changes from its maximum value L v tmax max� (17) to its minimum value L v tmin min� . (18) At time t, the cross section of the pulse which in- cludes the z axis (see Fig. 1) is bounded by the front, back, and side surfaces. These are formed by quasi- paricles which were in the corresponding surfaces of the initial pulse. According to (6), the front and the side surfaces of the pulse at time t are formed by the points z and r, which are the centers of segments of spherical surfaces with radius Lmax and solid angle p with axes di- rected oppositely to the z axis. The spherical surface of Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 581 a segment should not intersect the front surface of ini- tial pulse but just have a common point with the ini- tial front surface. A plane section through the symme- try axis is shown in Fig. 1. In other words, the points that form the front and side surfaces are at the maxi- mum distance Lmax from the initial pulse. As a result, when r varies between the limits 0 � � �r L (19) the front surface is the vertical plane with coordinate z L LF � �| | max . (20) As r increases in the range: L r r tF� � � ( ), (21) where r t L LF p( ) sinmax� �� � , (22) the front surface becomes an arc, which is described by the equation ( ) ( )| | maxz L r L L� � � �� 2 2 2 . (23) This equation determines the r and z coordinates of the front surface of the pulse. At r rF� the side surface is a cylinder of radius rF , with front surface at z L L p� �| | max cos� and back surface at z L p� max cos� . At points r rF� the signal is equal to zero. Here it should be noted that the second term in Eq. (22) de- scribes the widening of the pulse caused by its finite angular width. According to (6), the back surface of the pulse at time t is defined by the points z and r, which are the centers of segments of spherical surfaces with radius Lmin and solid angle p with axes directed oppositely to the z axis. The spherical surface of this segment should intersect the back surface of the initial pulse and the segment should include the maximum volume of the initial pulse. In this case, the points that form the back surface are at the minimum distance from the back surface of the initial pulse, for a given Lmin; see Fig. 1. As a result, at t t z v p � �� � �min cos , (24) where z L p� �� � tan (25) (see at Fig. 1 pulse at t t� 1), the back surface is the vertical plane with the coordinates z LB p� min cos� , (26) and the coordinates r are confined by the following inequality: 0 � �r r tB( ), (27) where r t L LB p( ) sinmin� �� � . (28) At the time t t� � (see pulse at t t� 2 in Fig. 1) the back surface is the vertical plane with coordinates zB , if the coordinates r are confined by the following in- equality r t r r ta B( ) ( )� � , (29) where r t L La p( ) sinmin� � �� . (30) For smaller values of r, in the range 0 � �r r ta ( ), (31) the back surface becomes an arc and is described by the equation z r L L2 2 2� � ��( ) min . (32) This surface connects the points r ra� , z zB� with the points r z L L� � � �0 2 2; min . (33) The points of the back vertical surface with coordi- nates r rB� , z zB� are connected with the end of the 582 Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt 0 1 2 3 4 5 6 7 8 9 4 2 0 2 4 z, mm r, m m Fig. 1. The shape and form of the structure of the phonon pulse with t sp � 4� and �p � 18�, and L� � 1mm at times t s t1 83� �. � � and t s t2 19� �� � in the plane of the cross section which includes the axis of the initial pulse. Solid lines represent phonons which move with maximum velo- city vmax � 363 m s, and dotted lines correspond to the minimum velocity vmin � 265 m s. At time t1, the rectangu- lar region, in the center of the pulse, has the initial value of the energy density for phonons with a given velocity. The intersection of all these regions gives the region with the initial value of the energy density. At time t2 the region in the center of the pulse has a relatively high energy den- sity, but the dispersion has already led to its decrease. side face, which has the following coordinates: r r L LM p� � �� max sin� , z z LM p� � max cos� by a cone-like surface formed by the dispersion. The length of the pulse L t rz ( , ), at different t at various r, is given by the relation L t r z t r z t rz ( , ) ( , ) ( , )� �front back . (34) At t t� � when r is changing between the limits (19), using the relations (20), (26), and (34) we get L t t r L L v v t ct v tz p p( , ) ( ) .| |� � � � � � � ��� � �0 1 2 2 (35) The second term on the right-hand side of this equa- tion describes the widening of the pulse caused by dispersion, the third term is caused by the finite an- gular width of the pulse, and the fourth term is due to the product of dispersion and finite angular width. At t t� � and r � 0 from relations (20), (33) and (34) we obtain L t t r L v v tz ( , ) ( )| |� � � � � �� 0 1 2 � � � � � � � � � � �L L L min min 1 2 2 . (36) At L L� ��2 2 1min in the third term of the right-hand side of Eq. (36) one can extract the contribution due to the finite angular width L t t r L v v tz ( , ) ( )| |� � � � � �� 0 1 2 � � � � � �� � � �� � �L ct L ct c v 2 2 2 2 1 min . (37) Comparing the results of (35) with (36) and (37) it follows that at any time, the widening of the pulse caused by dispersion, grows. The influence of a finite angular width has the opposite effect (near the axis of the pulse). At t t� � the length of the pulse grows linearly, and at t t� � the pulse width, caused by the finite angular width, begins to decrease as 1 t. 4. The energy density of a phonon pulse The energy density of a phonon pulse at different points in space and time is determined by the initial phonon distribution function (4). Let us consider the case when the initial distribution (5) is given in the following form: n n p g p p z p0 0( , ) ( ) ( cosp r r� � � � �� � � ��)� � . (38) In order to simplify this problem we will not take dis- persion into account. Then the energy density at dif- ferent points in space and time can be calculated by Eq. (6) taking into account (38), can be rewritten as follows: E t E g ct p d p p ( , )r r p � � � � �� � � ��� 0 , (39) where E n p p dp p p p 0 0 2 3 0 1 2 � � � � ( ) ( ) � (40) is the initial energy density. According to Eq. (39), in order to find the energy density at the point r at time t one should build a cir- cular cone with apex at the point r and radius L ctc � , (41) solid angle p and with a spherical end cap of radius Lc. The axis is directed oppositely to the z axis. Be- low, we refer to this circular cone as simply the cone. The energy density at the center of the pulse is de- termined by that part of the spherical surface which lies inside the initial pulse (38). If it lies completely inside the initial pulse, then, according to Eq. (39), the energy density at the point of the apex of this seg- ment is equal to the initial value of the energy density (40). If the spherical surface includes no points of the initial pulse, then the energy density at the apex is zero. This graphical method of explanation helps us to consider the distribution of the energy. At times t t� � (see pulse at t t� 1 on Fig. 1) the energy density has a maximum value inside the circle with coordinates z L LR c p� �| | cos� , and 0 � �r r za R( ), (42) where r z L z L za ( ) | | � � �� � � � � � � �� 1 � (43) for the case L Lc p| | � � . (44) The spherical surfaces of the cones whose apexes are situated on the circle (42) are found to lie completely inside the initial pulse, and so the energy density on the circle (42) is equal to the initial value. It should be noted that according to Eq. (35), when the inequality (44) is valid, the length of the initial pulse is greater than its longitudinal broaden- ing caused by the finite angular width. In this case, according to inequality (44), there exists a region in- side the pulse, in which the energy density is maximal and equal to the value of the energy density in the ini- Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 583 tial pulse (40) (see cylindrical region at t t� 1 on Fig. 1). This region is closed by two circles at each end, one of which is given by Eq. (42) and another is given by the following expressions: z LL c� , and 0 � �r r za R( ). (45) The side surface is described by the equation r r za R� ( ). The length of this region is z z L LR L c p� � �| | � , (46) and the length of the pulse at this time is equal to L L Lz c p� �| | � . (47) As a result, the region with the maximum value of energy density, E0, is separated from the front and back surfaces by the distance Lc p� , which is equal to the widening of the pulse. During the motion of the pulse along the z axis, the size of this region with E0 shows a monotonic decrease and disappears at the point z L z� �| | � , which is the point of intersection of the lines r za ( ) with the z axis. Here we should note that according to [14] a region of the pulse with density E0 exists in the opposite li- miting case, when the pulse is formed by strongly in- teracting phonons. However, the physical reasons for this region and its decrease are quite different in the two cases. Moreover, a phonon pulse of strongly inter- acting phonons, with the initial distribution function (38), moves along the axis z with much less transverse broadening than for noninteracting particles, in spite of the angular width and dispersion (see [14]). At t t� � , the integral (39) always has a region which gives no contribution to the energy density, i.e., where the initial pulse is absent (see pulse at t t� 2 in Fig. 1). That is why at times t t� � , the en- ergy density, at all points in the pulse, is less than the initial energy density (40). However, in this case the pulse has a region where the energy density is greater than at other points in the pulse. This region has coordinates z L L LR c � �� � �| | 2 2 and r � 0, (48) since the spherical surface of the cone, with apex at the point (48), intersects the edges of the front sur- face of the initial pulse at the point with coordinates z L� | | and r L� � . Clearly this region contributes most to the integral (39). The situation corresponding to inequality (44) in the case considered case takes place at L L L L c c | | � � � � � � � � � � � �1 1 2 2 , (49) when the length of initial pulse is greater than the widening caused by finite angular width, described by the third term in right-hand side of Eq. (36). When the inequality (49) is satisfied, the pulse at t t� � has a relatively high energy density, separated by the surfaces described by Eq. (48) and the point z L rL c � � �and 0, (50) which is the apex of the cone whose spherical surface intersects the side face of the initial pulse and touches its back face. 5. The beginning, end, and duration of the quasiparticle pulse at different spatial points Consider the pulse at t � 0, described at the begin- ning of Sec. 3. Taking into account the symmetry of the initial pulse, it is sufficient to consider the solu- tion of the problem in the plane of the cross section which includes the z axis, and coincides with the axis of the cylinder (14), which describes the surface boundaries of the initial pulse. In experiments, the signals are detected with detec- tors which are oriented so that all the quasiparticles in the pulse have v Nagr 0. (51) So, below we will consider the condition (51) to be valid. The time of the start of the signal t z rbeg ( , ) detected at the point with coordinates z and r can be found from the geometric interpretation of formula (6) in the way we now describe. One should again construct a cone with apex at the point z, r and axis directed op- positely to axis z, with solid angle p . The surface of the cone’s spherical end cap should have just one com- mon point with the front or side surfaces of the initial pulse. Then the radius of this segment is calculated and this value gives us the minimum value R z rmin( , ), determined by position of point z, r. The time of the start of the signal is given by the following formula: t z r R z r vbeg ( , ) ( , )min max � . (52) Using this interpretation, the dependence on coordi- nates z at any r is given by t z r L z L vbeg ( , ) | | max 0 � � � � � . (53) The signals at time (53) is formed by phonons moving parallel to z axis with velocity vmax from the front surface of the initial pulse. On Fig. 1 at t t� 1 and t t� 2 these phonons create the front surface. 584 Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt Figure 2 shows the dependence on time of the sig- nals at points with different values of r with a fixed value of z. In Fig. 2 the time (53) at z z� 1 and z z� 2 is shown by the vertical lines which are determined by the interval (19). On further increasing r, when L r r zF� � � ( ) , (54) where r z L z L L z z L zF p( ) ( )| | | | � � � � � � � �tan� � � . (55) The time of the start of the signal is given by t z L r r zFbeg ( , ( ))� � � � � � � �� 1 2 2 v r L z L max | |( ) ( ) . (56) The signal at time (56) is formed by quasiparticles that move at a definite angle 0 � �� �p with respect to the z axis, with velocity vmax from the point z L� | | and r L� � of the initial pulse. Figure 1 shows that at t t� 1 and t t� 2 these quasiparticles create the arc on the front surface of the pulse (23). In Fig. 2 at z z� 1 and z z� 2, the time (56) is shown by the curve (56). And, finally, when r lies between the limits r z r r zF B( ) ( )� � , (57) where r z L z L z z zB p( ) � � � � � �tan� � � , (58) and the time of the signal start is given by t z r z r r z r L vF B p beg ( , ( ) ( )) sinmax � � � � � � . (59) The signals at time (59) are formed by quasi- particles, which move at angle �p with respect to z axis with the velocity vmax from the side surface of the initial pulse. In Fig. 1 at t t� 1 and t t� 2 these quasiparticles form the cylindrical side surface r r tF� ( ) (see Eq. (22) and the second sentence after Eq. (23)). Figure 2 at z z� 1 and z z� 2 shows the time (59) by the segment of line inside the interval (59). At r r zB� ( ) there is no quasiparticle signal because �p is the maximum angle of quasiparticle motion in the initial pulse. Using the geometric interpretation of expression (6), let us find the time t z rend( , ) of the end of the sig- nal at a definite point z, r. In order to do this, it is suf- ficient to create the same cone as before but now this segment must intersect the back face of the initial pulse and include as many points of the initial pulse as possible. Then we find the radius of this segment R z rmax( , ) and calculate the desired time: t z r R z r vend( , ) ( , )max min � . (60) Using this approach, for points with z z� � , when the radius changes inside the interval 0 � �r r zB( ), (61) we find the maximum value for the time of the signal end, as will be shown below, t z v p end max( ) min( ) � �1 � . (62) The signal at time (62) is created by quasiparticles which move at the angle �p with respect to the z axis, with velocity vmin, from the back surface of the ini- tial pulse. In Fig. 1, at t t� 1, these phonons form the back surface of the pulse. Figure 2, at z z� 1, shows the time (62) by the segment of the vertical line in- side the interval (61). In Fig. 2 at z z� 1 and z z� 2 an imaginary horizon- tal line connects the end of the line describing the be- ginning of the signal (59) with the end of the vertical line, which describes the time of the end of the signal. These correspond to v v� max and v v� min. In Fig. 1 at t t� 1 and t t� 2 the corresponding construction con- nects the ends of the back vertical surfaces. Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 585 0 5 10 15 20 25 30 4 2 0 2 4 t,�s r, m m Fig. 2. The beginning, end, and duration of the same phonon signal as in Fig. 1 at points with coordinates z z1 � � and z z2 � � at different values of the radius r. Solid lines represent phonons, moving with maximum velocities, and dotted lines, with minimum velocities. The cone-like region at z z� �1 3 mm represents the maximum possible ampli- tude of the signal. At point z z� �2 7mm there is no region of constant signal. At z z� � , when the radius changes inside the inter- val r z r r zA B( ) ( )� � , (63) where r z L z z zA( ) � � � � � , (64) the time of the signal end is given by formula (62) and this signal is formed by the same quasiparticles as those at the end of the pulse at z z� 1. In Fig. 1 at t t� 2 these quasiparticles form the back surface of the pulse. In Fig. 2 at z z� 2, the time (62) is shown by the segment of the right vertical line, which is limited by inequalities (63). At z z� � , when the radius changes inside the interval 0 � �r r zA( ), (65) the time of the end of the signal is given by t z z r r z v z r LAend( , ( )) ( ) min � � � � � � �� 0 1 2 2 . (66) At this time, signals are created by quasiparticles that move at a definite angle 0 � �� �p with respect to the z axis with velocity vmin, from the point z � 0, r L� � of the initial pulse. In Fig. 1 at t t� 2 these quasi- particles form the arc of the back surface, and in Fig. 2 at z z� 2, the time (66) is described by the right arc of a hyperbola. The duration of the pulse at the point z, r is deter- mined by the equality �t z r t z r t z r( , ) ( , ) ( , )� �end beg . (67) According to (53), (56), (59), and (60) at z z� � the duration of the pulse has a maximum value, when r is inside the interval 0 � � �r L and equal to �t t z v v v v z vp p p max max min min( ) � � � � � 1 2 1 � � , (68) where t L vp � | | max is the initial length of the pulse. In the right-hand side of Eq. (68) the second term de- scribes the increase in pulse length due to dispersion, and the third term is caused by the finite angular width. Using relations (56), (62), and (67) at z z� � when L r rF� � � the pulse duration is given by � �t t v z LF � � � �max max | |( 1 � � � � �( ) ( ) )| |z L r L2 2 . (69) At r L� � the increase in the pulse length, according to (69), reaches its maximum value at �tmax. In- creasing r causes the signal length to decrease so, that at r rF� we have �t r r t z v v v vF F p p p ( ) max min � � � � � � 1 1 1 2 � � . (70) According to Eqs. (59), (62), and (67) at z z� � when r r rF B� � , the pulse length is given by �t z v r L vB p p � � � � min maxcos sin� � . (71) From this equation it follows that the length of the pulse decreases from the maximum value (70) to the minimum value at r rB� and is given by �t r r z v v v vB B p ( ) max min � � � � 1 1 2 � , (72) where the length of the signal does not depend on tp and is determined by the dispersion together with an additional contribution from the finite value of the angular width. At r rB� there is no signal. Expres- sions (68)–(72) give the length of the pulse at z z� � for all r. At z z� � , from Eqs. (53), (66), and (67), for a value of r inside the interval 0 � � �r L rAmin( , ), (73) it follows that �t t z v v v v z v r L zp� �� � � � � �� � � � � � � � � � �� � � 1 2 2 1 1 max min min � �� . (74) At r LA � � expression (74) has a maximum value (68), when r is changing inside the interval r r LA � � � . At z z� � and r LA � � , when L r rF� � � , the length of the pulse can be found from expression (69), and for greater r (r r rF B� � ) by expression (71). Using the equations (56), (62), and (66) at z z� 2 �, for r inside the interval L r rA� � � for the pulse length we get �t t z v v v v r L zv r L z LA p� � � � � � � � � �1 2 2 2 2 2max min min | | ( ) ( ) ( ) . maxv (75) On the right-hand side of this equation, the third and the fourth terms describe the contribution of the finite angular width, which decreases with increasing z. 586 Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt At z z� 2 � and L r r rA F� � � � the length of the signal is given by Eq. (69), and at r r rF B� � , by the relation (71). At r rB� there is no signal. The result (68)–(75) describes the length of the pulse at any value of r and z. From an examination of these results, it follows that the length of the pulse al- ways grows linearly due to dispersion. As to the con- tribution of the finite angular width, at first (when z z� � and 0 � � �r L ) the signal length grows lin- early with z. At z z� � such a dependence takes place only at z z z� �� � 2 and r r LA � � � . In all other cases, the contribution of the finite angular width to the increase of the pulse length is either constant (or grows when � p � 0 — see (70)), or decreases with in- creasing z. For large z the contribution of the finite angular width to the increase of pulse length can be neglected (see (74),(75)). However, it should be emphasized that according to Eq. (58) it is the finite angular width that always leads to the linear increase of the transverse width of the pulse with increasing z. As energy is conserved, this results in the signal decreasing monotonically along the z axis, during the motion of the signal. 6. Temporal dependence of the signal amplitude in space Let us consider the case when the initial distribu- tion is given by relation (38) and omit dispersion, which is not so important for the signal amplitude. In this case relation (7) can be written as I ,t c E g ct p p d p p p (r r p p Na) � � � � �� � � ��� 0 . (76) The temporal dependence of the signal amplitude at the point z r, can be obtained from the geometrical in- terpretation of the formula (76) which we used above. To find the amplitude at the point z r, we cre- ate a cone with its apex at this point, with solid angle p and radius R, and with its axis directed oppositely to the z axis. The signal amplitude at the apex is deter- mined by the averaged value of (76) using the initial energy density (40). The time t R z r( , , ) of the signal, detected at the considered point z r, , is determined by t R z r R c( , , ) � . (77) If all of the spherical-end surface of the cone lies inside the initial pulse, then, according to (76) at time (77) the signal amplitude is equal to the initial value which would be detected by the bolometer at t � 0, with the given Na , situated inside the initial pulse. If the curved surface of the cone has no points within the ini- tial pulse, the signal at time (77) is equal to zero. Using this approach we can find, without integra- tion, the point z r, at which the amplitude will have maximum value and find the time interval when the maximum value will be observed. During the initial stage of the evolution, there ex- ists, for a definite time, a region (denoted by the sub- script «spot») where the amplitude is equal to its ini- tial value. It follows from (76) at time t that this region is limited by two planes z ct z z L ctL R p� � � � �| | cos� (78) and the cylindrical surface which describes the right-hand side of the inequality 0 � � � ��r r z L cta R p( ) sin� . (79) From inequalities (78) and (79), we see that the time duration that this region exists is determined as follows: t t L c L cp p � � � � � � � � � � � � spot min sin , ( cos ) | | � �1 . (80) At t t� spot the region lies within the pulse. In the case of a sufficiently long pulses, for which L L p p | | cos sin1 � � � � � , (81) the maximum values of z where one can observe the maximum amplitude according to inequalities (78) and (80) is equal to z z t L L L zR pmax | | | |( ) cos� � � � ��spot � � . (82) In the case of sufficiently short pulses L L p p | | cos sin1 � � � � � (83) from the inequalities (78) and (80), we get z z t L R p max | | ( ) cos � � �spot 1 � . The evolution of the maximum-region for the case (81) is presented in the Figs. 1, 2. In the case (81), at fixed value z z z� �1 max in the plane r t, (see Fig. 2) the high energy density region, with the initial value of the signal amplitude, is deter- mined by the inequalities (see inequalities (78)) z L c t z cp 1 1� � � | | cos� (84) and limited by the conical side surface (see inequa- lity (79)) Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 587 r L ct p� �� sin� . (85) When z z� max the high energy density region is absent and there is no region where amplitude of the signal is constant (see Fig. 2). 7. Discussion The results obtained above allow us to discuss the experimental data presented in [16]. In these experi- ments the amplitude of phonon signal was measured by a bolometer, situated at the distance zB � 8 mm. The helium had such a low temperature (THe K� 01. ) that the thermal excitations can be neglected. A heater injects a pulse of phonons into this pure and isotropic superfluid helium («the superfluid vacu- um»). The duration of the pulse was 5�s and it moved in the direction normal to the surface of the heater. The helium was under 24 atm pressure at which phonon interactions are determined by the relatively slow four-phonon processes. This fact justified the au- thors using a model of noninteracting phonons, mo- ving ballistically. We start from this model, as it al- lows us to use the results obtained in this paper. Here we present the numerical values of all para- meters that describe experiments [16]. The sound ve- locity in helium at pressure of 24 atm is equal to c � 363 m s. The heater was made from a gold film and had dimensions L� � 0 5. mm. The velocity of longitu- dinal sound in the gold is equal to cl � 3240 m s and transverse sound, ct � 1200 m s. Using the acoustical theory of transmission [18,19], we find the maximum value �0 for a phonon radiated by a perfectly plane heater into helium: �0 0 307� �arcsin rad.c ct . (86) Below, we use this numerical value of the parameter � �p � 0 even though the evaporated gold film is rough and the conditions for a critical cone are not met. However, the phonon emission from a gold film has been measured and is found to be strongly peaked in the direction of the surface normal [20]. Then for typical distance z� , according to (25) we have z� � 158. mm, and for the length of the initial pulse L t cp| | � in this case we have L| | .� 182 mm. The size of bolometer was 1 1� mm. The amplitude of the phonon signal was measured at distance zB , which is five times more than z� . In this case we will use expression (74) for calculating the pulse length. To calculate the velocitiesv1 andv2, included to (74), we use the energy—momentum relation for phonons � �� �cp p( )1 2 , (87) where, according to [16], at a pressure of 24 atm, the dispersion parameter is equal to � � � � �8 96 10 4 2. K . According to (87) one can consider that v1 0� see (15), (16), and v2 can be found from the relation v c v2 � �( )gr , (88) where the line denotes averaging with respect to mo- mentum with an equilibrium function for the initial phonons in the pulse. Integration of (88) with the Bose—Einstein distribution function gives v c T2 0 20 028� . , (89) where T0 is the temperature that describes the distri- bution of the initial phonon pulse. Comparing the numerical values of the second and the third terms in right-hand side of Eq. (74), which describe the broadening of the pulse caused respec- tively by dispersion and the finite angular width of the initial pulse, we conclude that the contribution of the latter is negligibly small. The numerical value of the contribution of the dispersion is equal to � �t z c v v T T B disp s� � � 2 0 2 0 2 0 62 1 0 028min . . . (90) According to the result presented in Sec. 6, when zB exceeds z� by a factor of five, even in the region of maximum phonon energy density, the energy density is much less than its initial value. Under the condi- tions of the experiments considered, the intensity was 0.055 of its initial value. Figure 3 shows the temporal dependence of signal amplitude for a definite value of temperature T0 10� . K in the initial Bose—Einstein distribution function. The plots in the Fig. 3 are obtained from the formula (7) with the numerical values of the parameters given above. Figure 3 shows that the half width of the pulse broadening due to dispersion is close to (90), and the contribution of the angular width is negligible, in ac- cordance with the result obtained from Eq. (74). Estimations of the temperatures of the heater in the experiments [16] was � 1K for heater powers � 0 3. mW. However, the theoretical length of a pulse at T = 1 K ex- ceed the observed ones for heater powers � 0 3. mW. The theoretical values of the increase in pulse length, at dif- ferent initial temperatures (see [16]) agree with those calculated here, for initial temperatures in the interval 0.2–0.8 K. The same results were found in [16] using a computer calculation. These calculations were made us- ing the expression in [16], which describes the ballistic motion of phonons moving strictly parallel to the z axis, with group velocities given by the dispersion law (87). 588 Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, and A.F.G. Wyatt This model phonon pulse changes its form just because of dispersion. The conclusions from both analyses is that the tem- perature of the initial pulse of phonons is much lower than the temperature of the heater which creates them. According to the experimental conditions in [16] angular spreading does not contribute to the increase of the length of the pulse. However, the finite value of angle �p determines the transverse spreading of the pulse. Under the experimental conditions [16], this spreading results in a considerable decrease of the am- plitude of the signal along the axis of the pulse. Thus, the amplitude of the pulse at the detector is a factor 0.055 smaller than the amplitude of the initial pulse. Undoubtedly, a series of experiments at high pres- sure similar to experiments at the saturated vapor pressure [13] are of great interest. This would give the temporal and spatial dependence of the phonon energy density for different pulse lengths. The comparison of results of such experiments at high pressure and satu- rated vapor pressure would make it possible to com- pare the observed values with those calculated with the model of strongly interacting phonons (see [14]) and the results obtained here for a system of phonons without interactions. This comparison will give us im- portant information about anisotropic systems of phonons in superfluid helium. We express our gratitude to EPSRC of the UK (grant GR/S24855 and GR/N20225), and to GFFI of Ukraine (grant N02.07/000372) for support for this work. 1. A.F.G. Wyatt, N.A. Lockerbie, and R.A. Sherlock, Phys. Rev. Lett. 33, 1425 (1974). 2. M.J. Baird, F.R. Hope, and A.F.G. Wyatt, Nature 304, 325 (1982). 3. F.R. Hope, M.J. Baird, and A.F.G. Wyatt, Phys. Rev. Lett. 52, 1528 (1984). 4. M. Brown and A.F.G. Wyatt, J. Phys.: Condens. Matter 2, 5025 (1990). 5. M.A.H. Tucker and A.F.G. Wyatt, Science 283, 1150 (1999). 6. H.J. Marris and W.E. Massey, Phys. Rev. Lett. 25, 220 (1970). 7. M.A.H. Tucker and A.F.G. Wyatt, J. Phys.: Condens. Matter 6, 2813 (1994). 8. M.A.H. Tucker and A.F.G. Wyatt, J. Phys.: Condens. Matter 6, 2825 (1994). 9. I.N. Adamenko, K.E. Nemchenko, A.V. Zhukov, M.A.H. Tucker, and A.F.G. Wyatt, Phys. Rev. Lett. 82, 1482 (1999). 10. A.F.G. Wyatt, M.A.H. Tucker, I.N. Adamenko, K.E. Nemchenko, and A.V. Zhukov, Phys. Rev. B62, 9402 (2000). 11. I.N. Adamenko, K.E. Nemchenko, and A.F.G. Wyatt, Fiz. Nizk. Temp. 28, 123 (2002) [Low Temp. Phys. 28, 85 (2002)]. 12. M.A.H. Tucker and A.F.G. Wyatt, J. Low Temp. Phys. 113, 621 (1998). 13. R.V. Vovk, C.D.H. Williams, and A.F.G. Wyatt, Phys. Rev. B68, 134508 (2003). 14. I.N. Adamenko, K.E. Nemchenko, V.A. Slipko, A.F.G. Wyatt, Phys. Rev. B68, 134507 (2003). 15. R.C. Dynes and V. Narayanamurti, Phys. Rev. Lett. 33, 1195 (1974). 16. A.F.G. Wyatt, R.A. Sherlock, and D.R. Allum, J. Phys.: Condens. Matter 15, 1897 (1982). 17. A.F.G. Wyatt and D.R. Allum, J. Phys.: Condens. Matter 15, 1917 (1982). 18. I.M. Khalatnikov, An Introduction to the Theory of Superfluidity, Addison-Wesley Publishing, New York (1988). 19. I.M. Khalatnikov, Zh. Eksp. i Teor. Fiz. 22, 687 (1952) (in Russian). 20. N.A. Lockerbie, Ph. D. Thesis, Nottingham University (1975). Evolution of a pulse of noninteracting quasiparticles with dispersion and initial angular width Fizika Nizkikh Temperatur, 2004, v. 30, No. 6 589 16 18 20 22 24 26 28 0 0,2 0,4 0,6 0,8 1,0 R el at iv e en er gy de ns ity flu x t,�s 10� 10� Fig. 3. The temporal dependence of the energy density flux at z � 8 mm, t sp � 5 � , �p 176. �, and L� � 05. mm, relative to that of the initial pulse. The initial temperature of the pulse is T0 10� . K. The dotted line represents the signal when there is only dispersion, the dashed line repre- sents the signal without dispersion but with a finite angu- lar width in the initial pulse, and the solid line represents the signal with both dispersion and finite angular width. Note that the dashed and solid lines have had their y values multiplied by 10.