SunnySuite / Sunny.jl

Spin dynamics and generalization to SU(N) coherent states
Other
65 stars 19 forks source link

Restore scalar biquadratic in `mode = :dipole_large_S` #251

Closed kbarros closed 4 months ago

kbarros commented 4 months ago

This feature is not yet supported in set_pair_coupling!. This full solution, for general biquadratic interactions, would require to introduce new symbolic spin variables for pairs of sites (extending our use of DynamicPolynomials). As a workaround, we should get the special case of scalar biquadratic (temporarily) working again as a biquad keyword argument to set_exchange!. A low level function call is:

set_pair_coupling_aux!(sys, 0.0, J - biquad*I/2, biquad, zero(TensorDecomposition), bond)