Black rings

Assuming black holes exist, it is true that every one of them must spin, since the chance of one forming with exactly 0 angular velocity is 0. Black holes are thought to create a singularity in their centre, which is hidden behind their Schwartzchild radius, or event horizon. It seems to me (without thinking too hard about the maths) that a spinning black hole should produce an infinitely dense ring rather than a single point. A ring spinning around its axis should hold itself open.

What is interesting about this is that a singularity in the form of a ring is an allowable kind of analytic function of 3d space with two or more layers of volumes. By this I mean that you could pass through the ring and into a different version of 3d space without I think violating general relativity. This is a bit like ideas about space-time wormholes but it is different, it doesn't connect two distance areas of space by a shortcut, it is just a particular geometry of space where everywhere is continuous (apart from the ring singularity) but the world through the ring is different from the world if you pass by the ring without going through it. Moreover, there are potentially an unbounded number of layers each time you circle around the ring and through it.

Would make a nice concept for a sci-fi movie.

I wonder if such a physical setup could be simulated... the ring doesn't have to be massive, it could be the size of a door, so long as the ring is infinitely dense (which doesn't mean massive).

Reaction Diffusion Fractals

Previous work with fractal automata has the disadvantage that it lacks continuous rotational symmetry (it also lacks continuous translational and scale symmetry). An interesting idea is to instead work with reaction-diffusion systems which are already continuous in rotation and translation symmetry, and try to add scale symmetry. The normal formula is as follows (and is explained in the link above):

Scaling Du and Dv scales the size of the patterns formed, it seems that the pattern size is proportional to the square root of this scale. Therefore we choose to scale by t² where t is our scaling factor.

Next we change all of our variables to be vectors rather than scalars. This represents a list of independent reaction diffusion systems. We then choose vector t to be (1,t,t²,....), in other words the difference between each reaction diffusion system is a geometric increase in scale. If we view the first three components of v in the red,green and blue channels respectively, then the result is simply a superimposed set of three reaction diffusion systems (using Du = 1, Dv = 0.5):

Here the red scale is the smallest, the green is noticeably twice the scale, and blue twice again.

The final ingredient is to have these separate scales interact with eachother, we do that in the reaction part of the formula, which is the uv²...

each u is reacting with two vs, we replace the v vector (v for each scale) with a weighted average of all the vs where the weighting is an exponential dropoff s^|x| around each component. Effectively the vector v has been convolved with an exponential dropoff function, or low pass filtered. This gives us an extra parameter s, when s=0 we have no interaction like the image above, s=1 gives equal interaction which prevents any separation of the different components. Values in between are interesting.

Here we choose a small value of s = 0.05. t = 0.5 so we double the scale of each component (red,green,blue). I plot the reaction diffusion system with varying parameters F and k on the horizontal and vertical axes respectively (increasing right and up). 

zoomed out

s=0.05. zoomed in slightly. Varying the parameters in x,y shows the different patterns possible.

Next are with s = 0.1. The three components are more correlated. Notice that the thin bridges are more red (small scale) and the large areas more blue. 

zoom around 0.027, 0.058, range 0.08

zoom around 0.026, 0.057, range 0.04

zoom around 0.023, 0.055, range 0.02

zoom around 0.023, 0.051, range 0. (unvarying in x,y)

similar area, unvarying parameters. Notice the similar worm shapes at each scale.

The system does not require the scales are doubled for each component, here we show a zoomed out image for t=0.7, so green is roughly 1.4 times the width of red.   

In fact the results are an approximation of a continuous scale symmetry as t→1 from below. Make sure to change the dropoff to s^{|x|log₂(1/t)}

Symmetric binary fractals

