1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
/* -*- c++ -*- ----------------------------------------------------------
LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
http://lammps.sandia.gov, Sandia National Laboratories
Steve Plimpton, sjplimp@sandia.gov
Copyright (2003) Sandia Corporation. Under the terms of Contract
DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
certain rights in this software. This software is distributed under
the GNU General Public License.
See the README file in the top-level LAMMPS directory.
------------------------------------------------------------------------- */
#ifdef COMMAND_CLASS
CommandStyle(create_bonds,CreateBonds)
#else
#ifndef LMP_CREATE_BONDS_H
#define LMP_CREATE_BONDS_H
#include "pointers.h"
namespace LAMMPS_NS {
class CreateBonds : protected Pointers {
public:
CreateBonds(class LAMMPS *);
void command(int, char **);
private:
int igroup,group1bit,group2bit;
int btype,atype,dtype;
tagint batom1,batom2,aatom1,aatom2,aatom3,datom1,datom2,datom3,datom4;
double rmin,rmax;
void many();
void single_bond();
void single_angle();
void single_dihedral();
};
}
#endif
#endif
/* ERROR/WARNING messages:
E: Create_bonds command before simulation box is defined
Self-explanatory.
E: Cannot use create_bonds unless atoms have IDs
This command requires a mapping from global atom IDs to local atoms,
but the atoms that have been defined have no IDs.
E: Cannot use create_bonds with non-molecular system
Self-explanatory.
E: Illegal ... command
Self-explanatory. Check the input script syntax and compare to the
documentation for the command. You can use -echo screen as a
command-line option when running LAMMPS to see the offending line.
E: Cannot find create_bonds group ID
Self-explanatory.
E: Invalid bond type in create_bonds command
Self-explanatory.
E: Cannot use special no with create_bonds many
UNDOCUMENTED
E: Invalid angle type in create_bonds command
UNDOCUMENTED
E: Invalid dihedral type in create_bonds command
UNDOCUMENTED
E: Create_bonds requires a pair style be defined
Self-explanatory.
E: Create_bonds max distance > neighbor cutoff
Can only create bonds for atom pairs that will be in neighbor list.
W: Create_bonds max distance > minimum neighbor cutoff
This means atom pairs for some atom types may not be in the neighbor
list and thus no bond can be created between them.
E: Create_bonds command requires special_bonds 1-2 weights be 0.0
This is so that atom pairs that are already bonded to not appear in
the neighbor list.
E: Create_bonds command requires no kspace_style be defined
This is so that atom pairs that are already bonded to not appear
in the neighbor list.
E: New bond exceeded bonds per atom in create_bonds
See the read_data command for info on setting the "extra bond per
atom" header value to allow for additional bonds to be formed.
E: Create_bonds single/bond atoms do not exist
UNDOCUMENTED
E: Create_bonds single/angle atoms do not exist
UNDOCUMENTED
E: New angle exceeded angles per atom in create_bonds
UNDOCUMENTED
E: Create_bonds single/dihedral atoms do not exist
UNDOCUMENTED
E: New dihedral exceeded dihedrals per atom in create_bonds
UNDOCUMENTED
*/