ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
morse.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 morse.hpp
24
*/
25
#include "
morse.hpp
"
26
27
#ifdef MORSE
28
#include "
nonbonded_interaction_data.hpp
"
29
30
#include <cmath>
31
#include <stdexcept>
32
33
Morse_Parameters::Morse_Parameters
(
double
eps,
double
alpha,
double
rmin,
34
double
cutoff)
35
: eps{eps}, alpha{alpha}, rmin{rmin}, cut{cutoff} {
36
if
(
eps
< 0.) {
37
throw
std::domain_error(
"Morse parameter 'eps' has to be >= 0"
);
38
}
39
if
(cutoff < 0.) {
40
throw
std::domain_error(
"Morse parameter 'cutoff' has to be >= 0"
);
41
}
42
auto
const
add1 = std::exp(-2.0 *
alpha
* (
cut
-
rmin
));
43
auto
const
add2 = 2.0 * std::exp(-
alpha
* (
cut
-
rmin
));
44
rest
=
eps
* (add1 - add2);
45
}
46
47
#endif
// MORSE
morse.hpp
Routines to calculate the Morse potential between particle pairs.
nonbonded_interaction_data.hpp
Various procedures concerning interactions between particles.
Morse_Parameters::rest
double rest
Definition
nonbonded_interaction_data.hpp:162
Morse_Parameters::cut
double cut
Definition
nonbonded_interaction_data.hpp:161
Morse_Parameters::rmin
double rmin
Definition
nonbonded_interaction_data.hpp:160
Morse_Parameters::Morse_Parameters
Morse_Parameters()=default
Morse_Parameters::eps
double eps
Definition
nonbonded_interaction_data.hpp:158
Morse_Parameters::alpha
double alpha
Definition
nonbonded_interaction_data.hpp:159
src
core
nonbonded_interactions
morse.cpp
Generated on Fri Nov 8 2024 02:12:53 for ESPResSo by
1.9.8