Inf<strong>in</strong>ite-range sp<strong>in</strong> glass3D <strong>Coulomb</strong> glassequilibriumglass phaseyesequilibriumglass phase?saturatedl<strong>in</strong>ear gap<strong>in</strong> 1-DOSTABLE I: Critical exponents for <strong>the</strong> fluid-COP transitionyesalong BC <strong>in</strong> Fig. 1, compared with <strong>the</strong> RFIM values [37].γ/ν ¯γ/ν β/ν ν<strong>Coulomb</strong> glass 1.69(17) 2.89(15) 0.06(3) 1.11(12)RFIM 1.44(12) 2.93(11) 0.011(3) 1.37(9)We thank A.Möbius for helpful correspondence. Thisresearch is supported by <strong>the</strong> Generalitat de Catalunya<strong>and</strong> <strong>the</strong> M<strong>in</strong>isterio de Ciencia e Innovación (contractsFIS-2006-13321-C02-01 <strong>and</strong> AP2007-01005). The computationswere performed on <strong>the</strong> Spanish Supercomput<strong>in</strong>gNetwork (RES) node at Universidad de Cantabria.saturatedMP thanks <strong>the</strong> Aspen Center for Physics where part ofthis work was carried out.scale-freequadratic gap<strong>avalanches</strong>∗ Electronic address: palass<strong>in</strong>i@ub.edu[1] M. Pollak, Discuss. Faraday Soc. 50, 13 (1970); G. Sr<strong>in</strong>ivasan,Phys. Rev. B 4, 2581 (1971).<strong>in</strong> 1-DOSFIG. 4: (Color onl<strong>in</strong>e) Charge-order correlation length along [2] A. L. Efros <strong>and</strong> B. I. Shklovskii, J. Phys. C 8, 49 (1975).path BC <strong>in</strong> Fig. 1. Top <strong>in</strong>set: order parameter Ms along paths [3] E. I. Lev<strong>in</strong>, V. L. Nguyen, B. I. Shklovskii, <strong>and</strong> A.AB, BC, <strong>and</strong> DE. Bottom <strong>in</strong>set: specific heat along path BC. L. Efros, Sov. Phys. JETP 65, 842 (1987).[4] A. Hunt, Philos. Mag. Lett. 62, 371 (1990).[5] T. Vojta, W. John, <strong>and</strong> M. Schreiber, J. Phys. Condens.GP does not exist above <strong>the</strong> dashed l<strong>in</strong>e <strong>in</strong> Fig. 1.Matter 5, 4989 (1993).Critical behavior – S<strong>in</strong>ce at W = 0 <strong>the</strong> fluid-COP transitionhas a positive specific-heat exponent [23], disorder [7] J. G. Massey <strong>and</strong> M. Lee, Phys. Rev. Lett. 75, 4266[6] See e.g. Ref.[7] <strong>and</strong> references <strong>the</strong>re<strong>in</strong>.is relevant [34] <strong>and</strong> <strong>the</strong> W ≠ 0 transition will be governedby a r<strong>and</strong>om fixed po<strong>in</strong>t which, by analogy with <strong>the</strong> [8] B. S<strong>and</strong>ow, K. Gloos, R. Rentzsch, A. N. Ionov, <strong>and</strong> W.(1995).RFIM [12], we expect to be at zero temperature [35]. Assum<strong>in</strong>g<strong>the</strong> transition is second order, as we ascerta<strong>in</strong>edSchirmacher, Phys. Rev. Lett. 86, 1845 (2001).[9] Q. Li <strong>and</strong> P. Phillips, Phys. Rev. B 49, 10269 (1994).[10] M.H. Overl<strong>in</strong>, L.A. Wong, <strong>and</strong> C. C. Yu, Phys. Rev. Bby <strong>in</strong>spect<strong>in</strong>g <strong>the</strong> distribution of ms for <strong>in</strong>dividual samples,we obta<strong>in</strong> <strong>the</strong> critical exponents <strong>in</strong> Table I. We es-[11] M. Sarvestani, M. Schreiber, <strong>and</strong> T. Vojta, Phys. Rev. B70, 214203 (2004).timated β/ν <strong>and</strong> ¯γ/ν with <strong>the</strong> quotient method [36] for 52, R3820 (1995).<strong>the</strong> observables Ms <strong>and</strong> ¯χ L = N[〈ms 2 〉]av respectively [12] A. Möbius, M. Richter, <strong>and</strong> B. Drittler, Phys. Rev. B 45,[<strong>the</strong> quotient estimates from (L, L ′ ) = (6, 8), (6, 10) <strong>and</strong> 11568 (1992).(8, 10) agree with<strong>in</strong> <strong>the</strong> errors], while γ/ν was obta<strong>in</strong>ed [13] J. H. Davies, P. A. Lee, <strong>and</strong> T. M. Rice, Phys. Rev. Bby fitt<strong>in</strong>g aL γ/ν29, 4260 (1984).to <strong>the</strong> height of <strong>the</strong> peak of <strong>the</strong> susceptibilityN[〈ms 2 〉−〈|ms|〉 2 [14] See e.g. Ref.[21] <strong>and</strong> references <strong>the</strong>re<strong>in</strong>.]av (data not shown). Thespecific heat cL =1/(NT 2 )[〈H 2 〉 − 〈H〉 2 [15] A. Díaz-Sánchez, M. Ortuño, A. Pérez-Garrido, <strong>and</strong>]av shows a peak E. Cuevas, Phys. Stat. Sol. (b) 218, 11 (2000).that <strong>in</strong>creases slowly with L (Fig. 4, bottom <strong>in</strong>set), which [16] E. R. Grannan <strong>and</strong> C. C. Yu, Phys. Rev. Lett. 71, 3335suggests ei<strong>the</strong>r α < 0 or a logarithmic divergence (α = 0). (1993).We could not estimate ν directly <strong>in</strong> a reliable way, but [17] B. Surer, H. G. Katzgraber, G. T. Zimanyi,we obta<strong>in</strong> ν =1.11(12) from <strong>the</strong> modified hyperscal<strong>in</strong>grelation [35] (d − θ)ν =2− α, assum<strong>in</strong>g α = 0 <strong>and</strong> us<strong>in</strong>gθ =¯γ/ν −γ/ν =1.20(20). As shown <strong>in</strong> Table I, <strong>the</strong> criticalexponents agree fairly well with <strong>the</strong> known values formean-fieldB. A. Allgood,<strong>and</strong> G. Blatter, cond-mat/08054640 (2008).<strong>the</strong>ory[18] A. A. Pastor <strong>and</strong> V. Dobrosavljević, Phys. Rev. Lett. 83,4642 (1999).[19] S. Pankov <strong>and</strong> V. Dobrosavljević, Phys. Rev. Lett. 94,046402 (2005).<strong>the</strong> RFIM [37], which suggests that <strong>the</strong> system rema<strong>in</strong>s [20] M. Müller <strong>and</strong> L. B. Ioffe, Phys. Rev. Lett. 93, 256403effectively short-range near <strong>the</strong> transition.(2004).To conclude, our results show that, although meanfield<strong>the</strong>ory seems to capture correctly <strong>the</strong> DOS near <strong>the</strong>[21] M. Müller <strong>and</strong> S. Pankov, Phys. Rev. B 75, 144201(2007).[22] R. Dickman <strong>and</strong> G. Stell, cond-mat/9906364 (1999).<strong>Coulomb</strong> gap, <strong>the</strong> <strong>in</strong>teraction rema<strong>in</strong>s well screened <strong>and</strong>[23] A. Möbius <strong>and</strong> U. K. Rößler, cond-mat/0309001 (1999).<strong>the</strong> correlations rema<strong>in</strong> short-range down to ra<strong>the</strong>r low [24] V. Malik <strong>and</strong> D. Kumar, Phys. Rev. B 76, 125207 (2007).temperatures. If an equilibrium glass phase exists, its [25] A. Möbius, talk given at TIDS11 (2005).onset must be at exceptionally low temperatures.[26] M. Goe<strong>the</strong> <strong>and</strong> M. Palass<strong>in</strong>i, <strong>in</strong> preparation.4beyond mean fieldscale-free?<strong>avalanches</strong>?Pankov, Dobrosavljevic 2005Müller, Ioffe 2005Müller, Pankov 2007
“Equilibrium glass phase”Inf<strong>in</strong>ite-range sp<strong>in</strong> glass3D <strong>Coulomb</strong> glassTT goverlapP (q)q = 1 N∑Si a Sibparamagnetqi3D mean field <strong>the</strong>ory :(locator approximation)Pankov, Dobrosavljevic 2005Müller, Ioffe 2005Müller, Pankov 2007T < T g ∼ 1/ √ WCHARGE ORDERP (q)Almeida-Thouless l<strong>in</strong>eglassqGoe<strong>the</strong>, Palass<strong>in</strong>i 2009∼ e αN- Many (# ) equilibrium states- “Marg<strong>in</strong>al criticality”:power law correlations belowT g(〈S i S j 〉−〈S i 〉〈S j 〉) 2 ∼ 1/r θ ijhT0.150.10.050ACFluidBEChargeOrderD0 0.1 0.2 0.3 0.4 0.5WEquilibrium MC:no glass phasedown to verylow T(see alsoSurer et al. 2009)