arxiv:cond-mat/ v3 22 May 2003

Transcription

1 The Fokker-Planck operator at a contnuous phase transton Moshe Schwartz Rayond and Beverly Sackler Faculty of Exact Scences School of Physcs and Astronoy, Tel Avv Unversty, arxv:cond-at/315 v3 May 3 Raat Avv, Tel Avv 69978, Israel I consder a physcal syste descrbed by a contnuous feld theory and enclosed n a large but fnte cubcal box wth perodc boundary condtons. The syste s assued to undergo a contnuous phase transton at soe 4 crtcal pont. The theory that s a contnuous verson of the Isng odel s such a syste but there are any other exaples correspondng to hgher spn, hgher syetry etc. The egenfunctons of the correspondng Fokker-Planck operator can be chosen, of course, to be egenfunctons of the oentu operator. It s shown that the egenvalues of the FP operator, correspondng to each egenvalue of the oentu operator, evaluated at a transton pont of the fnte syste, accuulate at zero, when the sze of the syste tends to nfnty. There are any reasonable ways of defnng a crtcal teperature of a fnte syste, that tends to the crtcal teperature of the nfnte syste as the sze of the syste tends to nfnty. The accuulaton of egenvalues s nether affected by the specfc choce of crtcal teperature of the fnte syste nor by whether the syste s below or above ts upper crtcal denson.

2 The property of crtcal slowng down [1,] s known for a long te fro experental [3-7] and nuercal work [8-11]. The uanttatve descrpton s n ters of characterstc decay tes or alternatvely characterstc freuences, that govern the decay of a dsturbance of wave vector. It was found that at the z transton, behaves as soe postve power of,. Ths ples that the larger the scale of the dsturbance the longer t takes to decay and a dvergent scale results n a dvergent decay te. The evaluaton of the exponent z was the subject of theoretcal work usng a nuber of dfferent approaches but all based on stochastc feld euatons of the Langevn type to descrbe the dynacs of the syste [1-18]. More recently, t was suggested that not only do the characterstc freuences tend to zero wth but the decay functon tself becoes slower than an exponental (e.g. stretched exponental) for long tes, 1[19-1]. That property was shown to t hold for ute general Langevn feld euatons ncludng n addton to crtcal dynacs, euatons of the KPZ type. The dervaton of stretched exponental decay or any other for of slow decay has to rely on approxatons. As wll becoe evdent later (e.(1)) the te dependent structure factor needs evaluatons of the egenstates and egenvalues of a uantu feld Haltonan (e.(1)) obtaned by a standard slarty transforaton fro the correspondng Fokker-Planck operator. The accuracy wth whch these have to be evaluated has to ncrease wth the te arguent of the structure factor. Therefore the proble of very long te decay s extreely dffcult and t s portant to have as any exact results as possble n support of approxate dervatons. Fro e. (1) t wll becoe clear also that a necessary condton for slow decay of a dsturbance of wave vector s to have enough egenvalues correspondng to states carryng oentu n the vcnty of zero. It has been shown n a recent publcaton [], that there ust exst an egenvector of the FP operator, at the transton pont, that carres oentu, and has an egenvalue that s as close as we wsh to zero. In the present artcle, I go beyond that and show that for each the egenvalues of the FP operator accuulate at zero for the nfnte syste. The proof I present reles on the splest establshed lowest order results of fnte sze scalng eulbru theory and does not nvolve

3 even the next order correctons such as the dependence of a properly chosen crtcal teperature of the fnte syste on ts sze [, 3]. To y knowledge, those results have never been proven rgorously. Nevertheless ther use as conjectures n the present proof yelds nforaton about a non eulbru proble that s consderably ore dffcult and that nforaton s exact provded that those wdely accepted eulbru results are ndeed correct. Consder a syste descrbed n ters of a scalar feld and enclosed n a cube wth perodc boundary condtons. The statc statstcal propertes of the syste, are assued to be gven by the Gbbs dstrbuton, exp[ W ], where W s the classcal densonless Haltonan of the syste. I have n nd Haltonans that 4 are even functonals of the feld and are generalzatons of the theory that n ter of Fourer coponents of the feld has the for P e u l 1 l l 3 ( l1 l l3 ) 1 W r, (1) where the Fourer coponents of the feld defned by 1 ( r)exp( r) d r. () and s the volue of the syste. The oentu ndces are assued to be bound fro above by soe hgh oentu cut off and ther allowed values are gven by n, (3) L where L s the lnear sze of the syste and n a vector of ntegers. The followng dscusson and the proof I present s not lted, however, to the above typcal odel. The dscusson clearly covers other systes wth dfferent dependences ofw on the feld (e.g. systes descrbng hgher spns) or systes n whch the feld s not a scalar (e.g. systes wth O(n) syetry). As wll becoe evdent, the only property that I use s that there s a crtcal pont n the large volue lt and that at

