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