scieee AI-readable full text Open interactive document viewer

Implanted muon spin spectroscopy on 2-O-adamantane: a model system that mimics the liquid

Romanini, Michela,Tamarit Mur, José Luis,Pardo Soto, Luis Carlos,Bermejo, F. J.,Fernandez Perea, Ricardo,Pratt, F. L.

Abstract

The transition taking place between two metastable phases in 2-O-adamantane, namely the [Formula: see text] cubic, rotator phase and the lower temperature P21/c, Z = 4 substitutionally disordered crystal is studied by means of muon spin rotation and relaxation techniques. Measurements carried out under zero, weak transverse and longitudinal fields reveal a temperature dependence of the relaxation parameters strikingly similar to those exhibited by structural glass[Formula: see text]liquid transitions (Bermejo et al 2004 Phys. Rev. B 70 214202; Cabrillo et al 2003 Phys. Rev. B 67 184201). The observed behaviour manifests itself as a square root singularity in the relaxation rates pointing towards some critical temperature which for amorphous systems is located some tens of degrees above that shown as the characteristic transition temperature if studied by thermodynamic means. The implications of such findings in the context of current theoretical approaches concerning the canonical liquid-glass transition are discussed.

Full text

Implanted muon spin spectroscopy on 2-O-adamantane: A model system that mimics the liquid−→glasslike transitions M. Romanini, J.Ll Tamarit, L.C. Pardo Grup de Caracterizacio de Materials, Departament de Fisica, ETSEIB, Universitat Politecnica de Catalunya, Diagonal 647, 08028 Barcelona, Catalonia, Spain F.J. Bermejo, R. Fernandez-Perea Instituto de Estructura de la Materia,Consejo Superior de Investigaciones Cientificas, Serrano 123, E-20886 Madrid, Spain F.L. Pratt ISIS Neutron and Muon Source, STFC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire, OX11, 0QX, United Kingdom December 2, 2016 Abstract The transition taking place between two metastable phases in 2-O-adamantane, namely the Fm¯ 3mcubic, rotator phase and the lower temperature P21/c, Z = 4 substitutionally disordered crystal is studied by means of muon spin rotation and relaxation techniques. Measurements carried out under zero, weak transverse and longitudinal fields reveal a temperature dependence of the relaxation parameters strikingly similar to those exhibited by structural glass−→liquid transitions (Bermejo et al. Phys. Rev. B 70, 214202 (2004); Cabrillo et al. Phys. Rev. B 67, 184201 (2003)). The observed behaviour manifests itself as a square root singularity in the relaxation rates pointing towards some critical temperature which for amorphous systems is located some tens of degrees above that shown as the characteristic transition temperature if studied by thermodynamic means. The implications of such findings in the context of current theoretical approaches concerning the canonical liquid-glass transition are discussed. Keywords disordered systems, muon spin rotation, relaxation and resonance 1 Introduction The interest in the study of molecular materials showing some sort of disorder in one degree of freedom, while retaining ordered crystal structures for the rest, has reawakened in recent decades. This is mostly because of the possibilities such materials offer as model systems to extend our understanding of the properties exhibited by structurally disordered matter, but also due to some potential applications that were recently envisaged [1]. Materials prepared into such states are known to exhibit at low temperatures thermal, acoustic and dielectric properties which significantly differ from those shown by the same substances if prepared within their fully ordered, crystalline ground states. However, after some decades of extensive studies, our understanding of most features commonly ascribed to glassy dynamics mostly relies on phenomenological concepts [2], many of which have been developed as tools to quantitatively describe the physical state of simpler systems such as substitutional defects in alkali halides [3] as well as mixed halide crystals [4]. Such model systems display a full range of phenomena well understood microscopically, which go from tunneling motions of a single atom to interacting tuneling states in cases 1 with large concentrations of substitutional defects, which were shown to exhibit fully developed glassy behaviour . The phenomenology considered to be the distinctive signature of the glassy state is nowadays known to be shared by a wide number of systems [5] which go from partially ordered crystals to concentrated colloidal solutions, cold compressed crystals or idealized toy-models [6]. Several of such studies on partially ordered crystals have shown [7] these materials to exhibit phenomena remarkably close to those shown by the same material if prepared as a fully disordered substance. As a further step, recent additions to the literature [8] include incommensurately modulated crystals which exhibit at low temperatures similar features to those shown by amorphous materials, and are found to originate from the gapped phase and amplitude modes of the incommensurate crystal structure. Here we report on issues pertaining to a transition between disordered crystal states on the adamantane derivative 2-O-adamantane (C10H14O, also referred to as 2-adamantanone), which retains significant analogies with the canonical liquid −→ glass transition . The material, which can be prepared into different phases when cooling down the liquid from 557 K, exhibits a rich polymorphism at ambient pressure [9–13]. In fact, freezing the stable liquid leads to a rotator-phase (RP) or plastic-crystal state [9] where individual molecules having C2vpoint group symmetry sit at the nodes of a face-centered Fm¯ 3m cubic lattice and execute large amplitude, quasi-isotropic orientational excursions involving jumps between six possible orientations of the C=O group along the six fourfold crystal axes as well as π/2flops about the C=O axis. Further cooling, below 178 K leads to a metastable, P21c, Z = 4 monoclinic unit cell where statistical disorder is manifest by the site occupancy of the oxygen atom which may lie in three different crystallographic sites with fractional occupancies of 25%, 25% and 50%, here referred to as substitutionally-disordered phase (SDP). A true crystalline, fully-ordered ground state is only attained after repeated thermal cycles comprising the transition temperatures ≃100 K - 220 K leading to repeated RP SDP interconversions. Relatively fast molecular rotations were reported for adamantane derivatives and adamantanone in particular [14–16], which show that in contrast to rapid isotropic rotational reorientations found for the former in liquid or even within the solid at 294 K where rotational diffusion constants still are of the order of 1010rad2s−1, rather extreme anisotropic reorientations were found for the latter. Interestingly, while the isolated molecule closely behaves as a symmetric top with calculated rotational constants of 1678.6 MHz, 1201.5 MHz, and 1191.6 MHz, the principal values of the rotational diffusion tensor as determined from NMR spectroscopy in the liquid phase [14] unveil a marked motional anisotropy with a proportion 1:2:24, where the faster rate corresponds to reorientations along an axis passing through the C=O bond. The interest in carrying out a detailed study on the dynamics of the SDP −→ RP transition stems from findings previously reported by Brand et al. [10] and recently confirmed by Romanini at al. [11] which report on fairly well defined glassy dynamics features within the SDP phase by means of broadband dielectric spectroscopy, which were tentatively assigned to interconversions between the statistically disordered sites. Such findings contrast with previous spectroscopic results [17] which overlooked such microscopic motions most probably due to their available time windows. The presence of such motions thus suggests the non-ergodic nature of the SDP phase below a certain glass transition temperature [11] and therefore constitutes a valid analogue of a glass, since the former state is attained by virtue of the continuous slow down of the nearly free rotations present in the RP crystal. The present paper reports on results from measurements about the SDP −→ RP transition carried out by means of a spectroscopic technique able to detect single-particle motions taking place over a quite broad range of characteristic times (10−12 −10−4s) . In fact, here we show that spin relaxation spectroscopy with implanted positive muons provides a unique means to monitor the phase transformation referred to above. Previous studies on canonical glass −→ liquid transitions [18,19] have shown that measurements of the polarization function of the muon decay in weak transverse magnetic fields display unambiguous signatures of critical behaviour at temperatures somewhat above those signaling the thermodynamic glass transition, namely at the onset of stochastic molecular motions. The temperature dependence of transverse relaxation rates and amplitudes is found to provide sound evidence for the existence of a critical temperature Tcas predicted by kinetic theories of the glass transition. 2 2 Materials and methods The technique here employed consists in the measurement of the temperature and applied field dependences of the spin relaxation characteristics of implanted positive muons onto a sample of 2-O- adamantane previously prepared within the crystal state of interest. Because of the small muon mass, once the incoming muon reaches the end of the ionization track and gets fully thermalized, its motions once implanted at specific sites of the crystal structure will be driven by the excitations of its microscopic environment to which it is strongly coupled. The latter depends upon known radiochemical details pertaining to the fate of the implanted muon. The measurements were carried out using the ARGUS muon spin relaxation (µSR) spectrometer hosted by the ISIS pulsed source. The instrument was set up in a standard µSR configuration. The incoming beam was collimated down to a spot of 30 mm and spectra were left to accumulate some 10 million events. The sample was purchased from Aldrich Chemical with a stated purity of 99 % and used as received. Specimens for the muon experiments were prepared by enclosing some 0.8 g of the pristine powdered material within a 3 cm ×3 cm silver foil packet. Spectra under longitudinal (LF) and transverse (TF) fields were then measured from 20 K up to temperatures well above the SDP −→ RP transition, namely 300 K. The analysis of experimental data has been carried out employing software packages specially developed at ISIS, i.e. the WIMDA [20] and MantidPlot packages [21]. 3 Results 3.1 Characterization of muon sites From the chemical nature of the sample we are dealing with, several kinds of muon species are to be expected. Some muons may remain as bare µ+, whereas others will capture an electron, leading to formation of atomic Muonium Mu(≡µ+e−). Both these species will sit within the crystal interstices and may execute diffusive or jumping motions. Both µ+and Mu will stick to 2-O-adamantane molecules near the Oxygen attracted by its electronegativity [22]. The final state of such a muoniated adduct may either be as a charged molecule 2-O-adamantane-Mu+, closed shell in character and thus diamagnetic or rather a neutral 2-O-adamantane-Mu, which will probably have some distribution of unpaired spin density and is thus paramagnetic. Further implantation sites leading to the formation of adamantanyl radicals [23,24] are deemed to be far less abundant on account of the rather low chemical affinity of the known implantation sites if compared to that of the carbonyl group. A series of calculations employing quantum chemical tools have been carried out in order to clarify the chemical fate of the implanted muons as well as the properties of the muoniated molecules. The CASTEP [25] and DALTON [26] codes have been employed for such a purpose. The geometry and some properties of isolated molecules have been first calculated either at Hartree-Fock or DFT/B3LYP levels in order to get some indications on their relative stability as well as to guess some relevant spectroscopic properties such as the values of the hyperfine coupling parameters. The obtained results show that both neutral and charged muoniated species appear to be stable enough having a difference in energies of roughly 4 eV, depending upon the computational method employed for the purpose [27]. The obtained fully relaxed molecular geometries show the muon rest place as attached to the carbonyl Oxygen, separated by a distance of 0.97 Åmaking a C-O-Mu bond angle of 113.3 ◦with the muon staying within the C-C-O plane. Figure 1 displays the results of one calculation carried out using the CASTEP code showing the muon bonded to one of the carbonyl Oxygens. Detailed calculations show that the muoniated molecule is a symmetric top with a ratio of its principal moments of inertia of 0.7:1:1 with 72 molecular vibrational modes with frequencies comprising 19.3 meV - 463.8 meV. The lowest vibrational frequency corresponds to a normal coordinate in which several atoms including the stuck muon execute low amplitude motions. For the temperature range of interest however, the vibrational contributions of an isolated molecule to the total molecular partition functions are of the order of 2 ×10−6or below if compared to their rotational counterparts. 3.2 Properties of the muoniated radical Previous results on the chemically related molecule acetone [28] report on data where some 53 per cent of the implanted muons end up within diamagnetic molecular environments. Such a relatively large 3 Figure 1: (a) Molecular structure of the 2-O-adamantane molecule. Oxygen is shown as a red/white sphere, carbons as brown and hydrogens as pale brown spheres. (b) Structure of the molecule following addition of µ+or Mu. The bonded µ+/Mu is depicted as a pale blue/white atom-like sphere. (c) The arrangement of molecules within the P21c, Z = 4 monoclinic crystal structure showing multiple possible sites for the O (hydrogens are omitted for clarity). The colour code for the Oxygen and Mu is meant to give an indication of the fractional occupancy probabilities, that is the pale-blue and red depict occupational probabilities of a given species in such a crystal site. diamagnetic fraction is thus thought to arise from 2-O-adamantane-Mu+cations. Following the track of such studies which mostly focus on radiochemical details, we first attempted a measurement using the technique of avoided-level-crossing (ALC) resonance for the sought 2-O-adamantane paramagnetic radical fraction. The expectancy was to find some signature of the ∆1resonance which corresponds to muon spin flips induced through the coupling of the Zeeman states from the anisotropic part of the muon-electron hyperfine interaction as a dip in the signal amplitude. For such an avail we took the value for the muon hyperfine coupling constant (HFCC) of the Mu adduct of acetone of 25.4±0.6 MHz at 298 K [28] on the basis of the local chemical analogy of both compounds, leading to an expected resonance centered at about 930 G. If this were proven to be the case then the onset of rotation of the C=O bond would result in the muon HFCC changing significantly. Molecular rotations would then significantly reduce the dipolar coupling between the bonded muon and nearby protons and lead to a reduction in the width and amplitude of the resonance. Initial measurements of the relaxation parameters measured under longitudinal (LF) fields up to 4 kG and a full set of temperatures both within the SDP and the RP ranges of existence did not reveal any clear signature of the presence of such a resonance in the form of a sharp dip. To clarify this apparent discrepancy, a series of computations were carried out in order to derive estimates of the hyperfine coupling constants for the 2-O-adamantane-Mu adduct. Precise details about such calculations as well as the resulting estimates for the principal values of the symmetrized (i.e. diagonal) hyperfine coupling (HFC) tensor at the muon site and the axial and non-axial dipolar anisotropy terms are given as Supplementary Material. The results show that the estimation of the isotropic HFCCs are extremely sensitive to the local molecular geometry about the muon binding site and point towards a remarkable dependence of their values upon molecular internal motions which usually vibrate at frequencies which are comparable to those of the molecular lattice modes. In summary, our present findings suggest that the radical state shows a strong dependence upon molecular motions which may provide an interesting handle for further studies on the temperature dependence of the hyperfine coupling parameters thus leading to estimates for the relevant microscopic 4 correlation times. In what follows however, we will mostly rely upon data measured under transverse fields, whose interpretation is much more straightforward, taking data measured under zero field as a consistency test which is also reported on as Supplementary Material. 3.3 Transverse field studies Let us first discuss the results pertaining to the TF setup. Such data if taken at a pulsed muon source such as ISIS will provide access to details of the dynamics of diamagnetic species alone, due to the limited frequency window for rotation signals, whereas measurements carried out under longitudinal fields will also sense the paramagnetic muonium via the time-averaged polarisation. The spectra yield relaxation rates λT F (T)and their corresponding amplitudes which contain information on static or frozen-in spatial fluctuations of the local magnetic field since they induce dephasing among the spin precession of muons located at different sites. In other words, the relaxation is also driven by a zero-frequency term J(0) . Such spectral density is also sampled by experiments carried out under longitudinal fields (LF) which sample the J(ω0)spectral density of the local magnetic field fluctuations at the Larmor frequency for the muon states (or, for muonium, at various transition frequencies of the hyperfine-coupled spin states). Also, in the absence of any externally applied field, the implanted muons will precess at a frequency set by the internal fields. The physics sensed by such zero-field (ZF) experiments would be quite similar to the TF case although any contribution from paramagnetic species is to be expected. Signatures of such latter entities are in fact observed and a more detailed account is deferred to the relevant section within the Supplementary Material. Some representative histograms showing the time evolution of the muon polarization comprising two significant temperatures within the range of interest are shown in Fig.2. As expected, measurements carried out at the lower range of temperatures show damped oscillations in their decay curves whereas those measured for temperatures above the SDP −→ RP transition show small decay rates with a rather soft dependence with temperature. The transverse field spectra can be well approximated by a damped cosine signal namely, PT F (t) = aT F exp (−λT F t)cos(ω0t+φ) + bcgr (1) where ω0stands for the Larmor frequency corresponding to the effective field sensed at the muon site ω0=γµBeff ,γµ= 2π×13.55 kHzG−1stands for the muon gyromagnetic ratio and Beff is the field experienced at the muon site, which for a diamagnetic material is expected to be rather close to the applied external field, φis a phase shift related to the detector geometry and signal timing and λT F is a relaxation rate which accounts for both the distribution of local fields and their time-dependent fluctuations. A small background term bcgr was left as a free parameter in a first round of spectral fits. Its value was then set to the average taken over the whole temperature range and left as a constant in subsequent fits. A graph displaying the temperature dependence of the relaxation rate λT F and relative amplitude aT F as measured for increasing temperatures under an applied transverse field of 50 G is shown in Fig.3. Both quantities show an abrupt change when the transition SDP −→ RP is approached. Such features are here far more marked than those already reported for canonical glass −→ liquid transitions [18,19] and discussion of their physical relevance is deferred to the next section. As stated above, only diamagnetic species contribute to the measured intensity at a pulsed source such as ISIS under the applied fields here employed and therefore the interpretation should be straightforward. In fact, the technique is equivalent to the measurement of free induction decays in conventional 1H NMR and corresponds to precession of the muon spin at or near its Larmor frequency. In turn the much higher precession frequencies of atomic muonium can only be observed at continuous muon sources. Such measurements are sensitive to static, (i.e. time averaged ) spatial disorder and therefore can be compared to the second moment of the proton NMR line, M2(T)[29] which for our sample was studied by Amoreux and others in ref [17]. The reported rigid-lattice value for the SDP is some 21 G2larger than that for the RP. Such a difference is also sensed by our λT F data down to 10 K. However, contrary to the NMR results which show an abrupt jump at the crystal transition and basically no dependence with temperature down to a few kelvins, our data displayed in Fig.3 do show a marked dependence with decreasing temperature. The reason behind such different results coming out from the two techniques is well established and has to deal with the rather different rates of dynamic fluctuations explorable by 5 2 4 6 8 1 0 1 2 - 0 . 2 - 0 . 1 0 . 0 0 . 1 0 . 2 2 4 6 8 1 0 1 2 - 0 . 2 - 0 . 1 0 . 0 0 . 1 0 . 2 T i m e ( µs ) A s y m m e t r y 2 0 K T i m e ( µs ) A s y m m e t r y 220K Figure 2: A sample of spectra as measured under a transverse applied field of 50 G for two temperatures corresponding to the SDP at 20 K and for the RP at 220K. The solid lines depict fits to Eqn.1. magnetic resonance and muon spectroscopies. In fact, as it is now well established, µSR techniques are able to explore a wide range of fluctuation rates, namely 104Hz −1012Hz, whereas magnetic resonance techniques are generally restricted to explore phenomena with slower fluctuation rates within the range 10−2Hz −105H[30]. The reduction of λT F with increasing temperature is thought to arise from the onset of some molecular motions which will become rather fast within the RP state and consequently only a weak dependence of λT F with temperature is found there. As discussed below, the measured rates λT F monitor the development of a time-averaged component that parallels the behavior of the Lamb-M¨ossbauer (or Debye-Waller) factor which can be measured by various other spectroscopies. As shown in Figure 3 , above some 170 K, motional narrowing within the SDP crystal departs from the simple thermally activated behaviour expected for the Debye-Waller factor of a harmonic sample i.e., ln(λ(T)) significantly departs from a straight line. The drastic reduction in relaxation rate witnessed above some 160 K thus monitors the emergence of critical dynamics leading to the phase change. To compare the results with well defined physical quantities we have computed the rigid-lattice value 6 0 50 100 150 200 250 300 0 . 0 0 0 . 0 5 0 . 1 0 0 . 1 5 0 50 100 150 200 250 300 1 0 1 2 1 4 1 6 1 8 T F R e l a x a t i o n r a t e ( µs- 1 ) T ( K ) T F R e l a x . A m p l i t u d e T ( K ) Figure 3: Temperature dependence of the observed relaxation rates (upper frame) and amplitudes (lower frame) resulting from measurements carried out under a transverse applied field of 50 G. for a polycrystalline sample as calculated after an orientational average and reads, M2,poly =4I(I+ 1) 15 ~2γ2 µγ2 N N X j=1 1 r6 j (2) where γNstands for the gyromagnetic ratios of proton nuclei close to the muon site. For such a purpose a crystal structure having an implanted Mu has been evaluated by means of a DFT calculation which employed the crystallographic structure determined from diffraction experiments [9] and thus enabled the optimization of the Mu location within the 2-adamantanone-OMu adduct. The summation was run up to distances reaching 10 Å from the Mu location. The estimation was carried out using the programme developed by Goc [31] which enables the calculation of van Vleck second frequency moment [32] for materials composed of molecules which may have some internal degree of freedom. The result yielded a value of 24.9 G2which is to be compared with that of 22 G2derived from broadband NMR [17] and provides an estimation of the root-mean-squared field spread causing dephasing ∆within the static limit [33] of ≈10 G. Dynamic effects arising from temperature-induced fluctuations of the internal fields, resulting from motions of the molecular skeleton onto which muonium is stuck will lead to a reduction of the second moment through averaging of the local dipolar field due to the modulation of the dipole-dipole interactions [34]. Explicit expressions to carry out such calculations obviously require the specification of some dynamical model to represent the molecular motion and in most cases lead to rather cumbersome formulae the evaluation of which can only be performed by numerical means. Most practical situations however base their interpretation of results concerning the variation of second moment with the correlation time of the underlying motion on the basis of the classic Kubo and Tomita [35] formulation of motional narrowing which, under the assumption of the validity of the fast-fluctuation limit leads to simplified 7 expressions for the relaxation rate which is now given by, λT F =−γ2 µM2τc(3) where τcstands for the correlation time of the fluctuating fields. With values at hand we can now evaluate the order of magnitude of dynamic phenomena being sampled at the experiment. Lets first take the value for M2derived the calculation just referred to above as well as that resulting from NMR measurements [17]. Once plugged into the above equation and consideration of the order of magnitude of the relaxation rates it yields estimates for the correlation time for temperatures below the transition, that is below some 200 K of the order of 10−7s−1−10−6s−1depending upon the value for the second moment taken for the calculation. Above some 200 K the relaxation rate displayed in Fig.3 rapidly decreases with increasing temperature which, assuming the validity of Eq.3 within this range of temperatures, can be seen as a result of a concomitant reduction of the correlation time. The TF relaxation amplitude also displays a sudden jump at basically the same temperature as the relaxation rate does. The microscopic origin of such a steep step cannot be ascertained in full, but the data at hand suggests that well within the RP phase a free muon may combine with a nearly free rotating molecule to give rise to the diamagnetic 2-O- adamantane-Mu+adduct or, as another possibility, a neutral 2-O-adamantane-Mu radical may capture a mobile free electron to give the diamagnetic 2-O-adamantane-Mu−adduct. To wrap up this section it seems pertinent to carry out an order of magnitude comparison of the results just described with those previously obtained by means of dielectric spectroscopy. To such an avail here we take the measured dielectric relaxation characteristic times [11] τdand consider them as time constants of some underlying orientational motions, an assumption here justifiable since we are only pursuing a comparison of the magnitude of transverse relaxation rates. Under such assumptions, closed form expressions for the relaxation rate due to purely dipolar interactions under the assumption of rotation Brownian motion, are given in the literature pertaining to plastic crystals [36] and may be adapted for the case of relaxation from unlike spins with the help of results given in Ref. [37]. The formula used for such a calculation yields, 1 Tcalc 2 =M2×4τd+τd 1+(ωµ−ωp)2τ2 d +3τd 1 + ω2 µτ2 d +6τd 1 + ω2 µτ2 d +6τd 1+(ωµ+ωp)2τ2 d(4) where M2stands for the rigid lattice value of the spectral second moment given by Eq.2, ωµ,p stand for the muon and proton Larmor frequencies for precession within a 50 G field and the τdare assigned to the three different sets of relaxation times reported in [11]. The results tell that the expected relaxation rates for the times corresponding to the main αand βrelaxation processes detected within the SDP phases will be way out of the accessible window for µSR. A point of contact seems however to exist concerning rotational motions within the rotator phase where the calculated relaxation rates are within the interval 0.028µs−1≤1 Tcalc 2 ≤0.010µs−1. The temperature dependence of such quantities is however Arrhenius-like which is in sharp contrast with the weak dependence observed for the relaxation rates above 220 K. 3.4 Measurements under Zero Field As mentioned above, measurements were also carried out under zero applied field and varying temperature. The expectancy was that such measurements would encompass similar phenomena to those explored under TF as well as reveal some details pertaining the paramagnetic species which are to be expected. A detailed account of such measurements is given in the Supplementary Material section. The main findings which are shown in Figures S3 and S4 of such Material show that relaxation of the implanted muons towards thermal equilibrium proceeds following two steps with relatively well separated timescales. A fast process characterized by a rate λf ZF is clearly visible for times below the first microsecond of the experimental histogram. The measured rates are about one order of magnitude larger than those sampled under TF conditions and their temperature dependence shown if Fig. S4 appears to follow the general trends explored under TF, namely a relatively smooth decrease with temperature up to some 170 K followed by an abrupt decrease and a rather soft temperature dependence above 220 K. The results are however qualitative since the time scales of such phenomena are close to the limit of instrumental observation. In contrast, the slowest relaxation process named as λs ZF could be followed 8 in more detail. The data of Fig. S4 nicely follows the same trend explored under TF and its general shape could also accounted following the same steps described to parametrize the data given in the next section. Data concerning the temperature dependence of the fractional amplitude displayed in Fig. S4 of the Supplementary Material section show some significant difference with respect to those plotted in Fig.3 at low temperatures although also shows a steep riseup at the phase transition which is taken as a signature of conversion processes leading to a dominant diamagnetic fraction. 4 Discussion As known since the pioneering work of Adachi et al. [38], rapidly cooling some molecular crystals from their metastable rotator- or plastic-crystal phase can avoid the transition into their fully ordered ground states, driving those samples into some “glassy-crystal state”, where translational order is preserved but the strong rotational disorder characteristic of the plastic-crystal phases is partially quenched. Substituted adamantane compounds such as 1-cyanoadamantane (CN-a, C10H15CN) were found to conform to such a category of materials and dynamic and thermodynamic properties were soon investigated [39–41]. The results showed the persistence of uniaxial reorientational motions well within the “glassy-crystal state”. The material here investigated, which has some distinctive features with respect to CN-a, has shown the persistence of motions within its “glassy-crystal state” down to around 150 K, which are found to be fast enough to be explorable by means of dielectric broadband spectroscopy [11]. The thermal and dynamic characterization of such a metastable phase, as well as the comparison of its properties to those of the true crystal ground state, namely an orthorhombic crystal with Cmc21lattice structure, has recently been achieved [9,13]. Both crystals show a rather small difference in volume, which indicates that the transition into the truly ordered ground state involves a volume contraction, amounting a mere two per cent. In stark contrast, dynamical quantities such as the spectral frequency distribution, as well as the thermal conductivity of both phases, unveil significant differences attributable to substitutional disorder [13]. More specifically, the frequency distribution of the metastable phase shows a population of vibrational low energy states larger than that of the ground state, a difference most noticeable somewhat above 1 THz. Such frequencies which are found to correspond to well defined excitations of mostly librational character are strongly coupled to the acoustic modes [42,43] and become well resolved within the “glassy-crystal” as motions of molecules perpendicular to the molecular dipole axis. The results are also in agreement with previous results obtained from Raman spectroscopy reported by Bistricic et al. [17], who measured a well resolved, narrow band centered at some 1.1 THz, which shifts up to 1.26 THz in going down to the crystal ground state. Structurewise, one expects the SDP phase of the material here under consideration to share some common features with the corresponding glassy phase of the cyano-substituted analogue, which shows strong diffuse scattering in the diffraction patterns attributed to the presence of antiferroelectric domains [44] having a coherence length up to 20 Å. The expectancy of strong antiferroelectric interactions as the relevant processes driving the order-disorder transitions in crystalline adamantanone stems from the rather large value of its molecular dipole moment, 3.4 Debye (1 Debye ≈3.335 ×10−30 C. m), which compares with that of 3.83 Debye of CN-a. In fact, some hints of antiferroelectric behaviour can be inferred from its Kirkwood gkfactor which attains values well below 0.2 within the SDP and shows a sudden jump up to figures above 0.4 above some 200 K [11]. The experimental data for the relaxation rates can be well accounted for using a simple square-root law as previously applied to structural glasses that is [18,19], grounded on the kinetic theory of the glass transition [45,46] and should be valid for temperatures not too far below the glass-transition region, λT F =Ar|T−Tc Tc |+ξ(T), T ≤Tc λT F =ξ(T) = a+b|T−Tc Tc |, T ≥Tc(5) where a, b are parameters to account for the soft temperature dependence of data measured above the transition (i.e. above 220K). In the above formula, Tcstands for some critical temperature, Ais a global scaling constant and ξ(T)comprises a background term set by data measured above the transition. 9