Simulation of majority rule disturbed by power-law noise on directed and undirected Barabas

  • 格式:pdf
  • 大小:185.36 KB
  • 文档页数:6

a

r

X

i

v

:

0802.1476v1 [cond-mat.dis-nn] 11 Feb 2008Simulationofmajorityruledisturbedbypower-lawnoiseon

directedandundirectedBarab´asi-Albertnetworks

F.W.S.Lima

DepartamentodeF´ısica,UniversidadeFederaldoPiau´ı,64049-550,Teresina

-PI,Brazil

e-mail:wel@ufpi.br

Abstract:OndirectedandundirectedBarab´asi-AlbertnetworkstheIsingmodel

withspinS=1/2inthepresenceofakindofnoiseisnowstudiedthroughMonte

Carlosimulations.ThenoisespectrumP(n)followsapowerlaw,whereP(n)is

theprobabilityofflippingrandomlyselectnspinsateachtimestep.Thenoise

spectrumP(n)isintroducedtomimictheself-organizedcriticalityasamodel

influenceofacomplexenvironment.Inthismodel,differentfromthesquarelattice,

theorder-disorderphasetransitionoftheorderparameternotisobserved.For

directedBarab´asi-Albertnetworksthemagnetisationtendstozeroexponentially

andundirectedBarab´asi-Albertnetworks,itsremainconstant.Keywords:MonteCarlosimulation,spins,networks,Ising.

Introduction

ThispaperdealswithIsingspinondirectedandundirectedBarab´asi-

Albert(BA)networksinthepresenceofanoise.SumourandShabat[1,2]

investigatedIsingmodelswithspinS=1/2ondirectedBAnetworks[3]with

theusualGlauberdynamics.Nospontaneousmagnetisationwasfound,in

contrasttothecaseofundirectedBAnetworks[4,5,6]whereasponta-

neousmagnetisationwasfoundbelowacriticaltemperaturewhichincreases

logarithmicallywithsystemsize.LimaandStauffer[8]simulateddirected

square,cubicandhypercubiclatticesintwotofivedimensionswithheat

bathdynamicsinordertoseparatethenetworkeffectsofdirectedness.They

alsocompareddifferentspinflipalgorithms,includingclusterflips[9],for

Ising-BAnetworks.Theyfoundafreezing-inofthemagnetisationsimilarto

[1,2],followinganArrheniuslawatleastinlowdimensions.Thislackof

aspontaneousmagnetisation(intheusualsense)isconsistentwiththefact

thatifonadirectedlatticeaspinSjinfluencesspinSi,thenspinSiinturn

doesnotinfluenceSj,andtheremaybenowell-definedtotalenergy.Thus,

theyshowthatforthesamescale-freenetworks,differentalgorithmsgive

1differentresults.RecentlyStaufferandKu󰁍lakowski[10]simulatedtheIsing

two-dimensionalferromagnetinthepresenceofaspecialkindofnoisedif-

ferentfromthetraditionaltemperatureandBoltzmannprobabilities,where

thenoiseconsistofateachiterationofaL×Lsquarelatticewithfour

neighboursforeachsites,besidestheabovemajorityruletheyselectntimes

randomlyaspinandflipit.TheprobabilitydistributionfunctionP(n)of

thesenumbersnistakenasapowerlaw,

P(n)∝1/nα.(1)

02004006008001000time0.0107.0107.3107.5107.6107.7

magnetisationRelaxation, 500000 sites, 100 samples summed

Figure1:Summedmagnetisationsversustime,T=1,α=1,andm=2for

directed(BA)networks.

Inordertogetthisdistribution,theydeterminedrandomnumbersr,

homogeneouslydistributedbetweenzeroandone,andthentook

n=TL2r1/(1−α),(2)

2

02004006008001000time−107.30.0107.3107.6107.8

magnetisations100 samples, 500000,T: 0.001, 0.01, 0.05, 0.1, 0.5, 1.0 from top

02004006008001000time−107.30.0107.3107.6107.8

magnetization100 samples 500000, T: 0.1, 0.2, 0.3, 0.5, 1. from top

Figure2:Summedmagnetisationsversustimefor0.01≤T≤1.0,α=1

andm=2(left)and0.1≤T≤1.0,α=1andm=7(right)fordirected

(BA)networks.

forα<1,where1≤n≤L2and

n=Texp(rln(L2))(3)

forα=1andTdeterminestheamplitudeofthenoise.Herewesimu-

latedthesamemodeldescribedaboveondirectedandundirected(BA)net-

worksandourresultsaredifferentfromtheresultsobtainedbyStaufferand

Ku󰁍lakowski[10].

Results

InFig.1weshowthedependenceofthesummedmagnetisationM=󰀁iSiversustimeforα=1,T=0.01(eq.3)andconnectionnumberm=2

ondirected(BA)networkswith500000sites,whenwestartedwithallspins

up.Mrelaxesexponentiallytowardszero.InFig.2weplotsamepicture

ofFig.1,butforseveralvaluesofTandconnectionnumbersm=2(left)

andm=7(right),whereforsmallerTthemagnetisationtriestoremain

positive.IntheFig.3weplotthedependenceofthesummedmagnetisation

versustimeforα=1,severalvaluesofT(left),andforT=10(right)even

upto10000MonteCarlostepswithm=2onundirected(BA)networks

with10000sites.Thesummedmagnetisationremainspositiveanddoes

notdecaytowardszero.Itkeepsfluctuatingaroundapositivevalueofthe

magnetisationdependingonthevalueofT.

3

02004006008001000time0.00105.00105.30105.48105.60105.70105.78105.85105.90105.95

magnetisation100 samples 10000, T: 1, 10, 100, 1000 from top

0200040006000800010000time0.00105.00105.30105.48105.60105.70105.78105.85105.90105.95106.00

magnetisation100 samples 10000, T=10, m=2

Figure3:Summedmagnetisationsversustimefor1≤T≤1000.,α=1

andm=2(left)andforT=10,α=1andm=2(right)forundirected

(BA)networks.

InFig.4themagnetisationisshownasafunctionoftimeforonesample

forT=1,α=1andm=2ondirected(left)andundirected(right)(BA)

networks.Fordirected(BA)networksthesummedmagnetisationchanges

betweenpositiveandnegativevalueswithtime,andforundirected(BA)net-

worksthemagnetisationfluctuatesaroundapositivevalue.InFig.5we