You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
This example shows how to calculate the spin wave spectrum of the standard example bcc Fe and to understand the influence of the temperature on the spectra together with the influence of the uniaxial anisotropy. Files are found in the ``bccFeT1K`` and ``bccFeT300K`` folders.
12
-
13
-
Crystal & magnetic structure
14
-
^^^^^^^^^^^^^^^^^^^^^^^^^^^^
15
-
16
-
Using the lines below with the indicated files, the crystal and magnetic structure are readily available, so that a simulation of an Fe bcc system is set up.
17
-
18
-
::
19
-
20
-
simid bccFe100
21
-
ncell 20 20 20 System size
22
-
BC P P P Boundary conditions (0=vacuum, P=periodic)
23
-
cell -0.5000000000 0.5000000000 0.5000000000
24
-
0.5000000000 -0.5000000000 0.5000000000
25
-
0.5000000000 0.5000000000 -0.5000000000
26
-
Sym 1 Symmetry of lattice (0 for no, 1 for cubic, 2 for 2d cubic, 3 for hexagonal)
27
-
28
-
posfile ./posfile
29
-
momfile ./momfile
30
-
exchange ./jASD2S
31
-
anisotropy ./kfile
32
-
maptype 2
33
-
34
-
.. figure:: figures/tutorial1/fig1.png
35
-
36
-
Fig 1. Lattice and magnetic texture.
37
-
38
-
Thermalizing the system
39
-
^^^^^^^^^^^^^^^^^^^^^^^
40
-
41
-
Using the lines below, the system is brought to thermal equilibrium by means of Heat bath Monte Carlo.
42
-
43
-
::
44
-
45
-
ip_mode H Initial phase parameters
46
-
ip_mcanneal 1 --
47
-
10000 1.0 1.00e-16 0.3 --
48
-
49
-
Linear spin wave spectra
50
-
^^^^^^^^^^^^^^^^^^^^^^^^
51
-
52
-
Below the critical temperature bcc Fe has long range collinear ordering of spins. We calculate the adiabatic magnon spectra (AMS) using linear spin wave theory for collinear spin textures at the list of q points specified in the ``qfile.kpath``. Note that the spin wave is calculated for the T=0 K ground state as specified in the ``momfile``. The list of q points were calculated from the ``preQ.py`` script which analyses the space group symmetry of the crystal cell,
53
-
54
-
::
55
-
56
-
do_ams Y Collinear Adiabatic magnon spectra
57
-
do_magdos Y Calculate magnon density of states
58
-
qpoints D Direct coordinates
59
-
qfile ./qfile.kpath q points
60
-
61
-
**The first Brillouin zone of a body centered cubic lattice**
62
-
63
-
.. figure:: figures/tutorial1/fig3.png
64
-
65
-
Fig 3. Primitive and reciprocal lattice vectors in bcc.
66
-
67
-
.. figure:: figures/tutorial1/fig4.png
68
-
69
-
Fig 4. BCC 1st Brillouin zone.
70
-
71
-
.. figure:: figures/tutorial1/fig5.png
72
-
73
-
Fig 5. High symmetry points.
74
-
75
-
Plotting the spectrum
76
-
^^^^^^^^^^^^^^^^^^^^^
77
-
78
-
Use the UppASD graphical interface ``ASD_GUI`` or the ``postQ.py`` script to plot the linear spin wave spectra and the dynamical structure factor.
79
-
80
-
.. figure:: figures/tutorial1/fig6.png
81
-
82
-
Fig 6. Adiabatic magnon spectra.
83
-
84
-
Questions and exercises:
85
-
^^^^^^^^^^^^^^^^^^^^^^^^
86
-
87
-
1. Does the spectra follow the analytical expression?
88
-
2. Why the spectra is shift it up?
89
-
3. Plot the spectra without the gap around the center zone.
90
-
4. Why there are two branches, 1 acoustic and 1 optical?
91
-
5. Plot the spectrum for Fe fcc. Why now there is just 1 branch? Is it following the analytical expression?
92
-
93
-
.. figure:: figures/tutorial1/fig7.png
94
-
95
-
Fig 7. Adiabatic magnon spectra of Fe FCC.
96
-
97
5
Exercise 2: FM Heisenberg nearest-neighbour spin chain
This example shows how to calculate the spin wave spectrum of the standard example bcc Fe and to understand the influence of the temperature on the spectra together with the influence of the uniaxial anisotropy. Files are found in the ``bccFeT1K`` and ``bccFeT300K`` folders.
370
+
371
+
Crystal & magnetic structure
372
+
^^^^^^^^^^^^^^^^^^^^^^^^^^^^
373
+
374
+
Using the lines below with the indicated files, the crystal and magnetic structure are readily available, so that a simulation of an Fe bcc system is set up.
375
+
376
+
::
377
+
378
+
simid bccFe100
379
+
ncell 20 20 20 System size
380
+
BC P P P Boundary conditions (0=vacuum, P=periodic)
381
+
cell -0.5000000000 0.5000000000 0.5000000000
382
+
0.5000000000 -0.5000000000 0.5000000000
383
+
0.5000000000 0.5000000000 -0.5000000000
384
+
Sym 1 Symmetry of lattice (0 for no, 1 for cubic, 2 for 2d cubic, 3 for hexagonal)
385
+
386
+
posfile ./posfile
387
+
momfile ./momfile
388
+
exchange ./jASD2S
389
+
anisotropy ./kfile
390
+
maptype 2
391
+
392
+
.. figure:: figures/tutorial1/fig1.png
393
+
394
+
Fig 1. Lattice and magnetic texture.
455
395
396
+
Thermalizing the system
397
+
^^^^^^^^^^^^^^^^^^^^^^^
398
+
399
+
Using the lines below, the system is brought to thermal equilibrium by means of Heat bath Monte Carlo.
400
+
401
+
::
402
+
403
+
ip_mode H Initial phase parameters
404
+
ip_mcanneal 1 --
405
+
10000 1.0 1.00e-16 0.3 --
406
+
407
+
Linear spin wave spectra
408
+
^^^^^^^^^^^^^^^^^^^^^^^^
409
+
410
+
Below the critical temperature bcc Fe has long range collinear ordering of spins. We calculate the adiabatic magnon spectra (AMS) using linear spin wave theory for collinear spin textures at the list of q points specified in the ``qfile.kpath``. Note that the spin wave is calculated for the T=0 K ground state as specified in the ``momfile``. The list of q points were calculated from the ``preQ.py`` script which analyses the space group symmetry of the crystal cell,
411
+
412
+
::
413
+
414
+
do_ams Y Collinear Adiabatic magnon spectra
415
+
do_magdos Y Calculate magnon density of states
416
+
qpoints D Direct coordinates
417
+
qfile ./qfile.kpath q points
418
+
419
+
**The first Brillouin zone of a body centered cubic lattice**
420
+
421
+
.. figure:: figures/tutorial1/fig3.png
422
+
423
+
Fig 3. Primitive and reciprocal lattice vectors in bcc.
424
+
425
+
.. figure:: figures/tutorial1/fig4.png
426
+
427
+
Fig 4. BCC 1st Brillouin zone.
428
+
429
+
.. figure:: figures/tutorial1/fig5.png
430
+
431
+
Fig 5. High symmetry points.
432
+
433
+
Plotting the spectrum
434
+
^^^^^^^^^^^^^^^^^^^^^
435
+
436
+
Use the UppASD graphical interface ``ASD_GUI`` or the ``postQ.py`` script to plot the linear spin wave spectra and the dynamical structure factor.
437
+
438
+
.. figure:: figures/tutorial1/fig6.png
439
+
440
+
Fig 6. Adiabatic magnon spectra.
441
+
442
+
Questions and exercises:
443
+
^^^^^^^^^^^^^^^^^^^^^^^^
444
+
445
+
1. Does the spectra follow the analytical expression?
446
+
2. Why the spectra is shift it up?
447
+
3. Plot the spectra without the gap around the center zone.
448
+
4. Why there are two branches, 1 acoustic and 1 optical?
449
+
5. Plot the spectrum for Fe fcc. Why now there is just 1 branch? Is it following the analytical expression?
0 commit comments