Skip to content

Commit e12fb80

Browse files
committed
Edited the README.md file for the OPLSAA examples.
1 parent f2a1c33 commit e12fb80

File tree

1 file changed

+22
-20
lines changed
  • examples/all_atom/force_field_OPLSAA

1 file changed

+22
-20
lines changed

examples/all_atom/force_field_OPLSAA/README.md

Lines changed: 22 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -22,26 +22,6 @@ This will protect you from force-field parameter changes, and you
2222
will be able to continue using your existing atom and bonded type names safely.
2323

2424

25-
### Optional: Duplicate dihedrals, angles, and bonds
26-
27-
Sometimes, even after you have specified the (OPLSAA-specific) atom types
28-
for the atoms in your molecule, there may be multiple possible choices
29-
of dihedral, angle, or bond interactions between those atoms
30-
available in OPLSAA force field (stored in the "oplsaa.lt" file).
31-
When that happens, moltemplate.sh will *attempt to make a reasonable guess*,
32-
chosing the original (oldest, most common) version of the interaction between
33-
those atom types. However, you can override this choice:
34-
35-
- The new (2023) version of OPLSAA contains many additional choices for your dihedral, angle, and bond interactions. This gives you an opportunity to improve your simulation accuracy, but it also requires more effort on your part. To see the list of choices, you must now run moltemplate with the "-report-duplicates bytype __" arguments. For example:
36-
```
37-
moltemplate.sh system.lt -report-duplicates bytype __
38-
```
39-
- If you see a file named "warning_duplicate_dihedrals.txt", "warning_duplicate_angles.txt", "warning_duplicate_bonds.txt", or "warning_duplicate_impropers.txt" after running moltemplate, then it might be a good idea to read the first few warning messages
40-
in those files and modify your .lt files accordingly (for example, by adding a custom "Data Dihedrals" section). Several example .lt files demonstrate how to do that, including:
41-
- butane/moltemplate_files/butane.lt
42-
- benzene+benzoic_acid/moltemplate_files/benzoic_acid.lt
43-
44-
4525

4626

4727
### Atomic charges
@@ -55,6 +35,7 @@ LT files will be ignored.)*
5535
**These charges can be overridden.**
5636

5737

38+
5839
### Customizing atomic charges in OPLSAA molecules
5940

6041
#### Background information
@@ -108,6 +89,27 @@ of the [spce.lt](waterSPCE+methane/moltemplate_files/spce.lt) file.
10889

10990

11091

92+
### Optional: Duplicate dihedrals, angles, and bonds
93+
94+
Sometimes, even after you have specified the (OPLSAA-specific) atom types
95+
for the atoms in your molecule, there may be multiple possible choices
96+
of dihedral, angle, or bond interactions between those atoms
97+
available in OPLSAA force field (stored in the "oplsaa.lt" file).
98+
When that happens, moltemplate.sh will *attempt to make a reasonable guess*,
99+
chosing the original (oldest, most common) version of the interaction between
100+
those atom types. However, you can override this choice:
101+
102+
- The new (2023) version of OPLSAA contains many additional choices for your dihedral, angle, and bond interactions. This gives you an opportunity to improve your simulation accuracy, but it also requires more effort on your part. To see the list of choices, you must now run moltemplate with the "-report-duplicates bytype __" arguments. For example:
103+
```
104+
moltemplate.sh system.lt -report-duplicates bytype __
105+
```
106+
- If you see a file named "warning_duplicate_dihedrals.txt", "warning_duplicate_angles.txt", "warning_duplicate_bonds.txt", or "warning_duplicate_impropers.txt" after running moltemplate, then it might be a good idea to read the first few warning messages
107+
in those files and modify your .lt files accordingly (for example, by adding a custom "Data Dihedrals" section). Several example .lt files demonstrate how to do that, including:
108+
- butane/moltemplate_files/butane.lt
109+
- benzene+benzoic_acid/moltemplate_files/benzoic_acid.lt
110+
111+
112+
111113
### Minor issue: Improper angles
112114

113115
The style of improper interaction used by OPLS force fields depends on an angle which depends on the order of the atoms surrounding the central atom. When multiple atoms have the same type, this creates ambiguity in atom order. Since there is no guarantee that moltemplate will choose the same atom order as other molecule builders (such as VMD), this can lead to small unavoidable discrepancies in energies and forces computed by LAMMPS and NAMD. But their effect should be neglegible.

0 commit comments

Comments
 (0)