Implementing the Particle Mesh Ewald method

Click For Summary
SUMMARY

The discussion centers on implementing the Particle Mesh Ewald (PME) method for molecular dynamics simulations, highlighting its algorithmic complexity of O(n log(n)), which is significantly more efficient than the O(n^2) complexity of traditional methods. The original paper detailing PME is referenced, but participants express difficulty in translating the complex mathematics into code. There is a request for simple implementations or pseudocode, as well as resources that clearly explain the algorithm.

PREREQUISITES
  • Understanding of molecular dynamics simulations
  • Familiarity with the Particle Mesh Ewald method
  • Basic knowledge of algorithmic complexity
  • Proficiency in a programming language suitable for scientific computing (e.g., Python, C++)
NEXT STEPS
  • Research existing implementations of the Particle Mesh Ewald method in libraries like GROMACS or LAMMPS
  • Study the mathematical foundations of the Particle Mesh Ewald method through the original paper and supplementary materials
  • Explore pseudocode examples available in scientific literature or online repositories
  • Learn about optimizing molecular dynamics simulations using advanced algorithms
USEFUL FOR

This discussion is beneficial for computational chemists, molecular dynamics researchers, and software developers working on simulation tools who are looking to implement or optimize the Particle Mesh Ewald method.

cyberdash
Messages
2
Reaction score
0
I'm currently attempting to implement the "Particle mesh ewald" method for a molecular dynamics simulation. It has the advantage of having an algorithmic complexity of O(n log(n)), compared to a normal O(n^2) complexity for a regular approach.

I've been reading through the original papers for this ( http://dx.doi.org/10.1063/1.464397 ) But I am finding it to be quite dense, and it's difficult to translate the mathematics into actual code.

I'm wondering if anybody knows of a simple implementation that I could have a look at, even pseudocode would be helpful. If that doesn't exist, any resource explaning the actual algorithm clearly, or an explanation from yourself if you understand it, would be greatly appreciated.

Thanks!
 
Technology news on Phys.org
Or maybe point me to the right place / person to ask?
 

Similar threads

  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 13 ·
Replies
13
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
14K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K