Topas tools

These are set of small scripts and tools that were developed to help with structure refinement (of zeolites in particular) using the program TOPAS.

topasdiff

Topasdiff is a tool to generate nice looking difference maps. First, output the observed structure factors, and structure in cif format from Topas:

Out_fobs(fobs.hkl)
Out_CIF_STR(structure.cif)

Observed structure factors can be generated using this macro (add this to C:/topas5/local.inc):

macro Out_fobs(file)
{
   phase_out file load out_record out_fmt out_eqn
   {
       "%4.0f" = H;
       "%4.0f" = K;
       "%4.0f" = L;
       '"%4.0f" = M;
       "%12.4f\n" = (Iobs_no_scale_pks / M)^(0.5);
   }
}

Then to calculate the difference map, run the command:

topasdiff structure.cif --diff fobs.hkl

or use the GUI, accessible available via topasdiff.bat (Windows). The program will ask for the scale factor provided by TOPAS, and generates an input file for the program superflip (http://superflip.fzu.cz/). Superflip is then used to perform the fourier transform, and generates an XPLOR file that can be visualized using Chimera or VESTA.

Note: There is a bug in Topas where the cif files it outputs cannot be read using CCTBX, because they lack the data header. A work-around is to open C:/topas5/topas.inc and modify all instances of:

Out_String("\ndata_")

by:

Out_String("\ndata_topas_cif_out")

cif2topas

cif2topas transforms a crystal structure into cif format into the corresponding TOPAS code.

Usage:

cif2topas structure.cif

fh2topas

fh2topas is a script that converts from Fenske-Hall z-matrix to TOPAS code

Usage:

fh2topas zmatrix.fh [n]

Here, n is an optional parameter to give the number of molecules to generate. fh2topas will automatically number them to avoid naming conflicts.

topasrestraints

topasrestraints is a tool that can help with the generation of bond and angle restraints for structure refinement of zeolites with TOPAS. First, generate all bonds by using the TOPAS instruction:

append_bond_lengths

to generate all bonds and angles, including their symmetry codes for the current structure. Copy everything between the curly brackets to a new file called bonds.txt, and run:

topasrestraints bonds.txt

This generates a file called restraints.out that contains the restraints that can be added to TOPAS. The script checks for all distances of 1.61 +- 0.4 Angstrom to identify T-O bonds. The program checks the connectivity for every atom, and reports if there is a problem (too many / not enough distances per T-atom). There is no check for symmetrically equivalent connections, so some restraints may be redundant.

Requirements

• Python==2.7
• numpy==1.10
• matplotlib<3.0
• pandas==0.23
• CCTBX
• superflip (superflip.fzu.cz/)

Installation

1. Install miniconda (pick the one suitable for your platform, Python version does not matter here).

2. Download and extract: https://github.com/stefsmeets/topas_tools/archive/master.zip

3. Install using:

    conda env create -f environment.yml
    conda activate topas_tools
    python setup.py develop

(note that every time you want to use topas_tools, you must activate the environment using conda activate topas_tools)

4. Follow instructions for your platform to install CCTBX below.

Windows

Download and unzip cctbx_mini-0.1.0-x64.zip (patched version for windows)

Go into the directory and run python setup.py install

Linux/MacOS

Download and install the latest CCTBX build from here: cci.lbl.gov/cctbx_build/

Before running the programs listed here, you must run cctbx_env.sh / cctbx_env.csh to ensure CCTBX modules can be found.

If you use bash:

source /usr/local/cctbx-dev-xxx/cctbx_env.sh

If you use tcsh/csh:

source /usr/local/cctbx-dev-xxx/cctbx_env.sh