Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data
Different systems which could exhibit size-dependent second-order phase transitions have been studied experimentally during the last quarter century. The validity of the proposed theoretical results is verified by comparing high-resolution experimental data with an analytical evaluation of the he...
Saved in:
Date: | 2004 |
---|---|
Main Author: | |
Format: | Article |
Language: | English |
Published: |
Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
2004
|
Series: | Физика низких температур |
Subjects: | |
Online Access: | http://dspace.nbuv.gov.ua/handle/123456789/119845 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Cite this: | Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data / K.A. Chalyy // Физика низких температур. — 2004. — Т. 30, № 9. — С. 913–919. — Бібліогр.: 36 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-119845 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1198452017-06-11T03:03:41Z Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data Chalyy, K.A. Квантовые жидкости и квантовые кpисталлы Different systems which could exhibit size-dependent second-order phase transitions have been studied experimentally during the last quarter century. The validity of the proposed theoretical results is verified by comparing high-resolution experimental data with an analytical evaluation of the heat capacity of confined 4He. It is shown that the theoretical approach to the problem of the finite-size effect gives results that reasonably match the experimental data over a wide range of system sizes, from tens of nanometers up to a few micrometers for the cylindrical type of confinement geometry. The dependences of the shift of transition temperature on the cylinder size and boundary conditions are analyzed. The agreement of the results with finite-size scaling theory is confirmed. 2004 Article Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data / K.A. Chalyy // Физика низких температур. — 2004. — Т. 30, № 9. — С. 913–919. — Бібліогр.: 36 назв. — англ. 0132-6414 PACS: 65.20.+w, 67.40.Kh, 68.35.Rh http://dspace.nbuv.gov.ua/handle/123456789/119845 en Физика низких температур Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
topic |
Квантовые жидкости и квантовые кpисталлы Квантовые жидкости и квантовые кpисталлы |
spellingShingle |
Квантовые жидкости и квантовые кpисталлы Квантовые жидкости и квантовые кpисталлы Chalyy, K.A. Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data Физика низких температур |
description |
Different systems which could exhibit size-dependent second-order phase transitions have been
studied experimentally during the last quarter century. The validity of the proposed theoretical results
is verified by comparing high-resolution experimental data with an analytical evaluation of
the heat capacity of confined 4He. It is shown that the theoretical approach to the problem of the
finite-size effect gives results that reasonably match the experimental data over a wide range of
system sizes, from tens of nanometers up to a few micrometers for the cylindrical type of confinement
geometry. The dependences of the shift of transition temperature on the cylinder size and
boundary conditions are analyzed. The agreement of the results with finite-size scaling theory is
confirmed. |
format |
Article |
author |
Chalyy, K.A. |
author_facet |
Chalyy, K.A. |
author_sort |
Chalyy, K.A. |
title |
Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data |
title_short |
Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data |
title_full |
Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data |
title_fullStr |
Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data |
title_full_unstemmed |
Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data |
title_sort |
heat capacity of cylindrically confined helium: theoretical predictions versus experimental data |
publisher |
Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України |
publishDate |
2004 |
topic_facet |
Квантовые жидкости и квантовые кpисталлы |
url |
http://dspace.nbuv.gov.ua/handle/123456789/119845 |
citation_txt |
Heat capacity of cylindrically confined helium: theoretical predictions versus experimental data / K.A. Chalyy // Физика низких температур. — 2004. — Т. 30, № 9. — С. 913–919. — Бібліогр.: 36 назв. — англ. |
series |
Физика низких температур |
work_keys_str_mv |
AT chalyyka heatcapacityofcylindricallyconfinedheliumtheoreticalpredictionsversusexperimentaldata |
first_indexed |
2025-07-08T16:46:03Z |
last_indexed |
2025-07-08T16:46:03Z |
_version_ |
1837097988219994112 |
fulltext |
Fizika Nizkikh Temperatur, 2004, v. 30, No. 9, p. 913–919
Heat capacity of cylindrically confined helium: theoretical
predictions versus experimental data
K.A. Chalyy
Faculty of Physics, Kiev Taras Shevchenko National University, Acad. Glushkov Pr., 6, Kiev 03022, Ukraine
E-mail: kirchal@univ.kiev.ua
Received November 25, 2003, revised February 19, 2004
Different systems which could exhibit size-dependent second-order phase transitions have been
studied experimentally during the last quarter century. The validity of the proposed theoretical re-
sults is verified by comparing high-resolution experimental data with an analytical evaluation of
the heat capacity of confined 4He. It is shown that the theoretical approach to the problem of the
finite-size effect gives results that reasonably match the experimental data over a wide range of
system sizes, from tens of nanometers up to a few micrometers for the cylindrical type of confine-
ment geometry. The dependences of the shift of transition temperature on the cylinder size and
boundary conditions are analyzed. The agreement of the results with finite-size scaling theory is
confirmed.
PACS: 65.20.+w, 67.40.Kh, 68.35.Rh
1. Introduction: confined helium experiments
The physical properties of a liquid systems are
strongly dependent on its size. Spatial limitation of a
liquid matter system can cause significant changes in
its properties. Analysis of such microscopic systems
and their thermodynamic, static properties as well as
measurement of the optical and transport properties
has become a key task in various fields. The behavior
of condensed and soft matter systems at small length
scales is becoming a topic of technological importance
[1] in view of modern achievements in micro- and
nanofabrication techniques [2].
Most of the recent advances toward better under-
standing of the effect of spatial limitation and fi-
nite-size induced phenomena could not be realized
without basic research results of confined helium ex-
periments. Among the pioneering investigations in
this field was an experimental study by Chen and
Gasparini [3], which was conducted about quarter
century ago. Chen and Gasparini measured the spe-
cific heat of 4He near the superfluid transition for
films from about 10 � to 56 � thick and for samples
confined to a cylindrical geometry of 200 to 2000 � in
diameter.
Comparatively recently, the most reliable data
[1,4,5] on cylindrically confined helium have became
available from the Stanford research group lead by
J.A. Lipa. It was pointed out in [4] that the common
confining geometries are planes and cylinders, allow-
ing the study of crossover from bulk three-dimensional
behavior to two- and one-dimensional behavior, re-
spectively. By approaching the transition very closely,
it becomes possible to use well-defined geometries
amenable to detailed theoretical analysis. When the
correlation length � is on the order of the average dis-
tance to the walls, the sharp peak in the specific heat
of bulk helium is predicted to be rounded and shifted
to a lower temperature. However, in the region where
the confinement effects are strong, there does not ap-
pear to be a detailed theoretical prediction of the
shape of the heat capacity curve that is directly appli-
cable to cylindrical geometry [6].
In [6] Coleman and Lipa reported the results of the
heat capacity measurements as a function of tempera-
ture near the � point of 4He confined in 8-�m diameter
cylinders that have exceptionally uniform diameters
and smooth surfaces. Later in [1] the results of a more
specific analysis of the 8-�m data were presented and
the new data of an experimental study of a 0.26-�m
cylinder were discussed. In contrast to the much ear-
lier results [3] with smaller holes, they found good
agreement with scaling-law exponents derived by re-
normalization methods.
© K.A.Chalyy, 2004
For comparably large systems the size-dependent
shift of the critical temperature is extremely small,
which makes it difficult to detect due to the smearing
effect of Earth’s gravity. Up to this date most experi-
ments on the effects of confinement on helium in pla-
nar and cylindrical geometries have been restricted to
the submicron (or nano-) regime [3,7,8]. However,
Lipa et al. [1] stressed that using high-resolution ther-
mometry techniques [9], it has become possible to ex-
plore the region extending up to about 100 �m under
microgravity conditions. This circumstance provides
sufficient background to consider the comparably big
systems as «finite-size» one. This development has
dramatically increased the range over which length
scaling can be tested and has eased the problem of un-
controlled surface effects.
Experimental studies in finite-size liquids are quite
difficult to conduct due to the fact that blurring effect
of the Earth’s gravity often causes phenomena very
similar to the finite-size effect. This creates a high de-
mand for a theory that fits the existing experimental
results and which could be used for explicit calcula-
tions of the properties in a wide range of sizes and for
different types of system geometries. At present, for
the case of cylindrical geometry, the primary source of
predictions is the Monte Carlo studies (see, for exam-
ple, [10]).
In this paper, we present a comparative study of ex-
perimental data and theoretical approaches to the fi-
nite-size effect in the liquid helium heat capacity in
the close vicinity of the �-transition temperature T� .
Actual analytical calculations are conducted without
taking into account the gravity effect. Consequently,
it could be applicable for reliable comparison with the
results of microgravity experiments in Earth orbit or
with a number of Earth-based measurements of the
4He heat capacity in cylindrical confinement. In such
Earth-based measurements the characteristic size of
the system is supposed to be small enough to allow one
to neglect the influence of gravity on the precision of
new transition temperature detection.
2. Helium-in-cylinrical-confinment heat
capacity: theoretical predictions
A system can be considered as finite-sized near the
critical point or the phase transition point if its char-
acteristic linear size Z becomes comparable with the
correlation length � of the order parameter fluctua-
tion. As is well known from the fluctuation (scaling)
theory of phase transitions (see, e.g., [11–13]), the
correlation length is expressed by
� = �o�
–� , (1)
where the amplitude of the correlation length in clas-
sical liquids can reach the value �0 � 1–10 nm, and
� = |T–Tc|/Tc is the temperature variable. Nowadays,
it has become possible to approach the critical point
of classical liquids with respect to the tempera-
ture variable to within «distances» like |T–Tc|min�
� 10–2–10–3 K for Òñ � 102 K, or, using dimensionless
values, �min = |T–Tc|/Tc � 10–5 . As the critical expo-
nent of the correlation length is equal to � � 0.6705
[14] for the space dimensionality d = 3, the maxi-
mum values of the correlation length in classical
liquids can approach such values as � � 1–10 nm �
� (10–5)–0,67 � 103–104 nm.
Spatial insufficiency on an even bigger scale can be
realized for quantum liquids, where, because of the
very small values of the critical temperature, the cor-
relation length becomes two orders more; i.e., for in-
stance, in liquid helium � � 10–8 and �/�0 � 105.4.
So, classical liquid systems that have linear size Z up
to a micrometer (for quantum liquids, even to tens of
micrometers) could be considered as finite-sized near
the critical point (or � point). The physical properties
of finite-size systems preserve their singular behavior
only for the directions in which the system still has
semi-infinite size. However, taking into account the
limited temperature resolution of even extremely pre-
cise modern experiments [4,7], it should be noted that
over the limit of a system’s smallest confining size of
about a few hundreds micrometers the finite-size effect
seems to be almost impossible to detect even under
conditions of microgravitation. In the case of Earth-
based experiments due to the gravity effect that limit
will not exceed ten micrometers [15].
Using high-resolution thermometers [9] it has be-
come possible to determine the exact location of the
�-transition temperature with a precision of �5 nK [6].
Since the useful resolution of measurements at the
lambda point on Earth is limited in principle by the
pressure dependence of the transition temperature
(gravity effect), it is important to be sure that the ob-
served rounding of the specific heat is in fact produced
by finite-size effects. In spite of that, it has been sug-
gested [6] that gravity does not play a major role in
the actual Earth-based experiments with the typical
confined size less than 10 �m . This statement is sup-
ported by the fundamental studies of Ahlers (see, for
example, [16]) and could be explained in terms of di-
rect calculations of temperature variation due to the
gravity effect. In order to proceed, we have to con-
sider the limiting value of the cylinder length L over
which the results of a study of the finite-size effect be-
come independent of further increase of L. It was
found [17,18] that there is no need to take the actual
L
limit, because it turns out that for L = 5D and
914 Fizika Nizkikh Temperatur, 2004, v. 30, No. 9
K.A.Chalyy
larger (D is the diameter of a pore), the finite-size ef-
fect due to the finite value of L becomes insignificant.
The vertical variation of transition temperature �T�
over height L of the sample due to the hydrostatic
pressure induced by Earth’s gravity is determined [19]
by the expression �T� = �L with � = 1.273 �K/cm.
Thus, for cylinder of radius D = 8 �m [1,6] even the
enhanced limit of L = 80 �m corresponds to a maxi-
mum variation �T� � 10–8 K between the bottom and
top of the chamber. It is assumed [20] that data are
unaffected by gravity if |T – T�|
10�T� . In all cases
to be discussed below this condition is satisfied. This
makes it possible to neglect the gravity effect and si-
multaneously conduct a reliable study of the finite-
size effect at some range of confining sizes. The corre-
sponding theoretical predictions could be suitable for
comparison with actual Earth-based measurements
[3,5,6]. Furthermore, the same theoretical approach is
applicable for comparative study of the results of
microgravity Confined Helium Experiment (CHEX)
[1] conducted aboard the space shuttle and, hope-
fully, future experiments [17,20] which are scheduled
to be carried out on the Low Temperature Micro-
gravity Physics Facility (LTMPF) aboard the Inter-
national Space Station.
The major interest in the statistical physics ap-
proach to finite-size phase transitions is to find the
pair correlation function G2 of the order parameter
fluctuations and the associated correlation length �.
We consider here a cylindrical sample with diameter
D and length L>>D and study the region where the
value of the correlation length � becomes comparable
to or even larger than D but still smaller than L.
The correlation length could be defined [21] as � =
= (M2)
1/2 in terms of the pair correlation function G2
of the order parameter fluctuations for the reduced ge-
ometry of cylindrical form, where M2 denotes the nor-
malized second moment of G2 (see for details [22]).
The pair correlation function associated with the order
parameter fluctuations was deduced by applying
the Helmholtz operator. We have derived G2 from
the corresponding differential equation with Dirichlet
boundary conditions on the inner wall surface of a cyl-
inder and the condition of G2 decay at long distances,
i.e., zero value of G2 on the top and bottom ends of
the cylinder. In particular, this type of boundary con-
dition seems to represent the boundary conditions of
the real confined helium experiments more appropri-
ately [23,24]. As follows from the discussion in [22]
anomalous growth of the correlation length might oc-
cur at some new temperature Tc
* which is lower than
the bulk Tc and accordingly is located in the region of
negative � . For such a case the behavior of the correla-
tion length might obey a scaling law akin to (1),
namely
�* = �0�
*��, (2)
where �� = (T – Tc
*)/Tc
* is the new temperature vari-
able for a spatially limited system. It is necessary to
notice that the defining condition (2) used above
means physically that upon the achievement of the
new critical temperature Tc
* of a spatially limited liq-
uid in the sample of cylindrical geometry there might
exist anomalous growth of the longitudinal compo-
nent of the correlation length �* along the cylinder
axis as is described by
�* = �0 �
�
�
�
�
�
�
�
�
�
�
� �
�
�
�
�
�
�
�
�
�
1
1
1
K
/
( ) , (3)
where �0 is the amplitude of correlation length, K =
= r/�0 is the geometrical factor, r is the cylinder ra-
dius, � = (T – Tc)/Tc is the temperature variable, � is
the critical exponent of the correlation length, and
the quantity �1 is determined by the expression
J0(�1) = eA for the case of the constant boundary
condition G2 = A on the cylinder wall surface. Here
J0 is the zero-order Bessel function, and e is a tran-
scendental constant equal to 2.7182. In the case of
the zero boundary condition A = 0, the value of the
parameter �1 is equal to the first zero of the zero-or-
der Bessel function, �1 = 2.4048. Here, the boundary
effect to the pair correlation function G2 is translated
to the boundary condition for the corresponding dif-
ferential equation [22].
In order to verify the validity of the theoretical re-
sults proposed here, and before in [22,25], it is impor-
tant to have the possibility of comparing the results
with high-resolution experimental data. Since most of
the data [1,3,5,6] available up to this date are derived
from the study of the heat capacity of confined 4He, it
is necessary to obtain an analytical expression particu-
larly for this property. Such calculations could be con-
ducted on a background of the formulas for the tem-
perature dependence of the correlation length and
well-known scaling relations for the heat capacity.
While the correlation length � is proportional to �
–�,
the analogous formula for the temperature dependence
of the heat capacity C reads C � �
–� , where � is the
heat capacity critical exponent, which assumed posi-
tive. Combining it [26], one can get the scaling rela-
tion between the heat capacity and the correlation
length:
C � � �/�. (4)
The hyperscaling relation � = 2–3� [27] with � =
= 0.6705 yields � �/ = –0.0172. In fact, near the
Heat capacity of cylindrically confined helium
Fizika Nizkikh Temperatur, 2004, v. 30, No. 9 915
lambda point a small negative value of � � –0.026
was observed in early experiments [28]. Later, in the
most precise microgravity experiment in Earth orbit
[29], it was confirmed to be � = –0.01285. However,
it is known that near the liquid–gas critical point the
exponent � appears to be positive and close to 0.1 [30].
Using the general formula (4) it became possible to
obtain a formula for the heat capacity in a particular
type of confining geometry by assuming � = �*. The
correlation length �* expressed by Eq.(3) is directly
associated with the cylindrical geometry and corre-
sponds to the case of zero boundary condition for the
pair correlation function G2 on the inner surface of the
confining system. That is why the new expressions for
the heat capacity also will retain special parameters,
which are related to particular form of the system.
This makes possible an easy comparison with the data
of the corresponding experiments. For the case of cy-
lindrical geometry it reads
C
Kcyl
cyl
� �
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
2 11
1/
( ) . (5)
This equation (5) shows that the heat capacity in
the cylindrical geometry remains finite at the bulk T�
(� = 0) and grows up to its limit at some point below
T� in the region of negative �.
To carry out further comparison with experimental
data, here in Eq.(5) the geometrical factor was rede-
fined: Kcyl = D/�0 , with D being the cylinder diame-
ter. The amplitude of correlation length �0 for helium
is equal to 0.36 nm below and 0.143 nm above the
transition temperature [31]. The parameter �1 = 2.4048
is the first zero of the cylindrical Bessel function [32].
The critical exponent � is taken to be equal to 0.6705
[14] and � = –0.01285 [29].
Using Eq.(5) it is possible to evaluate the shift ��
of the transition temperature from the location of
the heat capacity maximum. The corresponding new
transition temperature T Kc
*( ) is determined by the ex-
pression
T Kc
*( ) = T
c
[1 + (2�1/Kcyl)
1/�]��. (6)
In Table the experimental data from Refs. 1,3,5,6 for
confined helium are combined with the results calcu-
lated from Eq. (5) for the same values of the system’s
sizes. Here, in Table ��E = (T
�
– Tm)/Tm represents
the shifts of the helium heat capacity maximum Tm
from its bulk value T
�
observed experimentally in re-
spect that T
�
> Tm , and ��T is the shift calculated ac-
cording to the proposed theoretical approach from the
expression:
��
T
= (2�1�0/D)1/�. (7)
3. Discussion
In Fig. 1 the dependence of the shift of the transi-
tion temperature ��T on the cylindrical pore diameter
D, which varies from 10 nm up to 10 �m, is presented
in log–log scale for the case of the zero boundary con-
dition. It shows that the value of ��T is decreasing
with respect to cylinder diameter growth. The slope
of the plot in Fig. 1 is 1/� � 1.49, in agreement with
the finite-size scaling theory predictions [13,33]: �� =
= aH –1/�, where a is a constant depending on the ge-
ometry. In the case considered above, the linear size of
the system H is treated as the cylinder diameter D,
and consequently the expression for the scaling coeffi-
cient a reads: a = (2�1�0)
1/�.
In Table the point set ��E C1–C4 is taken from the
graph, presented in the pioneering paper [3]. At that
916 Fizika Nizkikh Temperatur, 2004, v. 30, No. 9
K.A.Chalyy
Table
Transition temperature shift: Experimental data versus
correspondent theoretical calculation for the helium con-
fined in cylindrical samples of different diameters.
Diameter D
of the cylinder
Shift ��
E
,
Experiment
Shift ��
T
,
Theory
Point’s
number/[Ref.]
300 � 4.0·10–3 3.6·10–3 C1/[3]
800 � 6.5·10–4 8.3·10–4 C2/[3]
1000 � 4.3·10–4 5.9·10–4 C3/[3]
2000 � 1.2·10–4 2.1·10–4 C4/[3]
0.26 �m 1.53·10–4 1.44·10–4 C5/[1,5]
8 �m 1.20·10–6 0.87·10–6 C6/[6]
8.17 �m 1.24·10–6 0.85·10–6 C7/[1,5]
1·10
–4
1·10
–4
1·10
–3
1·10
–2
1·10
–1
1·10
–5
1·10
–5
1·10
–6
1·10
–6
1·10
–7
1·10
–7
1·10
–8
�
�
�
��
D , m
�
�
T
Fig. 1. The dependence of the shift of the transition tem-
perature ��T on the cylindrical pore diameter D in log–log
scale (according to Eq.(7)) — solid line. The slope of the
plot corresponds to � = 0.6705. � represents the experi-
mental data of Chen and Gasparini [3]. � represents the
experimental data of Lipa et al. [5].
time the cylindrical pores with diameters 300, 800,
1000, and 2000 � were investigated. Since it was diffi-
cult to determine ��E precisely from that graph in the
paper [3] (tabulated data were not included there),
the accuracy of these data is limited. A comparison
shows that the theoretical values ��T systematically
overestimate the shift of the new transition tempera-
ture by 30% on average, except C1, with a better
match. It should be noticed that points C1–C4 were
taken in experiments conducted about twenty five
years ago, but their results are still cited in the related
modern papers (see, for example, [6]). However, some
remarks [5] concerning the low accuracy and rough
agreement with scaling of these early results have to
be considered. Values ��T are calculated for the case
of zero boundary condition for the pair correlation
function G2 that provides the maximum possible value
of the shift �� for a particular size of the system. How-
ever, in a real experiment, depending on the specific
features of the surface interaction, the value of the
boundary constant could be nonzero — as is also ex-
pected to be the case for the binary-mixture experi-
ments — which lead to a decrease of the shift of the
new transition temperature. This could be one of the
possible reasons of this deviation ��E with ��T.
Judicious consideration of the boundary condition
is important. As we deal with a liquid in a reduced ge-
ometry of cylindrical form, here we introduce the
boundary condition to the pair correlation function G2
on the inner surface of the cylinder in the form G2 = A.
Here A is a certain constant, the exact value of which
depends on particular problem. This expression ap-
plies to the limiting case of completely hydrophobic
inner cylinder surface for A = 0 (i.e., the zero bound-
ary condition of the first type sometimes referred as
the «homogeneous condition») and hydrophilic one
for A = 1. In practice the situation could appear with
intermediate values of A, indicating a partial wetting
for which the numerical value in the function varies
from 0 to 1. In the general case of constant boundary
condition G2 = A on the cylinder surface the shifts ��T
will depend on A as
J0(D ���T /2�0) = eA. (8)
For example, for the helium sample in a cylindrical
chamber of diameter 1000 � (this value corresponds to
the point C3 in Table) in the case of the constant
boundary condition, the dependence of the shift of the
transition temperature ��T on the boundary constant
A is illustrated by Fig. 2. As can be seen from Fig. 2,
while the boundary constant A is decreasing in value
the shift ��T grows and reaches its maximum at A = 0.
Since A characterizes the border of the system, this be-
havior of the shift of the transition temperature may
be explained as being due to a competition between
the surface interaction near the walls and finite-size
effects.
Concerning point C6 for 8 �m cylinder [6], the the-
oretical value ��T appears to be about 30% smaller
than the results of that high-precision modern experi-
ment. Reference points C5 and C7 are taken from the
most recent measurements of the heat capacity of
8.17 �m and 0.26 �m cylindrically confined helium
that are available now [1,5]. This set of the points
shows better agreement between the theoretical calcu-
lations and experimental data in terms of the shift of
transition temperature, especially for the 0.26 �m di-
ameter pore. It was noticed [5] that the results of the
new measurement for 8.17 �m cylinders do not com-
pare well with previous observations [6] in the region
near the heat capacity peak. In contrast, the results
for 0.26 �m Anopore cylinders are in good agreement
with Monte Carlo predictions [10,18].
Heat capacity of cylindrically confined helium
Fizika Nizkikh Temperatur, 2004, v. 30, No. 9 917
1·10
–4
1·10
–3
1·10
–5
1·10
–6
1·10
–7
1·10
–8
�
�
T
0 0.1 0.2 0.3 0.4
A
Fig. 2. The dependence of the shift of the transition tem-
perature ��T on the boundary constant A for the sample in
a cylindrical pore of diameter 1000 � in semilog scale (ac-
cording to Eq.(8)).
1·10
–4
1·10
–5
1·10
–6
1·10
–7
80
75
70
C
,J
/m
o
l·K
�
Fig. 3. The dependence of the confined helium heat capac-
ity C on the reduced temperature � for the 8-�m diameter
cylindrical sample in semilog scale (according to Eq.(5)).
The heat capacity of 4He confined in an 8-�m diam-
eter cylindrical sample versus reduced temperature on
a logarithmic scale is presented in Fig. 3. This data
set, calculated from Eq.(5), obviously corresponds to
the region of negative values of reduced temperature
� 0 , where T < T�.
4. Conclusions
Finally, it should be noted that over the range of
temperature variable from 10–2 — with nanokelvin ac-
curacy — up to the new transition temperature of the
confined system the correlation length, calculated
from Eq. (3) in the same way at the heat capacity but
with the corresponding critical exponent � , demon-
strates growth at a rate of about 105 times. Conse-
quently, the correlation length reaches a value of the
order of tens of micrometers, which in most of the
cases considered cases is much greater than the diame-
ter of the cylinder. In fact, the same set of formulas
gave good agreement with the experimental results for
cylindrical geometry in terms of the shift of transition
temperature and also of the rate of heat capacity
growth near the new transition temperature. We be-
lieve that other static and dynamical physical proper-
ties of finite-size systems, which could be represented
near the criticality in terms of scaling relations with a
correlation length, could be described in the frame-
work of this study. At this point, the current theoreti-
cal approach gave results that are reasonably matched
to confined 4He heat capacity experimental data over
a wide range of system sizes, from tens of nanometers
up to about ten micrometers for the cylindrical type of
geometry.
In further development of this study we are going
to present the results [25] of the heat capacity calcula-
tions in the cases of planar and barlike geometry and
confirm their consistency with a number of available
experimental results [1,7,8,34–36] and Monte Carlo
simulations [10,18].
The results of the complex experimental, numerical
and theoretical research on confined systems are con-
sidered to be essential for fundamental studies. How-
ever, it could be even more important for providing
authentic and explicit information on the static and
dynamic properties of liquid matter on the micro- and
nano-scales to scientists and engineers for practical ap-
plications in the field of technical problem-solving
and the development of the novel technologies.
K. Ch. would like to express heartiest gratitude
to deceased Prof. K. Hamano, Japan and to thank
Prof. L. Bulavin at Kiev National University, Prof.
A. Chalyi at National Medical University, Ukraine,
professors K. Kubota and T. Yamamoto at Gunma
University, Japan for numerous useful discussions and
enlightening comments in the course of this research.
1. J.A. Lipa, M. Coleman, D.R. Swanson, J.A. Nissen,
Z.K. Geng, and K. Kim, Physica B280, 50 (2000).
2. G. Fasol, Science 280, 545 (1998).
3. Tar-pin Chen and F.M. Gasparini, Phys. Rev. Lett.
40, 331 (1978).
4. J.A. Lipa, D.R. Swanson, J.A. Nissen, Z.K. Geng,
P.R. Williamson, D.A. Stricker, T.C.P. Chui, U.E.
Israelsson, and M. Larson, Phys. Rev. Lett. 84, 4894
(2000).
5. J.A. Lipa, M. Coleman, and D.A. Stricker, J. Low
Temp. Phys. 124, 443 (2001).
6. M. Coleman and J.A. Lipa, Phys. Rev. Lett. 74, 286
(1995).
7. S. Mehta and F.M. Gasparini, Phys. Rev. Lett. 78,
2596 (1997).
8. S. Mehta, M.O. Kimball, and F.M. Gasparini, J. Low
Temp. Phys. 114, 467 (1999).
9. X. Qin, J.A. Nissen, D.R. Swanson, P.P. Williamson,
D.A. Stricker, J.A. Lipa, T.C.P. Chui, and U.E. Is-
raelsson, Cryogenics 36, 781 (1996).
10. N. Schultka and E. Manousakis, Phys. Rev. Lett. 75,
2710 (1995).
11. A.Z. Patashinskii and V.L.Pokrovskii, The Fluctu-
ation Theory of Phase Transitions, Pergamon Press,
Oxford, (1979).
12. M.E. Fisher, in: Critical Phenomena, Proc. Intern.
School of Physics «Enrico Fermi» Course 51, M.S.
Green (ed.), Academic Press, New York (1971).
13. M.E. Fisher, Rev. Mod. Phys. 46, 597 (1974).
14. L.S. Goldner and G. Ahlers, Phys. Rev. B45, 13129
(1992).
15. G. Ahlers, J. Low Temp. Phys. 113, 835 (1998).
16. G. Ahlers, J. Low Temp. Phys. 115, 143 (1999).
17. K. Nho and E. Manousakis, Phys. Rev. B64, 144513
(2001).
18. N. Schultka and E. Manousakis, J. Low Temp. Phys.
111, 783 (1998).
19. G. Ahlers, Phys. Rev. 171, 275 (1968).
20. P.B. Weichman, A.W. Harter, and D.L. Goodstein,
Rev. Mod. Phys. 73, 1 (2001).
21. J.L. Lebowitz and J.L. Percus, Phys. Rev. 122, 1675
(1961).
22. K.A. Chalyi, K. Hamano, and A.V. Chalyi, J. Mol.
Liquid. 92, 153 (2001).
23. V. Dohm, Phys. Script. 49, 46 (1993).
24. P. Sutterand and V. Dohm, Physica (Amsterdam)
B194–196, 613 (1994).
25. K.A. Chalyy, Doctor Thesis, Gunma Univ., Japan
(2001).
26. S. Bhattacharyya and J.K. Bhattacharjee, Phys. Rev.
B58, 15146 (1998).
27. M.E. Fisher, Rep. Prog. Phys. 30, 615 (1967).
28. K.H. Mueller, G. Ahlers, and F. Pobell, Phys. Rev.
B14, 2096 (1976).
29. J.A. Lipa, D.R. Swanson, J.A. Nissen, T.C.P. Chui,
and U.E. Israelsson, Phys. Rev. Lett. 76, 944 (1996).
918 Fizika Nizkikh Temperatur, 2004, v. 30, No. 9
K.A.Chalyy
30. J.A. Lipa, M.J. Buckingham, and C. Edwards, Phys.
Rev. Lett. 25, 1086 (1970).
31. A. Singsaas and G. Ahlers, Phys. Rev. B30, 5103
(1984).
32. I.S. Gradstein and I.M. Ryzhik, Tables of Integrals,
Sums, Rows, and Products, Nauka, Moscow (1971)
(in Russian).
33. K. Binder, Annu. Rev. Phys. Chem. 43, 33 (1992).
34. M.O. Kimball, S. Mehta, and F.M. Gasparini, J. Low
Temp. Phys. 121, 29 (2000).
35. S. Mehta, M.O. Kimball, and F.M. Gasparini, J. Low
Temp. Phys. 113, 435 (1998).
36. J.A. Nissen, T.C.P. Chui, and J.A. Lipa, J. Low
Temp. Phys. 92, 353 (1993).
Heat capacity of cylindrically confined helium
Fizika Nizkikh Temperatur, 2004, v. 30, No. 9 919
|