Structure and Dynamic Behavior of the Na–Crown Ether Complex in

1H NMR shows two kinds of molecular motions at room temperature: isotropic .... and small impurity signals around 0 ppm were observed in the range bet...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

C: Surfaces, Interfaces, Porous Materials, and Catalysis

Structure and Dynamic Behavior of Na-Crown Ether Complex in the Graphite Layers Studied by DFT and H NMR 1

Kazuma Gotoh, Shinya Kunimitsu, Hanyang Zhang, Michael M Lerner, Keisuke Miyakubo, Takahiro Ueda, Hyung-Jin Kim, Young-Kyu Han, and Hiroyuki Ishida J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b02965 • Publication Date (Web): 02 May 2018 Downloaded from http://pubs.acs.org on May 5, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Structure and Dynamic Behavior of Na-Crown Ether Complex in the Graphite Layers Studied by DFT and 1H NMR Kazuma Gotoh,*†‡ Shinya Kunimitsu,† Hanyang Zhang§, Michael M. Lerner§, Keisuke Miyakubo∥, Takahiro Ueda∥, Hyung-Jin Kim∇, Young-Kyu Han∇and Hiroyuki Ishida† †

Graduate School of Natural Science and Technology, Okayama University, 3-1-1

Tsushima-naka, Okayama 700-8530, Japan ‡

Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University,

Nishikyo-ku, Kyoto 615-8245, Japan §

Department of Chemistry, Oregon State University, Corvallis, OR 97331-4003, USA



Department of Chemistry, Graduate School of Science, Osaka University, Toyonaka, Osaka

560-0043, Japan ∇

Department of Energy and Materials Engineering and Advanced Energy and Electronic

Materials, Dongguk University-Seoul, Seoul 100-715, Korea

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 28

ABSTRACT Diffusion of alkali metals in graphite layers is significant for the chemical and electrochemical properties of graphite intercalation compounds (GICs). Crown ethers co-intercalate into graphite with alkali metal (Na and K) cations and form ternary GICs. The structures and molecular dynamics of 15-crown-5 and 18-crown-6 ether coordinating to Na+ or K+ in GICs were investigated by DFT calculations and 1H solid state NMR analyses. DFT calculations suggest a stacked structure of crown ether - metal complex with some offset. 1H NMR shows two kinds of molecular motions at room temperature; isotropic rotation with molecular diffusion, and axial rotation with fluctuation of the axis. The structure and dynamics of crown ether molecules in GIC galleries are strongly affected by the geometry of the crown ether molecules and the strength of the interaction between alkali metal and ligand molecules.

2 ACS Paragon Plus Environment

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1. INTRODUCTION Graphite is one of the prototypical hosts that can intercalate either cations or anions between its layers when it is reduced or oxidized, respectively.1 Because of their high electric conductivity, superconductivity,2-3 and thermoelectricity,4 graphite intercalation compounds (GICs) have attracted considerable attention from the science and engineering research communities. GICs are particularly significant for electrochemical energy storage in batteries because the binary GIC consisting of lithium and graphite (LiCx: x≥6) is employed as the negative electrode material in lithium ion batteries (LIBs).5-6 Furthermore, the potential use of GICs has been expanding in recent few years. For example, potassium GICs (KCx: x≥8) was evaluated as the anodes in potassium ion batteries,7-12 and GICs were explored in dual ion batteries, with graphite serving as both cathode and anode.13-14 For sodium-ion based batteries, binary GICs have not been considered because low-stage compounds Na (NaCx) do not form by electrochemical reduction.15-17 However, electrochemical co-intercalation and deintercalation of sodium and glyme molecules (symmetric linear glycol ethers) to form ternary GICs can be repeated over 6000 cycles at high rate, suggesting that these may provide useful electrode chemistry for sodium ion batteries.17-24 Additionally, these cells do not appear to generate a solid-electrolyte interphase (SEI) at the anode surface.25 The structure and the dynamic behavior of diglyme and sodium have been investigated using 2H solid state nuclear magnetic resonance (NMR),26 and also discussed by DFT calculation.27 Our group synthesized ternary GICs consisting of sodium, diglyme(d14), and graphite by a chemical reduction reaction, and investigated the dynamics and the coordination structure of molecules using 2H solid state NMR.

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 28

