Rule 146 average maximum run length vs time step for width 503

max run length for each system width and timestep
near peak around width ~500
noticed a gap of smaller run lengths around width ~503
what is this gap?

rule-146-monolith-width-vs-timestep-width-503-periodic_1.gif

zoom in on width 503

rule-146-monolith-width-vs-timestep-width-503-periodic_2.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_3.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_4.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_5.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_6.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_7.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_8.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_9.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_10.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_11.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_12.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_13.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_14.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_15.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_16.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_17.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_18.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_19.gif

Wow what’s happening here???:

rule-146-monolith-width-vs-timestep-width-503-periodic_20.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_21.gif

^^^^ Max run length falls short of system size for width 503, which indicated monoliths.  Widths 502 and 504 have complete or near-complete anihilations (run of 500 cells) after ~ 500 steps:

rule-146-monolith-width-vs-timestep-width-503-periodic_22.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_23.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_24.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_25.gif

Remember that the max run lengths for each time step are averages over 10 runs or so, so this periodic effect is a property of the average maximum run length for width 503.

evolutions for width 503

rule-146-monolith-width-vs-timestep-width-503-periodic_26.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_27.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_28.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_29.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_30.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_31.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_32.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_33.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_34.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_35.gif

rule-146-monolith-width-vs-timestep-width-503-periodic_36.gif

Lots of monoliths here.

But what is causing periodicity in the average max run length?

Controlling particle trajectories in rule 146

Can the trajectory of a particle in rule 146 be controlled?

Simplest case: 1 particle in the system, try to make it stay at the same position for all time.

Is there an initial condition that would achieve that?  Perhaps a perfectly symmetric one?

In an evolution where the particle stays at one position for all time, does the particle size distribution change vs an evolution where the particle meanders?  i.e. the run lengths of even runs, does the distribution of even run lengths change?

Is there some trade-off between the size of fluctuations in particle size and the deviation from a random walk?  I mean, if you try to control the particle trajectory, does it require messing with the size of the particle much?  I’m thinking of a Heisenberg uncertainty kind of thing, where localizing the particle in space means the momentum becomes uncertain.  Not sure what the right analogy is between momentum and particle size … or if there’s any anology there at all.  Maybe particle size is literally like the mass?

question-rule146-particle-control_1.gif

question-rule146-particle-control_2.gif

question-rule146-particle-control_3.gif

question-rule146-particle-control_4.gif

question-rule146-particle-control_5.gif

question-rule146-particle-control_6.gif

question-rule146-particle-control_7.gif

question-rule146-particle-control_8.gif

question-rule146-particle-control_9.gif

question-rule146-particle-control_10.gif

question-rule146-particle-control_11.gif

question-rule146-particle-control_12.gif

Hmmm, I thik it’s impossible to have a perfectly symmetric initial condition that has an odd run on the boundary, such that no particle appears there.

question-rule146-particle-control_13.gif

question-rule146-particle-control_14.gif

question-rule146-particle-control_15.gif

question-rule146-particle-control_16.gif

So I have to put the single particle on the boundary intentionally, and then RotateRight …

question-rule146-particle-control_17.gif

question-rule146-particle-control_18.gif

question-rule146-particle-control_19.gif

question-rule146-particle-control_20.gif

question-rule146-particle-control_21.gif

question-rule146-particle-control_22.gif

Right that works, but I happened to choose a width close to 2^n so there was a big annihilation:

question-rule146-particle-control_23.gif

question-rule146-particle-control_24.gif

Choose a different width:

question-rule146-particle-control_25.gif

question-rule146-particle-control_26.gif

question-rule146-particle-control_27.gif

question-rule146-particle-control_28.gif

question-rule146-particle-control_29.gif

question-rule146-particle-control_30.gif

Generalize to any width …
Symmetric random initial with a length of about n:

question-rule146-particle-control_31.gif

question-rule146-particle-control_32.gif

question-rule146-particle-control_33.gif

question-rule146-particle-control_34.gif

question-rule146-particle-control_35.gif

