@@ -23,6 +23,30 @@ differences.
23232 ) To avoid spurious coupling between the physical and ghost systems, which
2424can 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 for ring and
32+ sp2 bridges, where local geometry already constrains the ghost and 90°
33+ stiffening would introduce significant strain.
34+ - ** Residual term cleanup:** a post-processing pass removes mixed improper
35+ dihedrals and cross-bridge dihedrals missed by the per-bridge junction
36+ handlers, as well as angles where a ghost atom is the central atom and both
37+ terminal atoms are physical.
38+ - ** Mixed dihedral softening:** surviving mixed ghost/physical dihedrals can
39+ be softened via ` --soften-anchors ` to allow ghost groups to reorient and
40+ avoid steric clashes at small λ.
41+ - ** Rotamer stiffening:** ` --stiffen-rotamers ` replaces rotatable sp3 anchor
42+ dihedrals with a stiff single-well cosine to control ghost orientation
43+ through flexible bonds.
44+ - ** Ring-breaking perturbations:** adjacent bridges with independent ghost
45+ groups retain each other as physical neighbours; angles with a ghost central
46+ atom spanning two physical neighbours are replaced by a linear spacer
47+ (180°, soft force constant); and angles spanning the ring-making/breaking
48+ bond are removed in the state where that bond is absent.
49+
2650Ghostly is incorporated into the [ SOMD2] ( https://github.com/openbiosim/somd2 )
2751free-energy perturbation engine.
2852
0 commit comments