Out-of-Equilibrium Dissipative ac — Susceptibility in Quantum Ising Spin Glass

The imaginary part of the non-equilibrium magnetic susceptibility of Ising spin glass in a transverse field under timedependent longitudinal external magnetic field has been calculated at very low temperature on the basis of quantum droplet model and quantum linear response theory. Quantum and aging effects on the low temperature dynamics of the model are discussed. A comparison with recent theoretical and experimental data in spin glass is made.


Introduction
Over the last two decades there is a great deal of research on the experimental and theoretical description of disordered magnetic materials.The understanding of the interplay between disorder, quantum and thermal fluctuations remains among the most relevant problem of condensed matter physics.In strongly disordered systems the dynamics becomes very slow which is characteristic of glassy state and aging scenario [1][2][3][4][5][6][7][8].Aging phenomena and non-equilibrium slow dynamics have been investigated during last years in many materials with glassy properties such as spin glasses [1][2][3][4][5][6][7][8], polymer glasses [9,10], gels [11] and other areas like neural networks, information processing, optimization problems [12].Despite a great progress towards the understanding of nonequilibrium dynamics, some problems remain open.One of them is an investigation of a very low temperature nonequilibrium dynamics in quantum spin glasses, namely the nature of quantum channels of relaxation, the behavior of quantum glassy system subjected to periodic driving force, aging at very low temperatures.The natural basis for the interpretation of aging is based on coarsening ideas of a slow domain growth of a spin-glass type ordered phase [8,13,14].A large attention in the last decade was paid to the spin glasses representing a model systems for study of non-equilibrium dynamics providing a measure of processes causing the aging: the magnetic susceptibility [15][16][17][18][19][20][21][22].
In classical spin glasses in the ac susceptibility meas-urements the magnetic response of the system to a small oscillating magnetic field applied after quenching exhibits aging effects.This response depends on its thermal history and the time interval the system has been kept at a constant temperature in the glass phase [1][2][3][4][5][6][7][8][15][16][17].
It is assumed that isothermal aging of a d-dimensional spin glass is a coarsening process of domain walls, and the temporal ac susceptibility (real part   andimaginary part   ) at a given frequency of ac magnetic field ω at time t after the quenching scales as [23][24][25]   is size of the droplet being polarized by oscil- lating field and   R t is the typical domain size, L and R may scale according logarithmic growth law or algebraic one [23][24][25].
t , and 0 is a certainmicroscopic characteristic time [26].The logarithmic growth law (like the algebraic one) is supported by recent experiments [23][24][25].The expressions (1) and ( 2 In the ac susceptibility measurements, the ac field excites droplets of length scales L up to Because in aging experiments the time spent after quench is t   .These droplets have walls which partly coincide with walls of the domain of size R .
In this paper we investigate the real time non-equilibrium dynamics in d-dimensional Ising spin glass in a transverse field in terms of droplet model at very low temperatures [26].We calculate the dissipative component of the ac susceptibility as a function of the time elapsedsince a thermal quench and frequency of driven field.

Model Hamiltonian
The droplet model describing the low-dimensional shortrange Ising spin glass is based on renormalization group arguments [26,28].In dimensions above the lower critical dimension l (usually in spin glass l ) the droplet model finds a low temperature spin-glass phase in zero magnetic field.This phase differs essentially from the spin-glass phase in the mean-field approximation of the Sherrington-Kirkpatrick infinite-range spin-glass model [13].In the droplet model there are only two pure thermodynamical states related to each other by a global spin flip.In magnetic field there is no phase transition.A droplet is an excited cluster in an ordered state where all the spins are inverted.The natural scaling ansatz for droplet free energy ;  is the correlation length, L is the length scale of droplet and  is the zero temperature thermal exponent.The droplet excitations have a broad distribution of their free energies at scale L for large L in a scaling form [4] It is assumed that is a generalized temperature dependent stiffness modulus which is of order of character-  [26,27].
In this paper we use a phenomenological quantum droplet model of spin glass theory [26][27][28][29] (which does not use the mean-field approximation) in order to describe the non equilibrium behavior of the magnetic dynamical susceptibility at very low (but finite) temperatures .
We consider the following model hamiltonian of ddimensional Ising spin glass in a transverse field, x i z i  and where  are the Pauli matrices,  is the strength of the transverse field and the nearest neighbor interactions ij J are independent random variables of mean zero and the sum in Equation ( 4) is performed over nearest neighbors.
The quantum spin glass transition in a dilute dipole coupled magnet is described by this model hamiltonian [4,13,20].It is supposed that this model may also represent, for example, the physics of deuteron glass such as and mixed betaine phosphate-phosphite [28].The transverse field in [28] is interpreted as the frequency of the proton tunneling.

Rb ND D PO
In ref. [26,27], M. J. Thill and D. A. Huse have shown that for enough low T the quantum hamiltonian Equation (4) can be represented as independent quantum two-level systems (low energy droplets) with the hamiltonian,   where  is a microscopic tunneling rate and, finally, we assume L  is the same for all droplets L [26,27]. of scale In the quantum droplet model of Thill and Huse [26,27] the relative reduction of the Edwards-Anderson order parameter from its zero temperature value Here  is the stiffness modulus, and is the classical-to-quantum crossover length scale defined by .For droplets with (quantum regime [26][27][28]), the excitation energy is always greater than B and thermal fluctuations are irrelevant.These droplets behave quantum mechanically.The larger droplets and behave classically.In quantum regime the length growth is due to quantum fluctuations connected with droplet are quantum-mechanically active for is proportional to T L  .

Dynamic Nonequilibrium Magnetic Susceptibility
We consider the time dependent Hamiltonian Ĥ of the quantum system in the form [34]   where 0 H is the Hamiltonian of the unperturbed system and describes the equilibrium system.We suppose that the external perturbation is in some sense small.A is a linear operator which connects external timedependent force   h t   with the system.shall use the quantum-mechanical equations for the system dynamical response is the average value of in equilibrium; the sign + means conjugate value.It is necessary to approximate using the well-known perturbation expansion through first order in .We have where

   
We consider the functional of the dynamic response of the form [34] 0 0 0 is the time moment when the perturbating field is turned on, . Now we apply the aforementioned expressions for dynamic response to a magnetic droplet system.The response B t   M t is then the induced magnetization of the system.When one measures the ac susceptibility in spin glasses in the external magnetic oscillating field it is observed an aging effect: magnetic response of the system to the weak external field depends on the thermal history of the sample, on the time during which the system was kept in a spin glass phase.The sample is quenched in zero magnetic field from temperature The system is probed at the time after quench end ("the age").Using the linear response theory the magnetization of magnetic system is [35]

 
, e i t t t t will be equal to [15,35] t © 2013 We consider the behavior of the magnetic droplet system described by the Hamiltonian H and under ac . In our calculations, we suppose that There is a complicated crossover between classic and quantum behaviors of the droplets which depends on temperature, ac field frequency and length scale .According to [26,27] the dynamical crossover length is determined from the . The system behaves presumably classically or quantum mechanically when the dominant length scale is above or below dyn   .
for fixed frequency Following to aforementioned quantum droplet theory with model Hamiltonian (4) and domain growth ideas, we calculate the magnetic dynamic susceptibility using the dynamical response functional which includes first and second order linear response functions.The contribution of a single droplet to the ac susceptibility up to some factor is, is the Edwards-Anderson order parameter, is the static susceptibility of the droplet . The expression (11) is obtained for low frequencies under the condition   (and L    ), because this condition is used to observe non-stationary dynamics in susceptibility measurements [5].Now we have to average the susceptibility (11) over droplet energies L  and over droplet length scales .We use the droplet energy distribution given by Equation ( 3).Here we as-sume 0   .While integrating over , we note that the susceptibility is dominated by droplets of length scale is the natural length scale  of the problem when L   and it is the low limit ofthe integration over .The upper limit is energies and length scales we obtain the following expression for imaginary part of susceptibility of the droplet system: 1,  The Equation ( 12) is the main result of this paper.
As we see, the susceptibility depends on many parameters of the droplet system and the external ac magnetic field: on the form of droplet energies distribution  L L P  , on the droplet microscopic tunneling rate 0 , on the temperature, on the system "age" , on the ac field frequency and amplitude.Furthermore, we note that the expression (12) consists of terms which are time independent which describe oscillations with frequency, and terms which depend on and define nonstationary non-equilibrium dynamics of the droplet system.Thus the imaginary part of susceptibility can be represented as a sum of stationary part   , t   (Figure 3).The susceptibility quickly goes down and then slowly decays to some value with oscillations.Then we observe the stationary behavior of susceptibility.In particular, in Figure 2 we observe as, on longer times, the quantum fluctuations ( 0  dependence) becomes irrelevant. .On long times the susceptibility curve is a periodical function oscillating near some constant value (stationary process).We find temperature dependence of imaginary part of susceptibility and show how the quantum fluctuations influence the dynamic susceptibility of the droplet system at very low temperatures.If the ac field frequency increases then the nonequilibrium dynamics is suppressed.Thus the droplet system response to an external perturbing field depends on its thermal history.
In [17] it is shown that the behavior of response function       the dynamics is non-stationary.In [17] it is shown that quantum fluctuations in quantum glassy systems depress the phase transition temperature, in a glassy phase the aging effect survives the quantum fluctuations, and because of quantum fluctuations, the fluctuation-dissipation theorem is modified.In reference [19][20][21] it is shown that all terms in the dynamical equations governing the timeevolution of spin response and correlation function which are due to quantum effects, are irrelevant at long times.Quantum effects enter only through the renormalization of parameters in dynamical equations [19][20][21].The behavior of spin response as function of τ in [36] is similar to behavior of dynamical susceptibility in our paper.In [36] it is shown that quantum fluctuation slightly influences the aging regime and the quantum system behavior is approximately classic.
As far as we know there are no experiments on quantum spin glasses.In papers [17,23] there are experimental data for dynamic susceptibility in classic spin glasses.P. Svedlindh et al. [35] have investigated the behavior of and have found that decay is close to a logarithmic one.Shins et al. [22] also show that susceptibility decays with time in a nearly logarithmic way.
In our quantum system at very low temperature, we cannot find agreement with these data because classical and quantum spin glass has, in general, different behavior.We may compare our results with experimental data only approximately because these experimental data are on classical spin glasses and we consider here low-temperature dynamics of quantum spin glass.We observe a qualitatively similar behavior in the range of the small times elapsed since the quench.

.
The droplet model of classical Ising spin glass was considered by D.S. Fisher and D.A. Huse


are the Pauli matrices representing the two states of the droplet L D ; the sum is over all droplets L D at length scale L and over all length scales L The droplet length scale L is more or of order of the correlation length.The value   0 Lwhich regulates the strength of quantum fluctuations ( L corresponds to the classical limit) is the tunneling rate for a droplet of linear size L and  is a coefficient which is approximately the same for all droplets.0


be applied in z-direction.Here and h  are the amplitude and frequency of the ac field.
at time  .At this moment a very small external magnetic oscillating field   h t is applied to measure the ac susceptibility of the sample.The evolution continues in isothermal conditions, ac  is measured as a function of the time elapsed since the sample reached the temperature at

Figure 2
Figure 2. Imaginary part as function of time    t ,   t 1 50  0.01 at T  and 0.05    16 0 10   

Figure 3 3
Figure 3. Imaginary part as function of time    t ,  

Figure 1 4 .
Figure 1.Imaginary part as function of time at and