question-rule146-particle-control_36.gif question-rule146-particle-control_37.gif question-rule146-particle-control_38.gif question-rule146-particle-control_39.gif question-rule146-particle-control_40.gif
question-rule146-particle-control_41.gif question-rule146-particle-control_42.gif question-rule146-particle-control_43.gif question-rule146-particle-control_44.gif question-rule146-particle-control_45.gif
question-rule146-particle-control_46.gif question-rule146-particle-control_47.gif question-rule146-particle-control_48.gif question-rule146-particle-control_49.gif question-rule146-particle-control_50.gif
question-rule146-particle-control_51.gif question-rule146-particle-control_52.gif question-rule146-particle-control_53.gif question-rule146-particle-control_54.gif question-rule146-particle-control_55.gif
question-rule146-particle-control_56.gif question-rule146-particle-control_57.gif question-rule146-particle-control_58.gif question-rule146-particle-control_59.gif question-rule146-particle-control_60.gif
question-rule146-particle-control_61.gif question-rule146-particle-control_62.gif question-rule146-particle-control_63.gif question-rule146-particle-control_64.gif question-rule146-particle-control_65.gif
question-rule146-particle-control_66.gif question-rule146-particle-control_67.gif question-rule146-particle-control_68.gif question-rule146-particle-control_69.gif question-rule146-particle-control_70.gif
question-rule146-particle-control_71.gif question-rule146-particle-control_72.gif question-rule146-particle-control_73.gif question-rule146-particle-control_74.gif question-rule146-particle-control_75.gif
question-rule146-particle-control_76.gif question-rule146-particle-control_77.gif question-rule146-particle-control_78.gif question-rule146-particle-control_79.gif question-rule146-particle-control_80.gif
question-rule146-particle-control_81.gif question-rule146-particle-control_82.gif question-rule146-particle-control_83.gif question-rule146-particle-control_84.gif question-rule146-particle-control_85.gif
question-rule146-particle-control_86.gif question-rule146-particle-control_87.gif question-rule146-particle-control_88.gif question-rule146-particle-control_89.gif question-rule146-particle-control_90.gif
question-rule146-particle-control_91.gif question-rule146-particle-control_92.gif question-rule146-particle-control_93.gif question-rule146-particle-control_94.gif question-rule146-particle-control_95.gif
question-rule146-particle-control_96.gif question-rule146-particle-control_97.gif question-rule146-particle-control_98.gif question-rule146-particle-control_99.gif question-rule146-particle-control_100.gif
question-rule146-particle-control_101.gif question-rule146-particle-control_102.gif question-rule146-particle-control_103.gif question-rule146-particle-control_104.gif question-rule146-particle-control_105.gif
question-rule146-particle-control_106.gif question-rule146-particle-control_107.gif question-rule146-particle-control_108.gif question-rule146-particle-control_109.gif question-rule146-particle-control_110.gif
question-rule146-particle-control_111.gif question-rule146-particle-control_112.gif question-rule146-particle-control_113.gif question-rule146-particle-control_114.gif question-rule146-particle-control_115.gif
question-rule146-particle-control_116.gif question-rule146-particle-control_117.gif question-rule146-particle-control_118.gif question-rule146-particle-control_119.gif question-rule146-particle-control_120.gif
question-rule146-particle-control_121.gif question-rule146-particle-control_122.gif question-rule146-particle-control_123.gif question-rule146-particle-control_124.gif question-rule146-particle-control_125.gif
question-rule146-particle-control_126.gif question-rule146-particle-control_127.gif question-rule146-particle-control_128.gif question-rule146-particle-control_129.gif question-rule146-particle-control_130.gif
question-rule146-particle-control_131.gif question-rule146-particle-control_132.gif question-rule146-particle-control_133.gif question-rule146-particle-control_134.gif question-rule146-particle-control_135.gif

ECAM rules: dependence of final density on initial density

ECAM mean field density maps

From my thesis work many moons ago …

I made these really cool map of the number of iterations it took for mean field approximations of the density of ECAM rules to converge …

ECAM-T-density-vs-initial-density-110611_1.gif

ECAMProbs[] returns the probability vector for an ECAM, which is a modification of that for an ECA, where the centre cell in each 3-cell neighborhood has probability p[t-δ].

ECAM-T-density-vs-initial-density-110611_2.gif

ECAM-T-density-vs-initial-density-110611_3.gif

ECAM-T-density-vs-initial-density-110611_4.gif

ECAM-T-density-vs-initial-density-110611_5.gif

ECAM-T-density-vs-initial-density-110611_6.gif

ECAM-T-density-vs-initial-density-110611_7.gif