There are many reports on the preparation of ternary GICs containing Na or the other alkali metal ions with molecular organic co-intercalates forming a solvated intercalate structure.1 Most co-intercalates have been amines or ethers, for example, ethylenediamine (en) and tetrahydrofuran, due to their ability to strongly solvate alkali metal cations and their reductive stabilities.28-42 We recently described the preparation of the ternary GICs containing crown ethers (15-crown-5 ether (15c5), 18-crown-6 ether (18c6)) and Na or K cations.43 These GICs are obtained either by reductive intercalation of an alkali metal-en complex followed by co-intercalate exchange or by the direct reaction of graphite with a crown ether, alkali metal, and an electrocatalyst. Powder X-ray diffraction, thermal analysis, and GC/MS data were combined to determine a stage-1 bilayer gallery structure for these products. However, the dynamic behavior of the intercalated metal complex, which is a key to understanding the intercalate orientation, has not been evaluated. Comparison of the coordination structure with linear ether (glymes) GICs can be valuable in elucidating the different intercalate diffusion rates in these GIC galleries. We have been exploring the dynamic structure of organic amines and ethers using 2H solid state NMR.26, 44 However, this method can only be applied to GICs synthesized using deuterated solvents. We here investigate the dynamics of 15c5 and 18c6 in Na-15c5-GIC using 1

H static NMR instead of 2H NMR. The reported compositions of Na-15c5-GIC and

K-15c5-GIC are [Na1.7(15c5)]C21 and [K1.2(15c5)]C23, indicating the mole ratios of Na and K cations to 15c5 ether are 1.7 and 1.2, respectively.43 The energies of some structures models of sodium and crown ethers in graphene layers were estimated using DFT calculation, in order to discuss the coordination structure of crown ethers to alkali metal cations.

4 ACS Paragon Plus Environment

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

2. EXPERIMENTAL METHODS 2.1. Solid State NMR analyses Ternary GIC samples consisting of alkali metals and crown ethers were synthesized according to a previous report.43 We prepared Na-15c5-GIC, K-15c5-GIC, and Na-18c6-GIC by co-intercalate exchange of Na-en-GIC to form Na-15c5, K-15c5 or Na-18c6. The compositions of Na-15c5-GIC and K-15c5-GIC are estimated to be [Na1.7(15c5)]C21 and [K1.2(15c5)]C23, respectively, as has been reported previously,43 whereas the composition of Na-18c6-GIC has not been determined. For Na-18c6-GIC, two GIC phases having different gallery height (∆d) have been reported; ∆d = 0.86 nm (α), and 0.70 nm (β). We confirmed that our Na-18c6-GIC sample consists of phase α mainly, with a minor β phase component. All samples used in the present study were characterized using powder X-ray diffraction (PXRD) and elemental analyses. For 1H NMR measurement, each GIC sample was sealed into a 5.0 mmφ glass sample tube after drying under reduced pressure. 1

H NMR spectra (resonance frequency at 200 MHz) were obtained at the range between

145 K and ambient temperature using a DSX-200 spectrometer (Bruker). A solid echo pulse sequence with 4.0 µs pulse length was applied. For the calculation of the second moment (M2) value of 15c5 and 18c6 molecules, the most stable structures of 15c5 and 18c6 were estimated by DFT calculations at B3LYP/6-31G(d,p) level by Gaussian 09W. The theoretical M2 values about the rigid state and under rotation of 15c5 and 18c6 were estimated using the Van Vleck equation.45-46

23

Na magic-angle spinning (MAS) NMR spectra (resonance frequency at 132.2

MHz) were measured at room temperature using a DD2 spectrometer (Agilent Technologies) 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

