4.1. Using Amber/APBS module

The APBS module in sander offers an alternative to the built-in PB sander module. The APBS module can be used for example for implicit solvent minimization and dynamics, calculation and visualization of miscellaneous electrostatic biomolecular properties. Please see the Examples section of this User's Guide.

The following table lists all Amber/APBS keywords with their description. The left side of the table also lists corresponding APBS keywords which Amber/APBS keywords mimic very closely. For detailed discussion of APBS keywords please see APBS documentation.

4.1.1. Amber/APBS Module Keywords

Table 4-1. APBS and Amber/APBS keywords

APBS keywordAmber/APBS keywordDescription
 apbsEnters the APBS module
mg-auto/mg-paracalc_type [0]0: manual MG; 1: autoMG; 2: parallel MG
lpbe/nbpenonlin [0]Linear/full Poisson-Boltzmann equation: 0: linear; 1: non-linear; 4: size-dependent PBE
bcflbcfl [1]Boundary condition method: 0: zero; 1: sdh; 2: mdh; 4: focus
srfmsrfm [1]Surface calculation method: 0: mol; 1: smol; 2: spl2; 3: spl4
pdiepdie [2.0]Solute dielectric
sdiesdie [78.4]Solvent dielectric
sdenssdens [10.0]Vacc sphere density
sradsrad [1.4]Solvent radius
swinswin [0.3]Cubic spline window
temptemp [298.15]Temperature (in K)
gammagamma [0.105]Surface tension for apolar energies/forces (in kJ/mol/A^2)
chgmchgm [1]Charge discretization method: 0: spl0; 1: spl2; 2: spl4
volsmvolumeThe parameter smvolume controls the lattice size (in Angstroms^3) used in the SMPBE formalism.
sizesmsizeThe parameter smsize controls the relative size of the ions (in Angstroms) such that each lattice site can contain a single ion of volume radius^3 or size ions of volume radius^3/size.
calcenergycalcenergy [1]Energy calculation flag: 0: Do not perform energy calculation; 1: Calculate total energy only; 2: Calculate per-atom energy components
calcforcecalcforce [0]Atomic forces calculation: 0: Do not perform force calculation; 1: Calculate total force only; 2: Calculate per-atom force components
write potwpotWrites electrostatic potential data to iapbs-pot.dx in DX format
write chargewchgWrites charge data to iapbs-charge.dx in DX format
write smolwsmolWrites molecular surface data to iapbs-smol.dx in DX format
write kappawkappaWrites the ion-accessibility kappa map to iapbs-kappa.dx in DX format
write dielwdielWrites dielectric maps to iapbs-diel[x,y,z].dx in DX format
read chargerchgReads charge data from iapbs-charge.dx in DX format
read kapparkappaReads the ion-accessibility kappa map from iapbs-kappa.dx in DX format
read dielrdielReads dielectric maps from iapbs-diel[x,y,z].dx in DX format
 nionNumber of counterions
ionqionqCounterion charges (in e)
ioncioncCounterion concentrations (in M)
ionrionrrCounterion radii (in A)
dimedimeGrid dimensions (in x, y and z)
cmethcmeth [1]Centering method: 0: Center on a point 1: Center on a molecule
gcentcenterGrid center if cmeth=0
ccmethccmeth [1]Coarse grid centering method: 0: Center on a point 1: Center on a molecule
cgcentccenterCoarse grid center if ccmeth=0
fcmethfcmeth [1]Fine grid centering method: 0: Center on a point 1: Center on a molecule
fgcentfcenterFine grid center if fcmeth=0
gridgridGrid spacings
glenglenGrid side lengths
cglencglenCoarse grid side lengths
fglenfglenFine grid side lengths
pdimepdimeGrid of processors to be used in calculation
ofracofrac [0.1]Overlap fraction between processors
 radiopt [0]Radii optimization: 0: use prmtop; 1: read charge and radius information from 'pbparamsin' file (format: atom_label charge radius); 2: Read charge/radius information for a PQR file; 3: Read radius information from a PQR file.
 pqrPQR file name
 calcnpenergy [1]Calculate nonpolar energy [0, 1]. 0: Don't calculate; 1: SASA-based apolar energy caclulation.
 sp_apbs [.FALSE.]Perform SP APBS calculation

Note

Values in square brackets [] are defaults.

4.1.2. Amber/APBS Notes

See note about disabling creation of io.mc file below before attempting any extensive minimization or MD simulation with the APBS module.