ECAM-T-density-vs-initial-density-110611_8.gif

ECAM-T-density-vs-initial-density-110611_9.gif

ECAM-T-density-vs-initial-density-110611_10.gif

ECAM-T-density-vs-initial-density-110611_11.gif

ECAM-T-density-vs-initial-density-110611_12.gif

ECAM-T-density-vs-initial-density-110611_13.gif

ECAM-T-density-vs-initial-density-110611_14.gif

ECAM-T-density-vs-initial-density-110611_15.gif

ECAM-T-density-vs-initial-density-110611_16.gif

ECAM-T-density-vs-initial-density-110611_17.gif

ECAM-T-density-vs-initial-density-110611_18.gif

Number of iterations required for density convergence of rule 30 as a function of the initial densities p1 and p2 on the x and y axis, respectively.  p1 and p2 are the densities of the two rows comprising the initial condition for the ECAM-1 rule.

ECAM-T-density-vs-initial-density-110611_19.gif

ECAM-T-density-vs-initial-density-110611_20.gif

ECAM-T-density-vs-initial-density-110611_21.gif

ECAM-T-density-vs-initial-density-110611_22.gif

ECAM-T-density-vs-initial-density-110611_23.gif

ECAM-T-density-vs-initial-density-110611_24.gif

It’s a pretty interesting map!

now do something similar without the mean-field approximation …
actually run the rule
plot the final density after T time steps as a function of the inital density

original notebook from April 2006:

ECAM-T-density-vs-initial-density-110611_25.gif

ECAM-T-density-vs-initial-density-110611_26.gif

ECAM-T-density-vs-initial-density-110611_27.gif

ECAM-T-density-vs-initial-density-110611_28.gif

ECAM-T-density-vs-initial-density-110611_29.gif

ECAM-T-density-vs-initial-density-110611_30.gif

ECAM-T-density-vs-initial-density-110611_31.gif

ECAM-T-density-vs-initial-density-110611_32.gif

ECAM-T-density-vs-initial-density-110611_33.gif

ECAM-T-density-vs-initial-density-110611_34.gif

MapECAMDensity[rl_,T_,n_,initdensity_] returns the density of T steps of evolution of rule rl ECAM-T, with n spatial sites and memory T.

ECAM-T-density-vs-initial-density-110611_35.gif

ECAM-T-density-vs-initial-density-110611_36.gif

ECAM-T-density-vs-initial-density-110611_37.gif

Evolution of rule 18 ECAM-20 (i.e. memory T=20) for 20 time steps:

ECAM-T-density-vs-initial-density-110611_38.gif

ECAM-T-density-vs-initial-density-110611_39.gif

Varying the initial density:

ECAM-T-density-vs-initial-density-110611_40.gif

ECAM-T-density-vs-initial-density-110611_41.gif

Density after T time steps as a function of the density of the initial condition:

ECAM-T-density-vs-initial-density-110611_42.gif

Graphics:rule 18 ECAM, T=20

Same for rule 232 ECAM-20:

ECAM-T-density-vs-initial-density-110611_44.gif

ECAM-T-density-vs-initial-density-110611_45.gif

ECAM-T-density-vs-initial-density-110611_46.gif

Graphics:rule 232 ECAM, T=20

… and for rule 146 ECAM-20:

ECAM-T-density-vs-initial-density-110611_48.gif

ECAM-T-density-vs-initial-density-110611_49.gif

ECAM-T-density-vs-initial-density-110611_50.gif

Graphics:rule 146 ECAM, T=20

Take an average of 10 runs for each point (and show the point plus and minus one standard deviation from the mean above and below the average curve):

ECAM-T-density-vs-initial-density-110611_52.gif

Graphics:rule 146 ECAM, T=20

Rule 18 particle: correlation statistics (jellyfish plots)

single-particle correlation data

rule-18-correlations-032711_1.gif

rule-18-correlations-032711_2.gif

rule-18-correlations-032711_3.gif

rule-18-correlations-032711_4.gif

rule-18-correlations-032711_5.gif

rule-18-correlations-032711_6.gif

rule-18-correlations-032711_7.gif

rule-18-correlations-032711_8.gif

rule-18-correlations-032711_9.gif

Single-particle correlations from just one evolution, for width 251:

rule-18-correlations-032711_10.gif

