Spin-glass dynamics and the barrier model: extraction of the Parisi physical order parameter

The fact that the model of barrier heights monotonically increases with increasing metastable state Hamming distance (the difference between the spin orientations of any two states) on a hierarchical (ultrametric) tree can account for the ageing phenomenon in spin glasses. We make use of this model to extract spin-glass dynamics, including the instantaneous value of the zero-field-cooled magnetization M-ZFC, and the time dependences of the thermoremanent magnetization M-TRM(t(w),t) or, equivalently, of the zero-field magnetization M-ZFC(t(w), t), where t is the measurement time after waiting a time t(w) below the spin-glass transition temperature T-g. The pure states introduced by Parisi are thought to originate from metastable states because finite barriers diverge as the temperature is lowered. These infinite barriers encompass metastable states separated by a self-similar distribution of barrier heights, independent of temperature. The model for extraction of spin-glass dynamics which we have used relies on the assumption that the states separated by barrier heights O less than or equal to Delta less than or equal to E-Z empty instantaneously into the new ground states upon a corresponding change in magnetic field energy E-Z associated with a change in magnetic field H. Previous magnetic field cycling experiments have shown that the exchange of occupation 'respects' the barrier heights; the new instantaneous ground-slate occupation is bounded by the same set of barrier heights, O less than or equal to Delta less than or equal to E-Z, as the original state occupation. Subsequent diffusion within the initial manifold of states takes place over barriers with reduced magnitudes Delta - E-Z, and to a sink arising from the change in Zeeman energy overcoming the barriers Delta less than or equal to E-Z. This model can account for M-TRM(t(w), t) if the initial manifold is identified as the field-cooled state (with magnetization M-FC), and the ground state upon H --> O as the zero-magnetization state. We have solved the diffusion problem on the finite hierarchical tree in the presence of a sink at Delta less than or equal to E-Z enabling us to calculate M-ZFC(t(w), t). Our results are representative of experiment. We also show that we can account quantitatively for thermal cycling results by use of the measured temperature dependence of the barrier heights. We address the infinite network by relating M-TRM(t(w), t) to the Parisi physical order parameter P(D) given by P(D) = omega(t(w))[d/dD(E-Z)][M-TRM(t(w), t)/M-FC], where omega(t(w)) is the number of explored states during t(w), and D(E-Z) is the Hamming distance appropriate to a barrier height Delta(D) = E-Z. We have measured M-TRM(t(w), t)/M-FC over a wide range of temperatures and a very fine mesh of magnetic fields. Extracting P(D), we find agreement with the form of the mean-field results at different temperatures T. We relate the instantaneous (stationary) value of M-ZFC as H --> O to the strength of P(D) at D = O: a delta function with area 1-x(T) in the Parisi mean-field solution. The temperature dependence of 1-(x) over bar(T) generates the shape of M-ZFC(t(w),t --> O) against temperature. The mean-field solution leads to [M-ZFC(t(w), t --> O)/M-FC]\(T --> O) = 1/2 The ratios at T --> O found experimentally for a metallic spin glass Ag-6 at.% Mn and an insulating spin glass CdCr1.7In0.3S4 approach 0.4 and 0.2 or less respectively. This implies that, as 1-(x) over bar(T --> O) diminishes more from the mean-field result of 1/2, the more discrete the probability density for the exchange interactions.