arxiv: v1 [cond-mat.stat-mech] 1 Jan 2017

Transcription

1 Pis ma v ZhETF Full relica symmetry breaking in -sin-glass-like systems T. I. Schelkacheva 1, and N. M. Chtchelkatchev Institute for High Pressure Physics, Russian Academy of Sciences, 14190, Troitsk, Moscow, Russia L.D. Landau Institute for Theoretical Physics, Russian Academy of Sciences, 1443, Moscow Region, Chernogolovka, Russia 3 Ural Federal University, 6000 Ekaterinburg, Russia arxiv: v1 cond-mat.stat-mech 1 Jan Deartment of Theoretical Physics, Moscow Institute of Physics and Technology, Moscow, Russia 5 All-Russia Research Institute of Automatics, Suschevskaya, Moscow 17055, Russia 6 Institute of Metallurgy, Ural Division of Russian Academy of Sciences, Yekaterinburg 60016, Russia It is shown that continuously changing the effective number of interacting articles in -sin-glass-like model allows to describe the transition from the full relica symmetry breaking glass solution to stable first relica symmetry breaking glass solution in the case of non-reflective symmetry diagonal oerators used instead of Ising sins. As an examle, axial quadruole moments in lace of Ising sins are considered and the boundary value c1 =.5 is found. Introduction The basis of understanding glasses is the Sherrington Kirkatrick SK) model 1: the Ising model with random links. A stable solution for SK model was obtained by Parisi, 3 with a full relica symmetry breaking FB) scheme. Later it was realised that relica symmetry is not abstract and academic question but it corresonds to formation of the secific hierarchy of basins in the energy landscae of the glass forming system. A natural generalization of the SK model with air interaction of sins is a model with -sin interactions 3, 4. Unlike of the SK model, -sin model has a stable first relica symmetry breaking B) solution that arises abrutly. A very low-temerature boundary of the 1B stability region is given by so called Gardner transition temerature intensively discussed last time 4, 5 where a valley in configuration sace transforms to a multitude of searated basins. Now there is a reborn of interest to sin models showingglassybehaviour5, 6,7,8,9,10,11,1,13, 14. It turned out that these models can qualitatively exlain hysics of real glasses 15. On the other hand, there is limited number of analytically solvable glassy models and each such model is interesting itself. Here we roose analytical solution of -sin-glass-like system and discuss hysical alications of this model. For a long time there was a conjecture that there are moreorlesstwoclassesofmodels, deendingonhowthe relica symmetry breaking aears 18. In one class of models full relica symmetry breaking FB) occurs continuously at the transition oint from the aramagnetic state to the glass state like, e.g., in SK model). The second class of models can be called 1B-models -sin model, Potts models). In this case there is a finite range of temeratures where stable 1B glass solution occurs. What is imortant that this 1B solution mostly aears abrutly. 1B-models and esecially -sin glasses in recent years attract much interest in connection with the fact that there is a close relationshi between static relica aroach and dynamic consideration. For examle the Random First Order Transition theory for structural glasses is insired by the -sin glass model 5, 6, 7, 8, 9, 10, 11, 1, 13, 14, 15. Itshouldbealsonotedthatthetwoclassesofmodels mentioned above distinguish essentially by their energy landscae 16, which is an imortant concet in the dynamics of liquids and glasses. In the context of SK-like and 1B-models it is ossible to develo very advanced theoretical tools that can be reused in other contexts. These relatively simle models based on well-designed solutions allow to exlore qualitatively an extensive range of is- we investigate: full B stable 1B continuous 1B transition discontinuous 1B transition c1 c Fig.1 We have got FB glass solution of a -sin-like model when the effective number of interacting articles < < c1. We considered general diagonal oerators Û with TrÛ k+1 0, k = 1,,... instead of Ising sins. 1

