Orientational Glasses : NMR and Electric Susceptibility Studies

We review the results of a wide range of nuclear magnetic resonance (NMR) measurements of the local order parameters and the molecular dynamics of solid ortho-para hydrogen mixtures and solid nitrogen-argon mixtures that form novel molecular orientational glass states at low temperatures. From the NMR measurements, the distribution of the order parameters can be deduced and, in terms of simple models, used to analyze the thermodynamic measurements of the heat capacities of these systems. In addition, studies of the dielectric susceptibilities of the nitrogen-argon mixtures are reviewed in terms of replica symmetry breaking analogous to that observed for spin glass states. It is shown that this wide set of experimental results is consistent with orientation or quadrupolar glass ordering of the orientational degrees of freedom.


Introduction
Understanding the underlying physics and dynamics of glass systems has remained a challenge for contemporary physics despite the large number of experimental and theoretical studies that have shown that a wide variety of glass systems exhibit a number of universal features [1][2][3][4][5][6].Considerable progress was made with the discovery of spin glasses [7][8][9][10][11][12] and the introduction of the concepts of frustration [13][14][15][16][17] and replica symmetry breaking [13,[18][19][20].Beyond the simple spin glasses, which displayed a random orientation of dipole moments, it was recognized very early in the history of spin glasses that molecular systems, and in particular diatomic molecules that have short-range anisotropic electrostatic interactions, were highly frustrated and that they formed a special class of glass systems with geometrical frustration.These new glass systems are characterized by (i) a broad distribution of order parameters provided that sufficient disorder was introduced and (ii) a strong temperature dependence for the molecular dynamics, often resembling a Fulcher-Vogel dependence [21,22].
The interest in the molecular glasses is that the interactions are accurately known, the frustration can be described in clear terms and the local order parameters can be measured directly by nuclear magnetic resonance (NMR) methods.Only recently has there been success in relating the low temperature properties of the molecular glasses to replica symmetry breaking [23,24], and in this paper, we review the experimental NMR methods that have been used to measure the order parameters of these systems and relate the observations to electric susceptibility studies and the models of replica symmetry.

Orientational and Quadrupolar Glasses: Basics
While spin glasses arise from the combination of frustration (of dipolar magnetic interactions) and disorder, the molecular orientational and quadrupolar glasses arise from the frustration of short-range interactions that are dominated by intermolecular electric quadrupole-quadrupole interactions resulting from the non-spherical electrostatic charge distribution of the molecules.The classic examples are ortho-para hydrogen mixtures [21,[25][26][27] and solid N 2 -Ar alloys [28,29].The order parameters are given by two sets of parameters: (i) the local axes (x,y,z) for the mean orientation of the molecular axes and (ii) the degree of alignment or quadrupolarization about those axes.The latter in the classical case is given by the expectation values of the spherical harmonics, σ class = Y 20 (θ) and η class = Y 22 (θ, φ) .For solid hydrogen, the orbital angular momentum J is a good quantum number.There are two molecular species: ortho-H 2 (with orbital angular momentum J = 1 and total nuclear spin I = 1) and para-H 2 (with J = 0 and I = 0).Only the ortho-H 2 molecules have an electrostatic quadrupole moment and are subject to orientational or quadrupole ordering.The para-H 2 molecules provide the disorder.Although the para-H 2 molecules are the true ground state, the ortho-para conversion is a very slow process in pure solid hydrogen, and one can conduct experiments for a wide range of ortho-para mixtures (for a review of the properties of sold hydrogen, see Silvera [30]).For the ortho-H 2 molecules, the order parameters are are given by σ qu = 3J 2 x − J 2 and η qu = J 2 x − J 2 y , which are the expectation values of the tensor operator equivalents of the spherical harmonics to within a simple numerical factor.Above a critical concentration of 55% ortho-H 2 , the ortho molecules order in a four-sublattice anti-ferro-orientational Pa3 structure [30].Below 55%, the ortho molecules form a quadrupolar glass with no sudden phase transition [31].The onset of glass formation was signaled by the observation of a broad distribution of local order parameters at low temperatures.For a simple orientational glass, only the axes (x,y,z) are randomly distributed, while for the quadrupolar glass, both the axes and the quadrupolarizations vary randomly.The distinction between the two types of molecular glasses is illustrated in Figure 1 [31].In addition to the static observations of the order parameters, experiments also showed strong temperature dependencies for the dynamics of the molecules (principally their thermal fluctuations about the mean molecular alignments) as inferred from various studies of the nuclear spin relaxation times.Two independent measurements of the relaxation times in solid ortho-para H 2 mixtures are shown in Figure 2. The change in the static order parameter, σ (as defined above and represented by the variable S ≈ (1 − |σ|) in Figure 2b), is considerably less dramatic than the changes observed in the dynamical parameters that determine the nuclear spin-lattice relaxation times.The main difference in the two measurements shown in Figure 2 is the frequency of the measurements; 22 MHz for Figure 2a (Ishimoto et al. [25]) and 33 MHz for Figure 2b (Cochran et al. [26]).This difference was the first evidence of the characteristic frequency dependence of the response functions observed below the glass formation and will be discussed further in the next section.Studies of solid nitrogen-argon mixtures [33][34][35] showed similar features to those referenced above for solid ortho-para H 2 mixtures.The randomness for the N 2 -Ar solid mixtures is introduced by the replacement of nitrogen molecules with argon atoms, which are spherically symmetrical and carry no electrostatic quadrupole moment.At high nitrogen concentrations (more than 77%), long-range Pa3 ordering is observed similar to pure ortho-H 2 .Below that critical concentration, glass behavior is observed at low temperatures.As indicated in Figure 3a,b, a lattice change is observed when the long-rage ordering occurs with the hcp lattice transitioning to an fcc lattice as the orientational ordering occurs.The critical concentration depends on the concentration dependence of the anisotropic interactions.A detailed restricted trace calculation [36] showed that the orientational fluctuations for ortho-para hydrogen mixtures varied as (2x − 1) where x is the ortho concentration leading to a predicted critical concentration of 50%.No detailed calculation has been made for solid N 2 -Ar mixtures.For the latter, it is noteworthy that for N 2 concentrations less than 56%, a cubic lattice is maintained at low temperatures.This is important for any analysis of the glass formation because in the hcp structure, the symmetry of the electric quadrupole-quadrupole interaction is incompatible with the non-cubic lattice symmetry.As a consequence, replacing a quadrupole by a spherical diluent is equivalent to adding a local conjugate quadrupolar field so that for the hcp lattice, the glass formation is analogous to that of a spin glass in a static magnetic field, and any phase transition would be rounded out.The solid N 2 -Ar mixtures at low concentrations (Phase II of Figure 3b) therefore provide a better testbed for understanding quadrupolar or orientational glass formation.

NMR Methods
NMR methods are sensitive to long-range ordering or orientational or quadrupolar glass ordering because the nuclear spin-spin interactions depend on the local order parameters.Both the nuclear dipole-dipole interactions and the anisotropic chemical shift interaction vary with the values of the local order parameter and thus alter the NMR spectra.
The intra-molecular dipole-dipole interaction between the two nuclei of one molecule can be written as [27]: where we have introduced the irreducible orthonormal tensor operators T 2m and N 2m to describe the orientational and nuclear spin degrees of freedom, respectively.Under rotations, these operators transform equivalently to the spherical harmonics Y 2m (θφ).In the manifold J = 1, the orientational operators are [27,39]: and The strength of the intra-molecular dipole-dipole interaction for H 2 is given by D/2π = 173.06kHz.For a high magnetic field, the dipole interaction is only a weak perturbation of the nuclear Zeeman interaction, H z = γhI z B z , where γ is the nuclear gyromagnetic ratio and B z is the magnetic field aligned along the z-axis.For the Zeeman term alone, the energy levels for one ortho-H 2 molecule with total nuclear spin I= 1 consist of three equally-spaced levels corresponding to I z = 1, 0, −1 with energy separations ∆E 01 = ∆E −10 = hω L where ω L = γB z is the nuclear Larmor frequency.In the absence of any orientational ordering, the dipolar interaction averages to zero, and this energy separation is unperturbed.However, for a non-trivial orientational ordering, the expectation value of the secular component of H D (the part that commutes with the Zeeman interaction), H D0 , is non-zero, and the energy levels are perturbed as shown in Figure 4.
If the orientation of the magnetic field with respect to local axes of molecule i is defined by the polar angles (θ i , φ i ), we have for axial symmetry (J 2 x − J 2 y ) = 0 (which is usually assumed in the literature) [21,27], where the local order parameter From Equation (4) and Figure 4, we see that each molecule with a non-zero order parameter σ i contributes a doublet to the NMR absorption line shape at: The experiments are almost always carried out for powdered samples, and we need to sum the contributions to the spectrum over a powder distribution P(Ω) of the crystalline orientations Ω(θ, φ).If I(∆ω, σ) is the line intensity at frequency ∆ω for a fixed σ, we have The line shapes therefore consist of a doublet given by: These are the familiar Pake doublets and are shown in Figure 5 for σ = 1.In order to calculate the line shape expected for a quadrupolar glass, we need to consider an approximate form for the distribution of local order parameters.If P(σ) is the probability of finding the value σ in the glass state, then we find P(σ) from the distribution of the values of local entropies.For each molecule with local order parameter σ, the contribution to the entropy is: For small σ, S(σ) is linear in σ, and assuming the distribution of entropies, P (s) is a constant, we have as a first approximation: Using this approximation and the expressions for the line shapes for the doublets given by Equation ( 7), the glass line shape is expected to be given by: The calculated glass line shape using Equation ( 10) is shown in Figure 6.This shape is to be compared with the shape observed at low temperatures.An example is shown in Figure 7 for an ortho-H 2 concentration of x = 0.23 and a temperature of 85 mK.The general form is in good agreement with the predictions.
Frequency (∆ω)  Interestingly, a very small narrow peak is observed at the center that is the contribution from molecules with zero order parameter.Reproduced with permission from [21].Copyright American Physical Society, 1978.
At lower temperatures and for high magnetic fields, the line shape shows a distinct asymmetry as shown in Figure 8.This asymmetry is due to the finite nuclear spin polarization (16% in the case of Figure 8).For non-trivial nuclear polarizations, the ratio of the amplitudes of the two branches of the Pake doublet is not unity.It is given by: where ω L is the nuclear Larmor frequency and T is the absolute temperature of the sample, obtained by a polynomial fit to the order parameter distribution function using P (σ) = ∑ n a N σ n to find: using The order parameter distribution estimated from the best fits to the NMR line shapes is shown in Figure 9 for ortho-para H 2 samples with 23% ortho hydrogen [27].There is a smooth variation with temperature with no discontinuity or other evidence of a sudden transition.The relevant parameter for a spin glass transition, however, is the Edwards-Anderson [40] order parameter given by: where s i is the component of the spin at site i, is a thermodynamic average and [ ] c is an average over configurations c.While the average [ s i T ] c vanishes for paramagnetism and for a spin glass, q EA (T) is non-zero for the spin glass.For the quadrupolar glass, we therefore consider the glass order parameter: The temperature dependence of the quadrupolar glass order parameter q Quad EA is calculated from the order parameter distribution functions P (σ) determined from the NMR line shapes shown in Figure 10 [41].Although there is no evidence for an abrupt transition in the temperature variation of q Quad EA (T), the observed temperature dependence is much stronger than that predicted for a non-collective simple random field model [42][43][44] shown by the broken dashed line of Figure 10.The observed variation is in good qualitative agreement with the model for a quadrupolar glass proposed by Lutchinskaia et al. [45] consistent with a collective quadrupolar glass ordering at low temperatures.Temperature dependence of the quadrupolar glass order parameter q EA (T).Reproduced with permission from [41].Copyright Taylor & Francis, 1986.
Classical orientational glass ordering can be observed in solid N 2 -Ar mixtures using NMR techniques at low temperatures.Instead of the common isotope 14 N, experimenters use 15 N, which has nuclear spin 1/2, and the large (∼3 MHz) nuclear quadrupole interaction for 14 N vanishes.For the diatomic molecule 15 N 2 , one can treat the total molecular nuclear spin as I = 1 and use a similar analysis of the NMR spectra as was used for diatomic H 2 .Because of the larger mass and large moment of inertia compared to H 2 molecules, the angular momentum is not a good quantum number, and the local order parameters are purely classical and given by: where θ i is the angle between the instantaneous molecular orientation and the mean orientation.
The NMR line shapes for local orientational ordering in solid N 2 -Ar mixtures have one major difference compared to those observed for solid hydrogen because there is an extra term in the nuclear spin Hamiltonian that depends on the orientational order parameter.That term is the anisotropic chemical shift resulting from the interaction between the nuclear spin and the magnetic fields of the molecular electron distribution given by: K = (4.0 10 −4 )ω L .The 15 N 2 NMR line shapes in the presence of a finite orientational order parameter form doublets for a powered sample, but unlike those for solid H 2 , they are antisymmetric with Equation ( 6) replaced by: as illustrated in Figure 11.A typical line shape is shown in Figure 12 and compared with a calculated line shape for a fixed order parameter σ = 0.86 that includes a small broadening, as shown by the inset to Figure 12.The agreement is very satisfactory.Narrow lines consistent with long-range ordering in a Pa3 structure are only observed for N 2 concentrations above 76%, where the lattice structure remains fcc (Part III of Figure 3b).In the region where the lattice remains hcp for 56% < X N 2 < 77% (Part II of Figure 3b), only broad NMR line shapes characteristic of quadrupolar glass-like ordering are observed.The remarkable observation is that in Region IB of the phase diagram (see Figure 3b), one once again observes a relatively narrow line shape (see Figure 13) for x = 50% [38].Because of the presence of high substitutional disorder, this narrow line shape is believed to be a true orientational glass (analogous to that depicted in Figure 1a) with a broad distribution of alignment angles, but a relatively narrow distribution of order parameters about those angles and with a mean σ = 0.74.On further dilution of X N 2 to 40%, only very broad NMR line shapes are observed in Region IA of Figure 3b.While the NMR line shapes of solid ortho-para hydrogen mixtures and N 2 -Ar mixtures are consistent with quadrupolar and orientational glass states at low temperature and the growth of the quadrupolar order parameter as the temperature is lowered suggests that the local ordering is collective in nature and not driven by fixed local fields, it is important that the dynamics be investigated to determine whether the characteristic rapid freezing of the orientational fluctuations occurs at the onset of glassiness.This is discussed in the next section.

Dynamics
NMR methods can be used to investigate the dynamics of molecular fluctuations and in particular the orientational or quadrupolar fluctuations over a wide frequency range [25,26,37].The sensitivity of NMR methods to motions arises because fluctuations of the molecular alignments or local order parameters render the intra-molecular nuclear dipole-dipole interactions time dependent [46].Specifically, the nuclear spin-lattice relaxation time, T −1 1D , due to the modulation of the dipole-dipole interactions, is given by: D measures the strength of the intra-molecular dipole-dipole interaction, and ω 0 is the nuclear Larmor frequency.J(ω) is the spectral weight of the fluctuations at frequency ω.
The T lm are the tensorial operators defined in Equation (2).
The expectation value T 2m (t)T 2m (0) T must be calculated with a quantization axis parallel to the magnetic field.The order parameters, however, are calculated with respect to the local molecular symmetry axis.In order to carryout this evaluation, we therefore consider the rotation: The d mµ are the rotation matrix elements, and the angles (α, β) are polar angles defining the orientation of the magnetic field in a local reference frame aligned along the molecular equilibrium orientation.One then assumes [47] that the time dependence can be separated from the order parameter variations as: The original data shown in Figure 14 [48] were interpreted in terms a single correlation time τ Q .The result, a strong almost exponential temperature dependence, resembles, but is not identical to a Fulcher-Vogel temperature dependence.The data in Figure 14 along with the results observed for the very slow dynamics determined by NMR stimulated echo studies [49] are among the most convincing experimental evidence for traditional glassy dynamics for the quadrupolar glass state of solid ortho-para hydrogen mixtures.Further studies, however, showed that the nuclear spin relaxation was not a simple exponential, but rather a sum of different exponentials associated with different parts of the NMR line shape, as illustrated in Figure 15.This spectral inhomogeneity of the nuclear spin-lattice relaxation can be expected from Equation (22).If one assumes that the relaxation is dominated by the fluctuations of the axial order parameter, σ, then Equation ( 22) leads to: for the contribution to the relaxation by the dipole terms.The factor (2 − σ − σ 2 ) is the mean square variation of the operator T 20 .The angular terms are associated by different isochromats of the NMR line shape, and if one assumes Lorentzian forms with correlation time τ Q for the spectral densities g(mω 0 ), we have [47]: ∆ν is the frequency of an isochromat with respect to the center of the NMR absorption line.
The spin-rotational (SR) interaction also contributes to the relaxation rate: where C = 100 kHz for a hydrogen molecule.These two contributions to the relaxation rate from the intra-molecular dipole-dipole and spin-rotation interactions are additive, and the overall relaxation given by T −1 1SR is shown in Figure 16.The agreement with experiment is satisfactory.
Yu et al. [50] also noted that cross-relaxation, T −1 12 , between different isochromats will tend to smooth out the inhomogeneities in the relaxation rates, and the observed relaxation will be determined by The resulting variation of the relaxation as a function of ∆ν is shown in Figure 16.

Thermodynamic Measurements
Measurements of the heat capacities of solid ortho-para hydrogen [51,52] and N 2 -Ar mixtures [53,54] showed rather featureless temperature variations and were at first sight disappointing, but studies of the very slow dynamics of solid hydrogen with dilute ortho-H 2 concentrations were able to provide a uniform picture for the very slow motions and the thermodynamic behavior.The low frequency dynamics was investigated using NMR stimulated echo techniques employing three RF pulses [49].After a first pulse, the nuclear magnetization is rotated to a transverse axis (x-axis), where it evolves under the intra-molecular dipole interactions, which depend on the local order parameter.After a time τ, a second pulse stores the transverse magnetization along the original z-axis.The value stored depends on the original value of the order parameter.After a waiting time t w , a third pulse returns the stored magnetization back to the x-axis, and after a further time τ, a stimulated echo is formed, provided that the order parameter varies only very slowly in the waiting time t w .A variety of different pulse sequences can be considered but the simplest (π/2 − τ − π/2 − t w − π/2) about the same axes produces an echo signal: ..... designates a configurational average over sites i.The dependence on t w allows one to explore slow temporal fluctuations of the molecular order parameters over time τ.The results (Figure 17) were surprising, as they showed a unique logarithmic decay.The logarithmic decay can be understood in terms of the droplet model proposed by Fisher and Huse [55,56].Fisher and Huse considered the density of low energy states for connected clusters and argued that because the clusters are rare, they could be treated as "non-interacting two level systems".For a free energy barrier E B , the tunneling rate at temperature T is: where Γ 0 is an attempt frequency.For long times, Γ 0 is well defined because it is associated with a characteristic cluster size.The relaxation is hierarchical as low-energy barriers must be overcome before larger barriers can be crossed.The dynamics of of the faster degrees of freedom therefore constrain the slow motions.As a result, in time t, the barriers overcome will be those with energy less than E max (t) where: For all energy barriers traversed in time t, the order parameters change significantly, and the correlation functions cos[Dσ i (0)τ]cos[Dσ i (t)τ] will vanish after the barriers are crossed.If the density of low energy excitations is P(E B ), the stimulated echo amplitude will be given by: At low temperatures, we only need consider the density of states near E ≈ 0, and the echo amplitude becomes: t 0 is determined by the attempt frequency Γ 0 .Since the order parameter changes by the order of unity in time Γ −1 0 : The decay of the stimulated echos is therefore directly related to the density of states at low energy, P(0), which is expected to be constant and can be related to the thermodynamic measurements.From Figure 17, we find P(0) = 0.59 ± 0.05 K −1 .
Knowledge of the density of energy levels P(E) allows one to calculate the heat capacity and test the interpretation of the dynamics against the experimental results for heat capacity studies [51][52][53][54].Using the energy scheme of Figure 18, the heat capacity is calculated from the sum over the distribution of energy states S(∆) as: with z = 3∆/k B T. X is the ortho-H 2 concentration.The result of the calculation shown in Figure 19 gives a good description of the characteristic linear temperature dependence at low temperatures, and an excellent fit is obtained using small deviations from a constant density of energy states shown by the solid line of Figure 19.The diamonds represent the data of Ricketson [51] and the squares that of Haase et al. [54].Reproduced with permission from [39].Copyright Elsevier, 1969.

Dielectric Measurements
While the analysis of the NMR line shapes of solid ortho-para H 2 and solid N 2 -Ar mixtures with high disorder (that is, low concentrations of the active quadrupole bearing molecule (ortho-H 2 or N 2 )) show clear evidence of the formation of glassy states at low temperatures and the temperature dependence of the dynamics cannot be understood except in terms of local collective ordering, the characteristic hysteresis of the susceptibility (the parameter conjugate to the local order parameters) remained elusive until the studies of Pilla et al. [57,58] using electric susceptibility measurements.These measurements were carried out on dilute solid N 2 -Ar mixtures for concentrations where the lattice remained cubic.In this case, the concern that there exist local electric field gradients due to the introduction of substitutional disorder is greatly lessened, and the system is much closer to the class of the familiar spin glasses.
The use of a high sensitivity capacitance bridge [59] to measure dielectric susceptibilities showed distinct hystereses in the observed values of the dielectric susceptibility on thermal cycling [57].The break in the observed values of at 9.8 ± 0.1 K (Figure 20) clearly marks a well-defined dynamical transition for solid N 2 -Ar mixtures with 51% N 2 that was absent in the NMR measurements.The really significant consequence of the dielectric susceptibility measurements for the molecular glasses is that they can be used to test the applicability of replica symmetry breaking theories to this class of glass states [14,60].For simple replica symmetry breaking, Cugliandolo and Kurchan [18] developed a generalized fluctuation dissipation theorem for the generalized susceptibility: where q i is a generalized coordinate.The essential result is that below the characteristic glass freezing temperature of the dynamics, T dyn , the change in the response function depends on the waiting time t w because of the changes in the landscape of free energy barriers that occur in the time t w .The response function is defined by: R i (t, t w ) = 1 k B T χ(C)dC (34) where χ(C) is the relevant susceptibility.Pilla et al. [23] showed that as a consequence of this generalized fluctuation dissipation theorem, χ ∝ T/T dyn for large excursions compared to χ ∼ 1 (the classical result) for small δT.The results are shown in Figure 21.

Discussion
The experimental results for the NMR studies of the line shapes of solid ortho-para H 2 and solid N 2 -Ar mixtures using continuous wave techniques show that the quadrupole-bearing molecules in these mixtures have a broad distribution P(σ) of local order parameters σ at low temperatures.These parameters consist of two kinds: (i) quadrupolar, referring to the deviation from spherical symmetry of the average molecular alignments; and (ii) orientational, referring only to the orientation of the average alignment of the molecular axes.Ortho-para H 2 mixtures for ortho concentration less than 55% and N 2 -Ar mixtures for 56% < X(N 2 ) < 77% form quadrupolar glasses with both the quadrupolar and orientational parameters broadly distributed.Solid N 2 -Ar mixtures for 40% < X(N 2 ) < 56% form orientational glasses for which the local quadrupolar parameters are constant (but not necessarily zero), and only the orientational axes vary.The equivalent of the Edwards-Anderson order parameter, q(T) = σ 2 i , shows a rapid variation with temperature (stronger than non-collective para-orientational models would predict), but with no sharp transition for studies down to the lowest temperatures (35 mK for solid ortho-para H 2 mixtures and 1.5 K for solid N 2 -Ar mixtures).
The dynamics associated with the freezing of the molecular orientational and quadrupolar fluctuations in the quadrupolar glass states has been determined using pulsed NMR techniques, either using standard π 2 − τ − π 2 pulse sequences or π 2 − τ − π 2 − t wait − π 2 stimulated echo sequences.The latter is particularly useful for providing information about the very slow dynamics in the glass states.If the relaxation is interpreted in terms of a characteristic fluctuation time τ Q (which does give a good description of the data), τ Q increases rapidly with decreasing temperature below a temperature T dyn characteristic of a given ortho-H 2 or N 2 concentration.In addition, the unique logarithmic decay of the stimulated echos provides a measure of the density of energy states P(E) near E = 0.This density of states can be used to calculate the heat capacity C V as a function of temperature.Good agreement is found especially for the magnitude of the characteristic linear dependence on temperature for T < T dyn .
Finally, studies of the dielectric susceptibility that depends on the orientational parameters for the molecules have been made at very low audio frequencies, and the experiments show unmistakable hystereses on thermal cycling.Pilla and colleagues [57] have been able to show that this is consistent with the replica symmetry breaking models of Parisi [13], and these results put these disordered quadrupolar molecular solids in the same class as other frustrated systems that exhibit spin-glass-like behavior at low temperatures.Further studies of the dynamics for very low quadrupolar concentrations in solid N 2 -Ar mixtures as one varies the concentration through the lattice percolation concentration would provide a valuable test of the underlying physics since the belief is that it is the frustration of the short-range interactions that lie at the origin of the formation of these glass states.

Conclusions
NMR studies of the line shapes of solid ortho-para-H 2 and solid N 2 -Ar mixtures showed that at low temperatures the quadrupole bearing molecules (ortho-H 2 and N 2 ) were characterized by a broad range of local order parameters.This characterization lead to the interpretation of the low temperature behavior in terms of a quadrupolar glass analogous to a spin glass.The transition to the glass state as a function of temeprature was, however, smooth with no sudden phase transition.The temperature evolution of the Edwards-Anderson order parameter was nevertheless much stronger than that predicted for non-collective effects.Measurements of the nuclear spin relaxation times provided more insight into the local dynamics and implied a rapid freezing of the dynamical degrees of freedom.A deeper understanding of the underlying physics was provided by analysis of the decay of stimulated NMR echoes designed to explore ultra-slow motions that characterize the thermal behavior of spin glass materials as the system is governed by the crossing of energy barriers in a complex configuration landscape.This analysis provided an unambiguous measure of the density of low energy excitations and lead to a quantitative fit for the results of heat capacity measurements.Finally, the temperature of the quadrupolar susceptibility deduced from measurements of the dielectric susceptibility revealed hysteresis effects that are consistent with the replica symmetry breaking models of generalized spin glass theories.

Figure 2 .
Figure 2. Observed variation of the nuclear spin-lattice relaxation times of solid ortho-para hydrogen mixtures showing the sharp change in the temperature dependence attributed to the onset of quadrupolar glass formation.Figure 3b also shows the variation of the echo amplitudes S, indicating changes in the static order parameters below 0.55 K.In contrast to the variation of T 1 , S evolves smoothly with temperature.(a) Temperature dependence of relaxation times for 50% ortho H 2 in solid ortho-para H 2 mixtures.Reproduced with permission from Figure 6 of Ishimoto et al. [25].(b) Temperature dependence of relaxation times for 33% (diamonds), 40% (pluses) and 50% (crosses) ortho H 2 .Reproduced with permission from [26].Copyright American Physical Society, 1980.

Figure 4 .
Figure 4. Nuclear Zeeman energy levels for one ortho-H 2 molecule calculated for total molecular nuclear spin I = 1 and angular momentum J = 1.Each molecule contributes a doublet with absorption frequency ω = ω L ± Dσ i P 2 (cosθ i ).

Figure 5 .
Figure 5. Calculated Pake doublet line shape for the long-range ordered Pa3 structure of solid H 2 for a fixed value of the order parameter σ and a distribution of angles α for the the axes of a powdered sample.

Figure 6 .
Figure6.Calculated NMR line shape for quadrupolar glass state of solid H 2 for a broad distribution of order parameters σ.Reproduced with permission from[27].Copyright Springer Nature, 1998.

Figure 7 .
Figure 7. Observed NMR line shape for the quadrupolar glass state of solid H 2 for ortho-H 2 concentration x = 0.23 at a temperature of 85 mK.The experiments measured the derivative line shape using traditional lock-in methods, and the symmetrical dome shaped line is obtained by integration.Interestingly, a very small narrow peak is observed at the center that is the contribution from molecules with zero order parameter.Reproduced with permission from[21].Copyright American Physical Society, 1978.

Figure 8 .
Figure 8. Observed NMR line shape for quadrupolar glass state of solid H 2 for ortho-H 2 concentration x = 0.23 at a temperature of 45 mK in a magnetic field of 5 T, showing asymmetry due to a finite nuclear spin polarization.Reproduced with permission from [27].Copyright Springer Nature, 1998.

Figure 9 .
Figure 9. Temperature dependence of the probability distribution P (σ) for order parameters σ as deduced from observed NMR line shapes for solid hydrogen with ortho-H 2 concentrations of 23%.Reproduced with permission from[27].Copyright Springer Nature, 1998.

Figure 11 .
Figure 11.Calculated NMR absorption line shape for long-range Pa3 ordering for pure solid 15 N 2 .The line consists of two asymmetric Pake doublets.

Figure 12 .
Figure 12.Comparison of observed NMR absorption line shape at 4.2 K with the calculated NMR absorption line shape for long-range Pa3 ordering for pure solid 15 N 2 .The inset shows the order parameter distribution P (σ) that was used.Reproduced with permission from[38].Copyright Springer Nature, 1998.

Figure 13 .
Figure 13.Variation of observed low temperature NMR line shapes with nitrogen concentration in solid N 2 -Ar mixtures, showing inferred distributions of order parameters P (σ).Reproduced with permission from [38].Copyright Springer Nature, 1998.

Figure 15 .
Figure 15.Observed nuclear spin-lattice relaxation in the form of multi-exponentials in the glass state of solid ortho-para hydrogen mixtures.Reproduced with permission from [47].Copyright Springer Nature, 1986.

Figure 16 .
Figure 16.Variation of the nuclear spin-lattice relaxation with frequency for a glass ordering in solid ortho-para hydrogen mixtures.Reproduced with permission from [47].Copyright Springer Nature, 1986.

Figure 17 .
Figure 17.Decay of NMR stimulated echos for the glass state of ortho-para hydrogen mixtures.Reproduced with permission from [39].Copyright Elsevier, 1969.

Figure 18 .
Figure 18.Sketch of the energy levels of an individual ortho-H 2 molecule with angular momentum J = 1 with the degenerate states J z = ±1 separated from the state J Z = 0 by an energy gap 3∆.Reproduced with permission from [39].Copyright Elsevier, 1969.

Figure 19 .
Figure19.Heat capacity of the glass state of solid ortho-para hydrogen mixtures calculated from density of energy levels S(∆).The diamonds represent the data of Ricketson[51] and the squares that of Haase et al.[54].Reproduced with permission from[39].Copyright Elsevier, 1969.

Figure 21 .
Figure 21.Comparison of the hysteresis observed for the frustrated sold N 2 -Ar mixtures and solid O 2 as a function of the temperature swing of the hysteresis.H max is the maximum variation of the hysteresis.Reproduced with permission from [23].Copyright Institute of Physics, 2009.