[gmx-users] Long range electrostatic potential
G.N.Patargias at leeds.ac.uk
Wed Oct 17 16:28:00 CEST 2007
Thank you both Jeroen and Mark for your replies.
From: gmx-users-bounces at gromacs.org on behalf of Mark Abraham
Sent: Wed 10/17/2007 2:57 PM
To: Discussion list for GROMACS users
Subject: RE: [gmx-users] Long range electrostatic potential
> As far as I know it can't be done for individual groups. My best guess
> is that it has something to do with the 'M' in 'PME'. :-)
Well... sort of...
The real-space charges have to get gridded (there's that M), Fourier
transformed, then the long-range term in PME is formed from truncated
summation in reciprocal space, then back-transformed and interpolated off
the M. I'd guess you could do a group-wise decomposition of this
algorithm, but that would require P^2
transform-sum-backtransform-interpolation stages for P groups. This would
make the scaling of group-wise-LR PME something like P^2*N*log(N), which
you'd never implement as your default algorithm, since N*log(N) is
available. Then maybe you'd stop and ask yourself whether these group-wise
interaction energies mean much anyway...
PS Afterthought, maybe it would be less than P^2*N*log(N), since N might
be smaller. If you care, read the PME papers and go figure :-)
gmx-users mailing list gmx-users at gromacs.org
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 3647 bytes
Desc: not available
More information about the gromacs.org_gmx-users