Next: , Previous: PBE WRITE Command, Up: PBE Commands


8.3.6 PBE TEST Command

Syntax

     PBE TEST { CV cv-options }  repeat(pbe-test-options)
              { DH dh-options }
     
     cv-options ::= [CHARge real    ]
                    [RADIus real    ]
                    [CAVIty-eps real]
                    [QPOS real      ]
     
                    [CHARge real      ]
     dh-options ::= [RADIus real      ]
                    [CAVIty-eps real  ]
                    [IONStr real      ]
                    [TEMPERATURE real ]
     
                          [XDIM int    ]
                          [YDIM int    ]
     pbe-test-options ::= [ZDIM int    ]
                          [GRID real   ]
                          [SOLVent real]

Function

The PBE TEST command is used to generate test potentials for comparing the finite difference solution against two particular problems for which analytic solutions are available. The two problems are a charge in a spherical cavity of one dielectric with a different surrounding dielectric (the CV option), and a charged centered in a sphere surrounded by a Debye-Huckel fluid (the DH option). These calculations are done using analytic series representations, but they are not optimized for numerical accuracy, so beware of artefactual peaks. CONGEN will report the number of failures of the series convergence. Also, a setting of 3 for the PBE DEBUG variable will display the terms in each series calculation, see Debug Command. A debug setting greater than 10 will cause CONGEN to substitute the number of iterations in place of the potential.

The analytic potentials used in this code were derived by Malcolm Davis.

The command options have following meaning:

Option
Purpose
CV
Specifies a potential for a spherical cavity with an internal charge with an external dielectric be computed.
DH
Specifies a potential for a spherical cavity with a charge at the center surrounded by a fluid with Debye-Huckel charge screening.
CHARGE
Specifies the interior charge value. The default is 1.0.
RADIUS
Specifies the cavity radius. The default is 5.0 Angstroms.
CAVITY-EPS
Specifies the dielectric of the cavity. The default is 1.0.
QPOS
Specifies the position of the charge in the cavity test case along the X axis. The default is 0.0 (the origin).
IONSTR
Specifies the ionic strength in molarity units for the Debye-Huckel test case. The default is the current value in the pbe_info data structure.
TEMPERATURE
Specifies the temperature for the Debye-Huckel test case. The default is the current value in the pbe_info data structure.
XDIM
YDIM
ZDIM
Specifies the X, Y, and Z dimensions of the test grid, respectively. If there is no grid defined in the pbe_info data structure, then the default for all three dimensions is 20. Otherwise, the current grid dimensions provide defaults.
GRID
Specifies the grid dimension. The current pbe_info data structure provides the default.
SOLVENT
Specifies the solvent (external) dielectric constant. The current pbe_info data structure provides the default.