Dynamics in the Plastic Crystalline Phases of Cyclohexanol and

55 mins ago - User Resources. About Us · ACS Members · Librarians · Authors & Reviewers · Website Demos · Privacy Policy · Mobile Site ...
0 downloads 3 Views 708KB Size
Subscriber access provided by University of Winnipeg Library

B: Fluid Interfaces, Colloids, Polymers, Soft Matter, Surfactants, and Glassy Materials

Dynamics in the Plastic Crystalline Phases of Cyclohexanol and Cyclooctanol Studied by Quasielastic Neutron Scattering Eric Christopher Novak, Niina Jalarvo, Sudipta Gupta, Kunlun Hong, Stephan Förster, Takeshi Egami, and Michael Ohl J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.8b03448 • Publication Date (Web): 18 May 2018 Downloaded from http://pubs.acs.org on May 18, 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 27 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

Dynamics in the Plastic Crystalline Phases of Cyclohexanol and Cyclooctanol Studied by Quasielastic Neutron Scattering E. Novak,†,‡,§, N. Jalarvo*,‡,§, S. Gupta||, K. Hong¶,#, S. Förster§, T. Egami†,⊥,◊, and M. Ohl§ †

Department of Materials Science and Engineering, University of Tennessee, Knoxville, TN 37996, USA



Neutron Sciences Directorate, Oak Ridge National Laboratory (ORNL), Oak Ridge, TN 37831, USA

§

Jülich Centre for Neutron Science (JCNS), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany

Department of Chemistry and Macromolecular Studies Group, Louisiana State University, Baton

||

Rouge, LA 70803, USA ¶

Center for Nanophase Materials Sciences, Oak Ridge National Laboratory (ORNL), Oak Ridge, TN 37831, USA

#

Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN 37996, USA

ACS Paragon Plus Environment

1

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

⊥Materials

Page 2 of 27

Science and Technology Division, Oak Ridge National Laboratory (ORNL), Oak Ridge, TN 37831, USA

◊Department

of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996, USA

Notice of Copyright

This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC0500OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally

sponsored

research

in

accordance

with

the

DOE

Public

Access

Plan