2 T. I. Schelkacheva, and N. M. Chtchelkatchev sues, ranging from magnetic to structural glass transitions 5, 6, 7, 8, 9, 10, 11, 1, 13, 14, 15, 16, 17. Therefore, a natural generalizations of these basic models leads to a successful descrition of the various tyes of glasses, such as orientation glasses and cluster glasses. Relacing the simle Ising sins by more comlex oerators Û dramatically exands the range of solvable tasks. The oerator Û can be, for examle, the axial quadruole moment quadruole glass in molecular solid hydrogen at different ressures) or the role of Û can be layed by certain combinations of the cubic harmonics orientational glass of the clusters, for examle, in C 60 in a wide ressure range), see 18, 19, 0, 1,, 3, 4, 5, 6. It was shown recently that many statements related to the models based on Ising sins can be alied to -sin-like models that use instead of sins diagonal oerators Û when there is reflection symmetry : TrÛk+1 = 0, k = 1,... For such models we can build a stable 1B solution for > in a wide region of temeratures. It was shown that the oint = is secial for such models 6, 7, 8. For -sin-like model we reviously have got significantly different solutions when diagonal oerators Û have broken reflective symmetry, then TrÛk+1 0 like e.g., for quadruole oerators). And for these models 1B glass solution behavior have been recently investigated near the glass transition at different continuous) 6, 9. It turned out that there is a finite region of instability of 1B solutions for < c1 where c1 is determined by the secific form of Û. 1B solution is stable for > c1. Wherein the transition from arahase to 1B glass is continuous for c1 < < c. When > c 6, 9 1B glass occurs abrutly just as in the conventional -sin model. We should note that c is not universal, but it deends on the articular tye of Û. We have built 30, 31 FBSB solution for these models with a air interaction =. In this letter we investigate in detail the generalised -sin glass forming models in the region < c1 where instability of 1B glass solution is exected. We build a solution with full relica symmetry breaking. The very existence of the domain with full relica symmetry breaking is a surrising result esecially comared with the traditional -sin model of Ising sins. Continuously changing the effective number of interacting articles in -sin model allows us to describe the crossover from the full relica symmetry breaking glass solution to stable first relica symmetry breaking glass solution in our case of non-reflective symmetry diagonal oerators used instead of Ising sins). For illustrating examle we take oerators of axial quadruole moments in lace of Ising sins. For this model we find the boundary value c1 =.5. The model The staring oint is the -sin-glass-like Hamiltonian H = i 1 i... i J i1...i Û i1 Û i...ûi, 1) where the quenched interactions J i1...i are distributed with Gaussian robability: N 1 PJ i1...i ) = ex J i 1...i ) N 1!πJ!J. ) Here arbitrary non-reflective symmetry diagonal oerators Û are located on the lattice sites i instead of Ising sins, N is the number of sites. We remind that it imlies that TrÛk+1 0, k= 1,,... Here we investigate how order arameters of this model develo with the continuous arameter and the secific tye of the oerators Û. We write down in standard way the disorder averaged free energy using the relica aroach 3. We assume like it was done in 30, 31 that the order arameter deviations αβ are small from relica symmetry order arameterq. To describe the B-solution near the bifurcation oint temerature T 0, where transforms into B, we exand the exression for the free energy u to the fourth order of αβ. We write below in 3) the deviation F) of the free energy from its -art. It is imortant to note that the exression for Free energy have not been written before for arbitrary. We should also emhasize that it includes the terms with odd number of identical relica indices 33, 34, 35 unlike Ising-sin SK-models. All coefficients in 3) deend only on -solution at T 0. They are given in Sulementary Material. The rime on meansthat onlythe suerscritsbelonging to the same are necessarily different. The coefficients of the second and third orders, λ )rel,l,b 3,..B 4, have been obtained earlier 18, 6, 7, but we also write them in Sulementary Material for readability. This is the only overla with our revious ublications.

3 FB in -sin-glass-like systems... 3 F) NkT = lim n 0 { 1 t n 4 1) q ) λ))rel αβ ) t 4 L) αβ α t B 6 ), αβ αγ β + B ) αβ αγ α +B 3 ) αβ βγ γα +B 3) αβ ) αγ +B 4 ) αβ ) 3 + t 8 D ) αβ ) 4 +D31 ) αβ ) 3 αγ +D 3 ) αβ ) α) +D33 ) αβ ) αγ γβ +, D 4 ) αβ ) αγ α +D 43 ) αβ ) αγ β +D 45 ) αβ ) αγ γ + D 46 ) αβ αγ α βγ +D 47 ) D 54 ),µ αβ βγ γ α +D 53 ) αβ αγ α q βµ +D 55 ),µ,µ αβ αγ α q αµ + αβ αγ γ q µ }, 3) The order arameters and the relicon mode λ resonsible for B stability) we find as follows: λ) )rel = 1 t 1) q ) dz G TrÛ ex ˆθ ) TrÛ exˆθ Trexˆθ Trex ˆθ ; 4) ) TrÛ ex ˆθ q) = dz G ) ; 5) Tr exˆθ Û ˆθ) w) = dz GTr ex ) ; 6) Tr exˆθ dz π ex where dz G = t 0 + t; label n relicas and ˆθ) = zt q 1) z ), t = J/kT = Û +t w 1) q 1) Û. 4 FB Below we find answer of two questions: 1) if FB realises for given model we derive simle criterium), ) we rovide necessary tools for calculation of FB order arameter qx). To describe the FB function qx) of the variable x we include in the consideration the fourth-order terms intheexansionof F. We usethe standardformalized Parisi algebra rules, 3 to write the free energy as the functional of qx) and so to construct FB. The terms which are not formally described by the Parisi rules can be reduced to the standard form as well 30, 31. The equation for the order arameter qx) follows from the stationarity condition x) F) = 0. Therefore, we canis a similarwayfind the branching condition aearance of FB glass solution), which were derived in detail for = see 30, 31 and Refs. therein): λ) )rel T0 = 0, which roduces the temerature T 0 ). For clarity we write x) F) = 0 u to the terms of the second order. We get t 0 1) q ) d λ )rel dt t0 tqx) t 4 0L q +... = 0. 7) So if the oerators Û do not have the reflective symmetry therefore L 0 in that case) we need an additional branching condition, q 1 0 qx) = 0 + o t) 31, insuring the aearance of non-trivial new solutions. Integral equation, x) F = 0, that determines the function qx), as usually can be simlified using the differential oerator Ô = 1 d 1 d q q, where q = dqx) : t 6 {B 4 B 3 x}+t { D xD 47 xqx) 1 x } dyqy) + 4D +xd 33 qx) = 0. 8) The coefficients are given in Sulementary Material for.

