protected override double Spin_B(Parameters p, ERY.EMath.Vector3 k)
{
Complex Gn = G_n(k);
Complex Gm = G_m(k);
Complex result = 4 * (p.Jxy * Gm + p.Jxx * Gn);
return(result.Magnitude);
}
protected override double Orb_E(Parameters p, ERY.EMath.Vector3 k)
{
return(0);
}
protected override double Orb_D(Parameters p, ERY.EMath.Vector3 k)
{
return(0.5 * p.JI * Gamma_2(k));
}
protected override double Orb_B(Parameters p, ERY.EMath.Vector3 k)
{
return(Orb_A(p, k));
}
protected override double Orb_A(Parameters p, ERY.EMath.Vector3 k)
{
return(-p.Js);
}
protected override double Spin_Ad(Parameters p, ERY.EMath.Vector3 k)
{
return(0);
}
protected override double Spin_B(Parameters p, ERY.EMath.Vector3 k)
{
return(8 * p.Jxx * Gamma(k));
}
protected override double Spin_A(Parameters p, ERY.EMath.Vector3 k)
{
return(8 * p.Jxx + 0.5 * (p.Js - p.Jsigma * Math.Cos(k[0]) - p.Jpi * Math.Cos(k[1])));
}
protected override double Spin_B(Parameters p, ERY.EMath.Vector3 k)
{
return(4 * (p.Jxy * G_m(k) + p.Jxx * G_n(k)));
}
protected override double Spin_A(Parameters p, ERY.EMath.Vector3 k)
{
return(8 * p.Jbar);
}
