没有合适的资源?快使用搜索试试~ 我知道了~
量子力学或量子场论中的有限温度欧几里得两点函数的特征是与松原频率νk=2πk/β相关的离散傅立叶系数Gk,k∈Z。 我们表明,分析性意味着系数Gk必须满足无限个我们明确写下来的与模型无关的线性方程式。 尤其是,我们构建了“解析重整化组”线性映射Aμ,对于任何选择的截止μ,它都可以表示|νk| <μ的低能傅里叶系数(零模G0可能是例外), |νk|≥μ的高能量傅立叶系数,以及实时相关器和频谱函数。 通过运行一个简单的数值算法,我们证明了对Gk的精确通用线性约束可用于系统地改进任何随机近似数据集,例如,从蒙特卡洛模拟中获得的数据。 几个明确的例子说明了我们的结果。
资源推荐
资源详情
资源评论
Available online at www.sciencedirect.com
ScienceDirect
Nuclear Physics B 909 (2016) 880–920
www.elsevier.com/locate/nuclphysb
The analytic renormalization group
Frank Ferrari
a,b,∗
a
Service de Physique Théorique et Mathématique, Université libre de Bruxelles (ULB) and
International Solvay Institutes, Campus de la Plaine, CP 231, B-1050 Bruxelles, Belgique
b
Theoretical Physics Department, CERN, CH-1211 Genève, Suisse
Received 3
March 2016; received in revised form 1 June 2016; accepted 3 June 2016
Available
online 8 June 2016
Editor: Hubert
Saleur
Abstract
Finite
temperature Euclidean two-point functions in quantum mechanics or quantum field theory are
characterized by a discrete set of Fourier coefficients G
k
, k ∈ Z, associated with the Matsubara frequencies
ν
k
= 2πk/β. We show that analyticity implies that the coefficients G
k
must satisfy an infinite number of
model-independent linear equations that we write down explicitly. In particular, we construct “Analytic
Renormalization Group” linear maps A
μ
which, for any choice of cut-off μ, allow to express the low
energy Fourier coefficients for |ν
k
| <μ(with the possible exception of the zero mode G
0
), together with
the real-time correlators and spectral functions, in terms of the high energy Fourier coefficients for |ν
k
| ≥ μ.
Operating a simple numerical algorithm, we show that the exact universal linear constraints on G
k
can be
used to systematically improve any random approximate data set obtained, for example, from Monte-Carlo
simulations. Our results are illustrated on several explicit examples.
© 2016 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license
(http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP
3
.
*
Correspondence to: Service de Physique Théorique et Mathématique, Université libre de Bruxelles (ULB) and
International Solvay Institutes, Campus de la Plaine, CP 231, B-1050 Bruxelles, Belgique.
E-mail
address: frank.ferrari@cern.ch.
http://dx.doi.org/10.1016/j.nuclphysb.2016.06.003
0550-3213/© 2016
The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license
(http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP
3
.
F. Ferrari / Nuclear Physics B 909 (2016) 880–920 881
1. General presentation
Consider the space M of arbitrary two-point functions between bosonic operators A and B in
Quantum Mechanics or Quantum Field Theory, at finite temperature T = 1/β.
1
As is well-known
and will be reviewed in details in Section 2, the space M can be presented in several equivalent
ways. One can consider various real-time two-point functions (advanced, retarded, time-ordered,
etc.), which turn out to be all related to each other, since their Fourier transforms can be expressed
in terms of a unique spectral function ρ(
ω). Alternatively, one can work with the Euclidean-time
two-point function G(τ ). By the KMS condition, G is periodic and can be expanded in Fourier
series,
G(τ ) =
1
β
k∈Z
G
k
e
−iν
k
τ
, (1.1)
where the Matsubara frequencies are defined by
ν
k
= 2πk/β. (1.2)
We shall often refer to the set of Fourier coefficients (G
k
)
k∈Z
as the “data” which encodes the
two-point function. In a generic strongly coupled quantum mechanical model, this data can only
be computed numerically, using Monte-Carlo numerical simulations. Analytic non-perturbative
methods exist only in rare occasions.
2
By Carlson’s theorem [2], the real-time and Euclidean-time points of view are equivalent:
the continuous spectral function ρ(ω) can be expressed in terms of the discrete set of Fourier
coefficients G
k
and vice-versa, under some very general assumptions that are valid in all known
interesting physical theories.
3
The map between the real-time and the Euclidean-time formalism
is quite interesting and will be discussed very explicitly below.
The tw
o-point functions must satisfy general well-known constraints that follow straightfor-
wardly from the definitions and the spectral decomposition, see Section 2. For example, on top of
being β-periodic, G(τ) is analytic except at the points τ =kβ , k ∈ Z, where it is discontinuous if
A and B do not commute. This implies in particular that G
k
= O(1/k) at large |k|. The Fourier
coefficients also satisfy reality and positivity constraints depending on the reality properties of A
and B. We shall call F the real vector space of β-periodic functions satisfying all these standard
model-independent constraints.
One of the main goal of the present wo
rk is to show that M is a linear subspace of F of infi-
nite codimension. This may come as a surprise. It means that a typical set of Fourier coefficients
(G
k
)
k∈Z
satisfying all the usual constraints is actually inconsistent! Our central result is to show
that the Fourier coefficients must always obey an infinite set of universal, model-independent,
linear equations. For reasons that will become clear below, we call these equations “Analytic
Renormalization Group” (ARG) equations. We shall write down these equations ve
ry explicitly
in Section 3 and use them extensively in Sections 4 and 5.
1
We focus on the case of bosonic operators in the present paper. The case of fermionic operators can be discussed
along the same lines, with minor and straightforward modifications.
2
See e.g. [1] for a recent example from which the investigations presented in this paper originated.
3
In Quantum Field Theory, two-point functions of local operators do not in general satisfy the hypothesis of Carlson’s
theorem, due to the usual UV divergences at coinciding points. These divergences are governed by the Operator Product
Expansion which, in asymptotically free theories, can be reliably computed in perturbation theory. This problem is
handled in a standard way: one either considers smeared versions of the local operators or, more generally, one subtracts
explicitly the diverging piece in the correlator using the OPE.
882 F. Ferrari / Nuclear Physics B 909 (2016) 880–920
A spectacular concrete application of the existence of the ARG equations is as follows. Sup-
pose that we have at our disposal an approximate data set (G
a
k
)
k∈Z
, obtained, in a strongly
coupled model of interest, by using Monte-Carlo simulations, possibly combined with pertur-
bation theory at high energies. This data set corresponds to a point G
a
∈ F , randomly chosen
in a small neighborhood of the exact, but unknown, result G
e
∈ M . The point G
a
will never
belong to M , because a random approximate data set always violate the ARG equations (actu-
ally, this violation is always massive, even if the precision of the data is excellent; see below).
It is then possible to use the ARG equations to systematically improve the approximate data, by
suitably projecting G
a
onto M . An explicit algorithm implementing this idea will be presented
and tested in Section 5. In spite of its simplicity, our algorithm is able to improve the accuracy of
typical random approximate data by a factor of 2 to 4! The startling feature is that the procedure
is totally model-independent and can be applied straightforw
ardly to Monte-Carlo data in any
quantum mechanical system.
The fundamental ingredient at the basis of the ARG equations is analyticity, which is itself a
consequence of causality
. The fact that analyticity yields non-trivial constraints on real-time cor-
relation functions is of course well-known. For example, the famous Plemelj–Kramers–Kronig
identities relate the real and imaginary parts of the Fourier transform of the retarded tw
o-point
function,
Re ˜χ
r
(ω) =
1
π
P
+∞
−∞
Im ˜χ
r
(ω
)
ω
−ω
dω
, Im ˜χ
r
(ω) =−
1
π
P
+∞
−∞
Re ˜χ
r
(ω
)
ω
−ω
dω
. (1.3)
The starting point of our work is the very same analyticity property at the basis of (1.3), but we
use it in a more sophisticated way to derive the infinite set of linear constraints that the Fourier
coefficients G
k
must satisfy.
An interesting aspect of the construction is to mak
e an unexpected link, valid in any quantum
mechanical system, between the concept of Renormalization Group (RG) and analyticity. The
fundamental idea of the RG is to describe the physics below a certain RG scale μ in terms of
a Wilsonian action S
μ
that takes into account the physics at scales greater than μ. When μ is
lowered, the action S
μ
flows to a natural description of the physics at low energies. This flow is
constrained by the obvious fact that physics is independent of the arbitrary RG scale. This yields
the RG flow equations. The Analytic version of the RG that we find (the ARG) can be described
as follows. To an arbitrary RG scale μ, we associate the inte
ger k
μ
> 0 (that we also call the RG
scale by abuse of language) such that
ν
k
μ
−1
≤ μ<ν
k
μ
. (1.4)
We also introduce a strictly positive integer δ, that we call the “index” of the ARG. The ARG
states that there exist linear maps
A
+
k
μ
,δ
and A
−
k
μ
,δ
allowing to express the low energy Fourier
coefficients, for 1 ≤ k<k
μ
and −k
μ
<k≤−1 respectively, in terms of the high energy Fourier
coefficients G
k
μ
+δq
and G
−k
μ
−δq
, q ∈ N, respectively:
G
1
,G
2
,...,G
k
μ
−1
=
A
+
k
μ
,δ
G
k
μ
,G
k
μ
+δ
,G
k
μ
+2δ
,G
k
μ
+3δ
,...
, (1.5)
G
−1
,G
−2
,...,G
−k
μ
+1
=
A
−
k
μ
,δ
G
−k
μ
,G
−k
μ
−δ
,G
−k
μ
−2δ
,G
−k
μ
−3δ
,...
. (1.6)
Note that the zero mode G
0
plays a special role since it cannot be obtained, in general, from the
ARG. This subtlety is related to the phenomenon of Bose–Einstein condensation (see e.g. [1]).
F. Ferrari / Nuclear Physics B 909 (2016) 880–920 883
The above ARG maps A
±
k
μ
,δ
are not the most general one can build. The most general maps
allow to reconstruct the analytic continuations of the Fourier transforms of the retarded and
advanced two-point functions, ˜χ
r
(z) and ˜χ
a
(z), in the upper and lower half complex plane re-
spectively. This, in particular, automatically provides explicit maps between the real-time and
the Euclidean-time formalism. The general relations are of the form
˜χ
r
(ω +i) = A
+
k
μ
,δ
(ω, )
G
k
μ
,G
k
μ
+δ
,G
k
μ
+2δ
,G
k
μ
+3δ
,...
, (1.7)
˜χ
a
(ω −i) = A
−
k
μ
,δ
(ω, )
G
−k
μ
,G
−k
μ
−δ
,G
−k
μ
−2δ
,G
−k
μ
−3δ
,...
, (1.8)
for any 0 < ≤ μ. Explicit formulas for A
±
k
μ
,δ
(ω, ) are given in Section 3. From (1.7) and
(1.8), one can straightforwardly obtain similar ARG maps expressing the real-time correlators
χ
r
(t) and χ
a
(t), the spectral function, or any other real-time two-point function in terms of the
Euclidean Fourier coefficients above any RG scale k
μ
. The ARG maps (1.5) and (1.6) can also
be easily obtained from (1.7) and (1.8). Full details will be given below.
The lo
w energy Fourier coefficients G
k
, for |k| <k
μ
, obviously do not depend on the arbitrary
choice of RG scale k
μ
and index δ. On the other hand, the right-hand sides of the equations (1.5)
and (1.6) depend explicitly, and non-trivially, on k
μ
and δ. As usual, this yields renormalization
group equations. The same remark applies to the equations (1.7) and (1.8). More generally, by
abuse of language, we call “ARG equation” any universal linear relation between the Fourier
coefficients like (1.5) or (1.6).
Our results rely in an absolutely crucial wa
y on a generalization of mathematical tech-
niques first introduced in a remarkable paper by Cuniberti et al. [3].
4
The aim of [3] was to
provide explicit formulas for the reconstruction of the real-time correlators from the Euclidean-
time correlators, a notoriously difficult and important problem. One particular aspect of our
results is to provide a useful generalization and simplification of the reconstruction procedure
of [3].
The plan of the paper is as follo
ws. In Section 2, we review basic facts on two-point func-
tions in quantum mechanics: definitions of real-time and Euclidean-time correlators, spectral
decompositions and spectral function, the resolvent, analytic properties and Carlson theorem.
Section 3 is devoted to the derivation of the ARG maps and equations. We start in 3.1 e
xplaining
a simple idea at the basis of the ARG. We then present in 3.2 some useful mathematical results
on Laguerre polynomials, Pollaczek polynomials and the relation between them. These results
are used in 3.3 to build the general ARG maps
A
±
k
μ
,δ
(ω, ). We specialize to the case of the
maps
A
±
k
μ
,δ
in 3.4 and to the reconstruction of the real-time correlators from the Euclidean data
in 3.5. This eventually yields a very general multi-parameter continuous family of ARG equa-
tions. In Section 4, we discuss the numerical implementation of the ARG. We use in particular
the example of the damped harmonic oscillator to illustrate our results. In Section 5, we e
xplain
how the ARG equations can be used to systematically improve any given random approximate
data set. This is certainly the newest, most surprising and most central concrete application of
our work. We provide a simple numerical algorithm that we test successfully on several ex-
amples. Finally
, we briefly summarize our results and suggest future directions of research in
Section 6.
4
See also [4] for a concrete discussion of the construction in [3].
884 F. Ferrari / Nuclear Physics B 909 (2016) 880–920
2. Preliminaries
2.1. Basic definitions
Let H be the Hamiltonian. To simplify some formulas, we do as if the spectrum of H were
discrete. The generalization to the case of a continuous spectrum is completely straightforward
and the required modifications will be taken into account in our discussion.
5
Let {|p} be an
orthonormal basis of energy eigenstates, H|p = E
p
|p.
The partition function at temperature T = 1/β is
Z = tr e
−βH
=
p
e
−βE
p
. (2.1)
The expectation value of any operator O at temperature T is defined by
O
β
=
1
Z
tr
e
−βH
O
=
1
Z
p
e
−βE
p
p|O|p. (2.2)
The real-time and Euclidean-time evolutions are defined as usual by
O(t) = e
itH
Oe
−itH
, O
E
(τ ) = e
τH
Oe
−τH
(2.3)
respectively.
We consider tw
o particular bosonic operators A and B and denote their matrix elements as
A
pq
=p|A|q,B
pq
=p|B|q. (2.4)
The spectral function is defined by
ρ(ω) =
1
Z
p,q
E
p
=E
q
e
−βE
p
−e
−βE
q
A
pq
B
qp
δ(ω +E
p
−E
q
) +βn
0
ωδ(ω) , (2.5)
where the zero-frequency contribution
6
reads
n
0
=
1
Z
p,q
E
p
=E
q
e
−βE
p
A
pq
B
qp
. (2.6)
By using the δ-function constraint, we may also rewrite (2.5) as
ρ(ω) =
1 −e
−βω
Z
p,q
e
−βE
p
A
pq
B
qp
δ(ω +E
p
−E
q
), (2.7)
where the sum over energy eigenstates is now unconstrained and thus includes the terms with
E
p
= E
q
. Taking into account a possible continuous part in the spectrum of H , the spectral
function can be written as a sum
ρ(ω) = ρ
s
(ω) +ρ
d
(ω) +ρ
0
(ω) , (2.8)
5
Note also that most systems with a continuous spectrum can be obtained by taking the appropriate thermodynamic
limit of a compact system with a discrete spectrum.
6
This contribution is associated with the phenomenon of Bose–Einstein condensation. See [1] for a recent discussion.
剩余40页未读,继续阅读
资源评论
weixin_38728464
- 粉丝: 1
- 资源: 966
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功