4 4 T. I. Schelkacheva, and N. M. Chtchelkatchev x D a) b) T a) b) q 0 B4/B Fig. a).the deendence of the denominator x D = D +D 33x D 47x x of the function dqx) x on the effective number of interacting quadruoles obtained from Eq. 9). b) The deendence of the function q = dqx) x on the effective number of interacting quadruoles obtained from Eq. 9). Since qx) can only be a non-decreasing function of x we should consider how the sign of q = dqx) deends on the arameter. We obtain from Eq. 8): q ) = B 3 +4xqx) x 0 dyqy)d 47 +D 33 qx) t 4 D +D 33 x D 47 x. 9) This exression is one of the central results of our aer. Deending on the sign of q we can conclude if the system falls into FB ositive q ) state or not. Ifweconfine ourselvestothe termsofthe thirdorder over in F and 8) then we obtain t 6 {B 4 B 3 x}+... = 0. So the significantly deending on x art of qx) is concentrated in the neighborhood of x) = B 4 )/B 3 ). 10) Only in the case of oerators Û with TrÛk+1) 0 for k = 1,,.. we get x = B 4 /B 3 0. For Ising-like oerators with TrÛk+1) = 0 we get = see 6, 7, 9), B 4 = 0 and so x = 0 in accordance with usual standard case of the Parisi theory, 3. Examle As an illustrating examle, we consider quadruole glass in the sace, J = 1. Oerator Û = ˆQ is the axial quadruole moment and it takes values, 1, 1). So there is no reflection symmetry. Behavior of this model is significantly different from Fig. 3 a) Deendence of branching temerature T 0 of glass solution on the effective number of interacting quadruoles. b) Deendence of x) = B 4)/B 3) on the effective number of interacting quadruoles. what one usually has in the case of Ising-like oerators, 3, 18, 19, 0, 1,, 3, 4, 5, 6, 8. We have done earlier calculations of 1B solution for various values of that could change continuously 6. It has been found that for < < c1 c1 =.5) the solution is qualitatively similar to one which occurs in the case of air interaction = ). This behavior is natural for continuity reason. At high temeratures T > T 0 we remind that T 0 is the bifurcation temerature) there is a stable non-trivial solution for the order arameter q. But q 0 for any finite temerature T > T 0 unlike conventional Ising-sin model. For T < T 0 we have Almeida-Thoulles relicon mode λ )rel < 0 and -solution becomes unstable. Similarly, when T < T 0 we have λ 1B)rel < 0 determining the condition when the 1B glass is unstable 19, 4, 6, 7. So we investigated the stability of 1B solution and found that1bsolutionisunstablewhen <.5. The 1B solution becomes stable when >.5, see Fig. 1. The instability region of 1B solution is an area in which there may be FB solution. Using Eqs. 9)-10), we can describe FB solution for T < T 0 neart 0 for < < c1. We alsowant to find the value of c1 just numerically solving the equations Eqs. 9)-10). We intend to comare this result of c1 to c1 =.5 we had reviously received, considering the stability of the 1B solution 6. For air interaction, =, we earlier obtained 31

