1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
(*
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/.
*)
(* The B97 function g *)
b97_g := (gamma, cc, x) -> add(cc[i]*(gamma*x^2/(1 + gamma*x^2))^(i-1), i=1..5):
(* The parallel and perpendicular components of the energy *)
b97_fpar := (lda_func, mgamma, cc, rs, z, xs0, xs1) ->
+ lda_stoll_par(lda_func, rs, z) * b97_g(mgamma, cc, xs0)
+ lda_stoll_par(lda_func, rs, -z) * b97_g(mgamma, cc, xs1):
b97_fperp := (lda_func, mgamma, cc, rs, z, xs0, xs1) ->
lda_stoll_perp(lda_func, rs, z) * b97_g(mgamma, cc, sqrt(xs0^2 + xs1^2)/sqrt(2)):
b97_f := (lda_func, gamma_ss, cc_ss, gamma_ab, cc_ab, rs, z, xs0, xs1) ->
+ b97_fpar (lda_func, gamma_ss, cc_ss, rs, z, xs0, xs1)
+ b97_fperp(lda_func, gamma_ab, cc_ab, rs, z, xs0, xs1):
|