with a single pulse sequence of 2.0 µs pulse length. 1M NaCl aqueous solution was used as a reference for the chemical shift. Although differential scanning calorimetry (DSC) was performed using Pyris 1 analyzer (PerkinElmer) to observe possible phase transition in samples with the temperature range between 135 K and 298 K, no DSC signals were observed for these samples. 2.2. First-principles calculations Density functional theory calculations were performed using the Vienna ab initio simulation package (VASP).47-48 We employed the revised Perdew–Burke–Ernzerhof type exchange and correlation functional49-50 combined with the introduction of vdW-DF51 for the non-local correlation part to accurately account for the dispersion interactions. The projector augmented wave (PAW) method52 was used for the ion interaction. The Brillouin zone was sampled using a Γ-centered 3 × 3 × 3 k-point mesh, while the electronic states were smeared using the Methfessel–Paxton scheme with a broadening width of 0.2 eV. The electronic wave functions were expanded in a plane wave basis with a cutoff energy of 520 eV and the atomic relaxation was continued until the Hellmann–Feynman forces acting on the atoms were less than 0.2 eV nm–1. We treated 2p63s1 for Na, 3p64s1 for K, 2s22p2 for C, 2s22p4 for O, and 1s for H as the valence electron configuration. M (Na and K)-crown ether-GIC systems were simulated by 5 × 5 supercell including one graphene sheet consisting of 50 C atoms and M-crown ether complexes. The formation energy of M-crown ether intercalation were defined as Ef = [Etot(2M(15c5)C25) − Etot(C50) − 2Etot(M(15c5))]/2, where Etot(2M(15c5)C25), Etot(C50), and Etot(M(15c5)) are the total energy of M-(15c5)-GIC, graphene sheet, and M-crown ether complex, respectively. 6 ACS Paragon Plus Environment

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3. RESULTS 3.1. Solid State NMR Fig. 1 presents 1H NMR spectra of the Na-15c5-GIC, K-15c5-GIC, and Na-18c6-GIC samples taken between 293 K and 123 K. Although each spectrum doesn't correspond to a single Gaussian or Lorentzian curve, it was able to be fitted by two Gaussian components (Fig. 2) which correspond to two motional states of the crown ether molecule in each GIC sample. The crown ether molecules existing outside of layered structure should be quite few because of small specific surface area of graphite. The two components observed in 1H NMR are ascribed to molecules intercalated in the layered structure. The narrower components (C1) of Na-15c5-GIC (Fig. 1(a)) and K-15c5-GIC (Fig. 1(b)) having the FWHM (full width at half maximum) of ca. 10 ppm predominate at room temperature. The C1 components of Na-15c5-GIC and K-15c5-GIC decreased under 230 and 240 K, respectively, with a slight broadening, and almost disappeared at 157 K. On the other hand, the intensity of broader components (C2) of Na-15c5-GIC and K-15c5-GIC increased under 228 K with considerable broadening. Below 157 K, the spectra consisted of C2 only; the narrower component (C1) disappeared. For Na-18c6-GIC (Fig. 1(c)), C1 almost disappeared below 196 K. Only the broader component (C2) and small impurity signals around 0 ppm were observed in the range between 196 K and 146 K. The signal intensity ratio of C1/C2 and the second moment (M2) of C1 and C2, which are estimated by curve fitting, are shown in Figs. 3(a) and 3(b), respectively. The M2 values of C2 in Na-15c5-GIC and K-15c5-GIC increased from 2.0 Gauss2(G2: 10-2 mT2) to 20.8 G2 and 29.2 G2, respectively, with temperature reduction 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 28

to 149 K. The M2 values of C2 in Na-18c6-GIC also increased from 1.9 G2 to 29.9 G2. The consecutive signal-broadening of each component by cooling suggests no first-order structural phase transition in each GIC sample; this is confirmed by the absence of any DSC transition signal.

(a)

(b)

(c)

Fig. 1 1H static NMR spectra in Na-15c5-GIC (a), K-15c5-GIC (b), and Na-18c6-GIC (c). Some sharp components included in the spectra are noise caused in NMR hardware.

8 ACS Paragon Plus Environment

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Fig. 2 Deconvolution of 1H spectra taken at room temperature to two components; narrower component (C1) and broader component (C2).

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 28

(a)

(b)

Fig. 3 Signal intensity ratios of C1 / C2 for three GIC samples (a), and the estimated second moments (M2) of C1 and C2 for each sample.

10 ACS Paragon Plus Environment

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The 1H M2 values calculated for intramolecular 1H - 1H dipole - dipole interactions in the rigid and rotating state of 15c5 and 18c6 molecules are shown in Table 1. The evaluated M2 values of C1 and C2 at 293 K in Na-15c5-GIC (Fig. 2) are 0.04 G2 and 2.0 G2, respectively. Generally, the contribution to M2 by intermolecular 1H - 1H interactions can be estimated as 0.5∼1.5 G2. Therefore, the former M2 value (0.04 G2) should be assigned to the 15c5 molecule under isotropic rotation with liquid-like diffusion of intercalates because molecular diffusion can cancel the intermolecular contribution of M2 by averaging intermolecular interactions. This result means that most of Na-15c5 complexes are moving freely in the GIC galleries. The C1of K-15c5-GIC and Na-18c6-GIC can also be explained by the diffusion of crown ether molecules because the M2 values of C1 for K-15c5-GIC and Na-18c6-GIC at room temperature (0.01 G2 for K-15c5-GIC, and 0.04 G2 for Na-18c6-GIC) are similar to those obtained for Na-15c5-GIC.

Table 1 Calculated 1H second moments (M2) of 15c5 and 18c6 molecules assuming stable molecular structures. The values only consider intramolecular interaction without intermolecular

Motional mode

M2 / G2 (intramolecular) 15c5 18c6

static

15.6

16.2

rotation around an axis perpendicular to molecular plane

4.9

2.0

isotropic rotation

0.0

0.0

contribution.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 28