In a previous post I discussed fractal automata, which are like cellular automata (e.g. Conway's game of life) but include larger and smaller scales as neighbours as well as nearby locations. This is summarised here: https://sites.google.com/site/tomloweprojects/scale-symmetry/automataFinder

A subset of these systems are static fractals, generated by applying binary rules to lower resolution maps.
It is useful to try and find a setup with as much symmetry as possible and with as few degrees of freedom as possible. The default system operates on a 2d grid and so exhibits square symmetry. The 'type 7' fractals displayed in the above link use a semi-cheat to allow for octagonal symmetry. It works very well but has some drawbacks; it cannot be animated and isn't really octagonal in its symmetry, just close. It also doesn't extend to 3d which up to now has just cubic symmetry.

A new idea, which is the topic of this post, is to use extra dimensions to provide extra symmetry. If I use a 3d grid, and apply binary rules to convert low resolution voxels progressively to higher resolution ones, doubling the resolution each step, then take a long-diagonal cross section (the plane defined by x+y+z = 0) then the resulting 2d image has hexagonal symmetry. If I make the low resolution data just a 2x2x2 block of set voxels then the rules result in images such as these:

 

 

Of course, one can easily make a grid with hexagonal symmetry (a triangular lattice), but the problem is that we need to elegantly register each cell to its lower resolution cell on a coarser version of the grid. For triangular lattices the obvious choice is for 4 triangles to refer to one parent triangle, however these triangles have different geometric locations relative to their parent, it can be the centre triangle or one of the three corner triangles, consequently you need different rules for each case and the rule set becomes larger which defeats the goal of minimising degrees of freedom (the number of possible rules).

By contrast, for 3d voxels, each of 8 child voxels occupy an identical corner of the parent 2x2x2 scaled voxel in my implementation. The ruleset I chose for the above images is therefore quite succinct... I base it only on the seven closest parent voxels, which, assuming cubic symmetry, gives 40 possible configurations of these seven voxels. The ruleset provides an on/off of the child voxel for each of these combinations, therefore 2^40 possible rulesets.

I can reduce this large search space down further by enforcing bit symmetry, this specifies that on and off are arbitrary so if a particular low res shape produces a particular high res fractal, then the negative (on and off swapped) low res shape produces a negative of the fractal. This reduces the search space down to 2^20 rulesets.

A last improvement is to acknowledge that a ruleset and its negative are functionally no different if you have bit symmetry, so we can normalise by choosing the version which converts all seven parent voxels being off to an off state for the child. Resulting in 2^19 rulesets.

You can of course also view the 3d shape generated, but in 3d you miss the details behind the surface:

Optimisation

For the 2d images, we only need a single plane of high resolution voxels, and these trace back to a thin slice through the coarser resolution voxels. Therefore the processing required is much less than for viewing the 3d output and is almost definitely O(n^2) for resolution nxn unlike O(n^3) for the 3d case. This also means that the initial seeds are only needed in a thin slice of the 3d voxel grid.

Random seeds

I can also seed randomly, which gives a better overview of how each ruleset looks on arbitrary input:

 

 

 

Extensions

There are clearly variations on the rulesets that I'm using, one could look at a larger neighbourhood, or remove bit symmetry etc, I could also change the scaling factor (if the parent to child scaling is 3 rather than 2 then I am sure one can reconstruct the Koch snowflake for instance), but more degrees of freedom almost always gives less symmetric and less interesting (more noisy) results.

However, this method has several advantages, it probably can be made to animate (fractal automata), and it can be extended to higher dimensions.

3d fractals could be generated which have what is probably octahedral symmetry rather than cubic symmetry, by using a 3d long-diagonal cross section of a 4d binary fractal. 

I am also wondering whether one can make use of the highly symmetric D4 lattice in 4d to create even more symmetric 3d fractals... however I'm not sure how elegantly you can register one 24-cell to a parent larger 24-cell.

Simpler

I'll finish with some examples from a simpler system where I just look at the 4 closest parents. In this case I haven't used bit symmetry. You can sort of see there is more simplicity in the rules. The last image uses a random seed.

 

 

Mary Morris Diaries

This is quite unlike my previous posts but I wanted to keep a record of the events of the publication of my Grandma's diaries.
It is a fascinating story of a rebellious young Irish nurse who is thrown into the deep end as the war starts, having to care for gravely ill children with hopelessly inadequate supplies and in the midst of bombing raids (her sister hospital takes a direct strike). She follows the soldiers into France with the Normandy landings where she tends a field hospital on the frontline. Her patients are axis and allies alike and, in the face of daily tragedy and disease, she finds her strength in her patients with their varied and charming personalities. As a nurse and perhaps being Irish (who were neutral) she brings a culturally un-sided and stark account of the major events of WWII, dispairing at the poor treatment of German prisoners in British run PoW camps, and relaying British foot soldiers' dismay at not getting the air backup they were promised. Through the chaos she meets a man who is to become her husband as the war is ending, she tells of his role to protect the German civilians from Russian soldiers who are conducting revenge attacks.


Links:
Wikipedia: https://en.wikipedia.org/wiki/Mary_Morris_(Diarist)
Family blog: http://marymorrisauthor.blogspot.co.uk/2014/06/her-wartime-diaries.html
Old BBC entry by Grandpa: http://www.bbc.co.uk/history/ww2peopleswar/stories/53/a4106053.shtml
Irish times: http://www.irishtimes.com/culture/heritage/an-irishman-s-diary-on-an-irish-nurse-in-wartime-1.1817784
Daily Express: http://www.express.co.uk/news/uk/481145/The-forbidden-diary-of-a-wartime-nurse
Irish mail on Sunday, June 1st, 2014
Irish times book review: http://www.irishtimes.com/culture/books/discovering-an-irish-nurse-s-unique-wwii-diary-1.1829873?page=3
Connacht tribune: http://connachttribune.ie/caltra-nurses-d-day-diary-bestseller/
Daily mail: http://www.dailymail.co.uk/home/books/article-2680044/Nightingale-Dunkirk.html
Essay by Carol Acton: http://dx.doi.org/10.1080/09670880903533417
Irish times book review: http://www.irishtimes.com/culture/books/irish-nurse-s-account-of-the-war-a-very-private-diary-1.1874768

Operator before addition

I'm sure this idea is not new but, well it's a new idea to me. The question is, what comes before the plus operator on numbers? i.e. if we have + then * then power, then tetration... what do we get when going in the other direction?
These can be defined by the three valued Ackermann function phi(m,n, x) so phi(m,n, 0) = m+n, phi(m,n, 1) = m*n etc. The question may then be what is phi(m, n, -1)?
I'm going to call this operator # (or 'with').
I say 'may' because the extension of these numerical operators is perhaps ambiguous and so it may be that there are different correct answers depending on your interpretation of the existing operators. In a previous post I thought that # is the equality operator, it returns m if m==n, otherwise it returns void.

Today's idea is to make:
 m # n = log_2(2^m + 2^n)
So we get a sequence
a#a = a+1
a+a = a*2
a*a = a^2

The identity value for all four:
a#-inf = a
a+0 = a
a*1 = a
a^1 = a.

The operators are different because they consider the arguments in a different context... for addition the values m and n are considered as points along a line, for multiplication they are considered to be magnitudes, and for the # operator they are considered to be entropies. This is important because entropy is a very fundamental concept so having # as an operator as fundamental as + and * is useful in learning to deal with entropy and information.
For example, the information needed to store either a 16 bit colour or a 16 bit alpha is 16#16 = 17bits. Anyway, # is a sum 'as though the arguments are entropies'. It also has an inverse:
m -# n which is just log_2(2^m - 2^n).

Rolling fractals: Thistles and leaves

An interesting fractal construction comes from rolling circles inside a stack of larger circles. This has been considered here and here... however I don't think these designs are the correct maths for proper conformal rolling circles. The equation I uses gives a different result.

For a set of nested circles of power 2 diameters, rolling with a constant roll contact point speed, in alternating directions I get this:

It is part of a set of fractals which are defined by the curvature ratio from parent to child. Here we see the family for ratios -4,-3,-2,  2, 3, 4.  (-1,0 and 1 are not fractals).

None of these fractals appear to self intersect.

The golden egg

The value 0 is the set that is a distance 0 from its negative.
The unit diameter circle is the set that is a distance 1 from its negative.

The value 1 is the set that is a distance 0 from its inverse.
What is the set that is a distance 1 from its inverse?

One point in the set is the golden ratio (1.618..) since its inverse is 0.618... but it isn't the only point in the set, for instance 0.618 is also in the set. Also cos(30) + 0.5i is in the set since its inverse is cos(30) - 0.5i. The full set is drawn in blue:

The green curve is an ellipse, since the blue oval is 'fatter' at smaller real values and thinner at larger real values it is in egg shape; I'll call it the golden egg! In fact there are two of these, one being the negative of the other.

The left and right sides can be calculated iteratively and separately, for the left side (lesser real values):
a = 1/(c + a) where c is the unit radius circle.
For the right side:
a = 1/a - c

Its formula is:
|𝑧|⁴−|𝑧|²+1=2|𝑧|²cos(2arg(𝑧))  

Here's a better image:

We can extend this idea to the quaternions, giving a very similar looking shape:

but its a 3D surface in 4D, the above is a cross section at its largest point, but it is symmetric around the x (real) axis so it is simply a distorted 3-sphere in the manner shown.

Doesn't look much like an egg we're used to, but it is definitely an ovoid. If we scale down the non-x axes by two and colour it differently it looks more like your classic egg.

The width diameter is sqrt(5)/2 so with it scaled by half it is sqrt(5)/4. And its length is of course 1.