Skip to content

Commit 8dabe75

Browse files
committed
Add summary of extended features.
1 parent ea45db4 commit 8dabe75

1 file changed

Lines changed: 25 additions & 0 deletions

File tree

README.md

Lines changed: 25 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -23,6 +23,31 @@ differences.
2323
2) To avoid spurious coupling between the physical and ghost systems, which
2424
can affect the equilibrium geometry of the physical system.
2525

26+
Ghostly implements many extensions beyond the original modification scheme to
27+
handle the diversity of perturbations encountered in practice:
28+
29+
- **Anchor selection scoring:** physical anchor atoms are scored to avoid
30+
transmuting or bridge atoms, preventing geometrically inconsistent constraints.
31+
- **Ring and sp2 bridge handling:** angle stiffening is skipped by default for
32+
ring and sp2 bridges, where local geometry already constrains the ghost and
33+
90° stiffening would introduce significant strain. It can be re-enabled via
34+
`--stiffen-ring-bridges` and `--stiffen-sp2-bridges`.
35+
- **Residual term cleanup:** a post-processing pass removes mixed improper
36+
dihedrals and cross-bridge dihedrals missed by the per-bridge junction
37+
handlers, as well as angles where a ghost atom is the central atom and both
38+
terminal atoms are physical.
39+
- **Mixed dihedral softening:** surviving mixed ghost/physical dihedrals can
40+
be softened via `--soften-anchors` to allow ghost groups to reorient and
41+
avoid steric clashes at small λ.
42+
- **Rotamer stiffening:** `--stiffen-rotamers` replaces rotatable sp3 anchor
43+
dihedrals with a stiff single-well cosine to control ghost orientation
44+
through flexible bonds.
45+
- **Ring-breaking perturbations:** adjacent bridges with independent ghost
46+
groups retain each other as physical neighbours; angles with a ghost central
47+
atom spanning two physical neighbours are replaced by a linear spacer
48+
(180°, soft force constant); and angles spanning the ring-making/breaking
49+
bond are removed in the state where that bond is absent.
50+
2651
Ghostly is incorporated into the [SOMD2](https://github.com/openbiosim/somd2)
2752
free-energy perturbation engine.
2853

0 commit comments

Comments
 (0)