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.RecentlyStaufferandKulakowski[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
Kulakowski[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