From 40dbccc06a79fceaec4de7baa5d8298c35db2ae3 Mon Sep 17 00:00:00 2001 From: Kipton Barros Date: Tue, 18 Jul 2023 10:35:12 -0600 Subject: [PATCH] Update comment --- src/System/Interactions.jl | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/src/System/Interactions.jl b/src/System/Interactions.jl index ffc24faf3..87245551e 100644 --- a/src/System/Interactions.jl +++ b/src/System/Interactions.jl @@ -150,7 +150,7 @@ function local_energy_change(sys::System{N}, site, state::SpinState) where N if !iszero(coupling.biquad) J = coupling.biquad if sys.mode == :dipole - # Renormalization introduces a factor r and a Heisenberg term + # Renormalization defined in https://arxiv.org/abs/2304.03874. Sᵢ = (sys.Ns[site]-1)/2 Sⱼ = (sys.Ns[cellⱼ, bond.j]-1)/2 S = √(Sᵢ*Sⱼ) @@ -242,7 +242,7 @@ function energy_aux(sys::System{N}, ints::Interactions, i::Int, cells, foreachbo if !iszero(coupling.biquad) J = coupling.biquad if sys.mode == :dipole - # Renormalization introduces a factor r and a Heisenberg term + # Renormalization defined in https://arxiv.org/abs/2304.03874. Sᵢ = (sys.Ns[site1]-1)/2 Sⱼ = (sys.Ns[site2]-1)/2 S = √(Sᵢ*Sⱼ) @@ -316,8 +316,7 @@ function set_forces_aux!(B, dipoles::Array{Vec3, 4}, ints::Interactions, sys::Sy if !iszero(coupling.biquad) J = coupling.biquad if sys.mode == :dipole - # Renormalization procedure introduces a factor r and a - # Heisenberg term, https://arxiv.org/abs/2304.03874. + # Renormalization defined in https://arxiv.org/abs/2304.03874. Sᵢ = (sys.Ns[site1]-1)/2 Sⱼ = (sys.Ns[site2]-1)/2 S = √(Sᵢ*Sⱼ)