In contrast, the M2 values of C2 for Na-15c5-GIC (2.0 G2) and K-15c5-GIC (0.50 G2) observed at 293 K are comparable with the M2 values for isotropic rotation taking account of intermolecular interactions, but are too small for the axial rotation of 15c5. The assignment of molecular motion for 15c5 will be discussed after the calculation of the stable molecular structure of the Na+-15c5 complex in the GIC galleries, because an optimized stacked structure of these complexes is necessary to evaluate the molecular motion using the M2 values. The M2 of C2 for Na-18c6-GIC (1.9 G2) corresponds with axial rotation of 18c6 (2.0 G2) without intermolecular interactions (Table 1). The M2 of C2 for Na-15c5-GIC at 149 K (20.8 G2) might become larger at lower temperature because the M2 value is not constant even in the temperature range. The observed M2 values over 20 G2 (20.8, 29.2, and 29.9 G2 for Na-15c5-GIC, K-15c5-GIC and Na-18c6-GIC, respectively) at low temperature are significantly larger than the theoretical values for the rigid state, which means that the contribution of intermolecular interactions is not negligible, especially in K-15c5-GIC and Na-18c6-GIC. The crown ether molecules are densely packed in GIC galleries, with only a short separation between intercalates.

3.2. Structure of alkali metal and 15c5 between graphene layers estimated using DFT calculation and 23Na MAS NMR Figs. 4 (a) and (b) show the structures and total energy of two possible bilayer models for Na-15c5-GIC; (a) is a bilayer model with two stacked 15c5 molecules, whereas (b) shows an offset stacked structure. The total formation energy (Ef) of each system was estimated by the following equation: 12 ACS Paragon Plus Environment

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

