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