rule-18-correlations-032711_11.gif

statistical sampling of correlations at each correlation distance

system width: 251

rule-18-correlations-032711_12.gif

rule-18-correlations-032711_13.gif

rule-18-correlations-032711_14.gif

rule-18-correlations-032711_15.gif

rule-18-correlations-032711_16.gif

rule-18-correlations-032711_17.gif

rule-18-correlations-032711_18.gif

rule-18-correlations-032711_19.gif

Find a histogram of the correlations at a given correlation distance

rule-18-correlations-032711_20.gif

rule-18-correlations-032711_21.gif

rule-18-correlations-032711_22.gif

rule-18-correlations-032711_23.gif

rule-18-correlations-032711_24.gif

rule-18-correlations-032711_25.gif

rule-18-correlations-032711_26.gif

rule-18-correlations-032711_27.gif

rule-18-correlations-032711_28.gif

rule-18-correlations-032711_29.gif

rule-18-correlations-032711_30.gif

rule-18-correlations-032711_31.gif

rule-18-correlations-032711_32.gif

rule-18-correlations-032711_33.gif

rule-18-correlations-032711_34.gif

rule-18-correlations-032711_35.gif

rule-18-correlations-032711_36.gif

correlation histograms for various correlation distances

rule-18-correlations-032711_37.gif

rule-18-correlations-032711_38.gif

correlation envelopes

Now show the histogram as a function of correlation distance:

rule-18-correlations-032711_39.gif

rule-18-correlations-032711_40.gif

rule-18-correlations-032711_41.gif

rule-18-correlations-032711_42.gif

rule-18-correlations-032711_43.gif

rule-18-correlations-032711_44.gif

rule-18-correlations-032711_45.gif

The boundaries above and below in the plot define an envelope for the correlations at this system width:

rule-18-correlations-032711_46.gif

rule-18-correlations-032711_47.gif

How does the shape of the envelope change as the system width is increased?

width 398

rule-18-correlations-032711_48.gif

rule-18-correlations-032711_49.gif

rule-18-correlations-032711_50.gif

rule-18-correlations-032711_51.gif

rule-18-correlations-032711_52.gif

rule-18-correlations-032711_53.gif

rule-18-correlations-032711_54.gif

width 630

rule-18-correlations-032711_55.gif

rule-18-correlations-032711_56.gif

rule-18-correlations-032711_57.gif

rule-18-correlations-032711_58.gif

rule-18-correlations-032711_59.gif

rule-18-correlations-032711_60.gif

rule-18-correlations-032711_61.gif

width 1000

rule-18-correlations-032711_62.gif

rule-18-correlations-032711_63.gif

rule-18-correlations-032711_64.gif

rule-18-correlations-032711_65.gif

rule-18-correlations-032711_66.gif

rule-18-correlations-032711_67.gif

rule-18-correlations-032711_68.gif

rule-18-correlations-032711_69.gif

rule-18-correlations-032711_70.gif

width 1584

rule-18-correlations-032711_71.gif

rule-18-correlations-032711_72.gif

rule-18-correlations-032711_73.gif

rule-18-correlations-032711_74.gif

rule-18-correlations-032711_75.gif

rule-18-correlations-032711_76.gif

rule-18-correlations-032711_77.gif

rule-18-correlations-032711_78.gif

rule-18-correlations-032711_79.gif

width 2511

rule-18-correlations-032711_80.gif

rule-18-correlations-032711_81.gif

rule-18-correlations-032711_82.gif

rule-18-correlations-032711_83.gif

rule-18-correlations-032711_84.gif

rule-18-correlations-032711_85.gif

rule-18-correlations-032711_86.gif

rule-18-correlations-032711_87.gif

rule-18-correlations-032711_88.gif

width 3981

rule-18-correlations-032711_89.gif

rule-18-correlations-032711_90.gif

rule-18-correlations-032711_91.gif

rule-18-correlations-032711_92.gif

rule-18-correlations-032711_93.gif

rule-18-correlations-032711_94.gif

rule-18-correlations-032711_95.gif

rule-18-correlations-032711_96.gif

rule-18-correlations-032711_97.gif

width 6309

rule-18-correlations-032711_98.gif

rule-18-correlations-032711_99.gif

rule-18-correlations-032711_100.gif

rule-18-correlations-032711_101.gif

