on this, we determine for the first the relative per-
mittivity of Li-GICs, as a function of the SOC, for
the entire functional range of the material.
Beyond the previously outlined interest for
charge-kMC, we believe there to be many more
valuable applications for the relative permittivity
of Li-GICs also at higher levels of the multiscale
simulation hierarchy: one such motivation is un-
derstanding charge gradients, as they occur dur-
ing the fast charging of modern batteries in electric
vehicles, and the chemical pressure which leads to
plating and dendrite formation inside the batter-
ies under certain operational conditions. The latter
phenomenon is typically investigated by means of
continuum simulations like e.g. by Hein et al. [10],
which also rely at least implicitly on knowledge
of the dielectric response. Furthermore, a sim-
ple model of the charge carrier electrostatics could
be used as a physical baseline for otherwise short-
ranged machine learning models or cluster expan-
sions. Another related field is the development of
functional materials based on doped graphite [11].
The dielectric behaviour of Li-GICs (and solid
materials in general) is significantly more com-
plicated than the expression “dielectric constant”
would suggest, and is governed by vastly differ-
ent physics at different frequencies of a probing
field. In the static limit (which this work aims
to investigate – the probing field is essentially the
electrostatic field of the intercalated ion itself), no
periodic movement (beyond thermal fluctuations)
of the electrons is induced. However, this pic-
ture changes in the kHZ range, where the entire
charge carrier density oscillates with the probing
field, causing large polarization and large dielectric
screening. For example, Chung et al. [1] measured
an RP of r=2100 for highly oriented pyrolytic
graphite (HOPG) and even higher ones for other
carbon structures, at 2-10KHz. Moving on to the
GHz regime, a balance is reached where the field
oscillations are too fast for macroscopic bulk cur-
rents to build up, and a situation occurs that is
arguably similar to the static limit and may serve
as comparison for our research. Hotta et al. [4] put
forward a dielectric constant for graphite powder of
r≈15, at 6GHz. Finally, at even higher frequen-
cies beyond THz, the electric field becomes high
enough in energy to excite a significant number of
electrons, again creating a physically different sit-
uation with much lower dielectric screening, which
converges to transparency in the ω→ ∞ limit. A
study by Jellison et al. [12] in the frequency regime
of visible light finds an RP in xy-plane for HOPG
of r=4.21, which, for the previously mentioned
reasons, cannot be used as comparison either and
is expected to serve as a lower bound in the follow-
ing.
It is apparent that there is a glaring lack of
studies investigating the exact property of inter-
est to us, which is – again – the electrostatic di-
electric response of a perfect graphite crystal in
xy-plane, i.e. parallel to the graphene sheets, as
experienced by some internal charge carriers (in
this case Li-ions and vacancies). There are some
studies available on graphene, either on some sub-
strate or quasi-freestanding, with results ranging
from r=2.2−5.0by Elias et al. [13] to r=15.4
by Reed et al. [11], and another study by Bostwick
et al. finding r≈4.4[14], none of which can serve
as direct comparison to our research either. How-
ever, there is a study on bilayer graphene (which
according to our calculations can be compared with
graphite quite well) by Bessler et al. [15], putting
forward an RP of r=6±2. This is likely the most
reliable direct experimental comparison currently
available to us.
In terms of theoretical approaches to determin-
ing the dielectric response of materials, substantial
work has been done on water [16, 17, 18, 19]. There
is also some promising work by Gigli et al. [20] in
the development of an integrated machine learning
model predicting the dielectric response of BaT iO3.
3