f  tot 2Na155C  tot C

 2tot (Na155 /2 The Ef of system (b), -1.51 eV, was found to be lower than that for system (a), -1.34 eV, indicating that offset stacking is preferable for Na-15c5. The optimized gallery height (∆d) of 0.75 nm in these systems is slightly less than that observed experimentally (0.86 nm)43; however, offset stacking over the range of many molecules might result in a larger gallery height requirement. A similar result was also obtained for K-15c5-GIC (Fig. 5(a) and (b)). The formation energy of offset stacking (b) was 0.12 eV lower than for the stacked model (a). The

∆d in these systems, 0.80 nm, is also smaller than the experimental value (0.886 nm).

(a)

∆d = 0.83 nm Ef = -1.34 eV (b)

∆d = 0.75 nm Ef = -1.51 eV

Fig. 4 Structures and their total energy of two bilayer models of Na1.0(15c5)C25. (a) is a model of bilayer that two 15c5 molecules are stacked, whereas (b) considers an offset of stacking of ethers. (∆d : gallery height, Ef : total formation energy)

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

(a)

∆d = 0.82 nm Ef = -1.12 eV (b)

∆d = 0.80 nm Ef = -1.24 eV

Fig. 5 Structures and their total energy of two bilayer models of K1.0(15c5)C25. (a) is a model of bilayer that two 15c5 molecules are stacked, whereas (b) considers an offset of stacking of ethers.

The estimate of Ef also considers the state of the excess Na and K cations that do not coordinate to crown ethers. Fig. 6 shows two possible arrangements for Na1.5(15c5)C25; (a) with excess Na ions isolated on the graphene layers, or (b) with Na ions coordinated by forming a Na3(15c5)2 complex. Since the Ef of (a) is lower than that of (b), the formation of Na3(15c5)2 complexes appears unlikely. This result suggests that 1.0 of Na in the composition Na1.7(15c5)C25 is coordinated by 15c5, and the excess amount (0.7) of Na cations are on the surface of graphene layers without coordinating to 15c5.

23

Na MAS NMR spectra of

Na-15c5-GIC also show two components (Fig. S1(a)). Main peaks at 8 ppm and -23 ppm are ascribed to a second-order quadrupolar component (quadrupolar coupling constant; qcc= 4.2 14 ACS Paragon Plus Environment

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

MHz) of Na ion, whereas subordinate signal at +2 ∼ -5 ppm can be assigned to either Na ions at the graphene layer surface or to Na ions located in defects in the GIC. 23Na-18c6-GIC shows similar spectra (Fig. S1(b)), with a narrower qcc of the main component (qcc= 1.7 MHz) and the subordinate signal is indistinct. The anisotropy of the Na ion in Na-15c5 seems to be larger than in Na-18c6, although estimated Bader charges of Na ion in Na-15c5 and in Na-18c6 (Fig. S2) are almost the same.

(a)

∆d = 0.76 nm ∆E = 0 (b)

∆d = 0.80 nm ∆E = +0.35 eV per Fig. 6 Two possible structures and a comparison of their total energies of two bilayer models of Na1.5(15c5)C25. Na ion and 15C5 form 1:1 complex and the excess Na ion is isolate (a). Na and 15C5 form Na215C51 complex (b). The total formation energy (Ef) in (a) is 0.35 eV lower than (b).

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 28

4. DISCUSSION 4.1. Estimation of the motion of crown ethers in GIC galleries at room temperature According to the DFT calculations described in the previous section, an offset stacking model of two Na-15C5 complexes is the most stable structure for Na+-crown ether intercalates in these GICs. Based on the structure in Fig. 4(b), M2 values for the two offset stacked Na-15c5 complexes (Na2(15c5)2) were calculated and are shown in Table 2. The static M2 of 16.7 G2, which includes intramolecular interactions and intermolecular contribution of protons between the two 15c5 molecules, is larger than the 15.6 G2 obtained for static M2 of a 15c5 molecule (Table 1). The difference in these values, 1.1 G2, is the intermolecular contribution. The difference between experimental value at low temperature (20.8 G2) and 16.7 G2 can be ascribed to other intermolecular contributions between the neighboring molecules, or to 1H nuclei on the defect sites of the graphene surface.

Table 2 Calculated 1H second moments (M2) of Na2(15c5)2 complex (offset stacking) in Fig. 4(b). Motional mode of each molecule

M2 / G2 Na2(15c5)2

static

16.7

rotation around an axis perpendicular to molecular plane

5.3

rotation around an axis with 26° fluctuation of the axis

2.0

isotropic rotation

-(*)

(*)

Improbable because two 15c5 are too close.

16 ACS Paragon Plus Environment

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Since the estimated M2 for Na2(15c5)2 complexes under axial rotation of two 15c5 molecules, 5.3 G2, is still larger than the experimental value of 2.0 G2 for C2 in Na-15C5-GIC, we assumed a rotation of 15c5 molecules around an axis perpendicular to the molecular plane with some fluctuation of the rotation axis (z direction) in order to explain the M2 value. The experimental value (2.0 G2) of M2 coincides with the simulated M2 value when the fluctuation angle of the rotation axis to x and y directions is 26°. Considering the intermolecular contribution between the neighboring molecules, the fluctuation angles should be slightly larger than 26°. This molecular motion is illustrated in Fig. 7(a). For K-15c5-GIC and Na-18c6-GIC, the M2 values of C2 (0.51 G2 and 1.9 G2) were analyzed in the same way. We obtained a similar value of K-15c5-GIC by rotation of 15c5 around an axis perpendicular to molecular plane and another axis along to the molecular plane, which is almost equivalent to isotropic rotation of the molecule. In contrast, 1.9 G2 for 18c6 is almost explainable by rotation around the perpendicular axis with slight fluctuation, even though 0.3∼0.5 G2 in 1.9 G2 should be explained as an intermolecular contribution. These motions are illustrated in Figs. 7(b) and 7(c). The activation energies of the axial rotation for C2 of Na-15c5-GIC, K-15c5-GIC, Na-18c6-GIC are roughly estimated as 11.1, 7.0, and 12.2 kJ mol-1, respectively, by calculating ∆M2 between isotropic rotation and static state.46, 53

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(a)

Page 18 of 28

(b)

(c)

Fig. 7 Molecular motions of 15c5 and 18c6 in the Na-15c5-GIC (a), K-15c5-GIC (b), and Na-18c6-GIC (c).

4.2. Structure and motion of C2 crown ethers within GIC galleries The expected structures and their molecular motion of Na-15c5-GIC, K-15c5-GIC and Na-18c6-GIC estimated from the results of 1H NMR and DFT calculation are represented in Fig. 8. In Na-15c5-GIC (Fig. 8(a)), Na-15c5 complexes stack in pairs parallel to the encasing graphene layers with some offset. The C2 component of complexes are under axial rotation with a fluctuation of the rotation axis over 26°. In contrast, the motion of K-15c5 in K-15c5-GIC is more active than Na-15c5. The fluctuation is almost an isotropic rotation (Fig. 8(b)). In this case, a stacked structure of two complexes is unlikely because of their molecular motion. The nearly isotropic motion of K-15c5 is also supported by lower activation energy of the axial rotation for C2 of K-15c5-GIC (7.0 kJ mol-1). Indeed, our 2H NMR analyses suggest that K ions in K-piperazine-GIC and K-1,4-diazabicyclo[2.2.2]octane(dabco)-GIC are not strongly coordinated by organic molecules.44 In the present case, the 15c5 molecules under isotropic like rotation might be weakly anchored in the GIC galleries. These results imply that

18 ACS Paragon Plus Environment

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the structure and the molecular motion of intercalants are largely influenced by interactions from alkali metal cations. (a)

(b)

(c)

Fig. 8

The expected structure and the molecular motion for C2 of Na-15c5-GIC (a),

K-15c5-GIC (b) and Na-18c6-GIC (phase α) estimated from the results of 1H NMR and DFT calculation.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 28

For Na-18c6-GIC, although two GIC phases with different gallery heights (∆d) have been reported,43 the major component of our sample is phase α (∆d = 086 nm). In the GIC, Na-18c6 complexes can rotate around an axis perpendicular to the molecular plane with slight fluctuation at room temperature (Fig 8(c)). 18c6 is planar and symmetric, and has smaller distortion than 15c5 in the molecular structure. The structure might affect the fluctuation of the axial rotation of the molecules coordinating to alkali metal ions.

5. CONCLUSION The structures and molecular dynamics of 15c5 and 18c6 molecules coordinating to Na or K in ternary GICs were revealed by the combination of total energy estimation using DFT calculation and analyses of 1H solid state NMR spectra. DFT calculation suggested a more stable structure comprising Na-15c5 (1:1) stacked complexes stack with some offset. The excess Na or K can be ascribed to ions existing on the defect structure in graphene surface. By 1

H NMR line widths, two kinds of molecules undergoing different molecular motions were

observed at room temperature; isotropic rotation with molecular diffusion, and axial rotation with fluctuation of the axis. The rapid diffusion of crown ether - cation complex intercalates will be advantageous for electrochemical and chemical intercalation and deintercalation. The structure and dynamics of crown ether molecules in GIC galleries are strongly affected by the geometry of the crown ether molecules and the strength of the interaction between alkali metal and ligand molecules.

20 ACS Paragon Plus Environment

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: 23

Na MAS NMR spectra of Na-15c5-GIC and Na-18c6-GIC taken at room temperature,

Bader charges on Na ions in Na-15c5-GIC, Na-18c6-GIC, and Na on graphene estimated by DFT calculation.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] Notes The authors declare no competing financial interest. ACKNOWLEDGMENTS This work was supported by JSPS KAKENHI (Grant No. 17K06017), and Elements Strategy Initiative for Catalysts and Batteries (ESICB) in Kyoto University.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 28

