magres-format
Code for parsing the CCP-NC ab-initio magnetic resonance file format as used in the latest version of CASTEP, coming soon to other codes such as Quantum ESPRESSO. See more on this page
Documentation for the is available here.
A few IPython notebooks have been written using the library as examples. You can see them here.
IPython is an enhanced interpreter for Python and offers an excellent in-browser workbook experience, similar to Matlab or Mathematica. This is particularly useful when developing code using this library to process your magnetic resonance calculations. You can read instructions for installing it here.
A calculated J-coupling network in ethanol.
Installing
Clone the repository or download and extract the .zip file somewhere. From the command line, run:
sudo python setup.py install
to install it globally, or
python setup.py install --user
to install it just for your user account (read more here), and the Python module and associated scripts should now be installed.
If you have installed it locally with --user, you may have to add ~/.local/bin to your PATH. You can do this by adding
export PATH=$HOME/.local/bin:$PATH
to your ~/.bashrc and restarting your session or running "source ~/.bashrc". If you use tcsh you do
setenv PATH $HOME/.local/bin:$PATH
and then restart your session.
Extraction scripts
Some utility scripts for extracting values from a large number of calculation output files are provided. Look at their help information for detailed instructions
For magnetic shieldings (chemical shifts)
extract-ms.py --help
For electric field gradients (quadrupolar couplings)
extract-efg.py --help
For J-couplings (indirect spin-spin coupling)
extract-jc.py --help
These scripts can be called with an atom list to restrict which atoms or couplings are shown. These can select an entire species
H
will select all hydrogen atoms. They can select a single atom
H2
will select the second hydrogen atom. They select ranges of atoms
H2-5
will select the second to fifth hydrogen atoms. These can also be chained together with commas
H1-5,O
will select the first five hydrogen atoms and all oxygen atoms.
For example,
extract-ms.py . Zn
will print only the magnetic shieldings of zinc atoms in all the .magres
files found in the current directory.
The -N
flag optionally outputs in the first columns of the output an attempt at parsing out numbers in a path. This is useful for convergence tests. E.g. the path grid_scale=2/energy_cut_off=80/ethanol.magres
will output the numbers 2.0 and 80.0 in the first two output columns.
Conversion script usage
The convertoldmagres.py
script converts an old-style Castep magres file to the new-style CCP-NC format for use with the new tools. You use it from the command line like:
magres-convert.py sample.magres > sample.new.magres
and optionally with the associated job's .castep file, to capture the lattice information,
magres_convert.py sample.magres sample.castep > sample.new.magres
Python module usage
The magres.format
and magres.atoms
modules contain code for, respectively, a low level parser of the CCP-NC ab-initio magres format and a high-level collection of objects to represent its contents.
The module also include the useful magres.constants
module, which gives the best-known gamma constants and quadrupole moments for all isotopes, the most common isotopes used in experiments.
More documentation is available here. Also, see the IPython notebooks linked at the top of this document.
JSON schema
We use the JSONschema definition to provide a specification for the internal datastructure used by the parser and the format of the JSON emitted and consumed by .as_json() and .load_json() on MagresFile.
To dump the JSON representation of a .magres simply do
magresjson.py sample.magres > sample.magres.json
and sample.magres.json should now contain a schema-compliant JSON representation of sample.magres.