ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
bmhtf-nacl.cpp
Go to the documentation of this file.
1
/*
2
* Copyright (C) 2010-2022 The ESPResSo project
3
* Copyright (C) 2002,2003,2004,2005,2006,2007,2008,2009,2010
4
* Max-Planck-Institute for Polymer Research, Theory Group
5
*
6
* This file is part of ESPResSo.
7
*
8
* ESPResSo is free software: you can redistribute it and/or modify
9
* it under the terms of the GNU General Public License as published by
10
* the Free Software Foundation, either version 3 of the License, or
11
* (at your option) any later version.
12
*
13
* ESPResSo is distributed in the hope that it will be useful,
14
* but WITHOUT ANY WARRANTY; without even the implied warranty of
15
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16
* GNU General Public License for more details.
17
*
18
* You should have received a copy of the GNU General Public License
19
* along with this program. If not, see <http://www.gnu.org/licenses/>.
20
*/
21
/** \file
22
*
23
* Implementation of \ref bmhtf-nacl.hpp
24
*/
25
#include "
bmhtf-nacl.hpp
"
26
27
#ifdef BMHTF_NACL
28
#include "
nonbonded_interactions/nonbonded_interaction_data.hpp
"
29
30
#include <
utils/math/int_pow.hpp
>
31
32
#include <stdexcept>
33
34
BMHTF_Parameters::BMHTF_Parameters
(
double
a,
double
b,
double
c,
double
d,
35
double
sig,
double
cutoff)
36
: A{a}, B{b}, C{c}, D{d}, sig{sig}, cut{cutoff} {
37
if
(a < 0.) {
38
throw
std::domain_error(
"BMHTF parameter 'a' has to be >= 0"
);
39
}
40
if
(c < 0.) {
41
throw
std::domain_error(
"BMHTF parameter 'c' has to be >= 0"
);
42
}
43
if
(d < 0.) {
44
throw
std::domain_error(
"BMHTF parameter 'd' has to be >= 0"
);
45
}
46
if
(cutoff < 0.) {
47
throw
std::domain_error(
"BMHTF parameter 'cutoff' has to be >= 0"
);
48
}
49
computed_shift
=
C
/ Utils::int_pow<6>(
cut
) +
D
/ Utils::int_pow<8>(
cut
) -
50
A
* std::exp(
B
* (
sig
-
cut
));
51
}
52
53
#endif
// BMHTF_NACL
bmhtf-nacl.hpp
Routines to calculate the Born-Meyer-Huggins-Tosi-Fumi potential between particle pairs.
int_pow.hpp
nonbonded_interaction_data.hpp
Various procedures concerning interactions between particles.
BMHTF_Parameters::A
double A
Definition
nonbonded_interaction_data.hpp:143
BMHTF_Parameters::C
double C
Definition
nonbonded_interaction_data.hpp:145
BMHTF_Parameters::D
double D
Definition
nonbonded_interaction_data.hpp:146
BMHTF_Parameters::cut
double cut
Definition
nonbonded_interaction_data.hpp:148
BMHTF_Parameters::BMHTF_Parameters
BMHTF_Parameters()=default
BMHTF_Parameters::sig
double sig
Definition
nonbonded_interaction_data.hpp:147
BMHTF_Parameters::B
double B
Definition
nonbonded_interaction_data.hpp:144
BMHTF_Parameters::computed_shift
double computed_shift
Definition
nonbonded_interaction_data.hpp:149
src
core
nonbonded_interactions
bmhtf-nacl.cpp
Generated on Fri Nov 22 2024 02:13:54 for ESPResSo by
1.9.8