Quantum effects in an anharmonic crystal
A model of quantum particles performing D -dimensional anharmonic oscillations around their equilibrium positions which form the d -dimensional simple cubic lattice Zd is considered. The model undergoes a structural phase transition when the fluctuations of displacements of particles become macr...
Gespeichert in:
Datum: | 2002 |
---|---|
1. Verfasser: | |
Format: | Artikel |
Sprache: | English |
Veröffentlicht: |
Інститут фізики конденсованих систем НАН України
2002
|
Schriftenreihe: | Condensed Matter Physics |
Online Zugang: | http://dspace.nbuv.gov.ua/handle/123456789/120681 |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Zitieren: | Quantum effects in an anharmonic crystal / Yu. Kozitsky // Condensed Matter Physics. — 2002. — Т. 5, № 4(32). — С. 601-616. — Бібліогр.: 41 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-120681 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1206812017-06-13T03:03:43Z Quantum effects in an anharmonic crystal Kozitsky, Yu. A model of quantum particles performing D -dimensional anharmonic oscillations around their equilibrium positions which form the d -dimensional simple cubic lattice Zd is considered. The model undergoes a structural phase transition when the fluctuations of displacements of particles become macroscopic. This phenomenon is described by susceptibilities depending on Matsubara frequencies ωn , n ∈ Z . We prove two theorems concerning the thermodynamic limits of these susceptibilities. The first theorem states that the susceptibilities with nonzero ωn remain bounded at all temperatures, which means that the macroscopic fluctuations in the model are always non-quantum. The second theorem gives a sufficient condition for the static susceptibility (i.e. corresponding to ωn = 0 ) to be bounded at all temperatures. This condition involves the particle mass, the anharmonicity parameters and the interaction intensity. The physical meaning of this result is that, for all D and all values of the temperature, strong quantum effects suppress critical points and the long range order. The proof is performed in the approach where the susceptibilities are represented as functional integrals. A brief description of the main features of this approach is delivered. Розглядається модель квантових частинок, які виконують D -вимірні коливання довкола їх положень рівноваги, що утворюють d -вимірну просту кубічну ґратку Zd . Ця модель зазнає фазового переходу, коли флуктуації зміщень частинок стають макроскопічними. Таке явище описується сприйнятливостями, залежними від мацубарівських частот ωn , n є Z . Ми доводимо дві теореми, що описують термодинамічні властивості цих сприйнятливостей. Перша теорема стверджує, що сприйнятливості з ненульовими ωn залишаються обмеженими при всіх температурах, а це означає, що макроскопічні флуктуації в даній моделі є завжди неквантові. Друга теорема дає достатню умову на те, щоб і статична сприйнятливість (яка відповідає ωn = 0 ) теж була обмеженою при всіх температурах. Ця умова включає в себе масу частинки, параметри ангармонізму та інтенсивність взаємодії. Фізичний сенс цього результату полягає в тому, що для всіх D і для всіх значень температури сильні квантові ефекти унеможливлюють виникнення критичних точок і далекого порядку. Доведення проводиться в рамках підходу, у якому сприйнятливості представляються за допомогою функціональних інтегралів. Дається короткий опис головних аспектів цього підходу. 2002 Article Quantum effects in an anharmonic crystal / Yu. Kozitsky // Condensed Matter Physics. — 2002. — Т. 5, № 4(32). — С. 601-616. — Бібліогр.: 41 назв. — англ. 1607-324X PACS: 05.50.-d, 64.60.-i DOI:10.5488/CMP.5.4.601 http://dspace.nbuv.gov.ua/handle/123456789/120681 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
description |
A model of quantum particles performing D -dimensional anharmonic oscillations
around their equilibrium positions which form the d -dimensional
simple cubic lattice Zd is considered. The model undergoes a structural
phase transition when the fluctuations of displacements of particles become
macroscopic. This phenomenon is described by susceptibilities depending
on Matsubara frequencies ωn , n ∈ Z . We prove two theorems
concerning the thermodynamic limits of these susceptibilities. The first theorem
states that the susceptibilities with nonzero ωn remain bounded at all
temperatures, which means that the macroscopic fluctuations in the model
are always non-quantum. The second theorem gives a sufficient condition
for the static susceptibility (i.e. corresponding to ωn = 0 ) to be bounded
at all temperatures. This condition involves the particle mass, the anharmonicity
parameters and the interaction intensity. The physical meaning of
this result is that, for all D and all values of the temperature, strong quantum
effects suppress critical points and the long range order. The proof
is performed in the approach where the susceptibilities are represented as
functional integrals. A brief description of the main features of this approach
is delivered. |
format |
Article |
author |
Kozitsky, Yu. |
spellingShingle |
Kozitsky, Yu. Quantum effects in an anharmonic crystal Condensed Matter Physics |
author_facet |
Kozitsky, Yu. |
author_sort |
Kozitsky, Yu. |
title |
Quantum effects in an anharmonic crystal |
title_short |
Quantum effects in an anharmonic crystal |
title_full |
Quantum effects in an anharmonic crystal |
title_fullStr |
Quantum effects in an anharmonic crystal |
title_full_unstemmed |
Quantum effects in an anharmonic crystal |
title_sort |
quantum effects in an anharmonic crystal |
publisher |
Інститут фізики конденсованих систем НАН України |
publishDate |
2002 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/120681 |
citation_txt |
Quantum effects in an anharmonic
crystal / Yu. Kozitsky // Condensed Matter Physics. — 2002. — Т. 5, № 4(32). — С. 601-616. — Бібліогр.: 41 назв. — англ. |
series |
Condensed Matter Physics |
work_keys_str_mv |
AT kozitskyyu quantumeffectsinananharmoniccrystal |
first_indexed |
2025-07-08T18:23:22Z |
last_indexed |
2025-07-08T18:23:22Z |
_version_ |
1837104114629083136 |
fulltext |
Condensed Matter Physics, 2002, Vol. 5, No. 4(32), pp. 601–616
Quantum effects in an anharmonic
crystal
Yu.Kozitsky
Institute of Mathematics, Maria Curie-Skłodowska University,
20-031 Lublin, Poland
Received August 30, 2002
A model of quantum particles performing D -dimensional anharmonic os-
cillations around their equilibrium positions which form the d -dimensional
simple cubic lattice Z
d is considered. The model undergoes a structural
phase transition when the fluctuations of displacements of particles be-
come macroscopic. This phenomenon is described by susceptibilities de-
pending on Matsubara frequencies ωn , n ∈ Z . We prove two theorems
concerning the thermodynamic limits of these susceptibilities. The first the-
orem states that the susceptibilities with nonzero ωn remain bounded at all
temperatures, which means that the macroscopic fluctuations in the model
are always non-quantum. The second theorem gives a sufficient condition
for the static susceptibility (i.e. corresponding to ωn = 0 ) to be bounded
at all temperatures. This condition involves the particle mass, the anhar-
monicity parameters and the interaction intensity. The physical meaning of
this result is that, for all D and all values of the temperature, strong quan-
tum effects suppress critical points and the long range order. The proof
is performed in the approach where the susceptibilities are represented as
functional integrals. A brief description of the main features of this approach
is delivered.
Key words: displacements, structural phase transition, critical point
PACS: 05.50.-d, 64.60.-i
1. Introduction
We consider the following model of anharmonic crystal. To each point l of the
d-dimensional simple cubic lattice Zd there is attached a quantum particle of mass m
performing isotropic D-dimensional anharmonic oscillations in the crystalline field.
The Hamiltonian of our model has the form
H =
1
2
∑
l,l′∈Zd
dll′(ql,ql′) +
∑
l∈Zd
Hl , (1)
Hl =
1
2m
(pl,pl) +
k
2
(ql,ql) + V ((ql,ql)), k > 0,
c© Yu.Kozitsky 601
Yu.Kozitsky
where (. . . , . . .) stands for the scalar product in RD and dll′ is a dynamical matrix
of the model. We shall suppose that dll′ is:
• invariant under translations on Zd;
• ferroelectric, which means dll′ 6 0;
• such that
∑
l∈Zd |dll′| <∞.
These assumptions imply that there exists a function φ : R+ → R+
def
= [0,+∞),
such that
dll′ = −φ(|l − l′|), (2)
where |l− l′| =
√
(l1 − l′1)
2 + . . .+ (ld − l′d)
2 is the Euclidean distance on the lattice.
The anharmonicity potential V is supposed to be of the following form
V (ξ) = aξ + b2ξ
2 + . . .+ bpξ
p, a ∈ R, b2, . . . , bp−1 > 0, bp > 0. (3)
The simplest case of such a potential is where p = 2 [18,27,34]. If the parameter
a is negative and its absolute value is big enough, the equilibria at ql = 0, l ∈ Zd
become unstable (a double-well profile for D = 1) and macroscopic displacements of
particles from these unstable positions may occur. Such objects have been studied for
many years as quite realistic models of crystalline substance undergoing structural
phase transitions [14–16,26,27,34,38]. They, and similar models, are also employed
as parts of the models describing strong electron-electron correlations caused by the
interaction of electrons with vibrating ions [18,35,36].
Starting from the pioneering paper [31], many efforts were made to show that
“the more quantum is the model, the less possible is the phase transition”. To confirm
this physical common wisdom, different methods and approaches were applied. The
first fully rigorous proof of the suppression of the long range order in models of
this type with one-dimensional oscillations was done in [40]. This effect was also
demonstrated on certain exactly solvable models of anharmonic crystals [29,39].
In [4] it was shown (also for the scalar case) that not only the long range or-
der but any critical anomaly is suppressed if the model is strongly quantum. Later
in [22–24] his result was obtained for the vector case with arbitrary D ∈ N. The
present paper gives an extended and complete presentation of the results announced
in the latter papers. In section 2 below we formulate two theorems. The physical
meaning of the first theorem is that the large fluctuations of displacements of parti-
cles, which appear at the critical point, are always non-quantum. This means that
the only static susceptibility may diverge. The second theorem states that also the
static susceptibility is bounded (i.e. non-quantum large fluctuations are suppressed
as well) at all temperatures if the model is strongly quantum. This occurs if either
the particles are light and/or the energy levels of the one-particle Hamiltonian of the
corresponding scalar model (see below) are sufficiently separated from each other.
Both these mechanisms are observed experimentally as isotopic effects (light par-
ticles) [1,38] and external hydrostatic pressure effects (big separation) [13,37]. The
proof of the second theorem is performed in the approach which employs functional
602
Quantum effects in an anharmonic crystal
integrals, briefly described in section 3. This approach, initiated in [2], now is well
elaborated in [3,8–11,20,25]. A complete description of this approach in application
to the model considered in the present paper and an extended related bibliography
are given in the review article [6]. In section 4 we discuss some consequences of our
theorems.
2. The results
The description of the Gibbs states of lattice models is performed with the help
of local Gibbs states, which are constructed by means of local Hamiltonians (see
[8,9,19,20]). For our purposes, it is enough to consider periodic states only. Let Λ
stand for a finite subset of the lattice Zd, which we may choose as a cubic box of
finite length. For such a box, we define the periodic dynamical matrix
dΛ
ll′
= −φ(|l − l′|Λ), (4)
where the function φ is the same as in (2) and the periodic distance |l− l′|Λ is taken
on the torus which one obtains identifying the opposite walls of the cube Λ. The local
Hamiltonian HΛ corresponding to the periodic boundary conditions is constructed
according to (1) and has the following form
HΛ =
∑
l,l′∈Λ
dΛ
ll′
(ql,ql′) +
∑
l∈Λ
Hl . (5)
This Hamiltonian defines a density matrix
%β,Λ = exp (−βHΛ) /Zβ,Λ , Zβ,Λ = trace exp (−βHΛ) , (6)
where β−1 is the temperature. For the Hamiltonian (5), a general theory of the
Schrödinger operators [12] implies that the latter trace exists. The local periodic
Gibbs state γβ,Λ in the box Λ is defined as
γβ,Λ(A) = trace {A%β,Λ} , (7)
where A is an observable, which is a self-adjoint operator representing a certain
physical quantity. The state γβ,Λ maps the set of such observables into the real
numbers.
The physical properties of the system in Λ are described by Matsubara functions,
which, for the observables A1, . . . , An and τ1, . . . , τn ∈ [0, β], are defined as follows
Γ β,Λ
A1,...,An
(τ1, . . . , τn) = γβ,Λ
{
Tτh
Λ
τ1(A1) . . . h
Λ
τn
(An)
}
, (8)
where
hΛ
τk
(Ak) = exp (−τkHΛ)Ak exp (τkHΛ) , k = 1, 2, . . . , n. (9)
is the Heisenberg representation of the operator Ak taken at the imaginary time
tk = iτk and Tτ stands for the τ -ordering. By (8), each Matsubara function has the
property
Γ β,Λ
A1,...,An
(τ1 + θ, . . . , τn + θ) = Γ β,Λ
A1,...,An
(τ1, . . . , τn), θ ∈ [0, β], (10)
603
Yu.Kozitsky
where addition is understood to be of modulo β. A significant role in the Euclidean
approach is played by multiplication operators. Let |Λ| be the number of lattice
points in the box Λ. By xΛ we denote the vector xΛ = (xl)l∈Λ, which has |Λ| com-
ponents xl, each of which is a D-dimensional vector. For a wave function ψ = ψ(xΛ)
and for a continuous complex valued function A(xΛ), the corresponding multiplica-
tion operator, denoted also as A, acts
(Aψ)(xΛ) = A(xΛ)ψ(xΛ). (11)
Since such multiplication operators commute, the Matsubara functions constructed
with them are continuous as functions of τ1, . . . , τn on the whole cube [0, β]n.
In the model considered the phase transition is connected with the appearance
of the macroscopic displacements of particles breaking O(D)-symmetry, thus the
operator describing fluctuations of these displacements will play an important role
in our consideration. This operator is
Q
(α)
Λ
=
1√
|Λ|
∑
l∈Λ
q
(α)
l
, α = 1, 2, . . . , D. (12)
It is a multiplication operator. The Matsubara function (8) constructed with n = 2
and A1 = A2 = Q
(α)
Λ
will be in the center of our study. To simplify notations we set
G
(α)
β,Λ(τ, τ ′)
def
= Γ β,Λ
Q
(α)
Λ
,Q
(α)
Λ
(τ, τ ′), α = 1, 2, . . . , D. (13)
By (10), this function depends on the distance |τ−τ ′|β
def
= min{|τ−τ ′|, β−|τ−τ ′|},
which makes it possible to write its Fourier transformation in the following form
G
(α)
β,Λ(τ, τ ′) =
1
β
+∞∑
n=−∞
eiωn(τ−τ ′)Ĝ
(α)
β,Λ(ωn), (14)
Ĝ
(α)
β,Λ(ωn) =
∫ β
0
G
(α)
β,Λ(τ, τ ′)eiωn(τ−τ ′)dτ ′,
where ωn = (2π/β)n, n ∈ Z are Matsubara’s frequencies. The Fourier coefficients
Ĝ
(α)
β,Λ(ωn) have direct physical meaning – they are susceptibilities of the system in
the box Λ at the temperature β−1. For every box Λ and every β, these quantities
are well defined. Let {Λl}l∈N be a sequence of boxes such that, for every l ∈ N,
Λl ⊂ Λl+1 and ⋃
l∈N
Λl = Zd.
In the sequel we deal with such sequences of boxes only, thus we shall assume that
every sequence {Λl}l∈N has the above properties.
The phase transition in the model considered may be detected as follows. If, for a
given β and an arbitrary sequence of boxes {Λl}l∈N, all the sequences {Ĝ
(α)
β,Λl
(ωn)}l∈N,
n ∈ Z are bounded, then there is no symmetry breaking at this β. One can show
604
Quantum effects in an anharmonic crystal
(see below) that all G
(α)
Λ
(τ, τ ′) and all Ĝ
(α)
β,Λ(ωn) are nonnegative, thus by (14), one
has
Ĝ
(α)
β,Λ(ωn) 6 Ĝ
(α)
β,Λ(0), (15)
for all ωn. This yields that all the sequences {Ĝ
(α)
β,Λl
(ωn)}l∈N, n ∈ Z are bounded if
{Ĝ
(α)
β,Λl
(0)}l∈N is bounded. In this case one says that the fluctuations of displacements
of particles are normal. The abnormal fluctuations appear at the critical point of
the model or below it (with respect to the temperature, for more details see [39]),
they occur if some of the above sequences of susceptibilities become unbounded
(divergent). By (15), in all such cases the sequence {Ĝ
(α)
β,Λl
(0)}l∈N should diverge.
Such divergent sequences may be made bounded by multiplying by elements of an
infinitely small sequence of numbers, e.g. by λl = |Λl|
−δ/2, δ > 0. Suppose that, for
a divergent sequence {Ĝ
(α)
β,Λl
(ωn)}l∈N, n ∈ N, one finds δn > 0 and l0 ∈ N such that
{|Λl|
−δ/2Ĝ
(α)
β,Λl
(ωn)}l>l0 ⊂ [a, b], (16)
with certain 0 < a < b <∞. Then the exponent δn > 0 gives the rate of divergence
of the sequence {Ĝ
(α)
β,Λl
(ωn)}l∈N, and may be considered as a critical exponent. One
can show that the largest value of all such exponents is δ = 1. If δ0 = 1, the
spontaneous O(D)-symmetry breaking occurs at this β and the long range order
appears. The smallest limiting point of the sequence (16) with ωn = 0 is treated
as an order parameter. A detailed discussion of phase transitions, order parameters
etc, in quantum spin models may be found in [17]. If, for a given β, (16) holds with
n = 0 and a certain δ0 ∈ (0, 1), the fluctuations are not strong enough to destroy the
symmetry but, being abnormal, they correspond to the critical point of the model
and such β is the critical inverse temperature for this model. Thus, to prove that
the model considered undergoes a phase transition, connected with the appearance
of macroscopic displacements of particles, one has to show that (16) holds true (for
sufficiently large β) with n = 0 and δ0 = 1. This was done, for models of the type
of (1), in a number of papers [10,11,16,21,27]. Here one has to remark that, for the
isotropic ferromagnetic Heisenberg model, the appearance of a long range order has
not been proven so far. Another remark is that the existence of the long range order
((16) with ωn = 0 and δ0 = 1) does not mean the existence of the critical point ((16)
with a certain ωn and δn ∈ (0, 1)).
In this article we are going to show that, for the models considered: (a) (16)
never holds for ωn 6= 0 (quantum fluctuations are always normal); (b) (16) does not
hold for any β, if the model is “strongly quantum”. To make the latter condition
more precise we have to introduce the scalar version of the model (1). Namely, let
H̃ =
1
2
∑
l,l′∈Zd
dll′qlql′ +
∑
l∈Zd
H̃l , (17)
H̃l =
1
2
p2
l
+
k
2
q2
l
+ V (q2
l
),
and the anharmonicity potential V is given by (3). The one-particle Hamiltonian H̃l
has a pure point spectrum, each element of which εs, s ∈ N is its non-degenerate
605
Yu.Kozitsky
eigenvalue (see e.g. [12]). We set
∆ = min{εs+1 − εs | s ∈ N}. (18)
One can expect that, for the double-well potentials, the above ∆ is just the frequency
of the tunnelling motion between the wells. Since the model (17) in an evident way
is connected with the model (1), the parameter ∆ determines the properties of the
latter. In what follows, we will say that the model (1) is strongly quantum if the
following condition is satisfied
m∆2 > −
∑
l′∈Zd
dll′ =
∑
l∈Zd
φ(|l|)
def
= J. (19)
Its physical interpretation could be as follows: the total force acting on a given
quantum particle from all other particles is less than the energy of its quantum zero
oscillations (we use the units with ~ = 1).
Now we are at a position to formulate our results.
Theorem 1 For the model considered, the sequences {Ĝ
(α)
β,Λl
(ωn)}l∈N with ωn 6= 0
are bounded at all temperatures.
Theorem 2 For the model considered, let the condition (19) be satisfied. Then, for
all d,D ∈ N, for any β > 0, the sequence {Ĝ
(α)
β,Λl
(0)}l∈N is bounded, which means
that the model has no critical points and no symmetry breaking at all temperatures.
The proof of the first theorem is quite simple and we give it right now. The proof of
the second one will be performed in the framework of the Euclidean approach based
on the representation of the Matsubara functions (8) as functional integrals.
Proof of theorem 1. In fact we are going to prove that, for an arbitrary box
Λ and for any β,
0 6 Ĝ
(α)
β,Λ(ωn) 6
1
mω2
n
, ωn 6= 0, (20)
where m is the particle mass. One can show that also the Hamiltonian (9) has a
pure point spectrum, but this time the eigenvalues Es, s ∈ N are degenerate due
to the O(D)-symmetry (and the corresponding parameter ∆ would be zero in this
case). Let Ψs, s ∈ N be the corresponding eigenfunctions of HΛ. We set
Qss′ = 〈Ψs|Q
(α)
Λ
|Ψs′〉. (21)
Then by (8), (12)–(14), one has
Ĝ
(α)
β,Λ(ωn) =
1
Zβ,Λ
∑
s,s′∈N
Q2
ss′
Es − Es′
(Es − Es′)2 + ω2
n
[exp(−βEs′) − exp(−βEs)] .
This immediately yields Ĝ
(α)
β,Λ(ωn) > 0. Further, one may estimate the denominator
in the above expression
(Es − Es′)2 + ω2
n > ω2
n > 0,
606
Quantum effects in an anharmonic crystal
and obtain
Ĝ
(α)
β,Λ(ωn) 6
1
ω2
nZβ,Λ
∑
s,s′∈N
Q2
ss′ (Es − Es′)
× [exp(−βEs′) − exp(−βEs)]
=
1
ω2
n
γβ,Λ
{[
Q
(α)
Λ
,
[
HΛ, Q
(α)
Λ
]]}
=
1
mω2
n
,
where [. . . , . . .] stands for commutator.
3. Euclidean approach
The proof of our second theorem is not so elementary. To prove the boundedness
stated we have to have instruments for estimating Matsubara functions (8). These
instruments are provided by the Euclidean approach. Its basic element is the repre-
sentation of the functions (8) constructed with multiplication operators as functional
integrals. We begin the description of these integrals with introducing the sets of
functions on which we define them. Let Wβ stand for the set of continuous functions
w = w(τ) = (w(1)(τ), . . . , w(D)(τ)) defined on [0, β], taking values in RD and such
that w(0) = w(β). These functions are called periodic paths. The set of all vectors
wΛ = (wl)l∈Λ, each component of which is a periodic path at the point l ∈ Λ is
denoted by Wβ,Λ. Now we have to construct the measure on Wβ,Λ which will define
the functional integral. First we introduce the Gaussian measure dχβ which we shall
use as a reference measure. This measure in determined by the first two terms in Hl
(1) and describes a vector harmonic oscillator. Since our integration sets Wβ, Wβ,Λ
are infinite dimensional vector spaces, on which Lebesgue’s measure dw does not
exist, we cannot write our Gaussian measure in the way density × dw, as is usual
for the finite dimensional case. Nevertheless, the Gaussian measure is determined
by its Fourier transformation, which for our measure dχβ reads
∫
Wβ
exp
(
i
∫ β
0
(w(τ),x(τ))dτ
)
dχβ(w) =
= exp
{
−
1
2
∫ β
0
∫ β
0
(S(τ, τ ′)x(τ),x(τ ′))dτdτ ′
}
, (22)
where
(S(τ, τ ′)x(τ),x(τ ′)) =
D∑
α,α′=1
Sαα′
(τ, τ ′)x(α)(τ)x(α′)(τ ′), (23)
Sαα′
(τ, τ ′) = δαα′Ω0
exp[Ω0(β − |τ − τ ′|)] + exp[Ω0|τ − τ ′|]
exp(Ω0β) − 1
, (24)
and Ω0 =
√
k/m is the standard parameter of the harmonic oscillator. All proper-
ties of the measure dχβ are known (see e.g. subsection 2.2 in [6]). By means of it
607
Yu.Kozitsky
we construct the Gaussian measure which describes the system of non-interacting
harmonic oscillators located in Λ:
dχβ,Λ(wΛ) =
∏
l∈Λ
dχβ(wl). (25)
Now we take into account the anharmonicity potential and the interaction between
the particles. It is performed using the Feynman-Kac formula, the detailed descrip-
tion of which the reader may find in [30] and [32]. In what follows, we introduce the
measure
dνβ,Λ(wΛ) =
1
Zβ,Λ
exp
{
−
1
2
∑
l,l′∈Λ
dΛ
ll′
∫ β
0
(wl(τ),wl′(τ))dτ
−
∑
l∈Λ
∫ β
0
V ((wl(τ),wl(τ))dτ
}
dχβ,Λ(wΛ). (26)
This measure is called a periodic local Euclidean Gibbs measure. As usual, the factor
in front of the exp is the normalization constant, which provides that the normal-
ization condition ∫
Wβ,Λ
dνβ,Λ(w) = 1, (27)
is satisfied. On the other hand, Zβ,Λ is the partition function of the system in the
box Λ and has the representation following from (27)
Zβ,Λ =
∫
Wβ,Λ
exp
{
−
1
2
∑
l,l′∈Λ
dΛ
ll′
∫ β
0
(wl(τ),wl′(τ))dτ
−
∑
l∈Λ
∫ β
0
V ((wl(τ),wl(τ))dτ
}
dχβ,Λ(wΛ). (28)
Now we are at a position to write the main formula derived in [2]. It is
Γ β,Λ
A1,...,An
(τ1, . . . , τn) =
∫
Wβ,Λ
A1(wΛ(τ1)) · · ·A1(wΛ(τn))dνβ,Λ(wΛ), (29)
where A1, . . . , An are the multiplication operators by the functions A1(xΛ), . . . ,
An(xΛ) respectively. This representation is employed to obtain information regard-
ing mathematical properties of the Matsubara functions and their sequences. It is
worth noting that representations of this kind may be used (and are used) as a base
in the study of such models by means of computer simulations (see e.g. [33] and the
references therein).
Obviously, the same construction may be also performed for the scalar Hamilto-
nian (17). This will yield the local Euclidean Gibbs measure dν̃β,Λ(wΛ)
dν̃β,Λ(wΛ) =
1
Z̃β,Λ
exp
{
−
1
2
∑
l,l′∈Λ
dΛ
ll′
∫ β
0
wl(τ)wl′(τ)dτ
608
Quantum effects in an anharmonic crystal
−
∑
l∈Λ
∫ β
0
V (w2
l
(τ))dτ
}
dχ̃β,Λ(wΛ), (30)
where the Gaussian measure dχ̃β,Λ(wΛ) is defined by (25) with dχβ(wl) replaced by
dχ̃β(wl). The latter measure is defined by (22)–(24) with α = α′ = 1. The corre-
sponding space of scalar paths will be denoted W̃β,Λ. Accordingly, the Matsubara
functions given by the integrals (29) with ν replaced by ν̃ will be written Γ̃ .
Let us consider the set of continuous real valued functions F̃Λ defined on the
set R|Λ| = {xΛ = (xl)l∈Λ | xl ∈ R, l ∈ Λ} which have the following properties. A
function A : R|Λ| → R belongs to F̃Λ if: (a) there exist l ∈ Λ such that A depends on
xl only; (b) this dependence is either odd monotone or even positive and monotone
on R+. The basic element of the proof of theorem 2 is the following assertion, proved
in [23], see also [6], section 7.
Proposition 3 [Scalar Domination] Let the continuous functions A1, . . . , An :
RD|Λ| → R have the following property: there exist α ∈ {1, 2, . . . , D} and the func-
tions Ã1, . . . , Ãn ∈ F̃Λ such that Aj(xΛ) = Ãj(x
(α)
Λ
), j = 1, 2, . . . , n. Then for arbi-
trary D ∈ N and for all τ1, . . . , τn ∈ [0, β],
0 6 Γ β,Λ
A1,...,An
(τ1, . . . , τn) 6 Γ̃ β,Λ
Ã1,...,Ãn
(τ1, . . . , τn). (31)
Let us stress that all Aj, j = 1, . . . , n above are supposed to depend on their x
(α)
Λ
with one and the same α. As a corollary of this statement we have
0 6 G
(α)
β,Λ(τ, τ ′) 6 G̃β,Λ(τ, τ ′), (32)
where the function G̃β,Λ(τ, τ ′) is defined by (13) but with Q
(α)
Λ
(12) replaced by
Q̃Λ =
1√
|Λ|
∑
l∈Λ
ql ,
and with the Matsubara function Γ replaced by Γ̃ . Thus, in order to prove the
boundedness of {G
(α)
β,Λl
}l∈N it is enough to prove it for {G̃β,Λl
}l∈N. To this end we
introduce the following local Gibbs measure. For t ∈ [0, 1], we set
dν̃
(t)
β,Λ(wΛ) =
1
Z̃β,Λ(t)
exp
{
−
t
2
∑
l,l′∈Λ
dΛ
ll′
∫ β
0
wl(τ)wl′(τ)dτ
−
∑
l∈Λ
∫ β
0
V (w2
l
(τ))dτ
}
dχ̃β,Λ(wΛ), (33)
where Z̃β,Λ(t) is the corresponding normalization constant defined by (27), which
may be written similarly to (28). By means of this measure we define (see (29), (13),
(12), (9))
G̃β,Λ(τ, τ ′|t) =
∫
W̃β,Λ
(
1√
|Λ|
∑
l∈Λ
wl(τ)
)(
1√
|Λ|
∑
l∈Λ
wl(τ
′)
)
dν̃
(t)
β,Λ(wΛ). (34)
609
Yu.Kozitsky
Clearly
G̃β,Λ(τ, τ ′|1) = G̃β,Λ(τ, τ ′), (35)
and (see (8) and (7))
G̃β,Λ(τ, τ ′|0) =
1
|Λ|
∑
l∈Λ
1
Z̃β,l
trace
{
ql exp
(
−(τ − τ ′)H̃l
)
× ql exp
(
−(β − τ + τ ′)H̃l
)}
def
= U(τ, τ ′). (36)
Here we have taken into account that the traces in the latter formula do not depend
on l due to the translation invariance of the one-particle Hamiltonian (17). We also
set
Rβ,Λ
l1l2l3l4
(τ1, τ2, τ3, τ4|t) =
=
∫
W̃β,Λ
wl1(τ1)wl2(τ2)wl3(τ3)wl4(τ4)dν̃
(t)
β,Λ(wΛ)
−
∫
W̃β,Λ
wl1(τ1)wl2(τ2)dν̃
(t)
β,Λ(wΛ) ×
∫
W̃β,Λ
wl3(τ3)wl4(τ4)dν̃
(t)
β,Λ(wΛ)
−
∫
W̃β,Λ
wl1(τ1)wl3(τ3)dν̃
(t)
β,Λ(wΛ) ×
∫
W̃β,Λ
wl2(τ)wl4(τ4)dν̃
(t)
β,Λ(wΛ)
−
∫
W̃β,Λ
wl1(τ1)wl4(τ4)dν̃
(t)
β,Λ(wΛ) ×
∫
W̃β,Λ
wl2(τ2)wl3(τ3)dν̃
(t)
β,Λ(wΛ). (37)
One can show that all G̃β,Λ and Rβ,Λ, as functions of t ∈ [0, 1], are differentiable on
t ∈ (0, 1) and continuous on t ∈ [0, 1]. Furthermore,
Rβ,Λ
l1l2l3l4
(τ1, τ2, τ3, τ4|t) 6 0, (38)
for all values of its arguments. This estimate is known as Lebowitz’ inequality, its
validity for the measure (30) was proved in section 6 of [6]. Differentiating both sides
of (34) with respect to t after some calculations we obtain
∂
∂t
G̃β,Λ(τ1, τ2|t) = −
1
2|Λ|
∑
l1,...,l4∈Λ
dΛ
l1l2
∫ β
0
Rβ,Λ
l1l2l3l4
(τ1, τ2, τ, τ |t)dτ
+
(
−
∑
l∈Λ
dΛ
ll′
)∫ β
0
G̃β,Λ(τ1, τ |t)G̃β,Λ(τ, τ2|t)dτ. (39)
Here we have taken into account the translation invariance of the moments of the
measure (33), which is connected with the invariance of the dynamic matrix dΛ
ll′
. The
above formula may be considered as an integro-differential equation with respect to
G̃β,Λ subject to the initial condition (36). By (15), which holds also for the scalar
model, to prove the boundedness we need it to be sufficient to control the sequences
of the following functions
gβ,Λ(t)
def
=
∫ β
0
G̃β,Λ(τ1, τ2|t)dτ2 =
∫ β
0
G̃β,Λ(0, τ |t)dτ. (40)
610
Quantum effects in an anharmonic crystal
For this function, we obtain from (39) the following equation
∂
∂t
gβ,Λ(t) = −
1
2|Λ|
∑
l1,...,l4∈Λ
dΛ
l1l2
∫ β
0
∫ β
0
Rβ,Λ
l1,l2l3,l4
(τ1, τ2, τ, τ |t)dτdτ1
+
(
−
∑
l∈Λ
dΛ
ll′
)
[gβ,Λ(t)]2 , (41)
subject to the initial condition (see (36))
gβ,Λ(0) = u
def
=
∫ β
0
U(0, τ)dτ. (42)
Taking into account the assumed properties of the dynamical matrix dll′ one con-
cludes that, for every box Λ,
(
−
∑
l∈Λ
dΛ
ll′
)
6
(
−
∑
l∈Zd
dll′
)
= J > 0. (43)
Let us consider the following differential equation
dh(t)
dt
= J [h(t)]2 , t ∈ [0, 1], h(0) = u. (44)
If the condition
uJ < 1, (45)
is satisfied, it has a unique solution on t ∈ [0, 1], which has the form
h(t) =
u
1 − tuJ
. (46)
Now let us compare the equations (44) and (41). By (38), the first term in (41) is
nonpositive since all dΛ
ll′
are nonpositive. Both our equations are subject to the same
initial condition. Applying standard comparison methods of the theory of differential
equations (see e.g. [41]) one obtains
∂
∂t
gβ,Λ(t) 6
dh(t)
dt
,
for all t ∈ (0, 1). This immediately yields gβ,Λ(t) 6 h(t), for t ∈ [0, 1], in particular
∫ β
0
G̃β,Λ(0, τ)dτ = gβ,Λ(1) 6 h(1) =
u
1 − uJ
, (47)
which holds for all boxes Λ. Here we have also taken into account (35).
Proof of theorem 2. By (14), (32) and (47), the sequence {Ĝ
(α)
β,Λl
(0)}l∈N is
bounded provided (45) holds. Thus, to complete the proof one has to show that the
611
Yu.Kozitsky
condition (19) does imply (45). To this end we employ the method already used in
the proof of our first theorem. By (42), (36) and (18), one has
u =
1
Z̃β
∫ β
0
trace
{
qle
−τH̃lqle
−(β−τ)H̃l
}
dτ
=
1
Z̃β
∑
s,s′∈N
[〈ψs|ql|ψs′〉]2
e−βεs − e−βεs′
εs − εs′
6
1
∆2
1
Z̃β
∑
s,s′∈N
[〈ψs|ql|ψs′〉]2
[
e−βεs − e−βεs′
]
[εs − εs′ ]
=
1
∆2
·
trace
([
ql,
[
H̃l, ql
]]
e−βH̃l
)
trace e−βH̃l
=
1
m∆2
. (48)
Here ψs, s ∈ N are the eigenfunctions of H̃l corresponding to the eigenvalues εs.
Above, in the sums
∑
s,s′ the terms with coinciding s, s′ give zero contribution be-
cause of the corresponding matrix elements of ql vanishing. In deriving the above
estimate we have taken into account that the eigenvalues εs are simple, which is not
the case forD > 1. The latter fact gives the reason why we use our scalar domination
arguments.
Since the right-hand side in (47) is a monotone function of u, we apply there the
latter estimate and obtain
Ĝ
(α)
β,Λ(0) 6
∫ β
0
G̃β,Λ(0, τ)dτ 6
1
m∆2 − J
,
which completes the proof.
4. Discussion
First let us turn to the estimate (20). Since ωn = (2π/β)n, n ∈ Z, the series
∑
n∈Z\{0}
Ĝ
(α)
β,Λ(ωn), (49)
converges for any box Λ and for every β > 0. This fact has a deep mathematical
consequence – by Ascoli’s theorem (see e.g. [28], p. 72), the sequences of Matsub-
ara functions {Γ β,Λl
A1,...,An
}l∈N (for all multiplication operators A1, . . . , An, for which
they make sense) converge uniformly on [0, β]n provided the sequence {Ĝ
(α)
β,Λl
(0)}l∈N
is bounded. This convergence of the Matsubara functions imply in turn that the
limiting Gibbs measure, which describes the equilibrium state of the whole anhar-
monic crystal, exists in this case. It is called periodic Euclidean Gibbs measure (see
[9]). The boundedness of the sequence {Ĝ
(α)
β,Λl
(0)}l∈N means the absence of abnor-
mal fluctuations, it may be obtained by satisfying the condition (19), as follows
from theorem 2. The above fact, however, does not imply that there exists the only
612
Quantum effects in an anharmonic crystal
one limiting Gibbs measure (the latter would mean that there exists just one phase
and no phase transitions may occur). In fact, the uniqueness of the limiting Gibbs
measures may be proved based on the more sophisticated methods. Here one has to
mention that such a uniqueness was proved in [5] to hold for D = 1 and for the val-
ues of the particle mass from the interval (0, m∗(β)), where the bound m∗(β) tended
to zero as β → +∞. In a very recent paper [7] it is stated that the uniqueness of the
limiting Gibbs measures (again for D = 1 only) may be guaranteed by a condition
similar to (19), which holds for m ∈ (0, m∗), where m∗ is independent of β. In the
vector case D > 1, no uniqueness results have been obtained so far.
Now, let us say some words about the condition (19). In the harmonic case V = 0,
∆ = Ω0, thus m∆2 = k and this condition means k > J . For the potential V given
by (3), one can show that there exists C > 0 such that m∆2 ∼ Cm−(p−1)/(p+1) as
m→ 0, which means that the condition (19) is satisfied for m ∈ (0, m∗), where the
bound m∗ depends on the parameters J and a, b2, . . . , bp only.
Acknowledgement
This work was financially supported by the Deutsche Forschungsgemeinschaft
through the German-Polish project 436 POL 113/98/0–1 “Probability Measures”
which is gratefully acknowledged. The author is also grateful for kind hospitality
extended to him at the Forschungszentrum BiBoS, Universität Bielefeld in June-
July 2002.
References
1. Aksenov V.L., Plakida N.M., Stamenković S. Neutron Scattering by Ferroelectrics.
Singapure, World Scientific, 1990.
2. Albeverio S., Høegh–Krohn R. Homogeneous random fields and quantum statistical
mechanics. // J. Funct. Anal., 1975, vol. 19, p. 242–279.
3. Albeverio S., Høegh–Krohn R. Mathematical Theory of Feynman Path Integrals. Lec-
ture Notes in Math. 523, Berlin New York, Springer, 1976.
4. Albeverio S., Kondratiev Yu., Kozitsky Yu. Suppression of critical fluctuations by
strong quantum effects in quantum lattice systems. // Comm. Math. Phys., 1998,
vol. 194, p. 493–512.
5. Albeverio S., Kondratiev Yu., Kozitsky Yu., Röckner M. Uniqueness for Gibbs mea-
sures of quantum lattices in small mass regime. // Ann. Inst. H. Poincaré, Probab.
Statist., 2001, vol. 37, p. 43–69.
6. Albeverio S., Kondratiev Yu., Kozitsky Yu., Röckner M. Euclidean Gibbs states of
quantum lattice systems. Preprint BiBoS Nr. 01-03-034, Bielefeld, 2001 (to appear in
Rev. Math. Phys.).
7. Albeverio S., Kondratiev Yu., Kozitsky Yu., Röckner M. Gibbs states of a quantum
crystal: uniqueness by small particle mass. // C. R. Acad. Sci. Paris, Ser. I, 2002,
vol. 335, p. 693–698.
8. Barbulyak V.S., Kondratiev Yu.G. Functional integrals and quantum lattice systems:
I. existence of Gibbs states. // Rep. Nat. Acad. Sci of Ukraine, 1991, No. 9, p. 38–40.
613
Yu.Kozitsky
9. Barbulyak V.S., Kondratiev Yu.G. Functional integrals and quantum lattice systems:
II. periodic Gibbs states. // Rep. Nat. Acad. Sci of Ukraine, 1991, No. 8, p. 31–34.
10. Barbulyak V.S., Kondratiev Yu.G. Functional integrals and quantum lattice systems:
III. phase transitions // Rep. Nat. Acad. Sci of Ukraine, 1991, No. 10, p. 19–21.
11. Barbulyak V.S., Kondratiev Yu.G. The quasiclassical limit for the Schrödinger oper-
ator and phase transitions in quantum statistical physics. // Func. Anal. Appl. 1992,
vol. 26(2), p. 61–64.
12. Berezin F.A., Shubin M.A. The Schrödinger Equation. Dordrecht Boston London,
Kluwer Academic Publishers, 1991.
13. Blinc R., Žekš B. Soft Modes in Ferroelectrics and Antiferroelectrics. Amsterdam,
Noth-Holland Publishers, 1974.
14. Bruce A.D., Cowley R.A. Structural Phase Transitions. Taylor and Francis Ltd., 1981.
15. Choquard Ph. The Anharmonic Crystal. New York, W.A. Benjamin, 1967.
16. Driessler W., Landau L., Perez J.F. Estimates of critical lengths and critical tem-
peratures for classical and quantum lattice systems. // J. Stat. Phys. 1979, vol. 20,
p. 123–162.
17. Dyson F.J., Lieb E.H., Simon B. Phase transitions in quantum spin systems with
isotropic and nonisotropic interactions. // J. Stat. Phys. 1978, vol. 18, p. 335–383.
18. Freericks J.K., Jarrell M., Mahan G.D. The anharmonic electron-phonon problem. //
Phys. Rev. Lett., 1996, vol. 77, p. 4588–4591.
19. Georgii H.-O. Gibbs Measures and Phase Transitions. Vol. 9. New York Berlin, Walter
de Gruyter, Springer, 1988.
20. Globa S.A., Kondratiev Yu.G. The construction of Gibbs states of quantum lattice
systems. // Selecta Math. Sov., 1990, vol. 9, p. 297–307.
21. Kondratiev Ju.G. Phase transitions in quantum models of ferroelectrics. – In: Stochas-
tic Processes, Physics, and Geometry, vol. 2. Singapore New Jersey, World Scientific,
1994, p. 465–475.
22. Kozitsky Yu. Quantum effects in a lattice model of anharmonic vector oscillators. //
Lett. Math. Phys. 2000, vol. 51, p. 71–81.
23. Kozitsky Yu. Scalar domination and normal fluctuation in N -vector quantum anhar-
monic crystals. // Lett. Math. Phys. 2000, vol. 53, p. 289–303.
24. Kozitsky Yu. Quantum effects in lattice models of vector anharmonic oscillators. –
In: Stochastic Processes, Physics and Geometry: New Interplays, vol. 2, Leipzig, 1999,
p. 403–411, CMS Conf. Proc., 29, Amer. Math. Soc., Providence, RI, 2000.
25. Kozitsky Yu. Gibbs states of a lattice systems of quantum anharmonic oscillators. –
In: Noncommutative Structures in Mathematics and Physics, Kiev, 2000, p. 415–425,
(S. Duplij and J. Wess eds) NATO Sci Ser. II Math. Phys. Chem. 22, Kluwer Acad.
Publ. Dordrecht, 2001.
26. Kozitsky Yu.V., Levitskii R.R., Stasyuk, I.V. Polariton states in the order-disorder
type ferroelectrics. // Theor. Math. Phys. 1979, vol. 39, p. 106–117.
27. Pastur L.A., Khoruzhenko, B.A. Phase transitions in quantum models of rotators and
ferroelectrics. // Theor. Math.Phys. 1987, vol. 73, p. 111–124.
28. Mujica J. Complex Analysis in Banach Spaces. Amsterdam New York Oxford, North
Holland, 1986.
29. Plakida N.M., Tonchev M.S. Quantum effects in a d-dimensional exactly solvable
model for a structural phase transition. // Phys. A, 1986, vol. 136, p. 176–188.
30. Roepstorff G. Path Integral Approach to Quantum Physics. Berlin Heidelberg,
614
Quantum effects in an anharmonic crystal
Springer-Verlag, 1993.
31. Schneider T., Beck H., Stoll E. Quantum effects in an n-component vector model for
structural phase transitions. // Phys. Rev. 1976, vol. B13, p. 1123–1130.
32. Simon B. Functional Integration and Quantum Physics. New York San Francisco Lon-
don, Academic Press, 1979.
33. Sokal A. Monte Carlo methods in statistical mechanics: foundations and new algory-
thms. – In: Functional Integration, Cargèse, 1996, p. 131–192, NATO Adv. Sci. Inst.
Ser. B Phys., 361, New York, Plenum, 1997.
34. Stamenković S. Unified model description of order-disorder and displascive structural
phase transitions. // Condens. Matter Phys., 1998, vol. 1, No. 2(14), p. 257–309.
35. Stasyuk I.V. Local anharmonic effects in high-Tc superconductors. Pseudospin-
electron model. // Condens. Matter Phys., 1999, vol. 2(19), p. 435–446.
36. Stasyuk I.V. Approximate analitical dynamical mean-field approach to strongly cor-
related electron systems. // Condens. Matter Phys., 2000, vol. 3(22), p. 437–456.
37. Tibballs J.E., Nelmes R.J., McIntyre G.J. The crystal structure of tetragonal KH2PO4
and KD2PO4 as a function of temperature and pressure. // J. Phys. C: Solid State
Phys. 1982, vol. 15, p. 37–58.
38. Vaks V.G. Introduction to the Macroscopic Theory of Ferroelectrics. Moscow, Nauka,
1973 (in Russian).
39. Verbeure A., Zagrebnov V.A. Phase transitions and algebra of fluctuation operators
in exactly soluble model of a quantum anharmonic crystal. // J. Stat. Phys. 1992,
vol. 69, p. 37–55.
40. Verbeure A., Zagrebnov V.A. No-go theorem for quantum structural phase transition.
// J. Phys. A: Math. Gen. 1995, vol. 28, p. 5415–5421.
41. Walter W. Differential and Integral Inequalities. Berlin Heidelberg New York,
Springer-Verlag, 1970.
615
Yu.Kozitsky
���������
��
���������������������������� ������ � ��"!��#����$��#��%&�
')(+*�,�-�. / 0�12. 3
4 5�687:9:7<;=7�>
?:78@2>
?:7:92A=9CBD;E5GF HI@2JK6G9:78@:7�F >�LIMN?2J�F OEP�Q�J�F RTS�AVU2W�X
W�HI6GYIA<WGO
Z�WEU�YI[&?CB=\G]�RT]K^#_
`aQabEU�F 5
c�dKe
fhg�ihj2kmlhn"o�p#e
qCjKrtsCnhn
sme�u
v�WKwEx�UKy8Xh?GzI7:Y=6:y�> W�X
@:U�Y{A=H8?25878W�H�9I|{}=?G687:925�W�A:B~yKA=F�H�92A<W�58;EQ{7:Y
D
R�H�92>�F J�5GF
A<WEU�92H8?25G5:yaX
W�H�A<WEUK?aO |��EWEU2W���@25GYaJ�F H�5�W�H8?2xV9CB8[�W�;=7:HIW�J�Q�Q{7:Y
d
R�H�92>�F J�58;
�KJKW2687<;�AV;:b�F }:58;��VJG?:7:AV;
Z
d
L:��y�> W�X
@:U�Y�wI?EwE5:?Gz���?Ew:W�HIW�x�WN�E@2JK@I|:W�X�;IB:A<WCR
U�9���UE;EAV7<;I?2�2F OKwE>�F [�@25GY�}=?G687:925�W�A�687=?2Q{7:Y�>
?2A=JKW26GA<W��KF }:5G92>�9CL=��?2A<@�yKH�9�R
[�@mW��K9G68;:zI7:Y=6:y"6G�KJ�92�25:yI7VU�92HIW2687=yK>�9CB2wI?:U2@I��5G92>�9�H�F X�>
?2�:;:bK?2J�F HI6GYIA=9I|
}=?G6878WE7
ωn
B
n ∈ Z
LKMt9�X
W�HIW�X�92> W&X�H�FD78@GW�JK@2>�9CBG[�W�W��K9G68;EQ{7:Y�78@2J�> WCR
X�925:?2>�F }:5GFVH�UK?G687:92HIW2687:F��29I|{6G�KJ�92�25:yI7VU�92HIW26878@2�CL8Z�@2J���? 78@GW�JK@2>
?�687:HI@2J�R
X�� ;:z�BI[�W�6G�KJ�92�25:yI7VU�92HIW2687:FVw�5�@258;=U�Y=W�H�92>�9
ωn
wI?:U�92��?2Q{7:Y=6:y�W2b�> @I��@�R
5G92>�9&�KJ�9�HI6GF |{78@2>��E@2JG?:7<;EJG?I|2BD?��G@�WKwE5:?G}=?Gz�B=[�W&>
?2A=JKW26GA<W��KF }:5GF���UE;EAV7<;GR
?2�2F O:H�Xh?25GF ��> W�X
@:U�F�z�wI?2H<�
X�9�5�@2A=H8?25878W�H�F L���J�;Ex�?�78@GW�JK@2>
? Xh?Gz X
W2687=?:7:5GQ
;E> W�H�;�5:?�78@CBI[�W2bNFD687=?:7:9G}:5:?&6G�KJ�92�25:yI7VU�92H�F 687:Y���yKA~?�H�F X��EW�H�F Xh?Gz
ωn = 0 �78@I��bE;=UK?�W2b�> @I��@25�W�Q��KJ�9&HI6GF |{78@2>��E@2JG?:7<;EJG?I|2L8��ya;E> W�H8?�H�AVU�Qa}=?Gz�H�6�@�R
b2@�>
?G68;�}=?G687:925GA=9CB=��?2JG?2> @:7:J�9�?25Gx�?2J�> W�5GF wE>�;�7=?�F 5878@25�6G92H�5GF 687:Y�H8wI?GzK> WCR
X�F O�L:��F wE9G}:5G92��6�@25�6��2Y=W�x�W"JK@Ew=;=U�Y�7=?:7<;��EWEUKyKx�?Gz�Ha78W�>�;IB�[�W&X�UKy�HI6GF |
D
F
X�UKy&HI6GF |�wE5:?G}8@25GY
78@2>��E@2JG?:7<;EJ�9�6G9:U�YI5GF<A=H8?25878W�H�F�@E��@2AV7:9a;E5�@2> W�� U�92H�U�Q�R
Q{7:Y H�925G92A=5�@25G5:y�A=J�9:7:9G}:5G9I|
78W2}8W�A�F+Xh?:U2@2A<W�x�W&�EW�JGy8X�AV;ILT��W�HI@IX
@25G5:y��KJKWCR
HIW�X�9:7:Y=6:y�HaJG?2>�A~?I|a�KF X�|:W�X�;IB�;�yKA<W�>�;�6G�KJ�92�25:yI7VU�92HIW2687:F~�KJK@IX
687=?2H�UKyKQ{7:Y:R
6:y�wI?�X
W��EW�> W�x�W�Q���;E5GA=�2F W�5:?:U�YI5G9I|�F 5878@2xVJG?:U�F HKL~��?GzI7:Y=6:y�A<W�JKWE7:A=92��W��K9G6
x�WEU2W�H�5G9I|a?G6G�E@2AV7:F H��2Y=W�x�WN�KF X�|:W�X�;IL
�C���N�G�
�2�= G�
�
�K¡#¢�£I¤�¥ ¦t§E¨K¨E©Gª�«G¬TC®8¯I¬�®<h¨K°2±t²a³G£�´Gµ�°2±�¶�§I#§=·�¥ ¸
ª�¯Dh°C¬�°h¹~¨K³
¬D´�¹~¯�³
º2»{¼�½�¢E¾G¿2ÀÁ¿�¾KÀÃÂTÄ2ª�ÅEÆ�À�Å�¾KÀÃÂ�Ç
616
|