References 1.

Enoki, T.; Suzuki, M.; Endo, M. Graphite Intercalation Compounds and Applications; Oxford: New York, 2003.

2.

Matsumoto, R. Recent Advances in Superconductive Graphite Intercalation Compounds. Int. J. Condens. Matter Adv. Mater. Supercond. Res. 2013, 12, 75-83.

3.

Hannay, N. B.; Geballe, T. H.; Matthias, B. T.; Andres, K.; Schmidt, P.; MacNair, D. Superconductivity in Graphite Compounds. Phys. Rev. Lett. 1965, 14, 225-226.

4.

Matsumoto, R.; Okabe, Y.; Akuzawa, N. Thermoelectric Properties and Performance of N-Type and P-Type Graphite Intercalation Compounds. J. Electron. Mater. 2015, 44, 399-406.

5.

Yoshio, M.; Brodd, R. J.; Kozawa, A. Lithium-Ion Batteries; Science and Technologies; Springer: LaVergne, 2009.

6.

Larcher, D.; Tarascon, J. M. Towards Greener and More Sustainable Batteries for Electrical Energy Storage. Nat. Chem. 2015, 7, 19-29.

7.

Komaba, S.; Hasegawa, T.; Dahbi, M.; Kubota, K. Potassium Intercalation into Graphite to Realize High-Voltage/High-Power Potassium-Ion Batteries and Potassium-Ion Capacitors. Electrochem. Commun. 2015, 60, 172-175.

8.

Bie, X.; Kubota, K.; Hosaka, T.; Chihara, K.; Komaba, S. A Novel K-Ion Battery: Hexacyanoferrate(II)/Graphite Cell. J. Mater. Chem. A 2017, 5, 4325-4330.

9.

Zou, X.; Xiong, P.; Zhao, J.; Hu, J.; Liu, Z.; Xu, Y. Recent Research Progress in Non-Aqueous Potassium-Ion Batteries. Phys. Chem. Chem. Phys. 2017, 19, 26495-26506.

10.

Eftekhari, A.; Jian, Z.; Ji, X. Potassium Secondary Batteries. ACS Appl. Mater. Inter. 2017, 9, 4404-4419.

11.

Jian, Z.; Luo, W.; Ji, X. Carbon Electrodes for K-Ion Batteries. J. Am. Chem. Soc. 2015, 137, 11566-11569. 22 ACS Paragon Plus Environment

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

12.

Luo, W.; Wan, J.; Ozdemir, B.; Bao, W.; Chen, Y.; Dai, J.; Lin, H. Xu Y.; Gu, F.; Barone, V.; Hu L. Potassium Ion Batteries with Graphitic Materials. Nano Lett. 2015, 15, 7671-7677.

13.

Rodriguez-Pérez, I. A.; Ji, X. Anion Hosting Cathodes in Dual-Ion Batteries. ACS Energy Lett. 2017, 2, 1762-1770.

14.

Beltrop, K.; Beuker, S.; Heckmann, A.; Winter, M.; Placke, T. Alternative Electrochemical Energy Storage: Potassium-Based Dual-Graphite Batteries. Energy & Environ. Sci. 2017, 10, 2090-2094.

