Skip to content

Instantly share code, notes, and snippets.

@CatChenal
Last active November 17, 2024 22:47
Show Gist options
  • Save CatChenal/2b21a87423d8a8bfd1486f7d0705454c to your computer and use it in GitHub Desktop.
Save CatChenal/2b21a87423d8a8bfd1486f7d0705454c to your computer and use it in GitHub Desktop.

Tweaking the LJ Potential:

Unless I misunderstood the problem, this tweaking of the LJ parameters will be a tractable optimization problem once we know these two constraints:

  1. How much more favorable we want the new V to be?
  2. Over which additional tail range we want the favorable part to persist?

This is because there are numerous possible solutions.

Examples changing addtional with on sigma given n & m:

  • In these examples, the new V is still favorable at r=8:
  • Is the favorable difference sufficient, say at r=5, and do we have a handle on which extra energy is needed?

Fig.1 - New V(n.m) without additional width on sigma:
Fig.1

Fig.2 - New V(n.m) wit 0.15 additional width on sigma.
Fig.2

Fig.3 - New V(n.m) wit 0.25 additional width on sigma.
Fig.3

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment