Need to check implementation for calculating Q+,Q- in a triangular spin-lattice. The spin-texture I attach should have an homogenous topological charge density, while the actual implementation estimate Q-=-Q+, resulting in zero Q
n_Tsteps 100 #number of MC simulations at different T, uses Tini and Tfin
Total_MC_Steps 2000 #total number of MC steps per temperature
T_relax 50 #number of thermalization steps (will be T_relax*N)
T_auto 5 #measure steps occur every T_auto autocorrelation steps
tracker .F. #tracks skyrmions
ising .F. #Ising model: only up and down orientations
sphere_samp .T.
calc_topo .T.
do_fluct .F. #for fluctuations
do_fluct_spatial F
fluct_direction 1.0 0.0 0.0 #direction of first neighbour to consider when computing autocorrelations
Need to check implementation for calculating Q+,Q- in a triangular spin-lattice. The spin-texture I attach should have an homogenous topological charge density, while the actual implementation estimate Q-=-Q+, resulting in zero Q
git-version: '533aa548cced89b5dd0332b0a89a0958187ec5a9'
simulation metropolis
c_H_Exchten 1.0
H_ext 0.0 0.0 0.0 E_ext 0.0 0.0 0.00
Tini 100.0d0 Tfin 1.0d0
n_Tsteps 100 #number of MC simulations at different T, uses Tini and Tfin
Total_MC_Steps 2000 #total number of MC steps per temperature
T_relax 50 #number of thermalization steps (will be T_relax*N)
T_auto 5 #measure steps occur every T_auto autocorrelation steps
tracker .F. #tracks skyrmions
ising .F. #Ising model: only up and down orientations
sphere_samp .T.
calc_topo .T.
do_fluct .F. #for fluctuations
do_fluct_spatial F
fluct_direction 1.0 0.0 0.0 #direction of first neighbour to consider when computing autocorrelations
sym_mode 2
magnetic_r2_tensor 1 1 1 -4.0d-3 0.0 0.0 0.0 -2.8d-3 -0.7d-3 0.0 -0.7d-3 -3.65d-3 1.0 0.0 0.0 1 1 2 -0.15d-3 0.0 0.0 0.0 -0.2d-3 0.0 0.0 0.0 -0.025d-3 0.0 1.0 0.0 1 1 3 2.9d-3 0.0 0.0 0.0 2.9d-3 0.0 0.0 0.0 2.9d-3 1.0 0.0 0.0
magnetic_anisotropy 1 0 0 1 0.3d-3
Periodic_log .T. .T. .F. gra_fft .F. gra_log .F. gra_freq 10 # parameters of the lattice
Nsize 24 24 1
alat 0.396 0.396 1.0
lattice 1.0 0.0 0.0 -0.5 0.86602541212 0.0 0.0 0.0 1.0
atoms 3
Ni 0.0 0.0 0.0
I 0.333333335 0.666666669 -0.07349
I 0.666666620 0.333333300 0.07349
atomtypes 2
Ni 2.0 0.0 0.0 .F. #input for the dynamic
I 0.0 0.0 0.0 .F.