4 that pont and n that lt the structure factor or ts analog n the case of vector felds, dverges at sall s as rased to a power saller than -1. The dynacs s assued to be descrbed n ters of a set of Langevn feld euatons d W, (4) d where the nose obeys t and ' ' t. (5) A standard procedure leads to the Fokker-Planck euaton for the probablty dstrbuton of the felds, P, P t W P OP, (6) whereo s the FP lnear operator. The dstrbuton at eulbru, P e s the egenstate of the Fokker-Planck operator, O, wth egenvalue zero. It follows drectly fro euaton (6), that P exp[-w], (7) e s the Gbbs dstrbuton. I wll be nterested n the te dependent structure factor () ( ), where the eanng of the average s as follows: s easured at te t= at eulbru, the syste s then allowed to evolve freely and atheatcal for of the above stateent s s easured at te t. The ;, ( ) () (1) D D P P t, (8) e

5 () (1) where P e s noralzed and P { ;, } s the soluton of the Fokker-Planck () euaton (5) for a dstrbuton n the varables wth ntal condton () (1) ( ) (1), P. (9) ; l l l A standard transforaton, P P 1/ e, whch nduces a slarty transforaton on O, brngs the Fokker-Planck e. (6) and the defnton (8) nto fors ore falar fro uantu echancs. 1 W 1 W H, (1) and () ( ) exp[-h t]. (11) The Haltonan, H, obtaned by the slarty transforaton fro O, s Hertan and all ts egenvalues are postve apart fro the egenvalue zero, whch corresponds to the ground state. Snce the syste s translatonally nvarant, H coutes wth the oentu operator so that the egenstates of H can be chosen to carry defnte oenta. Let { n } be the set of egenstates of H carryng oentu and n the correspondng egenvalues. It s easly verfed that () ( ) n exp[ n ]. (1)

6 (Ths follows fro the fact that s an egenstate of the oentu operator wth oentu ). The spectru of H lted to the states { n }, n, has thus drect bearng on the decay of the te dependent structure factor. For exaple f the spectru n has a gap, the decay of the te dependent structure factor s faster than Defne next e t. n / n, (13) The te dependent correlaton functon can be wrtten as t exp[ ], (14) d It s clear, that for a fnte syste, t wll always be possble to fnd a regon ncludng the orgn and dependng on the sze of the syste, such that () vanshes dentcally n that regon, plyng a decay n te that s faster than soe exponental. Ths ay be characterzed by a lfete that depends on the sze of the syste. The ueston s how does () behave n the lt of an nfnte syste. If n that lt () does not vansh too fast as tends to zero, the decay wll be slower than exponental (e.g. stretched exponental, power law etc.) The defnton of () ples that n order for t not to vansh too fast as tends to zero two condtons have to be et. Frst, the egenvalues n have to accuulate at zero as the sze of the syste tends to nfnty. The second condton s that the atrx eleents n e.(13) do not vansh too fast as the correspondng egenvalues tend to zero. The actual long te dependence of the decay s strongly affected by the behavour of those atrx eleents. The accuulaton of egenvalues to be dscussed next s thus a necessary condton for slow decay but t does not deterne the for of the decay. In fact, t s not even suffcent. The structure factor n the lnear syste to be dscussed n the followng as a prelude to the general case decays exponentally n spte of the relevant egenvalues accuulatng at zero. The reason s that n (1) / 1/ s an exact noralzed egenstate of the Haltonan

