Skip to content

Commit 0adaca4

Browse files
committed
Added final example for thermalization session.
1 parent a9d856d commit 0adaca4

7 files changed

Lines changed: 90 additions & 11 deletions

File tree

Lines changed: 17 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,17 @@
1+
15
2+
2
3+
4
4+
9
5+
20
6+
44
7+
95
8+
206
9+
447
10+
968
11+
2097
12+
4543
13+
9841
14+
21316
15+
46169
16+
100000
17+
Lines changed: 34 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,34 @@
1+
simid bccFe100
2+
ncell 10 10 10 System size
3+
BC P P P Boundary conditions (0=vacuum, P=periodic)
4+
cell 1.00000 0.00000 0.00000
5+
0.00000 1.00000 0.00000
6+
0.00000 0.00000 1.00000
7+
Sym 0
8+
do_prnstruct 2
9+
10+
posfile ./posfile
11+
momfile ./momfile
12+
exchange ./jfile
13+
14+
Initmag 1 Initial config of moments (1=random, 2=cone, 3=spec., 4=file)
15+
16+
ip_mode N
17+
ip_mcanneal 1
18+
10000 100.0 1.00e-16 0.95
19+
100 100.0 1.00e-16 0.95
20+
100 100.0 1.00e-16 0.95
21+
22+
mode M
23+
Temp 0.001 K Temperature of the system
24+
hfield 0.00000 0.00000 0.00000 Static H field
25+
mcnstep 10000 Number of time-steps
26+
27+
do_avrg Y Measure averages
28+
plotenergy 1
29+
do_cumu Y
30+
31+
ea_model T
32+
33+
do_autocorr Y
34+
acfile ./acfile
Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
1 1 1.0 0.0 0.0 1.00
2+
1 1 -1.0 0.0 0.0 1.00
3+
1 1 0.0 1.0 0.0 1.00
4+
1 1 0.0 -1.0 0.0 1.00
5+
1 1 0.0 0.0 1.0 1.00
6+
1 1 0.0 0.0 -1.0 1.00
Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,2 @@
1+
1 7 -0.020 0.00000 0.0 1.0 0.0 0.1
2+
2 7 -0.020 0.00000 0.0 1.0 0.0 0.1
Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
1 1 1.00 1.0 0.0 0.0
Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
1 1 0.000000 0.000000 0.000000

content/finitetemp/finitetemp.rst

Lines changed: 29 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -160,14 +160,32 @@ Fig. 1 Spin spiral ground state for a DM system.
160160
* If you do not get the spiral state, try to obtain it with simulated annealing. How close can you come? Is the energy stabilizing?
161161

162162

163-
### Glassy systems
164-
### ^^^^^^^^^^^^^^
165-
###
166-
### As a final challenge, lets consider a spin glass system.
167-
### A simple yet illustrative model for a spin glass system is given by the Edwards-Anderson model, where a nearest neighbour Hamiltonian on a cubic
168-
### lattice, but with random exchange interactions are used.
169-
###
170-
### #.. math::
171-
### #H=-\sum_{i,j} J_{ij} \mathbf{m}_i \mathbf{m}_j
172-
###
173-
### The UppASD code can model Edwards-Anderson spin glasses by the keywords
163+
Glassy systems
164+
^^^^^^^^^^^^^^
165+
166+
As a final challenge, lets consider a spin glass system.
167+
A simple yet illustrative model for a spin glass system is given by the Edwards-Anderson model, where a nearest neighbour Hamiltonian on a cubic
168+
lattice, but with random exchange interactions are used.
169+
170+
The UppASD code can model Edwards-Anderson spin glasses by the keywords ``ea_model T`` and ``ea_sigma XX`` where ``ea_sigma`` contols
171+
the width of the Gaussian distribution of the randomized exchange interactions. Disregarding the ``ea_sigma`` keyword for now,
172+
we can define a simple cubic system with randomized exchange with the following ``inpsd.dat``
173+
174+
.. literalinclude:: EdwardsAnderson/inpsd.dat
175+
176+
While the ``posfile`` and ``momfile`` are made as easy as possible
177+
178+
.. literalinclude:: EdwardsAnderson/posfile
179+
.. literalinclude:: EdwardsAnderson/momfile
180+
181+
Even though the exchange interactions will be randomized, we still have to set up a Hamiltonian with existing couplings since only
182+
the magnitude and not the directions will be randomized. Thus we use the following ``jfile``
183+
184+
.. literalinclude:: EdwardsAnderson/jfile
185+
186+
With the glassy system setup. we end this exercise session with a friendly competition:
187+
188+
* Without changing the system size and Hamiltonian, apply your thermalization and minimization skills to get the lowest possible energy for this system.
189+
190+
The participant with the lowest energy will get a symbolic price during the conference dinner.
191+

0 commit comments

Comments
 (0)