of 5
Origin of low sodium capacity in graphite and generally
weak substrate binding of Na and Mg among alkali
and alkaline earth metals
Yuanyue Liu
a,b,1
, Boris V. Merinov
a
, and William A. Goddard III
a,1
a
Materials and Process Simulation Center, California Institute of Technology, Pasadena, CA 91125; and
b
The Resnick Sustainability Institute, California
Institute of Technology, Pasadena, CA 91125
Contributed by William A. Goddard III, February 19, 2016 (sent for review January 25, 2016; reviewed by Yi Cui and Michael L. Klein)
It is well known that graphite has a low capacity for Na but a high
capacity for other alkali metals. The growing interest in alternative
cation batteries beyond Li makes it particularly important to
elucidate the origin of this behavior, which is not well understood.
In examining this question, we find a quite general phenomenon:
among the alkali and alkaline earth metals, Na and Mg generally
have the weakest chemical binding to a given substrate, compared
with the other elements in the same column of the periodic table.
We demonstrate this with quantum mechanics calculations for a
wide range of substrate materials (not limited to C) covering a
variety of structures and chemical compositions. The phenomenon
arises from the competition between trends in the ionization en-
ergy and the ion
substrate coupling, down the columns of the
periodic table. Consequently, the cathodic voltage for Na and
Mg is expected to be lower than those for other metals in the
same column. This generality provides a basis for analyzing the
binding of alkali and alkaline earth metal atoms over a broad
range of systems.
energy storage
|
Na-ion battery
|
Mg-ion battery
|
quantum-mechanical calculations
D
evelopment of alternative cation batteries beyond Li could
solve issues related with cost, stability, and other perfor-
mance characteristics (1
3). Na is an obvious candidate, but its
storage in graphite (the commercial anode for Li-ion battery)
is rather poor, with an electrochemical capacity of less than
35 mAh/g (1, 4
7). Surprisingly, other alkali metals have a high
capacity (approximately hundreds) in graphite (8). To form a
basis for improving the battery performance, we seek to under-
stand the anomalously low capacity for Na.
One explanation in the literature for the low Na capacity is as
follows: Na intercalation expands graphite from its favorable
interlayer spacing, by a greater amount than Li, leading to a
higher strain energy for graphite and, therefore, a less favorable
formation energy for the Na-graphite compound compared with
the Li analog (9
11). However, this explanation would suggest
that graphite should have a low capacity for K, Rb, and Cs be-
cause of their even larger size, which is in stark contrast to the
experimentally observed capacity dramatically higher than that
for Na. This inconsistency calls for a revisit of the origin of the
low Na capacity in graphite.
In this work, we use quantum-mechanical methods to show
that the Na anomaly has its roots in a general phenomenon:
among alkali metals (denoted as M) and alkaline earth metals
(denoted as EM), Na and Mg generally have the weakest binding
to a given substrate, independent of variations in substrate
structure and chemistry. This phenomenon results from the
competition between the ionization of the metal atom and the ion
substrate coupling, which have opposite trends along the col-
umns. The universality of this phenomenon provides the basis
for analyzing trends in binding of alkali and alkaline earth
metals over a broad range of systems, and offers guidance for
designing improved systems.
Our density functional theory (DFT) calculations used the
Vienna ab initio simulation package (VASP) (12, 13) with pro-
jector augmented wave (PAW) pseudopotentials (14, 15). The
Perdew
Burke
Ernzerhof (PBE) exchange-correlation func-
tional (16) including van der Waals corrections (DFT-D3) (17) is
used to model the graphite and its compounds, whereas in other
cases where the van der Waals interaction is insignificant, we use
PBE only. The plane-wave cutoff energy is 400 eV, with suffi-
cient Monkhorst
Pack sampled
k
points (18) (for example, 15
×
15
×
7 for graphite). All structures are fully relaxed until the final
force on each atom becomes less than 0.01 eV/Å. We used a
single M or EM atom with 6
×
6 unit cells of each substrate
material to model the binding with graphene, its derivatives,
MoS
2
, SnS
2
, and TiS
2
, whereas a 4
×
4 cell is applied for V
2
O
5
.
For Pt(111), we used a slab consisting of 6
×
6 surface cells, and
three layers, with the bottom layer fixed in the same plane.
We consider first the case of graphite. Fig. 1 shows the cal-
culated formation energy (
E
f
) of M-graphite compounds, where
M
=
Li, Na, K, Rb, and Cs. The
E
f
is defined as follows:
E
f
=
½
E
ð
tot
Þ
n
M
E
ð
bulk
M
Þ
n
C
E
ð
graphite
Þ
=
ð
n
M
+
n
C
Þ
,
[1]
where
E
(tot) is the total energy of compound,
E
(bulk M) is the
energy per M atom for the bulk metal,
E
(graphite) is the energy
of C in graphite, and
n
M
and
n
C
are the number of M and C
atoms in the compound. Here, we focus on MC
6
and MC
8
be-
cause these stoichiometries are commonly found in non-Na com-
pounds. We find that the
E
f
follows the order: Na
>
Li
>
K
>
Rb
>
Cs, where all Ms except for Na have negative
E
f
with graphite, a
result consistent with calculat
ions in the literature that use
different method (19). Thus, Na-graphite compounds with high
Significance
The growing demand for energy storage urges the develop-
ment of alternative cation batteries, which calls for a system-
atic understanding of binding energetics. We discover a
general phenomenon for binding of alkali and alkaline earth
metal atoms with substrates, which is explained in a unified
picture of chemical bonding. This allows us to solve the long-
standing puzzle of low Na capacity in graphite and predict the
trends of battery voltages, and also forms a basis for analyzing
the binding of alkali and alkaline earth metal atoms over a
broad range of systems.
Author contributions: Y.L. and W.A.G. designed research; Y.L., B.V.M., and W.A.G. per-
formed research; Y.L., B.V.M., and W.A.G. analyzed data; and Y.L., B.V.M., and W.A.G.
wrote the paper.
Reviewers: Y.C., Stanford University; and M.L.K., Temple University.
The authors declare no conflict of interest.
1
To whom correspondence may be addressed. Email: wag@wag.caltech.edu or yuanyue.
liu.microman@gmail.com.
This article contains supporting information online at
www.pnas.org/lookup/suppl/doi:10.
1073/pnas.1602473113/-/DCSupplemental
.
www.pnas.org/cgi/doi/10.1073/pnas.1602473113
PNAS
|
April 5, 2016
|
vol. 113
|
no. 14
|
3735
3739
APPLIED PHYSICAL
SCIENCES
Na contents are not thermodynamically stable, in contrast to the
other four M-graphite compounds, in agreement with the exper-
imentally observed low Na capacity.
To understand why the
E
f
has a maximum at Na, we partition
the formation process for M-graphite compounds into three
steps, as illustrated in Fig. 2
A
. First, the bulk metal is evaporated
to form isolated atoms, with an energy cost of
E
d
(i.e., the cohesive
energy). Second, the graphite crystal is strained to the configura-
tion identical to that of the M-intercalated graphite, with an en-
ergy cost of
E
s
. This straining includes both interlayer expansion
and in-plane stretching. Third, the M atoms are intercalated into
the strained graphite, with an energy drop by
E
b
(to be consistent
with the other terms, we define
E
b
=
E
b
; a stronger binding
corresponds to a lower
E
b
). According to Hess
slaw:
E
f
=
E
d
+
E
s
+
E
b
.
[2]
This analysis helps identify the dominating contribution. Al-
though some terms are difficult to measure experimentally (such
as
E
s
and
E
b
), all are calculated easily using DFT.
As shown in Fig. 2
B
, neither
E
d
,
E
s
, nor their combination has
a trend similar to
E
f
, which suggests that they do not embody the
origin of the low Na capacity. For both MC
8
and MC
8
,
E
s
in-
creases monotonically as M moves down the periodic table due
to the increasing size of M atoms, whereas
E
d
decreases as a
result of weakened cohesion. The combination of
E
s
and
E
d
shows a monotonic drop, indicating that the change of
E
d
overwhelms that of
E
s
. Nevertheless, none of them explains the
maximum of
E
f
at Na.
On the other hand,
E
b
exhibits a maximum at Na, similar to
that of
E
f
. This suggests that the low Na capacity is related di-
rectly to
E
b
. In particular, compared with Li, the Na binding is
so weak that it exceeds the decrease of
E
s
+
E
d
, making
E
f
higher
for Na.
This weaker binding of Na compared with Li has been
reported for other intercalation compounds, which has been
proposed to account for the observed lower cathodic voltage
(20). Here we find that, of all of the five alkali metals, Na always
has the weakest binding for any given substrate. We first examine
dilute M adsorption on graphene (Fig. 3
A
), which shows a
maximum of
E
b
at Na. Then we modify the adsorption sites by
incorporating structural defects or foreign atoms. Remarkably,
Na always shows the weakest binding (Fig. 3
A
). We continue the
test by considering other 2D non-C materials that have been
tested for batteries, namely MoS
2
, TiS
2
, SnS
2
, and V
2
O
5
(Fig.
3
B
). We then extend this test to the surface of a typical bulk
materials Pt(111) (Fig. 3
B
; as commonly found in Pt-based ca-
talysis in alkaline solution). In all cases, Na has the weakest
binding among all five alkali metals, independent of the detailed
substrate chemistry/structure. This general phenomenon calls for
a unified explanation.
E
b
is the energy change when an M atom moves from the
vacuum to the binding site of the substrate. We consider this
process to first involve ionization of M by transferring its charge
Fig. 1.
Calculated formation energies (Eq.
1
) of alkali metal (M)-graphite
compounds. Note that, in contrast to other Ms, NaC
6
and NaC
8
have a pos-
itive formation energy.
Fig. 2.
(
A
) Partition of the formation process of M-graphite compound into
separate steps. C, brown; M, pink. (
B
) The energetics of each step, relative to
those of Li (see Eq.
2
and the related text).
Binding
means the reverse of
binding energy, i.e.,
E
b
in Eq.
2
. Note that only
E
b
shows a trend similar
with
E
f
.
3736
|
www.pnas.org/cgi/doi/10.1073/pnas.1602473113
Liu et al.
to the substrate, with an energy change by
E
ion
. This is followed
by the coupling of the cation to the substrate (negatively charged)
with an energy decrease of
E
cp
(which includes the electrostatic
and other quantum-mechanical interactions), as illustrated in Fig.
4
A
. Therefore,
E
b
=
E
ion
+
E
cp
.
[3]
As M moves down the periodic table, the ionization potential
decreases: 5.4 (Li), 5.1 (Na), 4.3 (K), 4.2 (Rb), and 3.9 (Cs),
which favors the binding and results in a decrease of
E
ion
. Note
that there is an abrupt drop in the ionization potential between
Na and K. However, at the same time, the distance between the
cation and the substrate becomes larger, which weakens their
coupling, leading to an increase of
E
cp
. To quantify this compe-
tition, Fig. 4
B
shows the
E
ion
and
E
cp
with respect to those of Li
on the same substrate. The relative
E
ion
is approximated by the
difference in the atomic ionization potential (IP), and the rela-
tive
E
cp
is then derived from Eq.
3
. Because the
E
cp
increases
smoothly from Li to Cs while the
E
ion
drops dramatically from
Na to K (in other words, the coupling strength increases smoothly
from Cs to Li while the ionization cost drops dramatically from K
to Na), we obtain maximum of
E
b
for Na.
The
E
cp
is a general term that includes various kinds of in-
teractions that are difficult to separate. However, for the cases of
M adsorption on metals [e.g., graphene and Pt(111)], one could
expect the difference in
E
cp
of different Ms is dominated by the
electrostatic contribution, which is approximately
14.38/(2*
d
),
where
d
is the distance between the cation and the substrate
according to the image charge method. Therefore, in these cases:
Δ
E
b
=
Δ
E
ion
+
Δ
E
cp
Δ
½
14.38
=
ð
2
p
d
Þ
+
Δ
IP.
[4]
As shown in Fig. 4
C
, the trend of
E
b,
given by the equation is
similar to that calculated by using DFT, which validates that the
competition between
E
ion
and
E
cp
is the reason for the maximum
of
E
b
at Na. Indeed, this weak Na binding is also found in
diatomic molecules M-X (where X
=
F, Cl, Br, I, and OH)
(21), which can be explained similarly.
It is interesting to consider whether the nonmonotonic trend
of the
E
b
is present in other columns. Based on the above ex-
planations, we anticipate that this can be observed in the col-
umns where
E
ion
and
E
cp
have a reverse trend when moving
down the periodic table. The
E
cp
perhaps always increases as the
atomic size gets larger; however,
E
ion
does not always decrease,
given the fact that the IP or electron affinity is nonmonotonic for
Fig. 3.
E
b
(the negative of the binding energy) for alkali metals binding to various substrate materials, relative to that of Li (the absolute values can be
found in
Supporting Information
). (
A
) Graphene and its derivatives (with defects or substitutional foreign atoms); (
B
) non-C materials. Note that Na always
has the weakest binding among alkali metals.
Liu et al.
PNAS
|
April 5, 2016
|
vol. 113
|
no. 14
|
3737
APPLIED PHYSICAL
SCIENCES
most groups. In fact, only the first two groups have a notably
decreasing IP. Therefore, one may expect a similar phenomenon
occurs for the EM group. Indeed, our calculations show that,
among the EM elements, Mg generally has the weakest binding
(Fig. 5) with a given substrate. This is consistent with the ex-
perimental fact that Mg has a low capacity in graphite, similar to
the case of Na (8, 22). Note that Be and Mg are only physisorbed
on pristine or nitrogen-doped graphene with the
E
b
approxi-
mately of
30 meV/EM, significantly weaker than other EMs.
This is because the work function of pristine or nitrogen-doped
graphene is too low to allow for a charge transfer from Be/Mg, as
shown by the band structures in
Supporting Information
. These
systems are perhaps not practically interesting as the physisorbed
adatoms could easily detach or cluster. It should also be noted
that, Be tends to have a stronger covalence than other EMs, due
to its high IP and small size. This might be the reason for the
significantly enhanced binding of Be with the monovacancy in
graphene, in which case the
E
cp
contributes a large energy drop.
For cathode materials, the weak binding with metal atoms
results in a low cathodic voltage. Therefore, we anticipate that
Fig. 4.
(
A
) Schematic illustration of the binding between the alkali metal (M) and the substrate. (
B
) Evolution of the ionization energy (
E
ion
) and the coupling
energy (
E
cp
) as M moves down in the periodic table. (
C
) Comparison of the calculated
E
b
(relative to that of Li) with that based on the model of Eq.
4
.
3738
|
www.pnas.org/cgi/doi/10.1073/pnas.1602473113
Liu et al.
Na and Mg have a low cathodic voltage compared with other
metals in the same columns. Indeed, the Na case has been ver-
ified by explicitly calculating the cathodic voltage for various
intercalation compounds (20). On the other hand, a weak binding
with anode is desired to enhance the voltage when connected with
cathode. However, this usually sacrifices capacity, as seen in the
case of Na in graphite. To improve the Na capacity, it is necessary
to reduce
E
f
. This can be achieved by using prestrained graphite (i.
e., reducing the
E
s
term in Eq.
2
) with expanded interlayer spacing
through intercalation of some other species, as has been demon-
strated experimentally (9). Our calculations show that the opti-
mum
E
f
for Na-graphite is reached when the graphite interlayer
distance is expanded to 4.3 Å, providing a target value for ex-
perimental design. Indeed, such an expanded graphite using or-
ganic pillars has been shown computationally to provide a very
promising hydrogen storage of 6.5 wt% at room temperature,
meeting the Department of Energy requirement (23). Alternatively,
the
E
f
might be reduced by enhancing the binding, which might be
achieved through incorporating defects (24, 25).
In summary, we use quantum-mechanical calculations, to find a
general phenomenon
among alkali and alkaline earth metals,
Na and Mg generally have the weakest binding for a given sub-
strate. We show that this results from the competition between the
ionization of the metal atom, and the ion
substrate coupling. This
finding elucidates the origin of the low Na capacity in graphite,
predicts the voltage trends for alkali and alkaline earth metal ion
batteries, and provides a basis for analyzing the binding of alkali
and alkaline earth metal atoms in a broad range of systems.
ACKNOWLEDGMENTS.
Y.L. thanks Drs. Brandon Wood, Suhuai Wei, and
Jiayu Wan for helpful discussions and Brandon Wood for providing access to
the Lawrence Livermore National Laboratory computational resources,
which were used for some of the computations (supported under the Lab-
oratory Directed Research and Development Program). Most of the calcula-
tions were performed on National Energy Research Scientific Computing
Center, a Department of Energy (DOE) Office of Science User Facility sup-
ported by the Office of Science of the US DOE under Contract DE-AC02-
05CH11231. Y.L. acknowledges the support from Resnick Prize Postdoctoral
Fellowship at Caltech. This research was funded by the Bosch Energy Re-
search Network and by NSF (CBET 1512759).
1. Yabuuchi N, Kubota K, Dahbi M, Komaba S (2014) Research development on sodium-
ion batteries.
Chem Rev
114(23):11636
11682.
2. Saha P, et al. (2014) Rechargeable magnesium battery: Current status and key chal-
lenges for the future.
Prog Mater Sci
66:1
86.
3. Lipson AL, et al. (2015) Rechargeable Ca-ion batteries: A new energy storage system.
Chem Mater
27(24):8442
8447.
4. Ge P, Fouletier M (1988) Electrochemical intercalation of sodium in graphite.
Solid
State Ion
28:1172
1175.
5. Stevens DA, Dahn JR (2001) The mechanisms of lithium and sodium insertion in car-
bon materials.
J Electrochem Soc
148(8):A803
A811.
6. Palomares V, et al. (2012) Na-ion batteries, recent advances and present challenges to
become low cost energy storage systems.
Energy Environ Sci
5(3):5884
5901.
7. Slater MD, Kim D, Lee E, Johnson CS (2013) Sodium-ion batteries.
Adv Funct Mater
23(8):947
958.
8. Dresselhaus MS, Dresselhaus G (1981) Intercalation compounds of graphite.
Adv Phys
30(2):139
326.
9. Wen Y, et al. (2014) Expanded graphite as superior anode for sodium-ion batteries.
Nat Commun
5:4033.
10. Cao Y, et al. (2012) Sodium ion insertion in hollow carbon nanowires for battery
applications.
Nano Lett
12(7):3783
3787.
11. Jache B, Adelhelm P (2014) Use of graphite as a highly reversible electrode with su-
perior cycle life for sodium-ion batteries by making use of co-intercalation phe-
nomena.
Angew Chem Int Ed Engl
53(38):10169
10173.
12. Kresse G, Hafner J (1993) Ab initio molecular dynamics for liquid metals.
Phys Rev B
Condens Matter
47(1):558
561.
13. Kresse G, Furthmüller J (1996) Efficient iterative schemes for ab initio total-energy
calculations using a plane-wave basis set.
Phys Rev B Condens Matter
54(16):
11169
11186.
14. Blöchl PE (1994) Projector augmented-wave method.
Phys Rev B Condens Matter
50(24):17953
17979.
15. Kresse G, Joubert D (1999) From ultrasoft pseudopotentials to the projector aug-
mented-wave method.
Phys Rev B
59(3):1758
1775.
16. Perdew JP, Burke K, Ernzerhof M (1996) Generalized gradient approximation made
simple.
Phys Rev Lett
77(18):3865
3868.
17. Grimme S, Antony J, Ehrlich S, Krieg H (2010) A consistent and accurate ab initio
parametrization of density functional dispersion correction (DFT-D) for the 94 ele-
ments H-Pu.
J Chem Phys
132(15):154104.
18. Monkhorst HJ, Pack JD (1976) Special points for Brillouin-zone integrations.
Phys Rev
B
13(12):5188
5192.
19. Nobuhara K, Nakayama H, Nose M, Nakanishi S, Iba H (2013) First-principles study of
alkali metal-graphite intercalation compounds.
J Power Sources
243(0):585
587.
20. Ong SP, et al. (2011) Voltage, stability and diffusion barrier differences be-
tween sodium-ion and lithium-ion intercalation materials.
Energy Environ Sci
4(9):
3680
3688.
21. Darwent Bd (1970)
Bond Dissociation Energies in Simple Molecules
(US National Bu-
reau of Standards, Washington, DC).
22. Kawaguchi M, Kurasaki A (2012) Intercalation of magnesium into a graphite-like
layered material of composition BC2N.
Chem Commun (Camb)
48(55):6897
6899.
23. Deng W-Q, Xu X, Goddard WA (2004) New alkali doped pillared carbon materials
designed to achieve practical reversible hydrogen storage for transportation.
Phys
Rev Lett
92(16):166103.
24. Liu Y, Wang YM, Yakobson BI, Wood BC (2014) Assessing carbon-based anodes for
lithium-ion batteries: A universal description of charge-transfer binding.
Phys Rev Lett
113(2):028304.
25. Liu Y, Artyukhov VI, Liu M, Harutyunyan AR, Yakobson BI (2013) Feasibility of lithium
storage on graphene and its derivatives.
J Phys Chem Lett
4(10):1737
1742.
Fig. 5.
E
b
(the negative of the binding energy) for alkaline earth metals binding to various substrate materials, relative to that of Be (the absolute values
can be found in
Supporting Information
). Note that Mg always has the weakest binding.
Liu et al.
PNAS
|
April 5, 2016
|
vol. 113
|
no. 14
|
3739
APPLIED PHYSICAL
SCIENCES