15.

Komaba, S.; Murata, W.; Ishikawa, T.; Yabuuchi, N.; Ozeki, T.; Nakayama, T.; Ogata, A.; Gotoh, K.; Fujiwara, K. Electrochemical Na Insertion and Solid Electrolyte Interphase for Hard-Carbon Electrodes and Application to Na-Ion Batteries. Adv. Funct. Mater. 2011, 21, 3859-3867.

16.

Yabuuchi, N.; Kubota, K.; Dahbi, M.; Komaba, S. Research Development on Sodium-Ion Batteries. Chem. Rev. 2014, 114, 11636-11682.

17.

Nayak, P. K.; Yang, L.; Brehm, W.; Adelhelm, P. From Lithium-Ion to Sodium-Ion Batteries: Advantages, Challenges, and Surprises. Angew. Chem. Int. Ed. 2018, 57, 102-120.

18.

Jache, B.; Adelhelm, P. Use of Graphite as a Highly Reversible Electrode with Superior Cycle Life for Sodium-Ion Batteries by Making Use of Co-Intercalation Phenomena. Angew. Chem. Int. Ed. 2014, 53, 10169-10173.

19.

Kim, H.; Hong, J.; Park, Y.; Kim, J.; Hwang, I.; Kang, K., Sodium Storage Behavior in Natural Graphite Using Ether-Based Electrolyte Systems. Adv. Funct. Mater. 2015, 25, 534-541.

20.

Kim, H.; Hong, J.; Yoon, G.; Kim, H.; Park, K.; Park, M.; Yoon, W.; Kang, K. Sodium Intercalation Chemistry in Graphite. Energy Environ. Sci. 2015, 8, 2963-2969.

21.

Zhu, Z.; Cheng, F.; Hu, Z.; Niu, Z.; Chen, J. Highly Stable and Ultrafast Electrode Reaction of Graphite for Sodium Ion Batteries. J. Power Sources 2015, 293, 626-634. 23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

22.

Page 24 of 28

Jache, B.; Binder, J. O.; Abe, T.; Adelhelm, P. A Comparative Study on the Impact of Different Glymes and Their Derivatives as Electrolyte Solvents for Graphite Co-Intercalation Electrodes in Lithium-Ion and Sodium-Ion Batteries. Phys. Chem. Chem. Phys. 2016, 18, 14299-14316.

23.

Seidl, L.; Bucher, N.; Chu, E.; Hartung, S.; Martens, S.; Schneider, O.; Stimming, U. Intercalation of Solvated Na-Ions into Graphite. Energy Environ. Sci. 2017, 10, 1631-1642.

24.

Kim, J.; Jung, M.; Kim, M.; Lee, Y. Electrochemical Performances of Lithium and Sodium Ion Batteries Based on Carbon Materials. J. Ind. Eng. Chem. 2018, 61, 368-380.

25.

Goktas, M.; Bolli, C.; Berg, E. J.; Novák, P.; Pollok, K.; Langenhorst, F.; Roeder, M. v.; Lenchuk, O.; Mollenhauer, D.; Adelhelm, P. Graphite as Cointercalation Electrode for Sodium-Ion Batteries: Electrode Dynamics and the Missing Solid Electrolyte Interphase (SEI). Adv. Energy Mater. 2018, 1702724.

26.

Gotoh, K.; Maruyama, H.; Miyatou, T.; Mizuno, M.; Urita, K.; Ishida, H. Structure and Dynamic Behavior of Sodium-Diglyme Complex in the Graphite Anode of Sodium Ion Battery by 2H Nuclear Magnetic Resonance. J. Phys. Chem. C 2016, 120, 28152-28156.

27.

Jung, S. C.; Kang, Y.; Han, Y. Origin of Excellent Rate and Cycle Performance of Na+ -Solvent Cointercalated Graphite vs. Poor Performance of Li+ -Solvent Case. Nano Energy 2017, 34, 456-462.

28.

Co Minh, D.; Mai, C.; Riviere, R.; Gole, J. Insertion of Solvated Alkali Metal Ions into Graphite. Structural Study of the Products Obtained in DME [Dimethoxyethane]. J. Chim. Phys. Physicochim. Biol. 1972, 69, 991-995.

29.

Beguin, F.; Estrade-Szwarckopf, H.; Conard, J.; Lauginie, P.; Marceau, P.; Guerard, D. Ternary Graphite-Lithium-Tetrahydrofuran Compounds: Synthesis and X-Ray and Magnetic Resonance Studies. Synth. Met. 1983, 7, 77-84.

30.

Schlögl, R.; Boehm, H. P. The Reaction of Potassium-Graphite Intercalation Compounds with Tetrahydrofuran. Carbon 1984, 22, 341-349. 24 ACS Paragon Plus Environment

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

