Coherent forward scattering as a robust probe of multifractality in critical disordered media Maxime Martinez1Gabriel Lemarié123Bertrand Georgeot1Christian Miniatura23456and Olivier Giraud7

2025-04-27 0 0 3.72MB 24 页 10玖币
侵权投诉
Coherent forward scattering as a robust probe of multifractality in critical disordered
media
Maxime Martinez,
1
Gabriel Lemarié,
1, 2, 3
Bertrand Georgeot,
1
Christian Miniatura,
2, 3, 4, 5, 6
and Olivier Giraud
7
1
Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France
2
MajuLab, CNRS-UCA-SU-NUS-NTU International Joint Research Unit,Singapore
3
Centre for Quantum Technologies, National University of Singapore, Singapore
4
Université Côte d’Azur, CNRS, INPHYNI, Nice, France
5
Department of Physics, National University of Singapore, Singapore
6
School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore
7
Université Paris Saclay, CNRS, LPTMS, 91405 Orsay, France
(Dated: October 10, 2022)
We study coherent forward scattering (CFS) in critical disordered systems, whose eigenstates are
multifractals. We give general and simple arguments that make it possible to fully characterize
the dynamics of the shape and height of the CFS peak. We show that the dynamics is governed
by multifractal dimensions
D1
and
D2
, which suggests that CFS could be used as an experimental
probe for quantum multifractality. Our predictions are universal and numerically verified in three
paradigmatic models of quantum multifractality: Power-law Random Banded Matrices (PRBM),
the Ruijsenaars-Schneider ensembles (RS), and the three-dimensional kicked-rotor (3DKR). In the
strong multifractal regime, we show analytically that these universal predictions exactly coincide
with results from standard perturbation theory applied to the PRBM and RS models.
PACS numbers: 05.45.Df, 05.45.Mt, 71.30.+h, 05.40.-a
I. INTRODUCTION
Wave transport in disordered systems is a long-standing
topic of interest in mesoscopic physics. In particular, wave
interference can have dramatic consequences on quantum
transport properties. The most celebrated example is
probably Anderson localization (AL) [
1
], that is, the
suppression of quantum diffusion and the exponential
localization of quantum states. AL is ubiquitous in wave
physics and has been observed in many experimental
situations: with acoustic waves [
2
,
3
], light [
4
8
], matter
waves [915].
Appearance of AL depends on several characteristics,
in particular dimensionality, disorder strength and corre-
lations. For instance, it is well established that 3d disor-
dered lattices undergo a genuine disorder-driven metal-
insulator transition (MIT), associated with a mobility
edge in the spectrum, separating the insulating phase
with localized eigenstates from the conducting phase with
extended eigenstates. Near the critical point of such
disorder driven transitions, eigenstates
φα
(with energy
ωα
) can display multifractal behavior, for instance at the
MIT in Anderson model [
16
18
] and graphs [
19
21
], but
also for Weyl-semimetal–diffusive transition [
22
]. They
are extended but non-ergodic, and characterized by the
anomalous scaling of their moments Iq(E):
Iq(E) = hPn|φα(n)|2qδ(Eωα)i
hPαδ(Eωα)iNDq(q1),(1)
where
Dq
are the multifractal dimensions, forming a con-
tinuous set with
q
real (
h. . .i
represents an average over
disorder configurations). Extreme cases
Dq
= 0 and
Dq
=
d
(the dimension of the system) for all
q
, cor-
respond respectively to localized and extended ergodic
eigenstates.
While Anderson MIT has been observed directly in
atomic matter waves [
13
], experimental observation of
multifractality remains challenging [
23
26
]. In particular,
there exists to our knowledge no direct experimental
observation of dynamical multifractality, i.e. manifestation
of multifractality through transport properties (e.g. power-
law decay of the return probability [27,28]).
Another celebrated wave interference effect is the co-
herent backscattering (CBS). It describes the doubling
of the scattering probability (with respect to incoherent
classical contribution) of an incident plane wave with
wave vector
k0
, in the backward direction
k0
. Coherent
backscattering has been observed in many experimental
situations: with light [
29
33
], acoustic waves [
34
,
35
], seis-
mic waves [
36
] and cold atoms [
37
,
38
]. Recently, it was
demonstrated that in the presence of AL a new robust scat-
tering effect emerges [
39
46
], namely the doubling of the
scattering probability in the forward direction +
k0
. This
phenomenon, which appears at long times, was dubbed
coherent forward scattering (CFS). CBS and CFS actually
have a distinct origin: CBS comes from pair interference
of time-reversed paths (and thus requires time-reversal
symmetry), while CFS is present even in the absence of
time-reversal symmetry [
39
,
40
]. From an experimental
point of view, CFS has recently been observed with cold
atoms [38].
In this work, we discuss the fate of CFS at the critical
point of a disorder-driven transition with multifractal
eigenstates. This problem was first addressed for a bulk
3d Anderson lattice [
44
], for which it was shown that CFS
survives at the transition, with however a scattering prob-
ability smaller than in the localized phase. More precisely,
arXiv:2210.04796v1 [cond-mat.dis-nn] 10 Oct 2022
2
Figure 1. CFS contrast Λ
N
(
k, t
;
E
)defined by
(40)
in critical
disordered systems.
k0
is the wave vector of the incident plane
wave and
Dq
are the multifractal dimensions of the eigenstates.
(a) In systems of infinite size
N→ ∞
, the emergence of the
CFS peak as a function of time is governed by the nonergodic
properties of multifractal eigenstates. The CFS wings decay
asymptotically like (
|kk0|t1/d
)
D2
, see Eqs.
(62)
and
(63)
,
while the CFS peak height grows algebraically in time like
tD2/d
and finally reaches the compressibility value
χ
= 1
D1
d
in the long-time limit
t→ ∞
, see Eq.
(58)
. (b) For systems
of finite size
N
, the long-time dynamics of the CFS peak
is governed by the box boundaries. The CFS peak height
reaches 1
αND2
for
t→ ∞
with
α
some numerical factor,
see Eq.
(56)
. The wings of the CFS peak are then described
by Eq. (54).
it was conjectured from numerical evidence that, instead
of a doubling of the classical incoherent contribution, the
forward scattering probability corresponds to a multipli-
cation by a factor (2
D1/d
), with
d
the dimension of
the system and
D1
the information dimension. In our
previous study [
46
], we gave scaling arguments that cor-
roborate this conjecture, backed by numerical simulations
on the Ruijsenaars-Schneider ensemble, a Floquet system
with critical disorder and tunable multifractal dimensions.
We also studied CFS at the transition in finite-size sys-
tems, unveiling a new regime, where CFS properties have
finite-size scaling related to the multifractal dimension
D2[46].
This article is based on the approach developed in our
previous work [
46
] and, somehow, in the spirit of the
random matrix theory point of view discussed in [
41
]. In
particular, we give a complete description of the dynam-
ics of CFS peak in critical disordered systems, including
height and shape of the scattering probability, in two
distinct dynamical regimes. Our findings are summa-
rized in the sketch in Fig. 1. In particular, we present
new links between CFS dynamics and the multifractal
dimension
D2
, that are relevant for most experimental
situations. Our analytical predictions are verified on three
different critical disordered models with multifractal eigen-
states: Power-Law Random Banded Matrices (PRBM),
Ruijsenaars-Schneider ensemble (RS) and unitary three-
dimensional random Kicked Rotor (3DKR). Our predic-
tions are also corroborated by perturbative expansions
for RS and PRBM models in the strong multifractality
regime. These results pave the way to a direct observa-
tion of a dynamical manifestation of multifractality in a
critical disordered system.
II. CRITICAL DISORDERED MODELS
As explained, in the following, our predictions will
be compared to numerical simulations on three different
models. All of them can be mapped onto the generalized
d
-dimensional Anderson model, defined by the following
tight-binding Hamiltonian
ˆ
H=X
n
εn|nihn|+X
n6=m
tnm |nihm|,(2)
where
|ni
are the lattice site states,
εn
the on-site ener-
gies and
tnm
the hopping between two sites at distance
|nm|
. Both
εn
and
tnm
can be considered arbitrary
random variables, whose exact properties will depend on
the system considered (see Table I). We will be interested
in finite-size effects, and will consider a system with linear
size N, i.e. with a total number of sites equal to Nd.
MIT in the generalized Anderson model
(2)
has been
intensively studied (see [
18
,
47
] and references therein).
The three relevant parameters are the spatial dimension
d
of the lattice, the range of the hopping
tnm
, and existence
of correlations in the random entries of the Hamiltonian.
We recall here some well established facts: (i) in the
absence of disorder correlations and if
h|tnm|i
decay faster
than 1
/|nm|d
, Anderson transition only occurs for
d >
2; (ii) in the absence of disorder correlations, critical
eigenstates can appear if
h|tnm|i
decay as fast as 1
/|n
m|d
; (iii) correlations in diagonal disorder
εn
weaken
localization while correlations in off-diagonal disorder
tnm can favour localization.
We now discuss the characteristics and properties of
the different models we used, as well as their link with
the Anderson model (2). A summary is given in Table I.
A. Power-Law Random Banded Matrices (PRBM)
Power-law random banded matrices were first intro-
duced in [
48
]. They were inspired from earlier random
banded matrix ensembles with exponential decay describ-
ing the transition from integrability to chaos [
49
]. The
PRBM model is defined by symmetric or Hermitian matri-
ces whose elements are identical independently distributed
(i.i.d.) Gaussian random variables with zero mean and
variance decreasing as a power law with the distance
from the diagonal. The critical PRBM model corresponds
to an Anderson model
(2)
with random long-range hop-
ping whose variance decays as the inverse of the distance
between sites.
More precisely, let
N
(
µ, σ
)be a Gaussian distribution
of mean
µ
and standard deviation
σ
. In the following
we use the version of PRBM considered in [
17
,
50
], with
periodic boundary conditions, where for
N×N
matrices
3
Model PRBM RS 3DKR
Tunable multifractal dimensions DqYes with b[0,[Yes with a[0,[No
Type Hamiltonian Floquet Floquet
Energy dependent properties Yes No No
Hopping range tnLong-range 1/n Long-range 1/n Short-range (exponential decay)
Dimension d= 1 d= 1 d= 3
Direct (disorder) space Position Momentum Momentum
Table I. Summary of some of the main properties of three models considered in this article (see text pour more details).
diagonal entries
εn
are i.i.d. with distribution
N
(0
,
1),
and real and imaginary parts of the off-diagonal entries
tmn are i.i.d. with distribution N0, σnm/2,
σ2
nm =1 + sin2(π|nm|/N )
(/N)21
.(3)
In particular we have
ph|tnm|2i
=
σnm
, which scales as
1/|nm|for b |nm|  N.
The density of states is defined as
ρ(E) = h1
NdX
α
δ(Eωα)i,(4)
which for this model gives
ρPRBM(E) = (1
2πexpE2
2b1,
1
224E2b1.(5)
Eigenvectors are multifractal, and their multifractal di-
mensions
Dq
, which depend on both
E
and parameter
b
, can be analytically computed [
18
,
50
]. Parameter
b
makes it possible to explore the whole range of multifrac-
tality regime: the weak multifractality regime
Dq
1is
reached for
b→ ∞
and the strong multifractality regime
Dq
0is reached for
b
0. All numerical data pre-
sented in this work are performed at the center of the
band E= 0.
B. Ruijsenaars-Schneider model
Let us consider the following deterministic kicked rotor
model [51,52]
ˆ
H=τˆp2
2+V(ˆx)X
n
δ(tn),(6)
with a 2
π
-periodic sawtooth potential
V
(
x
) =
ax
for
π < x < π
, and where
τ
is a constant parameter. As
a direct consequence of spatial periodicity of
V
(
x
), mo-
menta only take quantized values
pn
= 0
,±
1
,±
2
, . . .
(here
~
= 1). Additionally, we consider a truncated basis
in
p
space, with periodic boundary conditions, so that
the total number of momenta states
|pni
accessible is
N
.
This implies that position basis is also discretized (
xk
are
separated by intervals 2π/N, with kan integer).
It is well-known that kicked Hamiltonians such as
(6)
can be mapped onto the Anderson models
(2)
[
53
,
54
].
The
N
quantized plane waves
|pni
then play the role
of lattice site states
|ni
. The mapping is given (for an
eigenvector of the Floquet operator with eigenphase e
)
by
εn= tanω/2τn2/4,(7)
tnm =Zπ
π
dx
2πtan[V(x)/2]eix(mn),(8)
where the on-site energy
εn
takes evenly distributed
pseudo-random value, provided
τ
is sufficiently irrational.
As a consequence of the Fourier transform relation in
Eq.
(8)
, discontinuity of the sawtooth potential
V
(
x
)cre-
ates a long-range decay of the couplings
tnm
1
/|nm|
and actually induces multifractal eigenstates.
The Ruijsenaars-Schneider (RS) model was introduced
in the context of classical mechanics [
55
57
]. Its quantum
properties were studied in [
58
60
]. It is defined (for an
arbitrary real parameter
a
) by the Floquet operator of
the Hamiltonian (6) (with truncated basis in pspace)
ˆ
U= eˆpeiaˆx,(9)
where the deterministic kinetic phase has been replaced by
random phases
ϕˆp
(consequently the on-site energies
εn
in Eq.
(8)
are truly uncorrelated), and
x
is taken modulo
2π[61].
Importantly, unlike for PRBM, eigenstate properties
of the RS matrix ensemble do not depend on their quasi-
energy. In particular it has a flat density of states
ρRS(E) = 1
2π.(10)
Eigenvectors are multifractal; the multifractal dimensions
can be derived in certain perturbation regimes, and only
depend on the parameter
a
[
62
65
]. This parameter
a
allows us to explore the whole range of multifractality
regimes : the weak multifractality regime
Dq
1is
reached for
a
1and the strong multifractality regime
Dq0is reached for a0.
4
C. 3d Random Kicked Rotor (3DKR)
Our three-dimensional (3d) model is the deterministic
kicked rotor, defined by the following Hamiltonian [66]
ˆ
H=τxp2
x
2+τyp2
y
2+τzp2
z
2+V(q)X
n
δ(tn),(11)
where
τi
are constant parameters and the spatial potential
writes
V
(
q
) =
KV
(
x
)
V
(
y
)
V
(
z
),
K
the kick strength, with
V(x) = 2
2cos x+1
2sin 2x,(12)
so that the system breaks the time-reversal symmetry
[45].
As previously stated, the Hamiltonian
(11)
can be
mapped onto the 3d Anderson model
(2)
. For a given
eigenstate of the system with eigenphase e
, this mapping
writes
εn= tanω/2τxn2
x/4τyn2
y/4τzn2
z/4,(13)
tnm =ZZ π
π
dq
(2π)3tan[V(q)/2]eiq·(mn),(14)
where energies
εn
take pseudo-random values (provided
that (
τx, τy, τz
)are incommensurate numbers), and where
hopping terms
tnm
decay exponentially fast with distance
between sites |nm|[66].
The 3d random Kicked Rotor (3DKR) that we consider
in the following corresponds to the Floquet operator of
Hamiltonian (11)
ˆ
U= eˆ
peiV (ˆ
q),(15)
where deterministic kinetic phases are replaced by uni-
formly distributed random phases
φp
(this implies in
particular that energies εnin Eq. (13) are uncorrelated).
The 3DKR can be seen as the Floquet counterpart
of the usual 3d unitary Anderson Model. In particular
it undergoes an Anderson transition monitored by the
parameter
K
(that is related to the hopping intensity).
Using techniques inspired by [
66
,
67
], we found that the
critical value is
Kc
1
.
58 (see Appendix A). However,
unlike the 3d Anderson model, this unitary counterpart
has a flat density of states
ρ3DKR(E) = 1
2π,(16)
and no mobility edge.
Furthermore, we assumed that 3DKR has the same mul-
tifractal dimensions than the corresponding unitary 3d
Anderson model, because it belongs to the same universal-
ity class. The values that were determined in [
68
] (using
the same techniques as in [
69
,
70
]) are
D1
= 1
.
912
±
0
.
007
and D2= 1.165 ±0.015.
III. GENERAL FRAMEWORK FOR THE
STUDY OF CFS IN CRITICALLY DISORDERED
SYSTEMS
A. Eigenstates and time propagator
In the following, we will analytically and numerically
address CFS in critical disordered systems within a very
general framework, including both Floquet and Hamil-
tonian cases. The numerical methods are presented in
Appendix B.
For the sake of clarity, we use a common notation:
|φαi
refer to eigenstates (or Floquet modes) with energy (or
quasienergy)
ωα
. The time propagator of the system then
writes
ˆ
U(t) = X
α
eαt|φαihφα|,(17)
where time will be considered a continuous variable. In
particular, we use the following convention and notation
for the temporal Fourier transform:
f(ω) = Z
−∞
dt f(t)et, f(t) = Z
−∞
dω
2πf(ω)et.
(18)
B. Direct and reciprocal spaces
As illustrated by the models introduced above, in
generic critical disordered systems disorder can be present
either in position space (e.g. PRBM, Anderson model)
or momentum space (e.g. 3DKR, RS). From now on, we
refer to the basis where disorder is present (labeled with
kets
|ni
) as the direct space and to its Fourier-conjugated
basis (labeled with kets
|ki
) as the reciprocal space. This
distinction is particularly important because multifractal-
ity of eigenstates is a basis-dependent property that only
appears in direct space, where disorder is present, while
CFS is an interference effect taking place in reciprocal
space.
Importantly, we choose to use standard notations
of spatially disordered lattice systems, as in Eq.
(2)
.
For a
d
-dimensional system, direct space is spanned
by discrete lattice sites states
|ni
=
|n1, . . . , ndi
(
ni
=
N/
2 + 1
, . . . N/
2) (
N
will be considered even). The
dimension of the associated Hilbert space is
Nd
. Con-
sequently, the reciprocal space is spanned by a basis
|ki
=
|k1, . . . , kdi
(where
ki
=
±π
N,±3π
N··· ± (N1)π
N
).
We also choose the following convention for the change of
basis (see Appendix Cfor details)
φα(k) = X
n]N/2,N/2]d
φα(n)eik·n,(19)
φα(n) = 1
NdX
k]π]d
φα(k)eik·n,(20)
5
so that in the limit
N→ ∞
the system tends to a infinite-
size discrete lattice, that is,
φα(k)
N→∞
X
n1=−∞ ···
X
nd=−∞
φα(n)eik·n,(21)
φα(n)
N→∞ ZZ π
π
ddk
(2π)dφα(k)eik·n.(22)
We insist that for 3DKR and RS models, direct space
is the momentum space. For instance for the RS model
the basis
|ni
corresponds to plane waves with discrete
momenta
p
=
n~
(with
~
= 1) because of spatial 2
π
-
periodicity of kicked Hamiltonians. Consequently, the
reciprocal space corresponds to position space, so that
|ki
corresponds to discrete positions
xk
=
±π
N,±3π
N···±
(N1)π
N
. Spatial discretization comes from the imposed
periodic boundary conditions in the truncated momentum
basis, so that the linear system size in direct space is N.
C. Form factor and level compressibility
Previous studies [
39
46
] found that CFS dynamics
could be related to the form factor. We will show that it
is the same in critical disordered systems. We recall some
definitions that will be useful in forthcoming calculations.
1. Form factor
The form factor is the Fourier transform of the two-
point energy correlator; it is usually defined as
KN(t) = 1
NdhX
α,β
eαβ ti,(23)
with ωα,β =ωβωα. It can be rewritten as
KN(t) = ZdE ρ(E)KN(t;E)(24)
with
KN(t;E) = 1
Ndρ(E)*X
αβ
eαβ tδEωα+ωβ
2+.
(25)
The component
KN
(
t
;
E
)of the form factor can be in-
terpreted as coming from contributions of all interfering
pairs of states whose average energy is
E
. In order to
lighten forthcoming calculations, we introduce the follow-
ing implicit notation
hX
α,β
. . .iE*1
ρ(E)X
α,β
δEωα+ωβ
2. . .+,(26)
hf(ωα)iE*1
ρ(E)X
α
δ(Eωα)f(ωα)+,(27)
so that KN(t;E)writes
KN(t;E) = 1
NdhX
α,β
eαβ tiE.(28)
2. Compressibility and link to multifractal dimensions
The level compressibility χis defined as
χ= lim
t/N d0KN(t;E).(29)
It is a measure of long-range correlations in the spectrum.
It estimates how much the variance of the number of
states in a given energy window scales with the size of
the window. For usual random matrices (GOE, GUE. . . )
χ= 0, while for Poisson statistics χ= 1.
For critical systems that have intermediate statistics,
the level compressibility lies in between 0
<χ<
1[
27
].
It was proposed that
χ
could actually be related to mul-
tifractal dimension
D2
via
χ
= 1
D2/2d
[
27
,
71
], but
it was later observed that this relation fails in the weak
multifractal regime. Another relation was then conjec-
tured [
60
], relating
χ
to the information dimension
D1
χ= 1 D1
d,(30)
and has since been verified in many different systems [
63
,
7274] (see also Appendix B).
The information dimension
D1
appearing in Eq.
(30)
is defined trough the asymptotic expansion of Eq.
(1)
in
the limit q1
hPnδ(Eωα)|φα(n)|2ln |φα(n)|2i
hPαδ(Eωα)iD1ln N, (31)
and can be seen as the Shannon entropy of eigenstates
|φα(n)|2.
D. Energy decomposition and contrast definition
CFS is an interference effect that appears when the
system is initially prepared in a state localized in recip-
rocal space,
|ψ(t= 0)i
=
1
Nd/2|k0i
(our Fourier trans-
form and normalization conventions are listed in Ap-
pendix C). The observable of interest is the disorder
averaged scattering probability in direction
k
, defined
as
n
(
k, t
) =
1
Ndh|hk|ˆ
U(t)|k0i|2i
. Using
(17)
, it can be
expanded over eigenstates as
n(k, t) = 1
NdhX
α,β
eαβ tφα(k)φ?
α(k0)φβ(k0)φ?
β(k)i.
(32)
摘要:

CoherentforwardscatteringasarobustprobeofmultifractalityincriticaldisorderedmediaMaximeMartinez,1GabrielLemarié,1,2,3BertrandGeorgeot,1ChristianMiniatura,2,3,4,5,6andOlivierGiraud71LaboratoiredePhysiqueThéorique,UniversitédeToulouse,CNRS,UPS,France2MajuLab,CNRS-UCA-SU-NUS-NTUInternationalJointResear...

展开>> 收起<<
Coherent forward scattering as a robust probe of multifractality in critical disordered media Maxime Martinez1Gabriel Lemarié123Bertrand Georgeot1Christian Miniatura23456and Olivier Giraud7.pdf

共24页,预览5页

还剩页未读, 继续阅读

声明:本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。玖贝云文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知玖贝云文库,我们立即给予删除!
分类:图书资源 价格:10玖币 属性:24 页 大小:3.72MB 格式:PDF 时间:2025-04-27

开通VIP享超值会员特权

  • 多端同步记录
  • 高速下载文档
  • 免费文档工具
  • 分享文档赚钱
  • 每日登录抽奖
  • 优质衍生服务
/ 24
客服
关注