7 (1) so that apart fro the atrx eleent n all other atrx eleents n e. (1) vansh. The proof I ll present concerns only the frst necessary condton that as the sze of the syste tends to nfnty the egenvalues of O accuulate at zero for each. Ths has no bearng on the actual for of long te decay. The strategy of the proof s as follows: (a) The set of any exctaton states carryng oentu s defned and the correspondng expectaton values of the Haltonan H are obtaned and shown to accuulate at zero for any. (b) The relaton between those expectaton values and the expectaton values n an orthogonal set of states correspondng to the orgnal alost orthogonal set of any exctaton states s outlned. (c) It s shown how the above ples that the egenvalues of the Haltonan carryng oentu accuulate at zero. Before proceedng to the proof n the general case, I ll descrbe frst a nuber of coon results fro eulbru theory that wll be used as conjectures n the proof and defne soe useful defntons. Then I ll show how the accuulaton of egenvalues of H works for the soluble, lnear case. Ths wll gve soe hnts and otvatons for the general non-lnear case. (1) I keep the value of fxed (e. 3) as the sze of the syste s ncreased by consderng only ultplcaton of the orgnal sze, L, by an nteger. ()There are any reasonable ways to defne the crtcal teperature of the fnte syste. It ay be taken to be the teperature at whch the eulbrtu structure factor at the sallest possble nonzero s axal. It ay be chosen as the teperature where the specfc heat or one of ts dervatves wth respect to teperature s axal. For our purpose t s enough even to choose t as the transton teperature of the nfnte syste. (4) Regardless of the choce the fnte sze crtcal pont the structure factor has the scalng for. ( ) A for 1 L, (15) to leadng order n A and for sall copared to the cut-off oentu. Note that snce s kept fxed the condton L >>1 can always be attaned. Note that that s true also above the upper crtcal denson where = regardless of the fact that hyper scalng does not hold. Hgher order correlaton functons wll also be needed. I

8 wll need k... 1 k, where s even and k but no subset of the k s sus up to zero. For a large fnte syste, these correlatons have at the crtcal pont the scalng for 1 k... k A f k 1,..., k, (16) 1 d / L to leadng order n L. For k ' s sall copared to, that obey also k L >>1 the functon f s hoogeneous. Naely, f we scale all ts varables by, the functon of the scaled varables s related to the functon of the orgnal varables by f k k 1 / f k k,..., 1. (17) 1,..., The pre factor A() s of cobnatoral orgn and we wll need to assue that there exsts a fnte b such that ln A( ) ( b)ln( b). (4) A set of oenta, S l,..., l k 1 k, that can contan, n prncple, the sae oenta a nuber of tes, s sad to be rreducble f t does not contan real subsets, the oenta of whch su up to zero. (5) A set S k conjugate to Sk s the set l,..., l k 1. (6) The su S1 S of the sets S 1 and S, contans all oenta appearng n S1 and S, each oentu appearng exactly n the su of the nuber of tes t appears n each set separately. (7) The any exctaton set s the set of states of the for l,..., l ), (18) 1 l 1 where l 1 and where the set l 1,...,l s rreducble. A physcal proof would probably proceed along the followng lnes. Frst the l s are chosen to eual /. Then es.(1) and () are used to obtan the correspondng ( ) (1 ) expectaton value, /( A ). The thrd step would be to assert that because of

