Skip to content

NEB force: energy calculation #124

Open

Description

@sniblett402

In line 382 of _NEB.py, the spring energy appears to be divided by k instead of multiplied. Is this correct? It seems counterintuitive, and seems to contradict equation (1) of http://dx.doi.org/10.1063/1.1636455

Activity

  1. sniblett402 commented on Feb 11, 2015

    @sniblett402
    ContributorAuthor

    I now see that pele isn't actually using that line: it's using the Fortran version of the same function (_NEB_utils.f90) where dividing by k makes perfect sense. I still don't understand why the python version divides by k, however.

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions