Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Energy, etc. Tips + Needs #15

Closed
bkmi opened this issue Jan 17, 2018 · 2 comments
Closed

Energy, etc. Tips + Needs #15

bkmi opened this issue Jan 17, 2018 · 2 comments

Comments

@bkmi
Copy link
Owner

bkmi commented Jan 17, 2018

be able to switch energy from coulomb versus lenard jones energy calc

Don't use a quad deep look in k space. Build a grid of k vectors ahead of time and apply it. In addition, use symmetry to reduce the problem to only real space k vectors

And a more obtuse suggestion:
Take the equations from long range ewald, substitute structure factors, use Euler or symmetry... This process will yield coupled interactions between two particles (but still manage to avoid incorporating any images.) this notion of energy contributions from pairs of particles can be used for the forces. (Another small hint, it involves the complex square of the structure factor)

@zenpi-me
Copy link
Collaborator

zenpi-me commented Jan 30, 2018

force

First of all: Thank you for your suggestions and your time.

I am thinking that I can just take the energies and calc the negative gradient from them. I already did this and if you like you can see. Because I am not sure, if I can just do it like that. I actually just did "force = -np.diff(energy)/distance"...

The k vector grid is also implemented and I use symmetry. There is nothing imaginary anymore.

I don't understand your last suggestion :/ what images are you talking about?

@TracyMcBean
Copy link
Collaborator

I think the problem with the symmetry of the k vectors has been solved. In case this is not the case I would suggest opening a new issue which is up to date.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants