WolframAlpha.com
WolframCloud.com
All Sites & Public Resources...
Products & Services
WolframOne
Mathematica
WolframAlpha Notebook Edition
Programming Lab
Finance Platform
SystemModeler
Wolfram Player
Wolfram Engine
WolframScript
Enterprise Private Cloud
Enterprise Mathematica
WolframAlpha Appliance
Enterprise Solutions
Corporate Consulting
Technical Consulting
WolframAlpha Business Solutions
Resource System
Data Repository
Neural Net Repository
Function Repository
WolframAlpha
WolframAlpha Pro
Problem Generator
API
Data Drop
Products for Education
Mobile Apps
Wolfram Player
Wolfram Cloud App
WolframAlpha for Mobile
WolframAlphaPowered Apps
Services
Paid Project Support
Wolfram U
Summer Programs
All Products & Services »
Technologies
Wolfram Language
Revolutionary knowledgebased programming language.
Wolfram Cloud
Central infrastructure for Wolfram's cloud products & services.
Wolfram Science
Technologyenabling science of the computational universe.
Wolfram Notebooks
The preeminent environment for any technical workflows.
Wolfram Engine
Software engine implementing the Wolfram Language.
Wolfram Natural Language Understanding System
Knowledgebased broadly deployed natural language.
Wolfram Data Framework
Semantic framework for realworld data.
Wolfram Universal Deployment System
Instant deployment across cloud, desktop, mobile, and more.
Wolfram Knowledgebase
Curated computable knowledge powering WolframAlpha.
All Technologies »
Solutions
Engineering, R&D
Aerospace & Defense
Chemical Engineering
Control Systems
Electrical Engineering
Image Processing
Industrial Engineering
Mechanical Engineering
Operations Research
More...
Finance, Statistics & Business Analysis
Actuarial Sciences
Bioinformatics
Data Science
Econometrics
Financial Risk Management
Statistics
More...
Education
All Solutions for Education
Trends
Machine Learning
Multiparadigm Data Science
Internet of Things
HighPerformance Computing
Hackathons
Software & Web
Software Development
Authoring & Publishing
Interface Development
Web Development
Sciences
Astronomy
Biology
Chemistry
More...
All Solutions »
Learning & Support
Learning
Wolfram Language Documentation
Fast Introduction for Programmers
Wolfram U
Videos & Screencasts
Wolfram Language Introductory Book
Webinars & Training
Summer Programs
Books
Need Help?
Support FAQ
Wolfram Community
Contact Support
Premium Support
Premier Service
Technical Consulting
All Learning & Support »
Company
About
Company Background
Wolfram Blog
Events
Contact Us
Work with Us
Careers at Wolfram
Internships
Other Wolfram Language Jobs
Initiatives
Wolfram Foundation
MathWorld
ComputerBased Math
A New Kind of Science
Wolfram Technology for Hackathons
Student Ambassador Program
Wolfram for Startups
Demonstrations Project
Wolfram Innovator Awards
Wolfram + Raspberry Pi
Summer Programs
More...
All Company »
Search
Join
Sign In
Dashboard
Groups
People
Message Boards
Answer
(
Unmark
)
Mark as an Answer
GROUPS:
Mathematics
Physics
Graphs and Networks
Wolfram Language
Wolfram Summer School
4
Andoni Royo Abrego
[WSS20] Particles and fields from topological defects in spacetime mesh
Andoni Royo Abrego
Posted
3 months ago
566 Views

0 Replies

4 Total Likes
Follow this post

Goals of the Project:
Torsion of a connection has been studied as the substrate for matter gauge field since the early days of EinsteinCartan Theory. In this work we aim to study the minimal form of the ideas shared in such theory, on the setting of the graph rewrite system, where torsion is identified with topological defects. The analysis of different topological defects and their possible dynamics in a two dimensional flat background are aim to be described and hopefully find particlelike behaviour of defects and internal gauge degrees of freedom.
Main result:
a plausible interaction between
D
5
D
7
topological defects is defined in a centered hexagonal grid where these can be regarded as particles and an infinite range interaction field is analyzed. The computational tools for further study have been also created.
Future Work:
now that the setup of the interaction is established and the computational tools are created, many interesting experiments can be done to analyze the behaviour of the defects and see 1) how the interaction works in different scenarios and 2) how well does this behaviour match with nature. Before jumping into 3D grids, there are many things left to analyze like many defect problems, rotation of defects, defects with different nature (possibly symmetric), etc.
In most modern theories of physics nature is described by spacetime, a four dimensional smooth manifold, where concepts from differential geometry like curvature and torsion play a key role describing interactions like gravity. Wolfram models, however, use discrete structures (hypergraphs) to describe nature and build physical laws on it, and hence, this differential geometry objects have to be translated to discrete metric spaces. By doing this in a proper way, hopefully, particlelike behaviour, interactions and gauge fields might come out in a natural way.
There are different ways of detecting curvature on a graph, but let me work on this first approach, which I find very natural. There are only three ways of tiling a twodimensional flat space by regular polygons: four squares, three hexagons and six equilateral triangles in each vertex. Only the triangular grid (or centered hexagonal grid) has all the closest vertices at the same distance though, plus, it can be really easily constructed from the other two just by adding edges and vertices.
G
r
a
p
h
i
c
s
G
r
i
d
[
{
{
G
r
i
d
G
r
a
p
h
[
{
4
,
4
}
,
I
m
a
g
e
S
i
z
e
S
m
a
l
l
]
,
R
e
s
o
u
r
c
e
F
u
n
c
t
i
o
n
[
"
H
e
x
a
g
o
n
a
l
G
r
i
d
G
r
a
p
h
"
]
[
{
3
,
3
}
,
I
m
a
g
e
S
i
z
e
S
m
a
l
l
]
,
A
B
C
G
r
a
p
h
[
3
]
}
}
]
I
n
[
]
:
=
O
u
t
[
]
=
Pentagons, for example, do not tile the twodimensional flat space, since the angle of each vertex is of 108º and can not sum 360º. However, if the grid is positively curved the angles of a pentagon get wider and can indeed appear in the grid. Having a centered hexagonal grid, the existence of a pentagon is an indication of positive curvature (similarly, an heptagon is a sign of negative curvature).
s
p
h
=
D
i
s
c
r
e
t
i
z
e
G
r
a
p
h
i
c
s
[
G
r
a
p
h
i
c
s
3
D
[
S
p
h
e
r
e
[
]
]
]
;
D
i
s
c
r
e
t
e
S
p
h
e
r
e
=
M
e
s
h
C
o
n
n
e
c
t
i
v
i
t
y
G
r
a
p
h
[
s
p
h
,
V
e
r
t
e
x
S
h
a
p
e
F
u
n
c
t
i
o
n
"
P
o
i
n
t
"
,
E
d
g
e
S
h
a
p
e
F
u
n
c
t
i
o
n
"
L
i
n
e
"
]
;
H
i
g
h
l
i
g
h
t
G
r
a
p
h
[
D
i
s
c
r
e
t
e
S
p
h
e
r
e
,
{
N
e
i
g
h
b
o
r
h
o
o
d
G
r
a
p
h
[
D
i
s
c
r
e
t
e
S
p
h
e
r
e
,
{
0
,
1
2
}
,
1
]
,
{
S
t
y
l
e
[
N
e
i
g
h
b
o
r
h
o
o
d
G
r
a
p
h
[
D
i
s
c
r
e
t
e
S
p
h
e
r
e
,
{
0
,
1
3
3
}
,
1
]
,
P
u
r
p
l
e
]
}
}
]
I
n
[
]
:
=
O
u
t
[
]
=
D
5
and its Dynamics
At first it might look easy to define an update rule to describe the motion of a
D
5
in the hexgrid. However, it is pretty complicated actually. Let me make some definitions first:
NONPROPAGATING DYNAMIC RULE: it is a rule that defines the motion of a defect such that it does not generate more defects and it can be written by a finite set of rules even if the graph is infinite.
It is quite obvious that a nonpropagating dynamic rule has to be able to be executed within a closed subgraph, and hence the search of a nonpropagating dynamic rule can be reduced to the search of a rearrangement of the interior of the subgraph that also generates a single defect (in a different position), without changing the nature of the boundary. By “without changing the nature of the boundary”, I mean that if a vertex belonging to the boundary has a certain amount of edges towards the interior of the subgraph, after the arrangement this amount of edges has to be the same. This must be obeyed, otherwise, that vertex would generate a defect after the rearrangement.
D
D
5
=
G
r
a
p
h
[
E
d
g
e
L
i
s
t
[
N
e
i
g
h
b
o
r
h
o
o
d
G
r
a
p
h
[
D
i
s
c
r
e
t
e
S
p
h
e
r
e
,
{
0
,
1
2
}
,
6
]
]
,
V
e
r
t
e
x
L
a
b
e
l
s
A
u
t
o
m
a
t
i
c
,
I
m
a
g
e
S
i
z
e
L
a
r
g
e
]
;
l
i
=
V
e
r
t
e
x
L
i
s
t
[
D
D
5
]
[
[
A
l
l
,
2
]
]
;
D
5
=
V
e
r
t
e
x
R
e
p
l
a
c
e
[
D
D
5
,
T
a
b
l
e
[
{
0
,
l
i
[
[
i
]
]
}
i
,
{
i
,
1
,
L
e
n
g
t
h
[
l
i
]
}
]
]
;
H
i
g
h
l
i
g
h
t
G
r
a
p
h
[
D
5
,
{
S
t
y
l
e
[
{
2
5
2
7
,
2
5
5
1
,
5
1
4
5
,
4
5
1
1
,
1
1
4
4
,
4
4
6
6
,
6
6
1
5
,
1
5
6
7
,
6
7
6
0
,
6
0
3
9
,
3
9
5
9
,
5
9
5
6
,
5
6
3
2
,
3
2
3
0
,
3
0
8
9
,
8
9
2
7
}
,
B
l
u
e
]
,
S
t
y
l
e
[
{
6
,
2
,
3
,
4
,
5
,
6
2
,
2
3
,
3
4
,
4
5
,
5
6
}
,
R
e
d
]
}
]
I
n
[
]
:
=
O
u
t
[
]
=
It is not difficult to see that due to the presence of the defect the nature of the boundary is not symmetric (in the direction of the defect), which makes a nonpropagating dynamic rule very difficult to find (without loosing the nonplanarity). I think there might not even exist such a rule, no matter how big (noninfinity) the boundary is taken to be. Nevertheless, this suggests that there are nonpropagating dynamics rules for pairs of defects (which of course can be drawn inside a symmetric boundary). All the dynamic rules I have attempt generate inevitably more defects. Here there is one example:
O
u
t
[
]
=
Coordinate system for centered hexagonal grids
The coordinate system we are going to use in the centered hexagonal grid is the
Cubic Coordinate system,
which is obtained from a cubic grid by slicing out the diagonal plane x+y+z=0. This way there are three directions on the grid, and of course one of them is redundant. However, it is convenient to keep them all since we know that the sum of all three coordinates of any vertex in the hexgrid will sum zero. We have created some functions to play with this coordinate system, equipped with a very useful set of options. Here we have some examples of what we can do with them:
Create a hexgrid of desired radius with the Cubic Coordinate system ToolTip:
A
B
C
G
r
a
p
h
[
8
]
I
n
[
]
:
=
O
u
t
[
]
=
Use colors to indicate the values of the coordinates:
G
r
a
p
h
i
c
s
G
r
i
d
[
{
{
A
B
C
G
r
a
p
h
[
3
,
C
o
l
o
r
F
u
n
c
t
i
o
n
E
x
t
r
a
c
t
[
1
]
]
,
A
B
C
G
r
a
p
h
[
3
,
C
o
l
o
r
F
u
n
c
t
i
o
n
E
x
t
r
a
c
t
[
2
]
]
,
A
B
C
G
r
a
p
h
[
3
,
C
o
l
o
r
F
u
n
c
t
i
o
n
E
x
t
r
a
c
t
[
3
]
]
}
}
]
I
n
[
]
:
=
O
u
t
[
]
=
Label the coordinates directly:
A
B
C
G
r
a
p
h
[
2
,
L
a
b
e
l
i
n
g
F
u
n
c
t
i
o
n
F
u
n
c
t
i
o
n
[
#
]
,
I
m
a
g
e
S
i
z
e
M
e
d
i
u
m
]
I
n
[
]
:
=
O
u
t
[
]
=
Highlight a specific vertex or list of them:
A
B
C
G
r
a
p
h
[
2
,
G
r
a
p
h
H
i
g
h
l
i
g
h
t
{
{
0
,
0
,
0
}
,
{

2
,
1
,
1
}
,
{

1
,

1
,
2
}
}
]
I
n
[
]
:
=
O
u
t
[
]
=
A
B
C
G
r
a
p
h
[
2
,
G
r
a
p
h
H
i
g
h
l
i
g
h
t
{
F
u
n
c
t
i
o
n
[
F
i
r
s
t
[
#
]
>
1
]
R
e
d
,
F
u
n
c
t
i
o
n
[
L
a
s
t
[
#
]
>
1
]
B
l
u
e
}
]
I
n
[
]
:
=
O
u
t
[
]
=
D
5
D
7
Defect Pair
A
D
5
D
7
defect in a centered hexagonal grid, is a defect where two adjacent vertices have degree five and seven instead. This defect is very well known in Crystallography since it generates a line defect usually referred as disclination.
O
u
t
[
]
=
Since a
D
5
D
7
defect makes a whole new line appear in the grid, there is not a rule that simply generates a single
D
5
D
7
defect starting from a centered hexagonal grid without rearranging the whole grid. Nevertheless, instead of letting the “rift” (new line generated by the defect) propagate infinitely or until some boundary, it can terminate in another
D
5
D
7
defect oriented in the opposite direction, and then, the situation can be generated from a flat hexgrid in a finite and easily computable way. Now, taking a pair of
D
5
D
7
defects separated by a certain distance, there is a variety of rifts lying inside (or in contact with) the blue diamond in the picture that could connect them and all of them have the same trivially defined distance. Check that by stitching any of the possible rifts connecting the defects (gluing the pair of vertices in the both sides of the rift ), we can close or undo the rift, and the defects would come closer and closer until they annihilate each other and the flat hexgrid is recovered.
Zip dynamics?
The process of stitching a rift is very similar to closing a zipper. The defect acts like the head of a zipper in the sense that when it moves in one direction it opens the jacket (in our case it generates a new rift) and when it moves in the opposite direction it closes it, bringing both edges of the rift together and closing it. It is literally like stitching.
G
r
a
p
h
i
c
s
G
r
i
d
[
{
{
I
m
p
o
r
t
[
"
/
U
s
e
r
s
/
a
n
d
o
n
i
r
o
y
o
a
b
r
e
g
o
/
D
o
c
u
m
e
n
t
s
/
W
S
S
2
0
2
0
/
P
R
O
J
E
C
T
/
Z
i
p
p
e
r
1
.
p
n
g
"
]
,
I
m
p
o
r
t
[
"
/
U
s
e
r
s
/
a
n
d
o
n
i
r
o
y
o
a
b
r
e
g
o
/
D
o
c
u
m
e
n
t
s
/
W
S
S
2
0
2
0
/
P
R
O
J
E
C
T
/
Z
i
p
p
e
r
2
.
p
n
g
"
]
,
I
m
p
o
r
t
[
"
/
U
s
e
r
s
/
a
n
d
o
n
i
r
o
y
o
a
b
r
e
g
o
/
D
o
c
u
m
e
n
t
s
/
W
S
S
2
0
2
0
/
P
R
O
J
E
C
T
/
Z
i
p
p
e
r
3
.
p
n
g
"
]
}
}
,
I
m
a
g
e
S
i
z
e
L
a
r
g
e
]
I
n
[
]
:
=
O
u
t
[
]
=
In this illustrations the purple path indicates the chosen rift, the yellow point is the center of the pentagon, the green point is the center of the heptagon. Highlighted black pair of edges are the edges that have been stitched in the previous step and the shadowed parallelogram is the area is going to be stitched in the next step.
It is obvious that by closing/opening the zipper we are transporting the defect along the rift, and hence there is a simple update rule for the dynamics of
D
5
D
7
defects. I think there is another pretty obvious consequence here: there is an interaction between the defects. Imagine there is a
D
5
D
7
defect somewhere in an infinitely large hexgrid. Then, any other defect in the hexgrid could find a (finite) set of rifts that would connect with the first defect, and by the stitching the rift they would come closer and closer and eventually annihilate each other. These topological defects generate fields that provide interaction between defects, and these fields are “invisible” to simple observation. Let me explain this better. If a pair of
D
5
D
7
defects is generated and get separated one from the other by the “opening of the zipper” a long enough distance, one could simply look at the space in between them and see that this space is not different from the flat hexgrid, and he/she would not be able to tell wether there is some “extra space” generated by rifts or not.
Creation/annihilation of
D
5
D
7
defect pairs
Creation of a pair of
D
5
D
7
defects can be written in a simple rewriting rule by splitting the center of an hexagon, and by applying the same rule to the already existing heptagons, the rift keeps opening.
O
u
t
[
]
=
Strength of the interaction
H
o
w
c
o
u
l
d
o
n
e
n
o
w
c
h
a
r
a
c
t
e
r
i
z
e
t
h
e
s
t
r
e
n
g
t
h
o
f
a
n
i
n
t
e
r
a
c
t
i
o
n
i
n
t
h
i
s
s
e
t
u
p
?
W
e
l
l
,
t
h
e
r
e
i
s
s
o
m
e
t
h
i
n
g
t
h
a
t
w
e
k
n
o
w
:
t
h
e
l
o
n
g
e
r
t
h
e
d
i
s
t
a
n
c
e
b
e
t
w
e
e
n
t
w
o
d
e
f
e
c
t
s
,
t
h
e
l
a
r
g
e
r
t
h
e
a
m
o
u
n
t
o
f
p
o
s
s
i
b
l
e
r
i
f
t
s
c
o
n
n
e
c
t
i
n
g
t
h
e
m
.
L
e
t
’
s
a
n
a
l
y
z
e
t
h
e
n
,
h
o
w
d
o
e
s
t
h
e
n
u
m
b
e
r
o
f
p
o
s
s
i
b
l
e
r
i
f
t
s
c
h
a
n
g
e
w
i
t
h
t
h
e
d
i
s
t
a
n
c
e
.
I
n
t
h
e
e
n
d
o
f
t
h
e
d
a
y
,
t
h
e
n
u
m
b
e
r
o
f
p
o
s
s
i
b
l
e
r
i
f
t
s
i
s
t
h
e
n
u
m
b
e
r
o
f
p
a
t
h
s
i
n
s
i
d
e
t
h
e
d
i
a
m
o
n
d
c
o
n
n
e
c
t
i
n
g
t
h
e
t
w
o
p
e
n
t
a
g
o
n
s
,
w
h
i
c
h
i
s
j
u
s
t
t
h
e
c
o
m
b
i
n
a
t
o
r
i
a
l
n
u
m
b
e
r
i
n
t
h
e
c
e
n
t
r
a
l
l
i
n
e
o
f
t
h
e
P
a
s
c
a
l
T
r
i
a
n
g
l
e
.
T
h
i
s
n
u
m
b
e
r
i
s
c
a
l
l
e
d
t
h
e
c
e
n
t
r
a
l
b
i
n
o
m
i
a
l
c
o
e
f
f
i
c
i
e
n
t
a
n
d
i
s
g
i
v
e
n
b
y
N
(
d
)
=
2
d
d
,
w
h
e
r
e
d
i
s
t
h
e
t
r
i
v
i
a
l
d
i
s
t
a
n
c
e
b
e
t
w
e
e
n
t
h
e
p
e
n
t
a
g
o
n
s
.
What can we say about this set up for interactions between topological defects?
Well, this set up has many interesting properties one could analyze and compare with known physics in continuous spacetime, as well as some problems. Let me write down some of them.
 A
D
5
D
7
defect has an orientation, which is interesting in order to build an analog with spin 1/2 particles. The possible rifts coming out the defect, however, are delimited in an angular direction, which is a big problem, since they are the underlying structure that define the sense of the field. Therefore, the study of rotation of defects is essential in order to go further with this set up, and if analogies with fermions were to be made, spinorlike behaviour should be found.
The specific interaction proposed here has, in principle, an infinite range. Moreover, since the defect can take many different but equivalent paths it is not deterministic. The multiway graph though, is obviously confluent.
There are many simulations and experiments to be done yet, and now that the complicated tool for that have been created, this ideas can be computed and analyzed much further.
POSTED BY:
Andoni Royo Abrego
Answer
Mark as an Answer
Reply

Flag
Reply to this discussion
in reply to
Add Notebook
Community posts can be styled and formatted using the
Markdown syntax
.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Publish anyway
Cancel
Reply Preview
Attachments
Remove
Add a file to this post
Follow this discussion
or
Discard
Group Abstract
Be respectful. Review our
Community Guidelines
to understand your role and responsibilities.
Community Terms of Use
Feedback
Enable JavaScript to interact with content and submit forms on Wolfram websites.
Learn how »