Re: Effective damping program
Gavin Abo
Hi,
I don't believe the time-series is needed.
When you use the effective-damping, it does the equilibration, 30 degree rotation, and the time-series. At least, that is what the source code is showing:
I did try to reproduce FIG. 3 (a) but didn't have much luck in doing so as seen by the following plot:
In that plot, even my graphing of the analytic equations (9) for My and Mz don't look quite right.
It may be worth mentioning, though, that I used a VAMPIRE 7.0.0
develop version and haven't tried it yet using VAMPIRE 5.0 that
was mentioned in the paper [1].
One notable feature of FIG. 3 (a) in [1] is that it looks like it
was plotted to exclude the equilibration phase of the simulation,
whereas I left that in I my plot above. On the other hand, the
equilibration phase is shown for FIG. 5 of the same paper at [1],
where it looks like they must have shifted the outputted
simulation time so that t = 0 picoseconds corresponds with the end
of the equilibration phase.
I hope that will help you at least a little bit.
username@computername:~/effective-damping$ ls -l
total 12
-rw-r--r-- 1 username username 1198 Jul 2 10:08 CoFeB.mat
-rw-r--r-- 1 username username 1283 Jul 2 16:00 input
-rw-r--r-- 1 username username 617 Jul 2 16:07 plotm
username@computername:~/effective-damping$ mpirun -np 4
~/vampire/vampire-parallel
Processor 2 of 4 online on host computername
Processor 3 of 4 online on host computername
Processor 4 of 4 online on host computername
_
(_)
__ ____ _ _ __ ___ _ __ _ _ __ ___
\ \ / / _` | '_ ` _ \| '_ \| | '__/ _ \
\ V / (_| | | | | | | |_) | | | | __/
\_/ \__,_|_| |_| |_| .__/|_|_| \___|
| |
|_|
Version 7.0.0 Jun 21 2026 20:55:03
Git commit: 8707dedc7997cb60c06bb2efee2d64b9f6570a19
Licensed under the GNU Public License(v2). See licence file for
details.
Developers: Richard F L Evans, Sarah Jenkins, Andrea Meo,
Daniel Meilak, Andrew Naden, Matthew Ellis,
Oscar Arbelaez, Sam Morris, Rory Pond, Weijia Fan,
Phanwadee Chureemart, Pawel Sobieszczyk, Joe
Barker,
Thomas Ostler, Andreas Biternas, Roy W Chantrell,
Wu Hong-Ye, Razvan Ababei, Sam Westmoreland,
Milton Persson
Compiler Flags: MPI
Vampire includes a copy of the qhull library from C.B. Barber
and The
Geometry Center and may be obtained via http from www.qhull.org.
================================================================================
Thu Jul 2 16:00:56 2026
Processor 1 of 4 online on host computername
================================================================================
Initialising system variables
Creating system
System decomposed into 2 x 2 x 1 processors for parallel execution
Copying halo atoms to other processors...done!
Generating neighbour list..........done!
Initialising MPI communications...done!
Copying system data to optimised data structures.
Using generic/normalised form of exchange interaction with 16
total interactions.
Total number of atoms (all CPUs): 68600
Initialising dipole field calculation using macrocell solver
Dipole field calculation requires 12.7008 MB of RAM
Total memory for dipole calculation (all CPUs): 84.0158 MB of RAM
Precalculating rij matrix for dipole calculation using macrocell
solver..........done! [ 0.0188029 s ]
Time for initialisation: 0.435503
Starting Simulation with Program effective-damping...
Simulation run time [s]: 2790.9
Simulation ended gracefully.
username@computername:~/effective-damping$ gnuplot --persist plotm
[1] https://www-users.york.ac.uk/~rfle500/resources/paper_70.pdf
Kind Regards,
GavinVAMPIRE user
PHYSICAL REVIEW APPLIED 11, 044001 (2019)In this paper, it was mentioned that " The system is first equilibrated with a 10 T applied field along the z axis. A uniform rotation of the thermally equilibrated spin configuration is then made to initialize the magnetization to an angle of 30◦ in the y-z plane."
My doubt is1. Do we need to first run a time-series program with temperature 300, and 10T magnetic field applied in Z direction, and save the checkpoint using sim:save-checkpoint=end?
2. Now the spin configuration is thermally equilibrated, and we run the second program which is effective damping." A uniform rotation of the thermally equilibrated spin configuration is then made to initialize the magnetization to an angle of 30◦ in the y-z plane" how do we do this? just by changing the initial spin configuration in the material file? then what's the point of checkpoint?
--Vinod Naik Bhukya 在 2026年7月1日 星期三下午6:43:28 [UTC+8] 的信中寫道：
Dear vampire users,
Please help me with effective damping program. If someone has a input file on how to simulate for effective damping , it would be appreciable.My doubt is do we need to use both 'time-series' as well as 'effective-damping' programs and what are all the simulation attributes we need to mention in the input file?
Thank you.
You received this message because you are subscribed to the Google Groups "Vampire Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to vampire-user...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/vampire-users/8a3afe83-7009-4fe1-9e04-5a3307d14d9bn%40googlegroups.com.
Vinod Naik Bhukya
Best Regards,
Vinod (文诺德)
-------------------------------------------------------------------------------------
Vinod Naik Bhukya
Student Office: EF751-B
Dual Degree Ph.D. student (Indian Institute of Technology Kanpur)
International College of Semiconductor Technology (ICST)
National Yang-Ming Chiao-Tung University (NYCU), Hsinchu, Taiwan.
Mobile: +886 972126345
To view this discussion visit https://groups.google.com/d/msgid/vampire-users/f8e85e70-cc90-40b4-ad9e-1eca16204ad7%40gmail.com.