Cage 3 crystal

A visualization of how good of a model the minimum enclosing box is when applied to a real crystal structure.

Alpha shape applied to the Cage 3 molecule in the minimum enclosing box

Alpha shapes in a Cage 3 crystal

The alpha shapes do not overlap.  For some alpha.

Convex hulls in a Cage 3 crystal

The convex hulls of Cage 3 overlap.

Enclosing boxes in a Cage 3 crystal

 

First simulation

I started work with only with one enclosing box and then to continue to build upon that.

I took one enclosing box and constructed a unit cell by surrounding this central box by other boxes (as is visualized in previous post). Next I “optimized” lattice energy by Lennard-Jones potential using the centers of mass of the CC3 molecule. The central cube can by moved in 3-directions (x,y,z) but out of symmetry it’s sufficient to move only in direction in this particular case. The center of mass of CC3 is slightly off the center of mass of the small cube. That means there are 16 rotations of the central cube that can be potential minima of the energy function.

The result can by seen by below. The dots are centers of mass of CC3. The minimum was attained with no rotation of the central small box.

Reducing output matrix according to mean TE

The procedure for choosing output weights. First the training reservoir signals were collected and TE matrix was computed based on these signals. Next the mean TE of target unit was computed and only those units where chosen for output matrix computations which mean TE was higher then some Threshold. The threshold went from the minimum mean target TE to the maximum mean target TE in 100 steps. This was repeated 100 times. The network settings were same as in the IJCNN article for comparability reasons.

 

  

 

80: Sigma, Tau & Rho exploration

I have added spectral radius scaling of reservoir matrix to desired value \rho along \sigma and \tau parameters in W \sim N(0,\sigma^2) and W^{in} \sim Unif(-\tau,\tau) to the exploration. 100 neuron reservoir was used and 100 instances for every pair (\tau, \sigma) were generated.

The Z values are mean NRMSE from 100 instances.

 

NARMA

NARMA without unscaled

Mackey – Glass

 

Mackey-Glass unscaled

Lorenz

Lorenz unscaled

Minimum achieved mean:
Rendered by QuickLaTeX.com