5 FB in -sin-glass-like systems... 5 that T 0 = 1.37 and x = 0.43 for FB. Wherein q = qx) is a continuous increasing function of x in narrow neighborhood, e.g., x = near x for T T 0 ) = 0.. For other values of x, not in the neighborhood of x, qx) is very close to q 1B. So in our case < < =.5 we can with good accuracy ut q x) = q x) 1B = q0) 1B and xq x) x 0 dyqy) = 0 in 9). We also have got that q x) 1B is a small quantity and so we can neglect D 33 q x) comared with B 3. For other values of x not very close to x in exactly the same way as it was obtained for = 31 we find that qx) is very close to q 1B. The results of our calculations are resented in Figs. -3. At c1 =.5 the function dqx) x diverges and changes its sign when > c1 contrary to its conventional robabilistic interretation. It follows that FB is imossible for >.5. This result agrees with c1 obtained by us reviously for 1B solution stability border 6. For clarity we show the deendence the denominator x D = D +D 33 x D 47 x x on in Fig.. We show in Fig. 3 the deendence of the branching temerature T 0 and the deendence of x) = B 4 )/B 3 ) on the effective number of interacting quadruoles. Conclusions For the first time it is shown that continuously changing of one of the arameters of glass forming model with one tye of interaction allows moving from FB glass to stable 1B solution. As an examle the generalized -sin quadruole glass model is considered. Acknowledgments This work was suorted in art by the Russian Foundation for Basic Research No ) while numerical simulations were funded by Russian Scientific Foundation grant No ). 1. D. Sherrington, and S. Kirkatrick, Phys. Rev. Lett. 3, ); S. Kirkatrick, and D. Sherrington, Phys. Rev. B 17, ).. G. Parisi, J.Phys.A: Math. Theor. 13, L ). 3. M. Mezard, G. Parisi, and M. Virasoro, Sin Glass Theory and Beyond World Scientific, Singaore, 1987). 4. E. Gardner, Nuc. Phys. B57, ). 5. P. Charbonneau, Y. Jin, G. Parisi, C. Rainone, B. Seoane, F. Zamoni Phys. Rev. E 9, ). 6. A. Crisanti, H. Horner, and H. Sommers, Z. Phys. B 9, ). 7. T. R. Kirkatrick, D. Thirumalai and P. G. Wolynes, Phys. Rev. A ). 8. T. R.Kirkatrick, D. Thirumalai, Phys. Rev. Lett. 58, ), Phys. Rev. B 36, ). 9. S. Franz, H. Jacquin, G. Parisi, P. Urbani, F. Zamoni, J. Chem. Phys. 138, 1A ). 10. F. Caltagirone, U. Ferrari, L. Leuzzi, G. Parisi, F. Ricci-Tersenghi, T. Rizzo, Phys. Rev. Lett. 108, ). 11. T. Rizzo, Euro Phys. Lett. 106, ). 1. T. Rizzo, Phys. Rev. E 87, ). 13. G. Parisi and F. Zamoni, Rev. Mod. Phys. 8, ). 14. L. Berthier and G. Biroli, Rev. Mod. Phys. 83, ). 15. P.G. Wolynes and V. Lubchenko, Structural Glasses and Suercooled Liquids: Theory, Exeriment, and Alications Wiley, com, 01). 16. U. Buchenau, J. Phys. C 15, S ). 17. V. Alba, S. Inglis, L. Pollet, Phys. Rev. B 93, ). 18. T. I. Schelkacheva, E. E. Tareyeva, and N. M. Chtchelkatchev, Phys. Rev. E 79, ). 19. T.I. Schelkacheva, E.E. Tareyeva, and N.M. Chtchelkatchev, Phys. Rev. B 8, ). 0. T.I. Schelkacheva, E.E. Tareyeva, and N.M. Chtchelkatchev, Phys. Rev. B 76, ). 1. E.A.Lutchinskaia and E.E. Tareyeva, Phys. Rev. B 5, ).. N.M. Chtchelkatchev, V.N. Ryzhov, T.I. Schelkacheva, and E.E. Tareyeva, Phys. Lett. A 39, ). 3. K. Walasek, J. Phys. A 8, L ). 4. E.E.Tareyeva, T.I.Schelkacheva and N.M.Chtchelkatchev, Theor. Math. Phys. 160, ) 5. E.E. Tareyeva, T.I. Schelkacheva and N.M. Chtchelkatchev, Theor. Math. Phys. 155, ). 6. E.E. Tareyeva, T.I. Schelkacheva, N.M. Chtchelkatchev, J. Phys. A 47, ). 7. T.I. Schelkacheva, N.M.Chtchelkatchev, J. Phys. A 44, ). 8. T.I. Schelkacheva, E.E.Tareyeva, and N.M.Chtchelkatchev, Phys. Lett. A 358, 006). 9. E.E.Tareyeva, T.I.Schelkacheva and N.M.Chtchelkatchev, Theor. Math. Phys. 18, ). 30. E.E. Tareyeva, T.I. Schelkacheva, N.M. Chtchelkatchev, Phys. Lett. A 377, ). 31. T.I. Schelkacheva, E.E. Tareyeva, N.M. Chtchelkatchev, Phys. Rev. E 89, ). 3. A. Crisanti, L. Leuzzi, Phys. Rev. Lett. 93, ). 33. T. Temesvari, C. De Dominicis, I. R. Pimentel: Eur. Phys. J. B 5, ). 34. T. Temesvari, Nuclear Physics B, 89, ). 35. A. A. Crisanti, C. De Dominicis, J. Phys. A 43, ).