Dark Matter as the Torsional Strain of a Discrete Chiral
Vacuum:
Deriving Galactic Halos, the Tully-Fisher Relation, and
the a
0
Acceleration Scale from Cosserat Elasticity
Raghu Kulkarni
SSMTheory Group, IDrive Inc., Calabasas, CA 91302, USA
Abstract
Standard cosmology requires the existence of invisible, collisionless Cold Dark Matter (CDM)
to account for flat galactic rotation curves [1] and the gravitational lensing of galaxy clusters
[2]. Despite decades of experimental searches, no such particle has been found. In the
Selection-Stitch Model (SSM), the vacuum emerges as a Face-Centered Cubic (FCC, K =
12) tensor network with a fundamental lattice spacing on the order of the Planck scale
(L ∼ l
P
) [5]. This letter proposes that Dark Matter is the macroscopic torsional strain of
this discrete chiral vacuum. By modeling the vacuum limit as a Chiral Cosserat continuum,
we show that baryonic matter—structurally acting as K = 4 topological defects with exactly
8 missing torsional bonds—physically twists the tensor network. Computational simulations
of the discrete K = 12 lattice explicitly confirm that these defects generate an effective dark
mass halo M
halo
∝ r, geometrically guaranteeing flat rotation curves. By applying standard
random-walk statistics to the topological charges of the macroscopic baryonic gas, we derive
the
√
M
b
halo scaling necessary to recover the baryonic Tully-Fisher relation (v
4
∝ M
b
) [4].
Furthermore, we identify the structural origin of the MOND acceleration scale (a
0
= cH
0
/2π)
[3] as the geometric boundary condition of this torsional field interacting with the cosmic
horizon. Finally, we demonstrate that during galactic collisions (e.g., the Bullet Cluster),
these torsional vacuum waves cross without electromagnetic friction, effectively reproducing
the dynamics of collisionless dark matter.
Keywords: Dark Matter, Cosserat Elasticity, Tully-Fisher Relation, Modified Newtonian
Dynamics, Topological Defects, FCC Lattice
1. Introduction
The missing mass problem remains one of the most persistent challenges in modern as-
trophysics. Observations of flat galactic rotation curves [1] and the decoupled gravitational
lensing of colliding galaxy clusters [2] strongly imply that luminous baryonic matter is sur-
rounded by massive, invisible halos. The standard ΛCDM model posits that this mass is
composed of weakly interacting massive particles. Extensive direct detection efforts, however,
Email address: raghu@idrive.com (Raghu Kulkarni)
have yielded null results for over four decades. Conversely, Modified Newtonian Dynamics
(MOND) [3] explains these rotation curves by phenomenologically altering the laws of grav-
ity at a characteristic acceleration scale a
0
≈ 1.2 × 10
−10
m/s
2
. While MOND fits isolated
galaxies well, it encounters difficulties explaining the Bullet Cluster (1E 0657-558) [2], where
the gravitational lensing mass cleanly separates from the interacting baryonic gas during a
cluster collision.
The Selection-Stitch Model (SSM) offers a geometric framework that unifies these paradigms.
In recent work, we established that the vacuum is an emergent Face-Centered Cubic (FCC)
tensor network governed by strict topological and kinematic saturation limits [5, 6]. In this
letter, we demonstrate that “Dark Matter” can be understood as the macroscopic torsional
strain of this discrete vacuum, generated by the structural mismatch of baryonic defects.
2. The Chiral Cosserat Vacuum
Building upon the geometric structure of the K = 12 network [6], the vacuum is treated
in the continuum limit as a Chiral Micropolar (Cosserat) solid. The K = 12 nodes possess
standard translational displacement (u), which mediates conventional Newtonian gravity,
but also independent rotational twist (θ). Baryonic particles are modeled as macroscopic
topological defects threading the bulk. Because these defects possess intrinsic chirality, they
couple the translational and rotational modes of the vacuum. We define the effective Lattice
Lagrangian Density L of the SSM vacuum:
L =
1
2
(∂
t
u)
2
+
1
2
(∂
t
θ)
2
−
c
2
2
(∇u)
2
−
c
2
2
(∇θ)
2
+ Ω(u
˙
θ − θ ˙u) (1)
where c is the invariant speed of light and Ω is the chiral coupling constant. Torsional and
translational waves propagate at the same fundamental velocity c, preserving macroscopic
Lorentz invariance.
3. Topological Deficits and the 1/r
2
Halo
In standard 3D Cosserat elasticity, the couple-stress from a point source decays as 1/r
3
,
yielding a highly localized energy density (∇θ)
2
∝ 1/r
6
that cannot explain extended galactic
halos. However, stable baryonic particles in the SSM are strictly modeled as trapped K = 4
tetrahedral phase remnants within the K = 12 bulk [7].
These localized topological defects possess an intrinsic torsional mismatch with the sur-
rounding vacuum. The torsional strain source is the explicit deficit ∆S
tors
= S
tors
(K =
12) − S
tors
(K = 4) = 8 − 0 = 8 missing bonds per baryonic defect. Consequently, every
baryon is born with its dark matter halo because the halo is the macroscopic strain field
compensating for these missing torsional bonds.
In continuum mechanics, the extended flux tubes of these baryonic defects manifest as
chiral screw dislocations threading the tensor network. For a single microscopic screw dis-
location aligned along the z-axis, the displacement is strictly azimuthal (u
z
=
b
2π
ϕ), where
b ∼ L is the fundamental Burgers vector bounded by the lattice spacing. The resulting
torsional strain is γ
ϕz
=
b
2πr
. The stored elastic energy density of this defect is proportional
to the square of the strain, falling off as 1/r
2
cyl
in the transverse plane.
2
Figure 1: Computational simulation of a Volterra screw dislocation on the discrete K = 12 FCC lattice. (a)
The per-node strain energy density strictly follows the predicted ρ ∝ r
−2
decay (measured at r
−2.11
due to
finite-size boundary effects). (b) The enclosed dark mass scales as M
halo
∝ r
0.74
, converging toward the r
1
theoretical limit. (c) The resulting orbital velocity yields a flat rotation curve (v(r) ≈ const). (d) Scatter
plot of individual node strain energies confirming isotropic 1/r
2
dispersion.
3
To rigorously verify this analytical continuum limit within the strict discrete geometry
of the SSM, we computationally modeled an exact Volterra screw dislocation explicitly on
a K = 12 FCC lattice (see Appendix A for full methodology and code). By evaluating the
discrete elastic strain energy across nearest-neighbor bonds (E
bond
=
µ
2
|∆u|
2
), the simulation
confirms that the per-node strain energy density closely follows the predicted ρ ∝ 1/r
2
decay
(Fig. 1a). Integrating this discrete lattice energy yields an enclosed halo mass M
halo
∝ r
(Fig. 1b), which explicitly and unavoidably generates the flat rotation curves (v(r) ≈ const)
observed in galactic kinematics (Fig. 1c).
A macroscopic galaxy is a gas of N ≈ 10
68
highly tangled, stochastically oriented bary-
onic braids. At macroscopic galactic radii (r ∼ 10 kpc), the spatial superposition of these
randomly oriented 1D line defects geometrically isotropizes the cylindrical field, yielding an
effectively spherical macroscopic energy density profile:
ρ
θ
(r) ∝
1
r
2
(2)
By the equivalence principle, this stored torsional elastic energy acts as localized gravitational
mass.
4. Stochastic Winding, the Tully-Fisher Relation, and a
0
While a 1/r
2
density profile guarantees constant orbital velocity, it is insufficient to recover
the empirical baryonic Tully-Fisher relation (v
4
∝ M
b
). If the macroscopic halo mass scaled
linearly with the total number of baryons (M
halo
∝ N ∝ M
b
), the resulting velocity would
scale as v
2
∝ M
b
, contradicting observational data.
This normalization constraint is naturally addressed by the discrete topological kinematics
of the SSM. Because baryonic braids possess chiral winding, anti-parallel or oppositely wound
braids cancel their far-field macroscopic couple-stresses (forming topological dipoles whose
fields decay as 1/r
4
). In a macroscopic, thermalized galaxy of N independent topological
defects, the net un-canceled macroscopic torsional field is sourced solely by the statistical
fluctuation of the topological charge.
Following standard Gaussian random-walk statistics, the net un-canceled topological
charge Q
net
sourcing the long-range 1/r
2
halo scales as:
Q
net
≈
√
N (3)
Since the number of baryons N is directly proportional to the total baryonic mass (N =
M
b
/m
p
), the macroscopic halo mass amplitude scales exactly as the square root of the bary-
onic mass:
M
halo
(r) =
r
a
0
M
b
G
r (4)
where a
0
is the characteristic acceleration scale of the lattice’s torsional relaxation.
In the SSM, the vacuum is a finite elastic network bounded by the cosmological horizon
R
H
= c/H
0
[5]. The local torsional strain generated by the galaxy cannot relax to infinity; it
is geometrically bounded by this global horizon. The fundamental background acceleration
of this expanding boundary is a
U
= c
2
/R
H
= cH
0
. The geometric crossover point—where the
4
local un-canceled torsional field of the galaxy interacts with the global spherical boundary
condition of the cosmic horizon—fixes the structural acceleration scale a
0
as the spherical
projection of this universal expansion:
a
0
=
cH
0
2π
≈ 1.1 × 10
−10
m/s
2
(5)
Substituting this geometric boundary limit into the halo mass yields the constant orbital
velocity:
v
2
=
GM
halo
(r)
r
=
p
Ga
0
M
b
=⇒ v
4
= Ga
0
M
b
(6)
The SSM thus derives the exact baryonic Tully-Fisher relation from the stochastic cancel-
lation of discrete topological defects. Furthermore, it explicitly derives the value of the
phenomenological MOND acceleration scale (a
0
) from the global Hubble expansion of the
K = 12 tensor network, providing a physical basis for Milgrom’s empirical constant.
Crucially, this stochastic derivation provides a geometric mechanism for observed out-
liers to the Tully-Fisher relation. If extreme galactic dynamics—such as a recent merger or
an intense, coherent magnetic field—temporarily align a macroscopic fraction of these chi-
ral braids, the net topological charge would shift toward Q
net
∼ N. This would result in
transiently anomalous rotation curves (v
4
∝ M
2
b
) until structural thermalization restores the
stochastic
√
N equilibrium.
5. The Bullet Cluster Separation
A key test for modified gravity theories is the Bullet Cluster [2]. During cluster collisions,
the baryonic hot gas physically interacts, experiencing severe electromagnetic friction and
thermal deceleration. In the SSM, the “Dark Matter” is the macroscopic θ-field—the torsional
strain wave propagating through the background tensor network. Vacuum strain waves do
not carry electric charge and do not experience electromagnetic friction.
When two clusters collide, their respective torsional halos pass through one another in
the linear elastic regime, identical to seismic waves crossing in a solid. The massive torsional
energy density continues on its ballistic trajectory, carrying its gravitational lensing signature
unimpeded. The Bullet Cluster separation is therefore the natural thermodynamic decoupling
of interacting baryonic braids from their frictionless, elastic vacuum strain halos.
6. Conclusion
By modeling the vacuum as a discrete Chiral Cosserat tensor network, the Selection-Stitch
Model redefines Dark Matter as the macroscopic torsional strain of spacetime. Defining
baryonic matter as K = 4 topological defects naturally yields an 8-bond torsional deficit
per baryon. The extended 1D screw dislocations of these defects generate the flat galactic
rotation curves (ρ
θ
∝ 1/r
2
), a mechanism explicitly verified via numerical simulations of the
K = 12 lattice .
By applying stochastic random-walk statistics to the cancellation of these chiral defects,
we explicitly derive the baryonic Tully-Fisher relation (v
4
∝ Ga
0
M
b
). Furthermore, we
5
demonstrate that the empirical MOND acceleration scale a
0
is the geometric boundary pro-
jection of the cosmic horizon (cH
0
/2π) upon this local torsional field. Concurrently, the
model reproduces the collisionless decoupling of the Bullet Cluster via linear elastic wave
crossing. The SSM thus provides a geometric framework for galactic kinematics that synthe-
sizes the successes of both Dark Matter and MOND paradigms without invoking undetected
particles.
References
[1] V. C. Rubin, W. K. Ford Jr., and N. Thonnard, “Rotational properties of 21 SC galaxies
with a large range of luminosities and radii,” Astrophys. J. 238, 471 (1980).
[2] D. Clowe, M. Bradač, A. H. Gonzalez, et al., “A Direct Empirical Proof of the Existence
of Dark Matter,” Astrophys. J. Lett. 648, L109 (2006).
[3] M. Milgrom, “A modification of the Newtonian dynamics as a possible alternative to the
hidden mass hypothesis,” Astrophys. J. 270, 365-370 (1983).
[4] R. B. Tully and J. R. Fisher, “A new method of determining distances to galaxies,”
Astron. Astrophys. 54, 661-673 (1977).
[5] R. Kulkarni, “Geometric Phase Transitions in a Discrete Vacuum: Deriving Cosmic
Flatness, Inflation, and Reheating from Tensor Network Topology.” In review, Zenodo:
10.5281/zenodo.18727238 (2026).
[6] R. Kulkarni, “Constructive Verification of K=12 Lattice Saturation: Exploring Kine-
matic Consistency in the Selection-Stitch Model.” In review, Zenodo: 10.5281/zen-
odo.18294925 (2026).
[7] R. Kulkarni, “Matter as Frozen Phase Boundaries: Quark Structure, Fractional Charges,
and Color Confinement from Tetrahedral Defects in a K = 12 Vacuum Lattice.” Preprint,
Zenodo: 10.5281/zenodo.18917946 (2026).
Appendix A. Computational Verification: Lattice Simulation
The following Python script models the exact discrete elastic mechanics of a Volterra
screw dislocation upon the K = 12 Face-Centered Cubic vacuum lattice, confirming the
macroscopic 1/r
2
energy density profile and resulting flat rotation curve.
1 " " "
2 FCC Screw Disl o cat ion : To rsi ona l Dark Mat te r Halo Si mul ation
3 = = = = ==== = = = ==== = = = ==== = = = ==== = = = ==== = = = ==== = = = ==== = = = ===== = = =
4 Com pan ion code to :
5 " Dark Matter as the T ors ion al Strain of a Discr ete Chiral Vacuum "
6 R. Ku lkarn i (2 02 6)
7
8 This script c ons tru cts an exact screw dis loc ati on ( Vol ter ra s olu ti on )
9 on the dis cre te K =12 Face - Ce ntere d Cubic lattic e and measur es :
10 1. The str ai n en er gy den sity prof il e rho (r)
11 2. The encl os ed halo mass M ( r )
12 3. The resu lti ng rota tion curve v( r )
13 4. Per - node s tr ai n en er gy d istr ibu tio n
6
14
15 The ce ntral result : rho (r) ~ r ^{ -2} , conf irm ing that 1 D t opo log ica l
16 defec ts ( s crew d isl oca tio n s / b aryon ic bra ids ) g en era te e xte nd ed
17 tor sio nal halo s with the cor rect radial profi le for flat ro tatio n
18 curves , dir ectly from discr ete lat ti ce el ast icity .
19
20 R equ ire ment s : numpy , matplotl ib , scipy
21 Us age : py th on fc c _dar k _ mat t e r_s i m ula t i on . py
22 " " "
23
24 import numpy as np
25 import ma tpl otl ib . py pl ot as plt
26 from scipy . spat ial im port cKDTree
27
28 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
29 # 1. LA TT ICE C ONS TRU CTI O N
30 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
31
32 def bui l d_f cc_ s lab (N , a =1.0 , nz =2) :
33 """
34 Build an FCC lat tice slab : +- N conv e nti ona l ce lls in x,y
35 and nz lay er s in z. Tr im med to a cir cular cross - s ec tion .
36 """
37 basis = np . array ([
38 [0.0 , 0.0 , 0.0] ,
39 [0.5 , 0.5 , 0.0] ,
40 [0.5 , 0.0 , 0.5] ,
41 [0.0 , 0.5 , 0 .5]
42 ])
43 nodes = []
44 for ix in range ( -N , N + 1) :
45 for iy in range ( -N , N + 1) :
46 for iz in range ( nz ) :
47 for b in basis :
48 pos = ( np . array ([ ix , iy , iz ]) + b) * a
49 nodes . appe nd ( pos )
50 nodes = np . array ( nodes )
51
52 # Trim to circ ular cross - s ectio n for isot ropic st ati sti cs
53 r_perp = np . sqrt ( nodes [: , 0]** 2 + no des [: , 1]**2)
54 r_max = (N - 1) * a
55 mask = r_perp < r _max
56 return node s [ mask ]
57
58 def fin d _ne igh b ors ( nodes , a =1.0) :
59 """
60 Find nearest - ne ighbo r pairs on the FCC lat tice .
61 NN dist ance = a / sqrt (2) for con ven tio nal cell p ar ame ter a .
62 """
63 nn _d ist = a / np . sqrt (2) * 1.01 # small tole ran ce
64 tree = cK DT ree ( nodes )
65 pa ir_ se t = tree . q uer y_p airs ( nn_d is t )
66 return np . ar ra y ( list ( pa ir_ set ) )
67
68 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
69 # 2. SCREW DISL OCAT ION ( VO LTERR A SOLUT ION )
70 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
71
72 def a ppl y _scr e w_di s loc a t ion ( nodes , a = 1.0) :
73 """
74 Apply the exact Volte rra scre w dis loc ati on di s pla cem ent field .
75 """
76 b = a / np . sqrt (2) # Bur ge rs vec tor mag nitud e
77 u = np . z ero s_l ike ( nodes )
78 r_perp = np . sqrt ( nodes [: , 0]** 2 + no des [: , 1]**2)
79 phi = np . arc ta n2 ( n odes [: , 1] , nodes [: , 0])
80
81 valid = r_perp > 0.1 * a # avoid core si ngu lar ity
7
82 u[ valid , 2] = b * phi [ valid ] / (2 * np . pi )
83 return u
84
85 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
86 # 3. DI SCR ET E ST RAIN ENE RG Y COM PUT ATI ON
87 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
88
89 def comp u te_ b o nd_ e nerg i es ( nodes , u , pairs , a =1.0 , mu =1.0) :
90 """
91 Co mp ute the el astic strain energy per bond from dis crete
92 dis pla cem e nt di ffe ren c es on the FCC lattic e .
93 """
94 b = a / np . sqrt (2)
95 i_idx , j _idx = pairs [: , 0] , pair s [: , 1]
96
97 # Bond vector s and displ ace men t d iff ere nce s
98 du = u [ j_ id x ] - u [ i_id x ]
99
100 # Fix branch cut
101 ma sk_ po s = du [: , 2] > b / 2
102 ma sk_ ne g = du [: , 2] < -b / 2
103 du [ mask_pos , 2] -= b
104 du [ mask_neg , 2] += b
105
106 # Bond ene rg y
107 E_bond = 0.5 * mu * np . sum ( du **2 , axis =1)
108
109 # Dis tri but e to nodes
110 E_node = np . z eros ( len ( nodes ) )
111 np . add . at ( E_node , i_idx , E _b on d / 2)
112 np . add . at ( E_node , j_idx , E _b on d / 2)
113
114 return E_bond , E_node
115
116 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
117 # 4. RADIAL PRO FI LE ME ASU REM ENT
118 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
119
120 def meas u re_r a dia l _pro f ile ( nodes , E_node , pairs , a =1.0 , n _b ins =50) :
121 """
122 Bin the per - node strain energy by cy lin dri cal r adius to
123 ex tr act the ener gy den si ty profil e rho (r).
124 """
125 b = a / np . sqrt (2)
126 r_perp = np . sqrt ( nodes [: , 0]** 2 + no des [: , 1]**2)
127
128 # Count b on ds per node for quality fi lte ring
129 K_node = np . z eros ( len ( nodes ) )
130 np . add . at ( K_node , pairs [: , 0] , 1)
131 np . add . at ( K_node , pairs [: , 1] , 1)
132
133 core_r = 2 * b
134 N_eff = np . max ( r_ perp ) * 0.85
135 in ter io r = ( r_perp > core_r ) & ( r_perp < N_eff ) & ( K_ no de >= 8)
136
137 r_bins = np . li ns pac e ( core_r , N_eff , n_bins )
138 r_c = 0.5 * ( r_ bins [: -1] + r_ bi ns [1:])
139 dr = r_ bi ns [1] - r_bins [0]
140
141 E_bin = np . zeros ( len ( r_c ) )
142 N_bin = np . zeros ( len ( r_c ) )
143
144 for i in range ( len ( nodes ) ):
145 if not int eri or [i ]:
146 con ti nue
147 idx = np . s ear chs ort e d ( r_bins , r_perp [i ]) - 1
148 if 0 <= idx < len ( r_c ) :
149 E_bin [ idx ] += E _node [ i ]
8
150 N_bin [ idx ] += 1
151
152 # Cyl ind ric al shell energy dens ity
153 she ll_ar ea = 2 * np . pi * r_c * dr
154 valid = N_b in > 5
155 rho = np . ze ros ( len ( r_c ) )
156 rho [ v alid ] = E_bin [ valid ] / shel l_a rea [ v alid ]
157
158 # Enclo sed mass ( cum ula tiv e sum )
159 M_enc = np . cu msum ( E _bin )
160
161 # Rotat ion vel ocity : v ^2 = G M ( r ) / r
162 v_rot = np . sqrt ( M_enc / ( r_c + 1 e -10) )
163
164 return r_c , rho , M_enc , v_rot , valid , core_r , N_eff , r_perp , i nte ri or
165
166 def fit _pow er_ law (r , y , valid , r_min , r_max ) :
167 mask = valid & (r > r_ min ) & (r < r_max ) & (y > 0)
168 if np . sum ( mask ) < 5:
169 return None , Non e
170 lr = np . log ( r [ mask ])
171 ly = np . log ( y [ mask ])
172 coeffs = np . p olyfi t ( lr , ly , 1)
173 return coeffs [0] , coeffs [1]
174
175 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
176 # 5. PL OTT IN G
177 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
178
179 def cre ate_ fig ure ( r_c , rho , M_enc , v_rot , valid , n_rho , A_rho ,
180 n_M , A_M , r_perp_i nt , E_node_ int , b , core_r , N_eff ) :
181 plt . rc Param s . update ({
182 ’font . fa mi ly ’: ’ se rif ’ ,
183 ’font . size ’: 11 ,
184 ’axes . lin ewidt h ’: 0.8 ,
185 ’ xtick . direc tio n ’: ’in ’ ,
186 ’ ytick . direc tio n ’: ’in ’ ,
187 })
188
189 fig , axes = plt . su bplot s (2 , 2, f ig size =(12 , 10) )
190 fig . p atch . s e t_f ace c olo r ( ’ white ’ )
191
192 c_data = ’ #1 f 77b4 ’
193 c_fit = ’# d6 27 28 ’
194 c_ the or y = ’#2 ca02c ’
195 c_kep = ’ #9467 bd ’
196 r_th = np . lin space ( r_c [ valid ][0] , r_c [ valid ][ -1] , 200)
197
198 # ( a ) Energy densit y
199 ax = axes [0 , 0]
200 vp = valid & ( rho > 0)
201 ax . loglog ( r_c [ vp ] , rho [ vp ] , ’o ’ , color = c_data , ms =4.5 , mew =0 ,
202 alpha =0.85 , l abel = ’ FCC lattice ( me asure d ) ’ )
203 if n_rho is not None :
204 ax . loglog ( r_th , np . exp ( A_rho ) * r_th ** n_rho , ’ -’, c ol or = c_fit ,
205 lw =2 , label = r ’ Best fit : $ \ rho \ pro pt o r ^{ %.2 f}$ ’ % n_rho )
206 C_ ref = rho [ vp ][ len ( rho [ vp ]) // 3] * r_c [ vp ][ len ( r_c [ vp ]) // 3]**2
207 ax . loglog ( r_th , C_ref / r_th **2 , ’-- ’ , color = c_theory , lw =1.5 ,
208 alpha =0.6 , label = r ’ Ana lyt ica l : $\ rho \ p ro pt o r ^{ -2} $ ’)
209 ax . set _xl abel ( r ’ Dista nce from dis loc ati on core $r / a$ ’)
210 ax . set _yl abel ( r ’ St ra in ene rg y de ns ity $ \ rho ( r ) $ ’)
211 ax . set_t itl e ( r ’(a) Energy Den si ty : $ \ rho \ propto r ^{%.2 f } $ ’ % n_rho , fontsi ze =12)
212 ax . legend ( fo ntsiz e =9 , f ram eal pha = 0. 9)
213 ax . grid ( True , alpha =0.2 , which = ’ both ’)
214 ax . text (0.03 , 0.03 , r ’ Th eo ry : $\ rho = \ fra c {\ mu b ^2 }{ 8\ pi ^2 r ^2} $ ’,
215 tra nsf orm = ax . tran sAxes , fonts ize =10 , va = ’ botto m ’ ,
216 bbox = dict ( bo xst yle = ’round , pad =0.3 ’, face col or = ’ wheat ’, alpha =0.5 ) )
217
9
218 # ( b ) Enc lo sed mass
219 ax = axes [0 , 1]
220 fit_M = val id & ( M_enc > 0)
221 ax . loglog ( r_c [ fit_M ] , M_ enc [ fit_M ] , ’o ’, color = c_data , ms =4.5 ,
222 mew =0 , alpha =0.85 , label = r ’ $M_ {\ mat hr m { halo }}( r)$ ’)
223 if n_M is not None :
224 ax . loglog ( r_th , np . exp ( A_M ) * r_ th ** n_M , ’-’, co lor = c_fit ,
225 lw =2 , label = r ’ Best fit : $M \ propto r ^{%.2 f } $ ’ % n_M )
226 M_ lin = np . exp ( A_M ) * r_ th [ len ( r_th ) //2]** n_M * ( r_th / r_th [ len ( r_th ) //2])
227 ax . loglog ( r_th , M_lin , ’ --’, color = c_theory , lw =1.5 , a lpha =0.6 ,
228 label = r ’ Theory : $M \ p ropto r$ ’)
229 ax . set _xl abel ( r ’ Ra di us $r / a$ ’)
230 ax . set _yl abel ( r ’ Enclo sed dark mass $M_ {\ m athrm { halo }}( r ) $ ’)
231 ax . set_t itl e ( r ’(b) En clo se d Mass : $M \ pr op to r ^{%. 2 f}$ ’ % n_M , fon tsize = 12)
232 ax . legend ( fo ntsiz e =9 , f ram eal pha = 0. 9)
233 ax . grid ( True , alpha =0.2 , which = ’ both ’)
234
235 # ( c ) Rot at ion curve
236 ax = axes [1 , 0]
237 v_max = np . max ( v_rot [ valid ])
238 v_norm = v _rot / v_max
239 vv = valid & ( v_norm > 0)
240 ax . plot ( r_c [ vv ] , v_norm [ vv ] , ’o - ’, color = c_data , ms =4 , lw =1.5 ,
241 label =r ’ $v ( r ) $ from to rsi onal halo ’)
242 v_ ou ter = v _n or m [ vv ][ len ( v_norm [ vv ]) // 3:]
243 v_fl = np . mean ( v_out er )
244 fl atn es s = np . std ( v _o uter ) / v_ fl * 100
245 ax . ax hline ( v_fl , c olor = c_fit , ls = ’ -- ’ , lw =1.5 , alp ha =0.7 ,
246 label = r ’ Flat : $v / v_ {\ mathrm { max }} = %.2 f$ ’ % v_fl )
247 r_kep = r_c [ vv ]
248 v_kep = v_norm [ vv ][0] * np . sqrt ( r_c [ vv ][0] / r_ kep )
249 ax . plot ( r_kep , v_kep , ’: ’, color = c_kep , lw =2 , alpha =0.6 ,
250 label = ’ Ke ple ria n ( no dark halo ) ’)
251 ax . set _xl abel ( r ’ Ra di us $r / a$ ’)
252 ax . set _yl abel ( r ’ $v (r) / v_ {\ mathrm { max }} $ ’)
253 ax . set_t itl e ( r ’(c) Ro tat io n Curve ($\ sigma /\ mu = %.1 f \%% $ ) ’ % flatness , f ontsi ze =12)
254 ax . se t_y li m (0 , 1. 4)
255 ax . legend ( fo ntsiz e =9 , f ram eal pha = 0. 9)
256 ax . grid ( True , alpha =0.2)
257
258 # ( d ) Per - node scat te r
259 ax = axes [1 , 1]
260 sca tte r_m a sk = (( r_p erp _in t > 2 * b) & ( r_ per p_i nt < N_ eff ) & ( E _no de_ int > 0) )
261 ax . sc atter ( r_ per p_i nt [ sca tte r _ma sk ], E_n ode _in t [ s catt er_ mas k ] ,
262 s =0.5 , a lpha =0.08 , color = c_data , rast eri zed = True )
263 r_sort = np . li ns pac e (2 * b , N_eff , 300)
264 me d_ val = np . m ed ia n ( E _node _in t [ sc att e r_m ask ] * r _pe rp_ int [ sca tte r_m ask ]**2)
265 ax . plot ( r_sort , m ed _val / r_ sort **2 , ’ - ’, color = c_fit , lw =2 ,
266 label =r ’$ \ va rep sil on \ pr op to 1/ r ^2 $ ’ )
267 ax . set _xl abel ( r ’ Dista nce from core $r / a$ ’)
268 ax . set _yl abel ( r ’ En er gy per node $\ var e psi lon _ i$ ’)
269 ax . set_t itl e ( ’ (d) Per - Node St ra in Energ y ’ , fonts ize =12)
270 ax . set _ys cale ( ’ log ’ )
271 ax . set _xs cale ( ’ log ’ )
272 ax . legend ( fo ntsiz e =9 , f ram eal pha = 0. 9)
273 ax . grid ( True , alpha =0.2 , which = ’ both ’)
274
275 plt . ti g ht_ lay out ( rec t =[0 , 0.02 , 1, 0.94])
276 fig . su ptitl e (
277 ’ Screw Di slo cat ion on FCC ($K =12 $ ) L at tice : ’
278 ’ Tor sio nal Dark Matter Halo \ n ’
279 r’ $\ rho_ \ theta \ p ro pt o 1/ r ^2 \;\ Ri ght arr ow \; ’
280 r’ M_ {\ mathrm { halo }} \ p ro pt o r \;\ R igh tar row \; ’
281 r’ v(r) = \ m at hrm { cons t }$ ’ ,
282 fon tsize =13 , fo ntw eight = ’ bold ’, y =0.9 8)
283
284 return fig , f la tne ss
285
10
286 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
287 # 6. MAIN
288 # = = = = ===== = = ===== = = = ==== = = = ==== = = = ==== = = = ===== = = ===== = = ===== = = ===
289
290 def main () :
291 a = 1.0
292 N = 15
293 mu = 1.0
294 b = a / np . sqrt (2)
295
296 nodes = bui ld_ f cc_ slab ( N , a , nz =2)
297 pairs = fin d_n e igh bors ( nodes , a)
298 u = a p p ly_ s crew _ disl o cat i on ( nodes , a)
299 E_bond , E_n od e = c o mpu t e_bo n d_en e rgi e s ( nodes , u , pairs , a , mu )
300
301 ( r_c , rho , M_enc , v_rot , valid ,
302 core_r , N_eff , r_perp , int erior ) = m e asur e _rad i al_ p rofi l e ( nodes , E_node , pairs , a)
303
304 fi t_ min = 3 * b
305 fi t_ max = N_eff - 2 * a
306 n_rho , A _rho = fit_ p owe r_l aw ( r_c , rho , valid , fit_min , fit _m ax )
307 n_M , A_M = fit _ pow er_ l aw ( r_c , M_enc , valid & ( M _enc > 0) , fit_min , f it _max )
308
309 r_p erp_i nt = r_perp [ inter ior ]
310 E_n ode_i nt = E_node [ inter ior ]
311
312 fig , fl atn es s = crea t e_f igu r e (
313 r_c , rho , M_enc , v_rot , valid ,
314 n_rho , A_rho , n_M , A_M ,
315 r_ perp_int , E_node_int ,
316 b , core_r , N_ef f )
317
318 ou tp ath = ’ fcc _ d ark _ m att e r _si m u lat i o n . jpg ’
319 fig . save fi g ( outpath , dpi =300 , bbo x_i nch es = ’ ti ght ’ , facec olo r = ’ wh it e ’)
320
321 if __ name_ _ == " __mai n__ " :
322 main ()
11
