Da a o pa amagne ic and an i e omagne ic solu ions on he squa e
la ice
Pa amagne ic MIT
The mi and im di ec o ies con ain RAW gRISB ( ) da a o he pa amagne ic Mo ansi ion,
espec i ely o inc easing and dec easing in e ac ion s eng hs (me al- o-insula o and insula o - o-me al).
The ela ed Ma lab/Oc a e sc ip is non eeness_pa a.m . Running i in Ma lab/Oc a e would p oduce a
figu e esembling he published one:
The main diffe ences lie in he shape o he a ows (which was fine uned in pos -p oduc ion, o publica ion)
and he anspa en line ende ed as dashed (as anspa encies a e no well suppo ed in Oc a e).
AFM g ound s a e
The isb_a m and 2gh_a m di ec o ies con ain RAW da a o o he AFM g ound s a e o he model, as
ound espec i ely in RISB and gRISB ( ).
The ela ed Ma lab/Oc a e sc ip is non eeness_a m.m . Running i in Ma lab/Oc a e would p oduce a figu e
esembling he published one:
N =
ghos 2
N =
ghos 2
The Oc a e e sion o he sc ip does no suppo he b oken x-axis. You may wan o adjus he y- ange ( un
e.g. ylim([0,0.04]) in he command line) o be e isualize he RISB and gRISB alues o he non eeness
(a he cos o losing sigh o he Heisenbe g limi , o cou se). The Heisenbe g alues o he o de pa ame e
and he double occupancy a e ha d coded in he sc ip , ollowing Phys. Re . B 66, 024418 (2002), as de ailed
in he pape .
S uc u e o RAW (g)RISB da a
All RISB and gRISB calcula ions con ain he ollowing files:
├── 1bdm_1.da # one-body densi y ma ix o uni -cell si e 1
├── 1bdm_2.da # one-body densi y ma ix o uni -cell si e 2
├── dens.las # a ay con aining con e ged alues o <n_{i↑}> and <n_{i↓}>
├── docc.las # a ay con aining con e ged alues o <n_{i↑}n_{i↓}>
├── e o .SC.con g # sel -consis ency e o (||λ_old - λ_new|| + ||R_old - R_new||)
├── g isb.inpu # inpu ile o ou cus om code implemen ing he gRISB equa ions
├── k.da # 2D momen a used o sample he educed B illouin zone
├── lambda. agmen s. es a # con e ged alue o he λ enso in gRISB
├── maxdis ance.con g # check on gRISB oo condi ion (max(abs(λ-λ_c)))
├── o de _pa ame e s.da # alues o he Néel o de pa ame e
├── R. agmen s. es a # con e ged alue o he R enso in gRISB
├── Z_1.con g # con e ged alue o he quasi-pa icle weigh (si e 1)
└── Z_2.con g # con e ged alue o he quasi-pa icle weigh (si e 2)
To build ou o mula ion o he local non eeness we load he densi y and double occupancy alues om
dens.las and docc.las , espec i ely. Whene e he oo condi ion is iola ed o e a h eshold o we
disca d he poin (pu ing a NaN , which is no displayed in he plo s). This ne e happens in he published
da a, bu a p og amma ic check is ha d-coded in all sc ip s.
10−3