9 e.(17),the any exctaton states for an orthogonal set n the large volue lt (that s also orthogonal to the ground state at any fnte sze, because t carres a dfferent oentu), the egenvalues of the Haltonan n } accuulate at zero. { The atheatcal proof proceeds bascally along the sae lnes but has to take nto account two dffcultes. The frst s that f s an allowed oentu, / s genercally not allowed. The second s uch ore serous. The asserton that the set of representatve states s orthogonal n the large volue lt s probleatc. Many exaples exst n whch the fact that the scalar products of noralzed states tend to zero wth the sze of the syste does not ply conseuences expected fro orthogonal sets n fnte systes. The proof for the Gaussan case (u=) s straght forward. In that case, the states (18) are exact egenstates of the Haltonan H wth egenvalues. l1 l ( / ),..., l, (19) In the followng I wll ntroduce certan restrctons on the sze of, relatng t to.i wll use these restrctons for the Gaussan case as well for the general one. These do not lt the fnal result about the accuulaton of the egenvalues but enable an exact proof under the conjectures about eulbru correlaton functons at the transton. For any gven the sze of the syste can be chosen large enough so that n, the absolute value of the vector n correspondng to, s large enough so that after choosng a nuber n(1/, ), where wll be specfed later (e. (5)), the ntegers n the range ( n,1.5n ) are large enough, so that the followng wll hold: (a) It s possble to choose for each a representatve state of the for (18) such that the szes of the vectors l defnng that representatve are bound fro above by /(-1).(The dea s to choose all the l s as close as possble to choce / snce the l / for all yelds the lowest possble value for the egenvalue n the sector of exctatons states. The above obvous choce s possble, however, because f s an allowed oentu / s genercally not allowed. The condtons on ensure that t s possble to choose thel s n such a way that the correspondng egenvalue (e.(19)) departs n a controlled way fro the absolute

10 nu of the expresson ).b) Furtherore, the set S1 S s rreducble for any 1 and n the range. (Note that ths s not reured for the proof n the Gaussan case but wll splfy the proof n the general case.) Conseuently, the representatve egenvalues are bound fro above l 1,..., l ( / ) ( / ), to leadng order n n. (), rep 1 n Snce n L /, the eanng of e. () s that there are at least 1 L / states wth egenvalues saller than ( / ) L /. Naely, as the sze of the syste tends to nfnty, the egenvalues of H accuulate at zero for any. The proof for the general case wll proceed n two steps. The frst s to show that the expectaton values of H n the any exctatons states (18) accuulate at zero. The second s to show that ths accuulaton ples that the egenvalues of H also accuulates, at zero for any when the sze of the syste tends to nfnty. The representatve states are exactly those chosen before for the Gaussan case. The correspondng expectaton values of the H are 1, rep l1,..., l [ l,[ H,, ]] l / l l 1 1 /. (1) j1 j l l 1 l l 1 The above s an exact expresson that holds for a syste of any sze. The an ponts that lead to the frst eualty above are nvarance under reflecton ( H H and H =. The second eualty n (1) above s due to the fact that the double coutator [,[, ]] l l euals

11 .It ay see that the last expresson on the left hand sde of (1) does l l not depend on the non-lnear couplng u n the Haltonan. Ths s not true of course, because the ground state depends on u. To leadng order n the sze of the syste l l l j l j l l, () 1 j The correctons are ters of the sae order of agntude ultpled by a cobnatoral factor of the order of and dvded by the volue of the syste, the correcton s of the relatve order of L d, that tends to zero for any d (the densonalty of the syste) as the sze of the syste tends to nfnty. Thus, for a large enough syste, the approxaton () can be used to yeld, rep 1,..., [(3/ ) ] l l n. (3) A n ( 1 ) (Note that L here s stll large n spte of the large. In fact t s n and therefore uch larger than 1 so that eans, that there are H are saller than (3/ ) A e.(15) s used for the structure factor.) Ths 1 L / states, at least, for whch the expectaton values of ( L / ) (, L) 1. Snce s always saller than 1 ths ples that as the sze of the syste tends to nfnty, the expectaton values of H accuulate at zero for any gven. Note that so far I have used only the scalng for of the two pont functon (e.(15)) as the only nput conjecture n the proof. If the representatve states above were orthogonal, t would be easy to coplete the proof that the egenvalues of H correspondng to states carryng oentu accuulate at zero. Ths, however, s not the case. The purpose of the next part of the proof s to show that the set of representatve states, descrbed above, can be used to construct an orthogonal set of states, by a Grah-Schdt orthoganolzaton procedure, such that the expectaton values of H n the states of that set are dentcal to

12 leadng order n the sze of the syste to the correspondng expectaton values calculated above. Because of the rreducblty of any set, S1 S the scalar products of pars of the representatve states are gven by e. (16). Denote, for splcty the any exctaton representatve states by and the sallest and largest s between n and 1.5n by s and l respectvely. Further denote by ax and n the axal and nal values respectvely of, j for j n the range s, j. l The orthogonal set s obtaned by constructng lnear cobnatons of the ' s 1 j a, (4) s j j It s not dffcult to show that f n can be ade as large as we wsh copared to n A n n / L ax when the sze of the syste s ncreased, a j s gven to leadng order n the sze of the syste by, j for j>. Use of euatons (15) and (16), enables to show that the neualty relatng n and ax can be ade to be obeyed for large enough systes. It s easy to show now that the reured neualty relatng ax and n s obeyed f / b d. (5) Ths choce of the range of ensures that for j n the range,, j s onotoncally decreasng as a functon of j, for large enough L. It also ensures that,,, and H, (6) H, n where relatve correctons are at ost of the order of n A n ( ) ax / n L

13 (In fact, t ay see that a weaker condton on, d 1 / / b saller than ), suffces, to obtan e. (6) but snce the onotoncty of ( s always, j as a functon of ( j ) splfes atters consderably, the soewhat ore strngent ltaton s used.) Ths ples that to leadng order n L, the expectaton values of H n the orthogonal set are the sae as the expectaton values n the any exctaton states. Note that n the proof of the exstence of the orthogonal set that s close to the non orthogonal set of any exctatons I used the nput scalng for of the any pont correlaton functons (e.(17)). The fnal stage of the proof s to show that the accuulaton of the egenvalues of the Haltonan follows. Assue that ths s not the case. If so snce the egenvalues of H are postve, there exsts a fnte gap,, ndependent on the sze of the syste (for large enough L) and there are only a fnte nuber M (ndependent on the sze of the syste) of egenstates of H wth energes n the gap that tend to zero wth the sze of the syste. Denote these noralzed egenstates by and the correspondng egenvalue b and any noralzed vector vector belongng to the subspace spanned by egenstates of the Haltonan wth energes larger than the gap by {. Consder next the states } wth expectaton values of the Haltonan saller than (, L).I wll express each of these states by bj j b. (7) The expectaton value correspondng to s bounded by j bj b (8) On the other hand (, L). The concluson s that b (, L). Therefore, the projectons of the states on the M densonal space spanned by the s are alost noralzed and alost orthogonal to each other. The last result s possble

14 of course n vew of the enorous nuber of sates wth expectaton values less than (, L). Ths copletes the proof. The proof gven here concerns only the accuulaton of egenvalues of the FP operator at zero n the therodynac lt. Ths result s exact and general under the conjecture that the establshed results for eulbru correlaton functons at the transton are correct. It s thus a proof of the exstence of a necessary condton for slower than exponental decay at the transton pont. It has no bearng however on the actual for of decay whether exponental, stretched exponental or power law. The for of the decay depends on certan atrx eleents that are proble specfc. I hope to dscuss ths n future work. An nterestng ueston that has been rased a nuber of tes n connecton of the above descrpton s whether t s relevant to the glass transton. The frst tendency s to say that t s rrelevant, because n the glassy state translatonal syetry s broken and that was used heavly n the above. Furtherore, the glassy phase s supposed not to be n eulbru and the above treatent s based on eulbru dynacs. It sees, however, that there s soe chance of success by consderng the lud phase close to the transton. I hope to coe back to ths n the near future.

15 References [1] Van Hove L., 1954, Phys. Rev. 95, [] Landau L.D. and Khalatnkov I.M., 1954, Dokl. Akad. Novuk. SSSR 96, 469. [3] Palevsk H. and Hughes D.J., 1953, Phys. Rev. 9,Phys.. [4] Jacrot B. et. al.,1963, Inelastc Scatterng of Neutrons n Solds and Luds (Internatonal Atoc Energy Agency, Venna 1963) p [5] Powell L. et. al., 1963, Phys. Rev. 139 A1866. [6] Baratz M., 1968, Phys. Rev. 17, 4. [7] Cuns H.Z., 1971, Crtcal Phenoena LI Corso Acadec Press NY. [8] Stauffer D.,,1997, Inter. Jour. of Molecular Phys. C8, 163. [9] Stauffer D., 1997, Physca A44, 344. [1] Landau D.P., 1999, Jour. of Magn. and Magnetc Matter., 51. [11] Landau D.P., 1, Tsa S-H. and Bunker A., Materals 6-5, 55. [1] Halpern B.I. and Hohenberg P.C., 1969, Phys. Rev. 177, 95. [13] Halpern B.I., Hohenberg P.C. and Ma S.K., 1974, Phys. Rev. B1, 139. [14] Kawasak K.,1976, n Phase Transtons and Crtcal Phenoena Vol. eds. C. Dob and M.S. Green (Acadec Press, NY 1976). [15] Hohenberg P.C. and Halpern B.I., 197,Rev. Mod. Phys 49, 435 [16] Schneder T., et al.,1984, Phys. Rev. Lett, 53, 191. [17] Schneder T., Zanett M. and Bad R., 1985, Phys. Rev. B 31, 941. [18] Schneder T. and Schwartz M., 1985, Phys. Rev. B RC 31, [19] Calaor F. and Moore M.A. 1, Phys. Rev. E 63, [] Schwartz M. and Edwards S.F., 1, Europhys. Lett. 56, 499. [1] Edwards S.F. and Schwartz M.,, Physca A 33, 357. [] Cardy J. L. (Edtor),1988, n Current Physcs Sources and Coents Vol. Fnte Sze Scalng (North-Holland). [3] Prvan V. (Edtor),199,Fnte Sze Scalng and Nuercal Sulatons of Statstcal Systes( World Scentfc).