31.

Quinton, M. F.; Legrand, A. P.; Beguin, F. Structure of the Intercalated Layer in Graphite-Potassium-Tetrahydrofuran Compounds. Synth. Met. 1986, 14, 179-188.

32.

Beguin, F.; Guerard, D.; Goldman, M. Temperature Induced Collapse of the Lithium-THF Complexes in the Li(THF)2.13C12 Graphitide. Synth. Met. 1990, 34, 33-39.

33.

Inagaki, M.; Tanaike, O. Host Effect on the Formation of Sodium-Tetrahydrofuran-Graphite Intercalation Compounds. Synth. Met. 1995, 73, 77-81.

34.

Inagaki, M.; Tanaike, O. Formation of Na-THF-Graphite Intercalation Compounds. J. Phys. Chem. Solids 1996, 57, 795-798.

35.

Abe, T.; Mizutani, Y.; Inaba, M.; Ogumi, Z. Li-Graphite Intercalation Compounds Synthesized in Various Ether-Type Organic Solvents. Tanso 1998, 185, 290-295.

36.

Maluangnont, T.; Bui, G. T.; Huntington, B. A.; Lerner, M. M. Preparation of a Homologous Series of Graphite Alkylamine Intercalation Compounds Including an Unusual Parallel Bilayer Intercalate Arrangement. Chem. Mater. 2011, 23, 1091-1095.

37.

Maluangnont, T.; Gotoh, K.; Fujiwara, K.; Lerner, M. M. Cation-Directed Orientation of Amines in Ternary Graphite Intercalation Compounds. Carbon 2011, 49, 1040-1042.

38.

Maluangnont, T.; Lerner, M. M.; Gotoh, K. Synthesis of Ternary and Quaternary Graphite Intercalation Compounds Containing Alkali Metal Cations and Diamines. Inorg. Chem. 2011, 50, 11676-11682.

39.

Maluangnont, T.; Sirisaksoontorn, W.; Lerner, M. M. A Comparative Structural Study of Ternary Graphite Intercalation Compounds Containing Alkali Metals and Linear or Branched Amines. Carbon 2012, 50, 597-602.

40.

Sirisaksoontorn, W.; Lerner, M. M. Preparation of a Homologous Series of Tetraalkylammonium Graphite Intercalation Compounds. Inorg. Chem. 2013, 52, 7139-7144.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

41.

Page 26 of 28

Zhang, H.; Lerner, M. M. Preparation of Graphite Intercalation Compounds Containing Oligo and Polyethers. Nanoscale 2016, 8, 4608-4612.

42.

Zhang, H.; Wu, Y.; Sirisaksoontorn, W.; Remcho, V. T.; Lerner, M. M. Preparation, Characterization, and Structure Trends for Graphite Intercalation Compounds Containing Pyrrolidinium Cations. Chem. Mater. 2016, 28, 969-974.

43.

Zhang, H.; Lerner, M. M. Preparation of Graphite Intercalation Compounds Containing Crown Ethers. Inorg. Chem. 2016, 55, 8281-8284.

44.

Gotoh, K.; Sugimoto, C.; Morita, R.; Miyatou, T.; Mizuno, M.; Sirisaksoontorn, W.; Lerner, M. M.; Ishida, H. Arrangement and Dynamics of Diamine, Etheric, and Tetraalkylammonium Intercalates within Graphene or Graphite Oxide Galleries By 2H NMR. J. Phys. Chem. C 2015, 119, 11763-11770.

45.

Van Vleck, J. H. The Dipolar Broadening of Magnetic Resonance Lines in Crystals. Phys. Rev. 1948, 74, 1168-1183.

46.

Ueda, T.; Omichi, H.; Chen, Y.; Kobayashi, H.; Kubota, O.; Miyakubo, K.; Eguchi, T. Two-Dimensional Melting and Phase Change of Binary Mixtures of CCl4 and CHCl3 confined in ACF Nanospace Studied Using Solid-State 1H NMR. Bull. Chem. Soc. Jpn. 2010, 83, 1323-1332.

47.

Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169-11186.

48.

Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758-1775.

49.

Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868.

50.

Zhang, Y.; Yang, W. Comment on "Generalized Gradient Approximation Made Simple". Phys. Rev. Lett. 1998, 80, 890-890.

26 ACS Paragon Plus Environment

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

51.

Klimeš, J.; Bowler, D. R.; Michaelides, A. Van Der Waals Density Functionals Applied to Solids. Phys. Rev. B 2011, 83, 195131.

52.

Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979.

53.

Abragam, A., The Principles of Nuclear Magnetism; Oxford, New York, 1961.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 28

TOC GRAPHICS

28 ACS Paragon Plus Environment