rule-18-correlations-032711_102.gif

rule-18-correlations-032711_103.gif

rule-18-correlations-032711_104.gif

rule-18-correlations-032711_105.gif

rule-18-correlations-032711_106.gif

width 10000

rule-18-correlations-032711_107.gif

rule-18-correlations-032711_108.gif

rule-18-correlations-032711_109.gif

rule-18-correlations-032711_110.gif

rule-18-correlations-032711_111.gif

rule-18-correlations-032711_112.gif

rule-18-correlations-032711_113.gif

These plots totally look like jellyfish.

rule-18-correlations-032711_114.gif

rule-18-correlations-032711_115.gif

next: plot the raw particle trajectories in an interesting way …

rule-18-correlations-032711_116.gif

rule-18-correlations-032711_117.gif

rule-18-correlations-032711_118.gif

rule-18-correlations-032711_119.gif

rule-18-correlations-032711_120.gif

rule-18-correlations-032711_121.gif

rule-18-correlations-032711_122.gif

rule-18-correlations-032711_123.gif

Rule 18: correlations vs system size

autocorrelations in rule 18 particle position as a function of system size

It’s time to look at the autocorrelations in the rule 18 single-particle movement as a function of system size, for a larger range of system sizes (several orders of magnitude), with some statistics for each width.

Here’s the simulation code …

functions

rule-18-corr-vs-system-size-031211_1.gif

rule-18-corr-vs-system-size-031211_2.gif

rule-18-corr-vs-system-size-031211_3.gif

rule-18-corr-vs-system-size-031211_4.gif

rule-18-corr-vs-system-size-031211_5.gif

rule-18-corr-vs-system-size-031211_6.gif

rule-18-corr-vs-system-size-031211_7.gif

rule-18-corr-vs-system-size-031211_8.gif

simulation code:
correlation vs distance as a function of system size

rule-18-corr-vs-system-size-031211_9.gif

rule-18-corr-vs-system-size-031211_10.gif

rule-18-corr-vs-system-size-031211_11.gif

rule-18-corr-vs-system-size-031211_12.gif

data analysis

Now grab the resulting data files:

rule-18-corr-vs-system-size-031211_13.gif

rule-18-corr-vs-system-size-031211_14.gif

rule-18-corr-vs-system-size-031211_15.gif

Each index here is a different width:

rule-18-corr-vs-system-size-031211_16.gif

rule-18-corr-vs-system-size-031211_17.gif

rule-18-corr-vs-system-size-031211_18.gif

rule-18-corr-vs-system-size-031211_19.gif

I chose widths logarithmically spaced spanning 3 orders of magnitude, in the range 100 to 10,000:

rule-18-corr-vs-system-size-031211_20.gif

rule-18-corr-vs-system-size-031211_21.gif

These widths will be evenly spaced on a log plot.

Take a look at a single run at width 100:

rule-18-corr-vs-system-size-031211_22.gif

rule-18-corr-vs-system-size-031211_23.gif

rule-18-corr-vs-system-size-031211_24.gif

rule-18-corr-vs-system-size-031211_25.gif

rule-18-corr-vs-system-size-031211_26.gif

rule-18-corr-vs-system-size-031211_27.gif

Here’s the particle position vs. time, run for 50k steps:

rule-18-corr-vs-system-size-031211_28.gif

rule-18-corr-vs-system-size-031211_29.gif

The particle position went periodic.  Not surprising for such a small system width.  Here’s the correlation function, which also has a strong periodicity:

rule-18-corr-vs-system-size-031211_30.gif

rule-18-corr-vs-system-size-031211_31.gif

Let’s look at all of the runs for width 100:

rule-18-corr-vs-system-size-031211_32.gif

rule-18-corr-vs-system-size-031211_33.gif

rule-18-corr-vs-system-size-031211_34.gif

(Note that my particleposition code treats periodic boundary crossings as reflections.)

And all the correlations at this width:

rule-18-corr-vs-system-size-031211_35.gif

rule-18-corr-vs-system-size-031211_36.gif

rule-18-corr-vs-system-size-031211_37.gif

The correlations have a pretty regular structure, since all the evolutions seemed to have gone periodic within about 10k steps.  Given this reglar structure, it seems worthwhile to take the average autocorrelation at each distance:

rule-18-corr-vs-system-size-031211_38.gif

