Gaussian Cube file reader

User interface of the Gaussian Cube file reader, which appears as part of a pipeline’s file source.

This file format is used by the Gaussian simulation package and other ab initio simulation codes. It stores the simulation cell geometry, the atomic structure, and volumetric grid data.

Specifications of the format can be found here and here, for example.

If the imported file uses Bohr units, OVITO converts atomic coordinates and cell dimensions to Angstroms automatically.

Volumetric grid type

The user interface of the file reader provides an option that controls how the volumetric data should be interpreted by OVITO. The field values may either be attributed to the grid line intersections (the default) or the cell centers:

The selected grid type affects operations subsequently performed in OVITO, e.g. constructing an iso-surface from the volumetric data. In all cases, the file reader assumes 3d periodic boundary conditions for the volumetric grid and the atomic simulation cell.

The option Convert density values from Bohr units to Angstroms controls whether the field values loaded from the file are assumed to represent a density given in \(\text{bohr}^{-3}\) units and require conversion to \(\text{Å}^{-3}\) (OVITO’s internal units) to account for the volume change. This option is enabled by default. You can disable it if you know that the field values in your file are already given in \(\text{Å}^{-3}\) units or do not represent a density at all.

See also

• VoxelGrid.grid_type

• Create isosurface modifier

Python parameters

The file reader accepts the following optional keyword parameters in a call to the import_file() or load() Python functions.

import_file(location, grid_type=VoxelGrid.GridType.PointData, convert_field_bohr_to_angstrom=True, generate_bonds=False)

Parameters:

• grid_type – Selects how OVITO should interpret the volumetric data loaded from the Cube file. See VoxelGrid.grid_type for further information.

• convert_field_bohr_to_angstrom (bool) – Controls whether field values in the file are assumed to be density values given in \(\text{bohr}^{-3}\) and require conversion to \(\text{Å}^{-3}\) (OVITO’s internal units). You can disable it if you know that the field values in your file are already given in \(\text{Å}^{-3}\) units or do not represent a density at all.

• generate_bonds (bool) – Activates the generation of ad-hoc bonds connecting the atoms loaded from the file. Ad-hoc bond generation is based on the van der Waals radii of the chemical elements. Alternatively, you can apply the CreateBondsModifier to the system after import, which provides more control over the generation of pair-wise bonds.