## apbs-users

 [apbs-users] large potential values From: Nikolay Rodionov - 2012-09-09 05:58:58 Attachments: text/plain
 Re: [apbs-users] large potential values From: Baker, Nathan - 2012-09-09 20:49:14 ```Hi Nikolay -- I'm sorry but it is not possible to debug this based on the information you have provided. What is the solute? What is the basis for your by-hand calculation? Thanks, __________________________________________________ Nathan Baker Pacific Northwest National Laboratory Tel: 509-375-3997 http://nabaker.me; ________________________________ From: Nikolay Rodionov [nikolayrodionov387@...] Sent: Saturday, September 08, 2012 10:58 PM To: apbs-users@... Subject: [apbs-users] large potential values Could you look at my input file and see if there are any errors? I am trying to calculate the ESP of the system at 298K and .6M/10e-7 NaCl solution. The mean surface potential based on APBS is almost exactly an order of magnitude greater than what I get from calculating it by hand using: Charge_dens = sqrt(8 kB T Co Er Eo) * sinh( (e*zi*SurfPot) / (2 kB T) ), simple solution to PBE with a symmetric ion ....(I can solve for charge density because of the symmetric quaternary protein structure, and then use it to solve for the surface potential) APBS give me a values of -3.4e3 mV and -5.8e2 mV for the two molarities, while my hand calculations give me -3.9e2 mV and 6.1e1 mV I've checked my hand calculations and everything is correct so I am guessing that there is something wrong with my APBS input. Here is my .6M input file: read mol pqr cmplx.pqr end elec mg-auto dime 129 129 129 cglen 246.5325050354004 246.43500137329102 264.52200508117676 fglen 246.5325050354004 246.43500137329102 264.52200508117676 cgcent mol 1 fgcent mol 1 mol 1 lpbe bcfl sdh pdie 4 sdie 70.9 srfm smol chgm spl2 ion charge 1 conc .6 radius 1.16 ion charge -1 conc .6 radius 1.67 sdens 10.0 srad 1.4 swin 0.3 temp 298 calcenergy no calcforce no write pot dx SeaWater-pot write charge dx SeaWater-chargedist write atompot flat SeaWater-atompot end quit Thanks! Nikolay ```
 Re: [apbs-users] large potential values From: Nikolay Rodionov - 2012-09-10 14:11:46 Attachments: text/plain
 Re: [apbs-users] large potential values From: Baker, Nathan - 2012-09-12 02:33:55 Attachments: Message as HTML ```Hello – There are no blatant errors. However, the theory that you refer to in the link is an asymptotic approximation to the Poisson-Boltzmann model that APBS uses. There’s no reason to expect agreement between the analytic solution and the numerical solution when they come from different models. Thanks, __________________________________________________ Nathan Baker Pacific Northwest National Laboratory Tel: +1-509-375-3997 http://nabaker.me; From: Nikolay Rodionov [mailto:nikolayrodionov387@...] Sent: Monday, September 10, 2012 7:12 AM To: apbs-users@... Subject: Re: [apbs-users] large potential values More or less, I am curious if anything blatantly stands out as incorrect in the input file. My solute is a 14nm long particle of the tobacco mosaic virus. There is one thing really sparks my suspicions. When calculating for the surface potential on the external surface, the mean potential is positive at an extremely low molarity and negative at .6M. I don't understand how surface potential can be positive when surface charge density is negative. The basis of my calculations is eq127 of http://tfy.tkk.fi/kurssit/Tfy-3.363/lectures/lecture06a.pdf ... I solve for charge density in a region by summing the partial charges that make up the surface as according to Solvent Excluded Surface analysis and then dividing by an estimation of the surface area. Thanks, Nikolay On Sunday, September 9, 2012 1:58:51 AM UTC-4, Nikolay Rodionov wrote: Could you look at my input file and see if there are any errors? I am trying to calculate the ESP of the system at 298K and .6M/10e-7 NaCl solution. The mean surface potential based on APBS is almost exactly an order of magnitude greater than what I get from calculating it by hand using: Charge_dens = sqrt(8 kB T Co Er Eo) * sinh( (e*zi*SurfPot) / (2 kB T) ), simple solution to PBE with a symmetric ion ....(I can solve for charge density because of the symmetric quaternary protein structure, and then use it to solve for the surface potential) APBS give me a values of -3.4e3 mV and -5.8e2 mV for the two molarities, while my hand calculations give me -3.9e2 mV and 6.1e1 mV I've checked my hand calculations and everything is correct so I am guessing that there is something wrong with my APBS input. Here is my .6M input file: read mol pqr cmplx.pqr end elec mg-auto dime 129 129 129 cglen 246.5325050354004 246.43500137329102 264.52200508117676 fglen 246.5325050354004 246.43500137329102 264.52200508117676 cgcent mol 1 fgcent mol 1 mol 1 lpbe bcfl sdh pdie 4 sdie 70.9 srfm smol chgm spl2 ion charge 1 conc .6 radius 1.16 ion charge -1 conc .6 radius 1.67 sdens 10.0 srad 1.4 swin 0.3 temp 298 calcenergy no calcforce no write pot dx SeaWater-pot write charge dx SeaWater-chargedist write atompot flat SeaWater-atompot end quit Thanks! Nikolay ```