READ input file section¶
Note
This section has been ported to the new APBS syntax (see YAML- and JSON-format input files). See Data loading input file section (required) for more information.
The READ block of an APBS input file has the following general format:
READ
[ keywords... ]
END
where keywords
is or more of the keywords described below (the line breaks and indentation are for clarity; only whitespace is necessary).
Note
One of these sections must be present for every molecule involved in the APBS calculation. Molecule and “map” IDs are assigned implicitly assigned for each molecule/map read, based on order and starting at 1 and incremented independently for each input type. In other words, each input PQR file is assigned an ID 1, 2, 3, …; each input dielectric map is assigned an independent ID 1, 2, 3, …; etc.
charge¶
This command allows APBS to read the fixed (molecular) charge density function mapped to a mesh. The inputs are maps of charge densities; these values have units of ec Å-3, where ec is the electron charge. In general, this command will read charge-maps written by ELEC input file section write commands. The syntax of this command is:
READ charge {format} {path} END
format
Specify the format of the charge map. Acceptable values include:
dx
- OpenDX scalar data format
gz
- gzipped (zlib) compressed OpenDX scalar data format. Files can be read directly in compressed form.
path
- The location of the charge map file.
diel¶
This command allows APBS to read the dielectric function mapped to 3 meshes shifted by one-half grid spacing in the x, y, and z directions. The inputs are maps of dielectric variables between the solvent and biomolecular dielectric constants; these values are unitless. In general, this command will read dielectric maps written by by ELEC input file section write commands. The syntax of this command is:
READ diel {format} {path-x} {path-y} {path-z} END
format
The format of the dielectric map.
dx
- OpenDX scalar data format
gz
- gzipped (zlib) compressed OpenDX scalar data format. Files can be read directly in compressed form.
path-x
- The location of the x-shifted dielectric map file.
path-y
- The location of the y-shifted dielectric map file.
path-z
The location of the z-shifted dielectric map file.
Note
If you choose this option and have a non-zero ionic strength, you must also include a READ kappa statement.
kappa¶
This command allows APBS to read the ion-accessibility function mapped to a mesh. The inputs are maps of ion accessibility values which range between 0 and the build Debye-Hückel screening parameter; these values have units of Å-2. In general, this command will read kappa-maps written by by ELEC input file section write commands. The syntax of this command is:
READ kappa {format} {path} END
format
Specify the format of the charge map. Acceptable values include:
dx
- OpenDX scalar data format
gz
- gzipped (zlib) compressed OpenDX scalar data format. Files can be read directly in compressed form.
path
- The location of the map file.
Note
If you choose this option, you must also include a read diel statement.
mol¶
This command specifies the molecular data to be read into APBS. The syntax is
READ mol {format} {path} END
format
The format of the input data.
pqr
- Specify that molecular data is in PQR format.
pdb
- Specify that molecular data is in pseudo-PDB format. If this type of structure file is used, then a parameter file must also be specified with a READ parm statement to provide charge and radius parameters for the biomolecule’s atoms.
path
- The location of the molecular data file.
parm¶
This command specifies the charge and radius data to be used with pseudo-PDB-format molecule files. The syntax is:
READ parm {format} {path} END
format
The format of the parameter file.
flat
- Specify that the parameter file is in APBS flat-file parameter format.
xml
- Specify that the parameter file is in APBS XML parameter format
path
- The location of the parameter data file.
Note
APBS provides a few example files as part of the source code distribution. Currently, example files only contain the polar parameters that can also be assigned more easily through the PDB2PQR software.
pot¶
This command allows APBS to read the electrostatic potential mapped to a mesh. The inputs are maps of the electrostatic potential from a previous calculation. In general, this command will read potential-maps written by by ELEC input file section write commands. The syntax of this command is:
READ pot {format} {path} END
format
Specify the format of the charge map. Acceptable values include:
dx
- OpenDX scalar data format
gz
- gzipped (zlib) compressed OpenDX scalar data format. Files can be read directly in compressed form.
path
- The location of the map file.