rule-18-corr-vs-system-size-031211_39.gif

rule-18-corr-vs-system-size-031211_40.gif

rule-18-corr-vs-system-size-031211_41.gif

rule-18-corr-vs-system-size-031211_42.gif

rule-18-corr-vs-system-size-031211_43.gif

rule-18-corr-vs-system-size-031211_44.gif

Do this for all the other system widths:

rule-18-corr-vs-system-size-031211_45.gif

rule-18-corr-vs-system-size-031211_46.gif

rule-18-corr-vs-system-size-031211_47.gif

rule-18-corr-vs-system-size-031211_48.gif

rule-18-corr-vs-system-size-031211_49.gif

rule-18-corr-vs-system-size-031211_50.gif

rule-18-corr-vs-system-size-031211_51.gif

rule-18-corr-vs-system-size-031211_52.gif

Graphics:width 100
Graphics:width 158
Graphics:width 251
Graphics:width 398
Graphics:width 630
Graphics:width 1000
Graphics:width 1584
Graphics:width 2511
Graphics:width 3981
Graphics:width 6309
Graphics:width 10000

The correlation functions at each width have a prett similar structure.  It isn’t clear if I’m resolving any real differences it the correlations here at all actually.  It could be that the particle really isn’t much effected by the system size.

Vitaliy recommended that I look at the power spectrum, which is related to these correlation functions via the Weiner-Kinchin theorem.  That’s the next thing to do I think.

Rule 18: long-range particle correlations

I ran a few batch jobs with longer evolutions of rule 18 (100k steps, used to be 10k steps), and measuring much longer-range correlations than before (out to distances of 50k and 100k correlation distances).

My first task was to speed up the correlation distance function as much as possible.

speeding up the correlation function

Generate a random walk:

long-range-correlations-rule-18-walk_1.gif

long-range-correlations-rule-18-walk_2.gif

long-range-correlations-rule-18-walk_3.gif

My original version of the correlation function, which I want to speed up:

long-range-correlations-rule-18-walk_4.gif

Timing for a single correlation data point (lag 1):

long-range-correlations-rule-18-walk_5.gif

long-range-correlations-rule-18-walk_6.gif

I want to do this without using Partition, since I think it’s taking a lot of memory.

Try using RotateRight instead … already an order of magnitude faster:

long-range-correlations-rule-18-walk_7.gif

long-range-correlations-rule-18-walk_8.gif

long-range-correlations-rule-18-walk_9.gif

long-range-correlations-rule-18-walk_10.gif

long-range-correlations-rule-18-walk_11.gif

long-range-correlations-rule-18-walk_12.gif

A similar version which doesn’t wrap at the boundaries.  I don’t want to wrap around the boundaries in *time*, because there’s no reason to correlate the beginning of the timeseries with the end of the timeseries.  Normalization uses a few less points of the list, so the correlation comes out numerically slightly different:

long-range-correlations-rule-18-walk_13.gif

long-range-correlations-rule-18-walk_14.gif

… but it’s just as fast as the RotateRight method, and I think it’s more accurate for my purposes since it doesn’t do temporal wrapping.

So here’s my new correlation function:

long-range-correlations-rule-18-walk_15.gif

long-range-correlations-rule-18-walk_16.gif

long-range-correlations-rule-18-walk_17.gif

measuring long-range correlations in the rule 18 particle movement

Grab the particle correlation data for a rule 18 evolution of 100k timesteps (still using a width of 199 cells):

long-range-correlations-rule-18-walk_18.gif

And grab the actual position vs. time for the particle:

long-range-correlations-rule-18-walk_19.gif

long-range-correlations-rule-18-walk_20.gif

long-range-correlations-rule-18-walk_21.gif

long-range-correlations-rule-18-walk_22.gif

long-range-correlations-rule-18-walk_23.gif

Got about 50k correlation data points:

long-range-correlations-rule-18-walk_24.gif

long-range-correlations-rule-18-walk_25.gif

long-range-correlations-rule-18-walk_26.gif

long-range-correlations-rule-18-walk_27.gif

Here’s the beginning of the walk:

long-range-correlations-rule-18-walk_28.gif

long-range-correlations-rule-18-walk_29.gif

But it clearly went periodic:

long-range-correlations-rule-18-walk_30.gif

long-range-correlations-rule-18-walk_31.gif

