1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
(*
Copyright (C) 2017 M.A.L. Marques
This Source Code Form is subject to the terms of the Mozilla Public
License, v. 2.0. If a copy of the MPL was not distributed with this
file, You can obtain one at http://mozilla.org/MPL/2.0/.
*)
XX := (z, xs) -> xs*opz_pow_n(z,4/3)*2^(-4/3):
YY := (z, xt, xs0, xs1) -> 2*(XX(z, xs0)^2 + XX(-z, xs1)^2) - xt^2:
f_th := (rs, z, xt, xs0, xs1) -> add(params_a_omega[i]
* (n_spin(rs, z)^params_a_a[i] + n_spin(rs, -z)^params_a_a[i])
* z^(2*params_a_b[i])
* 1/2*(XX(z, xs0)^params_a_c[i] + XX(-z, xs1)^params_a_c[i])
* YY(z, xt, xs0, xs1)^params_a_d[i], i=1..params_a_n)/n_total(rs):
|