(http://energy.gov/downloads/doe-public-access-plan).

Abstract

Plastic crystals are a promising candidate for solid state ionic conductors. In this work, quasielastic neutron scattering is employed to investigate the center of mass diffusive motions in two types of plastic crystalline cyclic alcohols: cyclohexanol and cyclooctanol. Two separate motions are observed which are attributed to long range translational diffusion (α process) and cage rattling (fast β-process). Residence times and diffusion coefficients are calculated for both processes, along with the confinement distances for the cage rattling. In addition, a binary mixture of these two materials is measured to understand how the dynamics

ACS Paragon Plus Environment

2

Page 3 of 27 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

change when a second type of molecule is added to the matrix. It is observed that upon the addition of the larger cyclooctanol molecules into the cyclohexanol solution, the cage size decreases, which causes a decrease in the observed diffusion rates for both the α - and fast βprocesses.

ACS Paragon Plus Environment

3

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 27

Introduction In today’s highly mobile, information-based society, the demand for portable electronic devices powered by long lasting batteries is ever increasing. Currently, the energy market is dominated by rechargeable lithium-based liquid electrolyte batteries due to their high energy density, design flexibility, and high ionic conductivity1. However, these batteries have major safety concerns due to their dangerous liquid electrolyte, which can be flammable and are prone to leakage2. This concern has prompted significant research in the field of energy science recently, largely due to the emergence of electric vehicles. A safer alternative would be to use a solid-state ionic conductor, such as a glass or polymer, which would not leak if the battery is damaged. A promising candidate for solid state ionic conductors are plastic crystals (PCs), which are orientationally disordered crystals, also known as orientational glasses3. PCs are a type of molecular glass former that possesses properties of both liquids and solids. The center of mass of the molecules are positioned on lattice sites while the orientations of the molecules are dynamically disordered. Translational motion is largely confined but rotational degrees of freedom are available for motion. Ions diffuse through this system by the “revolving door” mechanism, where the molecules rotate to increase the free volume, creating space for the ions to diffuse efficiently through the opening of these doors4-7. Some PCs have been reported with technologically relevant ionic conductivities around ambient temperatures upon the addition of ions4-5, 8-9. PC systems tend to have complex phase diagrams with metastable phases, multiple phase transitions, and are highly dependent on cooling and heating rates. Investigations have shown that the PC phase can be stabilized and extended throughout a larger temperature range through the formation of a binary PC mixture6,

10

. This also allows for the tuning of the ionic conductivity by altering the

concentration of the constituents. However, the ionic conductivity of two different binary PC

ACS Paragon Plus Environment

4

Page 5 of 27 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

mixtures have been investigated recently with conflicting results4, 6. In one system, the addition of larger molecules to a solution of smaller molecules increased the ionic conductivity while in the second system the opposite occurred. In this paper, we investigate one of these systems, i.e. cyclooctanol (OCT) and cyclohexanol (HEX) and their mixture, using backscattering neutron spectroscopy to better understand the dynamics of their PC phases. The materials are measured in their pure forms and in a binary mixture to shed light on how the motion of the molecules change when the binary mixture is created. Ions are not added to the solutions to ensure that the measured signal only reflects the motion of the PC matrix. For the revolving door mechanism, the rotation of the molecules is the most important motion that determines the success of the ionic diffusion in PCs. However, translation of the center of mass of the molecules is also important and can influence the amount of free volume created. Within the context of Götze’s mode coupling theory (MCT), there exists two motions above the crossover temperature: a slower long range translational diffusion (α-process) and a faster cage rattling (fast β-process) motion11-12. A cage is formed by the nearest neighbors (and to a lesser extent second nearest neighbors) when their positions are fixed in space. When this occurs, a particle inside the cage can only undergo small amplitude, sub-diffusive motions that are less than the diameter of the molecule, known as cage rattling. This is representative of the local environment in PCs where the molecules are largely fixed in space on lattice sites. For a particle to move over longer distances, a “door” needs to be opened in the cage. This increase in free volume will allow the particle to break free from the cage and contribute to the structural relaxation (α-process). In the case of ionic diffusion, an ion that is stuck inside the cage will now be able to escape out through the door. However, the cage can only open when many particles in the system are able to cooperatively rearrange. This framework has been confirmed by molecular dynamics

ACS Paragon Plus Environment

5

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 27

simulations13. The free volume created by the α- and fast β-process thereby directly influences the success of the revolving door mechanism and the ionic diffusion. Quasielastic neutron scattering (QENS) can access the dynamics associated with these two processes and probe the changes in a binary mixture. This knowledge leads to a better understanding of how to tune the local environment so that the ionic conductivity can be optimized. Materials and Methods Three samples were investigated in this experiment: cyclooctanol (C8H16O), cyclohexanol (C6H12O), and 80% cyclohexanol – 20% deuterated cyclooctanol (C8D16O) in terms of mol% concentration. Protonated OCT was purchased from Tokyo Chemical Industry Co. (TCI) with a reported purity of >98.0%. Protonated HEX was purchased from Sigma-Aldrich with a purity of 99%. The fully deuterated OCT was synthesized at the Center for Nanophase Materials Sciences (CNMS) at Oak Ridge National Laboratory (ORNL). QENS measurements were performed using the backscattering spectrometer BASIS at the Spallation Neutron Source (SNS) at ORNL14. Annular aluminum sample cans were used with a thickness of 0.05 mm to minimize multiple scattering effects by ensuring the sample transmission was ≥ 95%. All samples were mixed in the appropriate amounts and loaded into the sample holders inside a glove box under an inert controlled atmosphere to reduce exposure to moisture. The standard instrumental setup was used with an energy resolution FWHM = 3.5 µeV yielding an accessible energy range of ± 120 µeV and Q range from 0.2 to 2.0 Å-1. Note that plastic crystals are known to have coherent structural peaks at this Q-range, usually above 1 Å-1. For this reason, a range of Q values had to be removed from the data analysis to avoid the effects of the coherent signal. The choppers were set at 60 Hz and the wavelength was centered at 6.4 Å. Samples were cooled slowly from ambient temperature to 20 K over the course of approximately 3 hours using a closed cycle refrigerator. The resolution

ACS Paragon Plus Environment

6

Page 7 of 27 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

function was measured at 20 K for each sample and measurements at higher temperatures were performed for approximately 3 hours/temperature. The data was normalized using an annular vanadium standard, reduced using Mantid, and model fit were performed with DAVE software package15. Data Analysis The QENS intensity is proportional to the dynamic incoherent structure factor Sinc(Q, E). This contains information about the correlation of the same scatterer at different times, allowing for accurate measurements of diffusive motions. Due to the large incoherent cross section of hydrogen and the high amounts of hydrogen contained in each molecule, the measured intensity will yield information about the motion of the hydrogen in the system. The data was binned to 18 averaged Q values and each of them was fit using the expression16-19: 𝐼(𝑄, 𝐸) = [𝑋(𝑄)𝛿(𝐸) + (1 − 𝑋(𝑄))𝑆𝑞𝑒 (𝑄, 𝐸)] ⊗ 𝑅(𝑄, 𝐸) + 𝐵(𝑄, 𝐸)

(1)

Here, δ(E) is a delta function representing the elastic signal and Sqe(Q,E) is the model representing the quasielastic scattering. R(Q,E) is the resolution function. X(Q) represents the fraction of elastic scattering in the signal, also known as the elastic incoherent structure factor (EISF). The symbol ⊗ represents the numerical convolution and B(Q,E) is the linear background term arising from the sample, sample environment and spectrometer. Two motions are clearly observable at the measured time scale, corresponding to a slower (narrower Lorentzian component) and a faster (broader Lorentzian component) dynamic process. We associate the narrower Lorentzian to the long range translational diffusion and the broader Lorentzian to cage rattling. Evidence supporting this reasoning will be presented in detail later in the discussion section. Therefore, a single Lorentzian model was insufficient and a two-Lorentzian model was needed for Sqe(Q,E) of the form:

ACS Paragon Plus Environment

7

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

𝑆𝑞𝑒 (𝑄, 𝐸) = (1 − 𝑃(𝑄))

Page 8 of 27

1 𝛤1 (𝑄) 1 𝛤2 (𝑄) + 𝑃(𝑄) 2 𝜋 𝐸 2 + 𝛤1 (𝑄) 𝜋 𝐸 2 + 𝛤22 (𝑄)

(2)

where (1–P(Q)) and P(Q) represent the spectral weights of the narrow and broad Lorentzian terms of width Γ1 and Γ2, respectively. The first term (1–P(Q)) is a measurable quantity called elastic incoherent structure factor (EISF), that is the fraction of total quasielastic intensity contained in the purely elastic peak20 𝐸𝐼𝑆𝐹 =

𝐼𝐸 𝐼𝐸 + 𝐼𝑄𝐸

(3)

where IE and IQE are the respective integrated elastic and quasielastic intensities. The total EISF contains the quasielastic intensity integrated from both Lorentzian functions. However, for the calculations in this paper, only the EISF with respect to the broad component intensity is investigated so that the length scales for the fast β-process can be calculated. Physically, the EISF is directly related to the region of space accessible to the scatterers and yields information about the geometry of the motion21. Diffusive motion inside the cage can be modeled using the diffusion within a sphere model16, 19-20, 22 2

3𝑗1 (𝑄𝑟𝐶𝑎𝑔𝑒 ) 𝐸𝐼𝑆𝐹 = 𝑝1 + (1 − 𝑝1 ) ( ) 𝑄𝑟𝐶𝑎𝑔𝑒

(4)

where rCage is the confinement radius of the sphere, j1 is the first-order Bessel function, and p1 represents the fraction of hydrogen in the system that is “immobile” and does not contribute to the diffusive motion. A relatively good fit is achieved using Eq. 4 for the higher temperature measurements, but this approach can be improved as follows. Since the molecules in PCs are positioned on lattice sites, they are expected to have small translational diffusion components compared to rotational motion. As the temperature is decreased, the molecules will also experience smaller amplitudes of cage rattling. This means that rotation of the molecule and segmental motion

ACS Paragon Plus Environment

8

Page 9 of 27 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

not associated with center of mass motion will have a larger effect in the EISF. This is worsened by the fact that there is a limiting momentum transfer value at higher Q-values that denotes the crossover from center of mass motion to intramolecular motion. Therefore, a second component has been added to Eq. 4 to account for this motion, known as the diffusion on a sphere model,16, 19-20, 22

to give the full expression: 2

3𝑗1 (𝑄𝑟𝐶𝑎𝑔𝑒 ) 𝐸𝐼𝑆𝐹 = [𝑝1 + (1 − 𝑝1 ) ( ) ] [𝑝2 + (1 − 𝑝2 )𝑗02 (𝑄𝑟𝐼𝑛𝑡𝑟𝑎 )] 𝑄𝑟𝐶𝑎𝑔𝑒

(5)

Here, p2 again represents the fraction of hydrogen that is immobile to the motion, j0 is the zeroorder Bessel function, and rIntra is the sphere radius. This additional term has been previously used to account for the skeletal motion of CH2 groups about the C-C bonds22. Let it be noted that the two models described here, i.e. diffusion “within” and “on” a sphere, are indeed distinct processes. The diffusion within a sphere model describes a diffusive motion that is free to move inside the confinement radius of a sphere. On the other hand, the diffusion on a sphere model describes a motion that is confined to move only along the surface of a sphere. This describes bonded atoms or side groups whose motion is confined to a distinct length scale associated with the bond length. The EISF is supposed to approach a value of 1 as Q approaches zero. However, the low Q values were too small, which is most likely the result of multiple scattering. These effects were minimized by using a thin sample with high transmission. Multiple scattering corrections are not necessary for the interpretation of our data and this rather complex procedure was therefore omitted. The Q dependency of the HWHM of the narrower Lorentzian exhibits the typical features associated with long range translational diffusion. This can be modeled using a jump diffusion model with a random distribution of jump lengths16, 18, 20, 23-24 𝛤(𝑄) =

ℏ 1 (1 − ) 𝜏 1 + 𝜏𝐷𝑄 2

(6)

ACS Paragon Plus Environment

9

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 27

where D is the diffusion coefficient and τ is the residence time between successive jumps of the molecule. The mean square jump length is related to the diffusion coefficient through the expression D = /6τ. The broader component shows broadening less than Q2 at higher momentum transfers indicating localized jump at the corresponding length scale. At the low Q regime, the HWHM does not approach zero at Q = 0 Å-1, and plateaus at a constant value instead. This shows evidence of the diffusion taking place in a restricted space since the widths do not continue to decay at longer length scales. For diffusion in a sphere, the confinement radius, RΓ, can be calculated using25 𝛤=

4.33296𝐷 𝑅𝛤2

(7)

where Γ is the HWHM plateau value at low Q and D is the diffusion coefficient extracted from Eq. 6. Results Cyclooctanol (OCT) The phases of OCT have been thoroughly investigated6, 10, 26-28. There are four main phases reported for OCT that are of interest: the liquid phase, PCs I and II, and the low temperature ordered crystal III. The phase transitions are reported as crystal III to PC II at 221 K (achieved by extrapolation), PC II-I at 261.7 K, and crystal I fusion at 295.3 K. For our experiment, a quasielastic signal was first observed at 190 K with very low intensity after slowly heating the sample from 20 K. Upon additional heating, it transitioned into the PC II phase for the measurements at 230 K and 240 K and subsequently the PC I phase at 270 K and 280 K. We report incoherent QENS data for protonated OCT in Fig. 1a for an average momentum transfer of Q = 1.0 Å-1 and for temperatures of 230 K, 240 K, 270 K and 280 K. The solid black lines in are the

ACS Paragon Plus Environment

10

Page 11 of 27 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

resulting fit of Eq. 1 with two Lorentzian functions. The dotted-dashed line is the resolution function that was measured at 20 K. The resolution function is strongly asymmetric as a result of the instrument design and the moderator emission spectrum14 . The EISF with respect to the broad component is displayed in Fig. 1(b). The solid lines are the fit given by Eq. 5. Data in the intermediate Q range of 1.2 and 1.3 Å-1 includes coherent structural features and was thus eliminated from the HWHM and EISF fits. An X-ray diffraction pattern for OCT at 280 K has been included as supplemental information to display this feature. The EISF for PC phases I and II look essentially very different, which indicates that a major change in the dynamics occurs upon the transition. PC phase I was satisfactorily fit with the model described in Eq. 5, whereas this model could not reproduce the EISF for PC phase II. The HWHM for the narrow component and broad component are displayed in Fig. 1c and 1d respectively and have been fit with Eq. 6. The fits for the narrow component are in rather good agreement with the experimental data. On the other hand, the broader component has some systematic deviation from the fits at the low Q regime, as discussed above. However, the diffusion coefficients and the residence times can be reliably extracted, and the values are listed in Table 1.

ACS Paragon Plus Environment

11

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 27

Figure 1. (a) QENS data for OCT at Q = 1.0 Å-1 and temperatures 230 K, 240 K, 270 K, and 280 K. The solid lines are a fit of Eq. 1 and the dotted-dashed line is the resolution function. The insert displays all the fit components for T = 280 K, with the narrow and broad Lorentzian shown in green and blue, respectively. (b) EISF for the broad Lorentzian where the solid lines are a fit of Eq. 5. Γ vs Q2 for the (c) narrow Lorentzian and (d) broad Lorentzian where the solid lines are fits using Eq. 6. Cyclohexanol (HEX) The phases of HEX have also been thoroughly investigated6, 29-30. The calorimetric study by Adachi et al30 highlights the complexity of PC phase diagrams and how temperature ramp rates effect phase formation. In general, HEX has a transition from ordered crystal III to PC II at 244 K,

ACS Paragon Plus Environment

12

Page 13 of 27 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

followed by PC II to PC I transition at 265K and PC I fusion at 299 K. The sample was measured at 20 K for the resolution function and then heated to 260 K and 280 K to measure the dynamics of PC II and PC I phases, respectively. The spectra and the fits of Eq. 4 are shown in figure 2a. The same process that was applied to treat the data for OCT was used for the HEX analysis and the HWHM and EISF are shown figures 2b-d. Again, Q values from 1.2 – 1.5 Å-1 were excluded due to coherent scattering effects. The calculated line shape given by Eq. 1 deviated from the experimental results in the high Q range for HEX. The error was due in part to the narrow Lorentzian attempting to fit the elastic intensity. Other sources for this error are presented in the discussion section. As a result, the jump diffusion model was only fit for Q ≤ 1.1 Å-1 for the narrow Lorentzian and we had to exclude the additional points of Q = 1.8 and 1.9 Å-1 from the EISF fitting procedure.

ACS Paragon Plus Environment

13

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 27

Figure 2. (a) QENS data for HEX at Q = 1.0 Å-1 and temperatures 260 K and 280 K. The solid lines are a fit of Eq. 1 and the dotted-dashed line is the resolution function. The insert displays all the fit components for T = 280 K, with the narrow and broad Lorentzian shown in green and blue, respectively. (b) EISF for the broad Lorentzian where the solid lines are a fit of Eq. 5. Γ vs Q2 for the (c) narrow Lorentzian and (d) broad Lorentzian where the solid lines are fits using Eq. 6. 80% Cyclohexanol – 20% Deuterated Cyclooctanol (80-20 Mixture) A calorimetric study on a wide range of concentrations of a binary mixture of OCT and HEX with 1 mol% lithium salt added has recently been reported by Reuter et al6. While the phase diagrams for the pure components are rather complex, the opposite was seen for the various mixtures. There was only one PC phase that was observed before melting. The authors comment that the salt seems to stabilize the PC phase. However, another investigation10 of OCT and cycloheptanol (HEPT) binary mixtures without salt also exhibited this increase in phase stability. Creating a binary mixture appears like a promising way to stabilize the PC phases. In our experiment, the sample transitioned into a PC state around 190 K and remained for the measurements at 230 K, 260 K, and 280 K. The same process that was applied to treat the data for OCT was used for the 80% protonated HEX – 20% deuterated OCT mixture, which will be referred to as the 80-20 mixture for the extent of this paper. The use of selective deuteration enables us to highlight the signal of the protonated molecules, since the incoherent neutron scattering cross section of hydrogen is about 40 times that of deuterium. In the mixture, we will observe the motion of the protonated HEX molecules and will neglect the signal from the deuterated OCT molecules. This allows us to focus on how the motion of the HEX molecules are affected by the addition of OCT to the system. The results of these measurements on the 80 – 20 mixture can be observed in Fig. 3. Note, the HWHM for 190 K was not added to Fig. 3c and 3d because the error bars were

ACS Paragon Plus Environment

14

Page 15 of 27 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

rather large and cluttered the figure. Diffusion coefficients were unable to be extracted from this data set so only the EISF will be displayed for 190 K.

Figure 3. (a) QENS data for the 80 – 20 mixture at Q = 1.0 Å-1 and temperatures 190 K, 230 K, 260 K and 280 K. The solid lines are a fit of Eq. 1 and the dotted-dashed line is the resolution function. The insert displays all the fit components for T = 280 K, with the narrow and broad Lorentzian shown in green and blue, respectively. (b) EISF for the broad Lorentzian where the solid lines are a fit of Eq. 5. Γ vs Q2 for the (c) narrow Lorentzian and (d) broad Lorentzian where the solid lines are fits using Eq. 6. Extracted Values

ACS Paragon Plus Environment

15

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 27

Table 1 displays the diffusion coefficients and residence times extracted from applying Eq. 5 to the HWHM for both the fast and slow components. The diffusion coefficients are in units of 10-11 m2/s for the α-process and 10-10 m2/s for the fast β-process. The characteristic radii for diffusion in a sphere, RCage, and diffusion on a sphere, RIntra, were extracted from fitting the EISF with Eq. 5 are listed in Table 2. These are also compared to the cage sizes, RΓ, extracted from HWHM fits. The meaning of the weight parameters, P1 and P2, have been explained in the text earlier. Note that a value of 0.25 for a weight parameter means that 75% of the hydrogen is participating in that motion and that 25% is essentially “immobile” on the measured time scale.

Table 1. Diffusion Coefficients and Residence Times Sample

T (K)

Dα (10-11 m2/s)

Dβ (10-10 m2/s)

τα (ns)

τβ (ns)

OCT (PC II)

230

0.54 ± 0.09

1.44 ± 0.17

2.02 ± 0.19

0.152 ± 0.007

240

0.70 ± 0.08

1.97 ± 0.21

1.62 ± 0.13

0.152 ± 0.006

270

4.52 ± 0.44

4.21 ± 0.45

1.43 ± 0.03

0.144 ± 0.003

280

6.87 ± 0.94

6.19 ± 0.81

1.32 ± 0.03

0.139 ± 0.003

HEX (PC II)

260

6.17 ± 0.72

7.08 ± 0.52

0.94 ± 0.04

0.139 ± 0.001

HEX (PC I)

280

6.12 ± 0.56

5.03 ± 0.89

1.28 ± 0.04

0.124 ± 0.003

80 – 20 Mix

260

(PC I)

280

OCT (PC I)

2.06 ± 0.39 3.27 ± 0.49

0.126 ± 0.004

2.08 ± 0.34

1.49 ± 0.06

0.114 ± 0.005

Table 2. EISF Fit Parameters and Calculated Cage Sizes Sample

T (K)

RCage (Å)

RIntra (Å)

P1

P2

RΓ (Å)

OCT (PC I)

270

4.52 ± 0.35

0.85 ± 0.06

0.88 ± 0.02

0.22 ± 0.05

6.44

280

3.67 ± 0.57

0.90 ± 0.10

0.85 ± 0.05

0.13 ± 0.08

7.81

ACS Paragon Plus Environment

16

Page 17 of 27 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

HEX (PC II)

260

5.11 ± 0.40

1.17 ± 0.02

0.95 ± 0.01

0.35 ± 0.01

7.53

HEX (PC I)

280

3.99 ± 0.90

1.13 ± 0.07

0.93 ± 0.04

0.17 ± 0.03

5.95

80 – 20 Mix

260

2.77 ± 0.44

0.70 ± 0.36

0.78 ± 0.13

0 ± 0.61

3.93

(PC I)

280

2.85 ± 1.23

1.05 ± 0.21

0.76 ± 0.25

0 ± 0.11

4.37

Discussion QENS investigations on PC materials typically focus on the rotation of the whole molecule or side groups, such as three-fold jumps observed in methyl groups31. We initially considered rotational motion in this study, but it became clear that we see center of mass diffusion instead. The DQ2 dependence of Γ(Q) is the characteristic feature of unconstrained long-range translational diffusion, rather than the Q-independent relationship typically seen with localized motion. While the slower motion is rather straightforward to analyze in terms of the jump diffusion model, the faster motion is more challenging. The plateau at low Q suggests that the motion becomes restricted at long length scales. Above the plateau, the broadening is described by the jump diffusion model, which suggests that this motion is a fast diffusion process within a confined geometry. A common relaxational process observed in glass formers is the fast β-process, which is often attributed to fast, sub-diffusive motions inside a cage formed by a particles nearest neighbors, also known as cage rattling11-12. This type of motion can be measured using dielectric spectroscopy at a frequency in the GHz – THz range, which appears as an excess in the minima between the main structural peaks (α and β) and the boson peak32-33. The frequency of this motion is perfectly aligned with the accessible energy range at BASIS, as illustrated by Fig. 1 in the paper by Gupta, et al17. Spherical symmetry can be used to model the cage in our system, where the nearest neighbors compose the cage boundaries. Therefore, the diffusion inside a sphere model

ACS Paragon Plus Environment

17

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 18 of 27

developed by Volino and Dianoux25 was chosen. This type of motion has been observed in a wide range of systems using QENS16, 22-23, 34. Comparing the diffusion coefficients displayed in Table 1, OCT is observed to diffuse much more efficiently in the PC I phase rather than the low temperature PC II phase, For the PC I phase, Dα is approximately an order of magnitude larger while Dβ is three times larger than in the PC II phase. In contrast, the lower temperature PC II phase in HEX has higher diffusion coefficients than the PC I phase for both processes. At first glance, this may seem surprising since diffusion rates typically increase with temperature. However, this is due to a change in the local environment between the two phases which enables the molecules to diffuse more efficiently in the PC II phase. This may be attributed to the larger cage size in the PC II phase, as seen in Table 2. The diffusion coefficients for OCT and HEX have very similar values for both processes in their respective PC I phases. To our knowledge, diffusion coefficients for the temperatures measured in this paper have not been reported for either molecule. However, self-diffusion coefficients have been reported for HEX using the pulse field gradient nuclear magnetic resonance (PFG-NMR) technique for the temperature range of 350 – 440 K35. Diffusivities measured by QENS are usually larger than those measured by PFG-NMR since these techniques measure different length scales3637

. Neutrons probe length scales from few Å to few dozen Å, while PFG-NMR probes much longer

length scales up to m order. Extrapolating this data to lower temperatures gives an estimated value of Dα = 2.85 ∙ 10-11 m2/s for 280 K. We report a diffusion coefficient of 6.12 ± 0.56 ∙ 10-11 m2/s at this temperature. Upon addition of OCT molecules into the HEX solution (80 – 20 mixture), both diffusion rates decreased by about half for 280K. Dα for the mixture is now in close agreement with the extrapolated value for pure HEX from ref. 17. At the temperatures below 280 K for the 80-20 mixture, the broadening for the α-process deviates from the DQ2 behavior and becomes Q-

ACS Paragon Plus Environment

18

Page 19 of 27 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

independent, which prevents diffusion coefficients to be calculated. As the temperature is decreased, it becomes increasingly more difficult for the molecule to overcome the jump barriers and the long-range diffusion becomes restricted. The molecule is confined to a lattice site where it can only participate in small scale motions, such as rotations. This feature can also be observed in the widths for the fast β-process in Fig. 3d. For the highest temperature, 280K, Γ(Q) covers a range from approximately 19 – 35 µeV with a clear diffusive relationship. However, this range becomes increasingly smaller as the temperature decreases and approaches a constant value at 230K. At this temperature, Γ(Q) only varies from 25 – 29.5 µeV and is considered constant. We would like to briefly address the poor fitting results at the high Q-values for the slow process in HEX, as seen in Fig. 2c. The jump diffusion model has a DQ2 dependence at low Q and plateaus to a value at high Q, as seen in OCT and the 80 – 20 mixture. However, the experimental observed Γ(Q) values decrease below the theoretical value given by the jump diffusion model at large momentum transfers for HEX. One cause for this is that the narrow Lorentzian was attempting to fit the elastic signal. The theoretical fit achieved using Eq. 1 was not perfect for these Q-values and it is possible that the experimental line shape is not sufficiently described by Lorentzian functions. Another reason is due to the limitations of the jump diffusion model, which in many liquid and glass forming systems can only accurately describe data up to Q ≈ 1 Å-1 22, 34. Above this point, measurements transition from observing the center of mass motion and focus on intramolecular motion instead. The assumption that the diffusive jumps are infinitesimally small becomes invalid. A more sophisticated model that includes the precise jump distances of the hydrogen atoms is necessary to model the diffusion in this momentum range. However, we are successful in explaining the low Q-values by the jump diffusion model. This enables us to accurately calculate diffusion coefficients by only using the low Q-values in the fitting procedure.

ACS Paragon Plus Environment

19

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 27

The EISF for the lower temperatures, i.e. T < 240 K, was modeled using Eq. 4 for both OCT and the 80 – 20 mixture. The overall quality of the fit was relatively good, but the extracted length scales of RCage and RIntra were very small. This suggests that the motion is not diffusive, and rather localized since the molecule is confined to a lattice site. Attempts were made to fit this data with other models, such as rotational jumps on a circle, all of which proved to be unsuccessful in describing the data. We believe that one of these models should adequately explain the data and that the failure is due to the very low quasielastic intensity measured at the low temperatures, which made the fitting procedure difficult. The values of RIntra displayed in Table 2 were extracted from Eq. 4. The values range from 0.7 – 1.17 Å, which is close to the C-H bond distance of 1.09 Å. The value for 0.7 Å is small but acceptable once the error of 0.356 Å is considered. This type of motion is attributed to the various small-scale motions of the hydrogen atoms about the carbon atoms. The effect of this motion in the data becomes more prominent at high Q, which is why we needed the extra correctional term in Eq. 4. Analyzing the weight parameters P1 and P2, there is the general trend that the parameter increases as the temperature is decreased. This is exactly what we would expect to see since the weight parameter represents the fraction of hydrogen in the system that is considered immobile in the energy window of the observed relaxational process. As the temperature is lowered, the motion becomes frozen in and more hydrogen become immobile. The small values for P2 suggests that a large majority of the hydrogen is participating in the motion, which supports the idea that these are just small amplitude motions of the hydrogens around the carbon atom. The fast β-process has been thought to be a prerequisite for the α-process since these small length scale relaxations can increase the free volume and open the cage. Confining the diffusive motion inside the cage will hinder the ability of the local environment to properly rearrange and cages will

ACS Paragon Plus Environment

20

Page 21 of 27 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

not sufficiently open to allow for long range translational diffusion. As seen in Table 2, the average cage sizes, RCage, for HEX and OCT are determined to be 4.5 and 4.1 Å, respectively. The diameters of the molecules are approximately 4.96 and 5.9 Å for HEX and OCT, respectively. Considering the sizes of the OCT molecule and the cage, the molecule is free to move about its center of mass in any direction with an amplitude around 1 Å inside the cage. The HEX molecule has even more freedom than OCT to move around the cage, approximately 2 Å, which explains the observed faster dynamics. These distances seem reasonable in comparison to the mean-squared displacements derived from neutron scattering experiments on various glass formers22, 38-39. In one of these studies, Mamontov measured the glass former glycerol at BASIS to investigate the diffusive dynamics in the liquid state at 360 K22. Identical to our investigation, two distinct motions were observed and attributed to a slower, long-range diffusion and a faster, cage rattling process. Using the diffusion within a sphere model, the calculated cage size, RCage, was determined to be 8.0 Å, which was close to the estimated a cage size of 5.7 – 7.8 Å. As expected, this value is much larger than the values we report for HEX and OCT since glycerol was measured in the liquid state. The main goal of this investigation is to determine how the dynamics change when a binary PC mixture is created. We have already reported the average cage sizes for the pure samples of HEX and OCT to be 4.5 and 4.1 Å, respectively. Upon addition of OCT into the HEX solution, the average cage size was reduced to only 2.8 Å. The HEX molecules are now largely confined to their lattice positions, with freedom to move approximately 0.3 Å in any direction. This increased confinement is most likely the reason why there is such a drastic decrease seen in the self-diffusion coefficients for the 80 – 20 mixture. The larger OCT molecules reduce the amount of free volume available for HEX movement inside the cage. This reduces the diffusion inside the cage and results in the less efficient formation of open cage doors for the particle to break free from the cage and

ACS Paragon Plus Environment

21

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 27

contribute to the α-process. In this framework, the fast β-process has a large effect on the success of the main structural relaxation in glass formers. This result supports the conclusion by Reuter et al. on the ionic conductivity in HEX-OCT binary PC mixtures6. They observed a decrease in the ionic conductivity of Li ions in the system with increasing concentrations of the larger OCT molecules in the solution. They argued that this is due to a decrease in the reorientational relaxation times of the PC matrix. In the current investigation, we observed a decrease in two different translational diffusion processes, which in addition to the slower rotation of the molecules, shows that the larger OCT molecules reduce the free volume and decrease the overall mobility of the smaller HEX molecules. This makes the PC matrix less efficient at conducting the Li ions. In addition, the features of increased phase stability and an extension of the temperature range for the PC phase, as observed in the HEX-OCT6 and HEPT-OCT10 systems, is probably a result of this increased confinement. While fitting the EISF with Eq. 5 is one way to estimate the cage size, there is a second method available that involves fitting Γ(Q) for the broad Lorentzian with Eq. 6 instead. Once Dβ is obtained, it can be inserted into Eq. 7 along with the plateau value of Γ(Q) at low Q. This will yield the confinement radius for the fast diffusion process. The values extracted from both methods are compared in Table 2. It is interesting that the distances from fitting Γ(Q), RΓ, are larger than those from fitting the EISF, RCage. However, both methods produce the same general result where the cage size for the 80 – 20 mixture is smaller than the pure samples. We believe that the RCage values are more realistic since the amplitude of cage rattling should be less than the diameter of the particle in accordance with MCT11-12. For PC materials, the ionic conductivity is dominated by the revolving door mechanism, especially at lower temperatures where translational movement is confined. Götz et al. argues that

ACS Paragon Plus Environment

22

Page 23 of 27 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

there should be no significant translational diffusion in PC systems since the centers of mass are completely fixed on lattice sites with strict symmetry40. This statement is true for low enough temperatures, such as the measurements at T < 240 K in this investigation. However, we have observed that there is significant translational movement in the PC phases at temperatures just below the liquid phase, i.e. T ≈ Tm – 40 K. The average position of the molecules is still considered to be positioned on lattice sites, but they can now make small amplitude motions inside the cage, as we have shown. The results in this paper show that the cage size influences both the fast, shortrange diffusion inside the cage and the slower, long-range translational diffusion between cages. Manipulating the cage size, which is possible through the formation of a binary mixture of varying concentration, can be a useful tool for controlling the local environment in PC systems. This can be used to tune the ionic conductivity for PC-based, solid state ionic conductors. It should be noted that rotation of the molecules is still very important here, but we have demonstrated that manipulating the cage size can be an additional way to control the diffusion of ions in a PC system. Our findings confirm a clear change in the observed dynamics when the mixture is formed, and a more in-depth QENS investigation of this concentration dependent behavior is needed over a wide range of binary mixtures to obtain a generalized behavior in such PCs. Conclusion In the present work, we have investigated the well-known PC systems of HEX and OCT in their pure forms and in an 80% HEX – 20% OCT mixture using QENS. Two distinct center of mass motions were observed, i.e. long-range translational diffusion (α-process) and cage rattling (fast β-process). A jump diffusion model was employed to determine diffusion coefficients and residence times for each sample. It was found that the diffusion coefficients for both processes decreased upon the formation of the 80 – 20 mixture. The diffusion within a sphere model was

ACS Paragon Plus Environment

23

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 24 of 27

applied to the EISF to determine the cage sizes for the fast β-process. The cage sizes were determined to be an average of 4.1, 4.5, and 2.8 Å for the OCT, HEX, and 80 – 20 mixture, respectively. The reduced diffusion coefficients for the 80 – 20 mixture are attributed to the increased confinement of the particle inside the cage, as observed by a reduction in cage size. A reduction in the rotational and translational mobility of the molecule hinders the ability of the PC matrix to efficiently conduct ions in this type of system. Therefore, manipulation of the cage size is possible by forming a binary PC mixture, which by altering the concentration of the constituents, can be a useful tool to control the ionic conductivity in PC-based, solid state ionic conductors. Acknowledgements Research conducted at ORNL's SNS and CNMS was sponsored by the Scientific User Facilities Division, Office of Basic Energy Sciences, US Department of Energy. Financial support was provided by Jülich Center for Neutron Sciences (JCNS). TE was supported by the US Department of Energy, Office of Science, Basic Energy Sciences, Materials Science and Engineering Division. We would like to thank Dr. Peter Lunkenheimer and Dr. Alois Loidl for their valuable discussions and to Dana Vieweg for the X-ray diffraction measurements.

Author Information Corresponding Author *E-mail: [email protected] Supporting Information X-ray diffraction spectra for cyclooctanol to demonstrate the large coherent signal in the intermediate Q-range.

ACS Paragon Plus Environment

24

Page 25 of 27 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

References 1. Tarascon, J. M.; Armand, M., Issues and challenges facing rechargeable lithium batteries. Nature 2001, 414 (6861), 359-367. 2. Motavalli, J., Technology: A solid future. Nature 2015, 526 (7575), S96-S97. 3. Loidl, A., Orientational Glasses. Annual Review of Physical Chemistry 1989, 40 (1), 29-60. 4. Geirhos, K.; Lunkenheimer, P.; Michl, M.; Reuter, D.; Loidl, A., Communication: Conductivity enhancement in plastic-crystalline solid-state electrolytes. The Journal of Chemical Physics 2015, 143 (8), 081101. 5. MacFarlane, D. R.; Forsyth, M., Plastic crystal electrolyte materials: new perspectives on solid state ionics. Advanced materials 2001, 13 (12-13), 957-966. 6. Reuter, D.; Geiß, C.; Lunkenheimer, P.; Loidl, A., Variation of ionic conductivity in a plastic-crystalline mixture. The Journal of Chemical Physics 2017, 147 (10), 104502. 7. Nilsson, L.; Thomas, J. O.; Tofield, B. C., The structure of the high-temperature solid electrolyte lithium sulphate at 908K. Journal of Physics C: Solid State Physics 1980, 13 (35), 6441. 8. Alarco, P.-J.; Abu-Lebdeh, Y.; Abouimrane, A.; Armand, M., The plastic-crystalline phase of succinonitrile as a universal matrix for solid-state ionic conductors. Nat Mater 2004, 3 (7), 476-481. 9. MacFarlane, D. R.; Huang, J.; Forsyth, M., Lithium-doped plastic crystal electrolytes exhibiting fast ion conduction for secondary batteries. Nature 1999, 402 (6763), 792-794. 10. Rute, M. A.; Salud, J.; Negrier, P.; Lopez, D. O.; Tamarit, J. L.; Puertas, R.; Barrio, M.; Mondieig, D., Twocomponent system cycloheptanol (C7) + cyclooctanol (C8): an extraordinary system. Journal of Physical Chemistry B 2003, 107 (24), 5914-5921. 11. Donth, E.-J., The Glass Transition : Relaxation Dynamics in Liquids and Disordered Materials. Berlin, Heidelberg : Springer Berlin Heidelberg : Imprint: Springer: Berlin ; New York, 2001. 12. Götze, W., The essentials of the mode-coupling theory for glassy dynamics. Condensed Matter Physics 1998, 1 (4), 873-904. 13. Barrat, J. L., The liquid-glass transition of the hard-sphere system. Journal of Physics: Condensed Matter 1989, 1 (39), 7163-7170. 14. Mamontov, E.; Herwig, K. W., A time-of-flight backscattering spectrometer at the Spallation Neutron Source, BASIS. Review of Scientific Instruments 2011, 82 (8), 085109. 15. R.T. Azuah, L. R. K., Y. Qiu, P.L.W. Tregenna-Piggott, C.M. Brown, J.R.D. Copley, and R.M. Dimeo, DAVE: A comprehensive software suite for the reduction, visualization, and analysis of low energy neutron spectroscopic data J. Res. Natl. Inst. Stan. Technol. 2009, 114 (341). 16. Mamontov, E.; Luo, H.; Dai, S., Proton dynamics in N,N,N',N'-Tetramethylguanidinium Bis(perfluoroethylsulfonyl)imide protic ionic liquid probed by quasielastic neutron scattering. J Phys Chem B 2009, 113 (1), 159-69. 17. Gupta, S.; Arend, N.; Lunkenheimer, P.; Loidl, A.; Stingaciu, L.; Jalarvo, N.; Mamontov, E.; Ohl, M., Excess wing in glass-forming glycerol and LiCl-glycerol mixtures detected by neutron scattering. Eur Phys J E Soft Matter 2015, 38 (1), 1. 18. Gupta, S.; Mamontov, E.; Jalarvo, N.; Stingaciu, L.; Ohl, M., Characteristic length scales of the secondary relaxations in glass-forming glycerol. Eur Phys J E Soft Matter 2016, 39 (3), 40. 19. Gupta, S.; Fischer, J. K.; Lunkenheimer, P.; Loidl, A.; Novak, E.; Jalarvo, N.; Ohl, M., Effect of adding nanometre-sized heterogeneities on the structural dynamics and the excess wing of a molecular glass former. Sci Rep 2016, 6, 35034. 20. Bée, M., Quasielastic Neutron Scattering. Hilger: Bristol, UK, 1988. 21. Bée, M., A physical insight into the elastic incoherent structure factor. Physica B: Physics of Condensed Matter 1992, 182 (4), 323-336. 22. Mamontov, E., Diffusion in confinement as a microscopic relaxation mechanism in glass-forming liquids. Chemical Physics Letters 2012, 530. 23. López-Cabarcos, E.; Batallán, F.; Frick, B.; Ezquerra, T.; Calleja Fj, B., Molecular dynamics of ferroelectric polymeric systems as studied by incoherent quasielastic neutron scattering. Physical Review B 1994, 50 (18), 1321413224. 24. Singwi, K. S.; Sjölander, A., Diffusive Motions in Water and Cold Neutron Scattering. Physical Review 1960, 119 (3), 863-871.

ACS Paragon Plus Environment

25

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 26 of 27

25. Volino, F.; Dianoux, A. J., Neutron incoherent scattering law for diffusion in a potential of spherical symmetry: general formalism and application to diffusion inside a sphere. Molecular Physics 1980, 41 (2), 271-279. 26. Andersson, O.; Ross, R. G., Thermal conductivity, heat capacity and phase diagram of cyclooctanol in liquid, solid and glassy crystal states under high pressure. Molecular Physics 1990, 71 (3), 523-539. 27. Osamu, Y., Calorimetric and neutron scattering studies of plastically crystalline cyclooctanol*. Journal of Physics: Condensed Matter 2003, 15 (32), 5439-5450. 28. Edelmann, R.; Würflinger, A., Phase Behaviour of Cyclo-octanol up to Pressures of 300 MPa and Temperatures between 100 and 370 K. Molecular Crystals and Liquid Crystals 1987, 148 (1), 249-254. 29. Adachi, K.; Suga, H.; Seki, S.; Kubota, S.; Yamaguchi, S.; Yano, O.; Wada, Y., Dielectric Relaxations in Various Crystal Modifications of Cyclohexanol and 2,3-Dimethylbutane. Molecular Crystals and Liquid Crystals 1972, 18 (3-4), 345-368. 30. Keiichiro Adachi, H. S., Syûzô Seki, Phase Changes in Crystalline and Glassy-Crystalline Cyclohexanol. Bulletin of the Chemical Society of Japan 1968, 41, 1073-1087. 31. Bermejo, F. J.; Jiménez-Ruiz, M.; Criado, A.; Cuello, G. J.; Cabrillo, C.; Trouw, F. R.; Fernández-Perea, R.; Löwen, H.; Fischer, H. E., Rotational freezing in plastic crystals: a model system for investigating the dynamics of the glass transition. Journal of Physics: Condensed Matter 2000, 12 (8A), A391. 32. Brand, R.; Lunkenheimer, P.; Loidl, A., Relaxations and fast dynamics of the plastic crystal cyclo-octanol investigated by broadband dielectric spectroscopy. Physical Review B 1997, 56 (10), R5713-R5716. 33. Lunkenheimer, P.; Brand, R.; Schneider, U.; Loidl, A., Excess wing and high frequency dynamics in plastic crystals. Philosophical Magazine B 1999, 79 (11-12), 1945-1951. 34. Volino, F.; Pineri, M.; Dianoux, A. J.; De Geyer, A., Water mobility in a water-soaked nafion® membrane: A high-resolution neutron quasielastic study. Journal of Polymer Science: Polymer Physics Edition 1982, 20 (3), 481496. 35. Reilly, D. E.; Peterson, E. M.; Hogenboom, D. L., Self‐Diffusion Coefficients and Rotational Correlation Times in Polar Liquids. V. Cyclohexane, Cyclohexanone, and Cyclohexanol. The Journal of Chemical Physics 1972, 57 (9), 3969-3976. 36. Jobic, H., QENS and the Benefit of Diffusion Measurements over Different Length Scales. Diffusion Fundamentals 2007, 6, 52.1-52.8. 37. Paoli, H.; Methivier, A.; Jobic, H.; Krause, K.; Pfeifer, H.; Stallmach, F.; Karger, J., Comparative QENS and PFG NMR diffusion studies of water in zeolite NaCaA. Microporous and Mesoporous Materials 2002, 55, 147-158. 38. Ngai, K. L., Why the fast relaxation in the picosecond to nanosecond time range can sense the glass transition. Philosophical Magazine 2004, 84 (13-16), 1341-1353. 39. Simmons, D. S.; Cicerone, M. T.; Zhong, Q.; Tyagi, M.; Douglas, J. F., Generalized localization model of relaxation in glass-forming liquids. Soft Matter 2012, 8 (45), 11455-11461. 40. Götz, M.; Bauer, T.; Lunkenheimer, P.; Loidl, A., Supercooled-liquid and plastic-crystalline state in succinonitrile-glutaronitrile mixtures. The Journal of Chemical Physics 2014, 140 (9), 094504.

Table of Contents Graphic

ACS Paragon Plus Environment

26

Page 27 of 27 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

ACS Paragon Plus Environment

27