Knowen - 2D: numerics for 4th order term's attachments

Now you are in the subtree of Energetic stability of quasicrystals project.

Attachments

Full vertex (Q, alpha). T = 0.01
GL4_2D_T_0.01.jpg

Full vertex (q,alpha). T = 0.05
GL4_2D_T_0.05.jpg

Full vertex (Q, alpha). T = 0.1
GL4_2D_T_0.1.jpg

Vertex as a fncn of angle. Q = 1.8. T = 0.1, 0.05, 0.02
GL4_2D__T_.jpg

Smooth limit of alpha ->0
4th_order_SE.jpg

Comparison of analytics and numerics for B term in self-interaction
B_term_num_vs_analyt.jpg

Matlab script for vertex calculation, mutual and self
boxd2d_new.m

Free energy denominator D(Q), electrons only for different symmetry states (F = -r^2/D)
F_denom_2D.jpg

Free energy denominator D(Q), electrons + extra repulsion C for different symmetry states (F = -r^2/D)
F_denom_2D___repulsion...

Matlab script for 2d variational calculation, including non-electronically mediated interaction C
boxd2d_var.m