Chemical Science Research Frontiers 2022 92 Phonon-assisted proton(deuteron) transfers leading the antiferro-electric ordering in superprotonic conductors Cs 3H(D)(SeO 4) 2 Proton conduction is a key functionality to improve a high-performance fuel cell that is a fundamental technology to reduce CO 2 emissions and to achieve a hydrogen society. In general, proton conduction is phenomenologically discussed with activation- type motions as a diffusion process. In crystalline proton conductors, a proton–phonon (PP) interaction is significant to understand microscopic-transfer mechanisms because the proton never transfers individually but interacts with a lattice system. Nevertheless, little knowledge on PP interactions has been obtained so far in protonic systems. Among a variety of proton conductors, we have focused on Cs 3 H(D)(SeO 4 ) 2 that exhibits superprotonic conduction above T spc and antiferroelectric (AF) ordering below T AF ( Fig. 1 (a)), because the good single crystal is obtained. The crystal structure on the ab plane of localized phases III and IV is shown in Fig. 1(b), where the lattice constant is summarized in the table. Two SeO 4 tetrahedra are connected through a hydrogen bond between apical O(2) and O’(2) in the form of dimer ( Fig. 1(c)). The proton and deuteron possess two equivalent sites in between O(2) and O’(2). A remarkable isotope effect relevant to the AF ordering was known about 40 years ago (Fig. 1(a)); T AF = 50 K (Cs 3 H(SeO 4 ) 2 ) and 168 K (Cs 3 D(SeO 4 ) 2 ). Similar isotope effects are observed in the corresponding systems such as Rb 3 H(D)(SeO 4 ) 2 , and ferroelectric materials like potassium dihydrogen phosphate (KDP). Employing a wide range of infrared spectroscopies combined with first principles calculations and 1 H-NMR experiments, we have clarified the difference of AF- ordering mechanisms and the role of PP interactions in Cs 3 H(SeO 4 ) 2 and Cs 3 D(SeO 4 ) 2 [1,2]. Absorbance spectra below 300 K were measured at 5 – 8000 cm −1 utilizing a terahertz time-domain spectrometer, and a FT-IR spectrometer equipped with a Cassegrain microscope. Especially in far- infrared range, the single-domain crystal with about 1 mm in diameter is quite difficult to measure with usual spectrometers. Thanks to the microscopic spectrometer facilitated at SPring-8 BL43IR providing strong far-infrared light, we have successfully detected the polarized absorbance spectra. In the far-infrared spectra of Cs 3 H(SeO 4 ) 2 ( Fig. 2(a)) and Cs 3 D(SeO 4 ) 2 ( Fig. 2(b)), a libration and SeO 4 -tetrahedral deformations accompanied with a O-Se-O bending (def bend mode) appear at 180 and 300 – 450 cm –1 , respectively [1,2]. The libration in Cs 3 H(SeO 4 ) 2 exhibits a spectral splitting at 30 – 70 K, while such anomaly is never observed in Cs 3 D(SeO 4 ) 2 . In Fig. 2(c), the splitting for a polarization at T AF in Cs 3 H(SeO 4 ) 2 is reproduced with two Lorentzian curves. The vibrational motif of libration illustrated in Fig. 2(d) induces a characteristic modulation of O(2)-O’(2) distance. From the mid-infrared spectrum, anomalous OH stretching mode proves that the proton anharmonically couples to the libration, and thereby the PP interaction plays the important role to generate the AF ordering [1]. In Cs 3 D(SeO 4 ) 2 , anomalous enhancements of the integrated absorbance ( Fig. 2(e) ) are observed at around T AF in the 310-cm –1 def bend mode that is painted in blue in Fig. 2(b), and the vibrational motif is depicted in Fig. 2(f). Taking the anharmonicity of OD stretching mode into account, the deuteron is found to couple to the 310-cm –1 def bend mode [2]. Therefore, the proton and deuteron interact to the different phonon modes. The proton and deuteron have to transfer between the equivalent sites to achieve the AF arrangement. According to the 1 H-NMR measurements, conventional thermal hopping in a rigid double-minimum potential ( Fig. 3(a)) diminishes below 150 K [1]. As the modulation of O(2)-O’(2) distance is huge enhanced through fluctuations towards the AF ordering in 456 K 448 K 372 K 168 K 369 K 50 K H D (a) (b) a b (c) O(2)-O’(2) = 0.254 nm O(2) O’(2) O(1) H + (D + ) Cs O(1) Se H + O(2) SPC AF T SPC T AF T II-III (I) (II) (III) (IV) a b 1.0888 nm Cs 3 H(SeO 4 ) 2 Cs 3 D(SeO 4 ) 2 0.6388 nm 1.0891 nm 0.6383 nm Fig. 1. (a) Transition temperatures in Cs 3 H(SeO 4 ) 2 and Cs 3 D(SeO 4 ) 2 . (b) Crystal structure on the ab plane and lattice constants. (c) Two equivalent sites in O(2)-O’(2) of the dimer. Research Frontiers 2022 93 Cs 3 H(SeO 4 ) 2 , the wave functions of two equivalent sites can overlap, and simultaneously a dynamical proton tunneling yields the splitting of the ground state as shown by red lines in Fig. 3(b). Since each level strongly couples to the libration, the split bands emerge in the far-infrared time scale. This phenomenon corresponds to a phonon-assisted proton tunneling (PAPT) associated with the libration, and hence quantum fluctuations conduct the AF ordering. A splitting is never observed in the spectra of Cs 3 D(SeO 4 ) 2 , though the anomalous integrated absorbance of 310-cm –1 def bend mode represents the increase of oscillator strength through huge modulations of the O(2)-O’(2) distance [2]. The modulation dynamically suppresses the potential- barrier height ( Fig. 3 (c)), and then a phonon-assisted deuteron hopping (PADH) associated with the 310-cm –1 def bend mode generates classical fluctuations to lead the AF ordering. The isotope effect is concluded to change the phonon that helps the transfer. Since the 310-cm –1 def bend mode is energetically larger than the libration, Cs 3 D(SeO 4 ) 2 possesses higher T AF than Cs 3 H(SeO 4 ) 2 . Finally, we would like to mention the 610-cm –1 band that exhibits the most remarkable isotope effect. The absorbance spectrum obtained by the first principles calculation has no bands at 500~700 cm –1 . The 610-cm –1 band is not attributed to phonons and molecular vibrations, but to a collective excitation relevant to a classical hopping, because the barrier height coincides with the result evaluated from the 1 H-NMR measurement [1]. The intensity in Cs 3 D(SeO 4 ) 2 ( Fig. 2(b)) is about ten times as large as one in Cs 3 H(SeO 4 ) 2 ( Fig.2(a)) [2]. In Cs 3 H(SeO 4 ) 2 , the PAPT associated with the 440-cm –1 def bend mode becomes dominant at 70– 250 K, so that very few protons contribute to the classical hopping. On the other hand, all the deuterons contribute to the classical hopping, and consequently the 610-cm –1 band in Cs 3 D(SeO 4 ) 2 has large intensity compared to Cs 3 H(SeO 4 ) 2 . Hiroshi Matsui Department of Physics, Tohoku University Email: hiroshi.matsui.b2@tohoku.ac.jp References [1] H. Matsui et al .: J. Chem. Phys. 152 (2020) 154502. [2] H. Matsui, K. Fukuda, S. Takano, Y. Ikemoto, T. Sasaki and Y. Matsuo: J. Chem. Phys. 156 (2022) 204504. Fig. 2. (a, b) Absorbance (A) of Cs 3 H(SeO 4 ) 2 and Cs 3 D(SeO 4 ) 2 in far-infrared range. (c, d) Splitting and vibrational motif of the libration in Cs 3H(SeO 4) 2. (e, f) Integrated absorbance (IA) and vibrational motif of the 310-cm –1 def str mode in Cs 3D(SeO 4) 2. Fig. 3. Schematic illustrates of double minimum potential for thermal hopping (a) , PAPT (b) and PADH (c) . 200 300 400 500 600 100 700 (a) (b) A (arb. unit) Wavenumber (cm –1 ) 300 K 200 K 100 K 6 K 300 K Cs 3 D(SeO 4 ) 2 Cs 3 H(SeO 4 ) 2 200 K 100 K 12 K 170 180 190 (c) 0 50 100 150 200 250 300 T (K) T AF (d) (e) (f) A (arb. unit) IA (arb. unit) Wavenumber (cm –1 ) PAPT Thermal hopping PADH (a) (b) (c)