And the correlations also went periodic:

long-range-correlations-rule-18-walk_32.gif

long-range-correlations-rule-18-walk_33.gif

So I did another evolution, this time using width of 999 cells:

long-range-correlations-rule-18-walk_34.gif

long-range-correlations-rule-18-walk_35.gif

long-range-correlations-rule-18-walk_36.gif

long-range-correlations-rule-18-walk_37.gif

long-range-correlations-rule-18-walk_38.gif

long-range-correlations-rule-18-walk_39.gif

long-range-correlations-rule-18-walk_40.gif

long-range-correlations-rule-18-walk_41.gif

Now the walk doesn’t go periodic, because the underlying evolution didn’t go periodic:

long-range-correlations-rule-18-walk_42.gif

long-range-correlations-rule-18-walk_43.gif

And behold the long-range correlations!

long-range-correlations-rule-18-walk_44.gif

long-range-correlations-rule-18-walk_45.gif

This is an interesting plot.  Correlations fall gradually here over about 50k lag, then hover around a negative correlation of around -0.4 from 50k to 100k lag.

Next steps:

1. Get a statistical distribution of correlations at each lag, from multiple runs of the rule 18 particle

Rule 18 single-particle correlations: larger correlation distances

Rule 18: Eating its own light cone (Part 5)

At the last ruleshack meetup, we discussed why the autocorrelation curve for the position of the rule 18 particle appears to be linear.  And we realized that I’m probably just not looking at large enough correlation distances (lags).

The particle position vs time looks like this:

rule-18-random-walk-022611_1.gif

rule-18-random-walk-022611_2.gif

rule-18-random-walk-022611_3.gif

rule-18-random-walk-022611_4.gif

rule-18-random-walk-022611_5.gif

rule-18-random-walk-022611_6.gif

And the autocorrelations out to lag 100 look like this:

rule-18-random-walk-022611_7.gif

rule-18-random-walk-022611_8.gif

rule-18-random-walk-022611_9.gif

rule-18-random-walk-022611_10.gif

During the meetup, we looked at the particle position with the mean subtracted (since that’s what the autocorrelation function is based on):

rule-18-random-walk-022611_11.gif

rule-18-random-walk-022611_12.gif

We noticed that the typical time between zero crossings (when the curve crosses the x-axis in this plot, or equivalently the instants in time when the particle passes through its average position), is greater than the maximum lag I was using of 100 time steps.  In fact, looking at the mean-subtracted position plot above, sometimes it’s more like 1000 time steps between crossings.  That means I need to sample much larger lags in order to capture more features of the random walk.

To be more precise, I can compute the crossing intervals, and the distribution of interval sizes:

rule-18-random-walk-022611_13.gif

rule-18-random-walk-022611_14.gif

rule-18-random-walk-022611_15.gif

rule-18-random-walk-022611_16.gif

rule-18-random-walk-022611_17.gif

rule-18-random-walk-022611_18.gif

rule-18-random-walk-022611_19.gif

The average interval size:

rule-18-random-walk-022611_20.gif

rule-18-random-walk-022611_21.gif

rule-18-random-walk-022611_22.gif

rule-18-random-walk-022611_23.gif

The interval distribution:

rule-18-random-walk-022611_24.gif

rule-18-random-walk-022611_25.gif

rule-18-random-walk-022611_26.gif

rule-18-random-walk-022611_27.gif

It looks like the average crossing interval size is of the order of a few hundred time steps, with a intervals of order 10 being most common, but with some rare intervals being on the order of 1000 time steps.  

Call the average crossing interval rule-18-random-walk-022611_28.gif.  Then the products p(t)*p(t+r) will tend to be positive for rule-18-random-walk-022611_29.gif.  But for lags rule-18-random-walk-022611_30.gif, there will be more products where p(t) and p(t+r) have different sign, and thus contributing negative terms to the sum over t.  So I’d expect the autocorrelations to start becoming negative somewhere around rule-18-random-walk-022611_31.gif.

Here’s what the correlations look like out to lag 2000 (which I ran on a server).  The curve is looking more interesting now, with the linearity going away at around lag 1000.

rule-18-random-walk-022611_32.gif

rule-18-random-walk-022611_33.gif

rule-18-random-walk-022611_34.gif

rule-18-random-walk-022611_35.gif