A Python tool to deconstruct MOFs into building units and compute porosity. The code removes guests from MOFs and all porous systems, computes SMILES strings and InChIKeys of all building units.
Project description
mofstructure
This is an elaborate python module that provides simple functions for manipulation metal-organic frameworks and other porous systems such as COFs and Zeolites. Some uses of the module involves
-
Computation of geometric properties of MOFs. It calls zeo++ in the background and enables a quick computation of all porosity information such pld, lcd, asa ...
-
Automated removal of unbound guest molecules
-
Deconstruction of metal-organic frameworks into building units. And for each buidling units (organic ligand, metal cluster, organic sbu and metal sbu) computes their cheminformatic identifiers such as SMILES strings, inchi and inchikey. It also identify the type of metal sbu and coordination number of central metal.
-
Wraps systems around unit cell so as to remove effect of pbc. This is often the case, when one tries to visualize the cif files or convert cif files into xyz, the system often appears to be uncoordinated
-
Seperation of building units into regions. This is essential when on wishes to subsitute a specific ligand or building unit.
Installation
Option 1
pip install mofstructure
Option 2
git clone https://github.com/bafgreat/mofstructure.git mofstructure
cd mofstructure
pip install .
Quick start
Run on the commandline
Building units
Simply run the following command on a cif file or any ase readable file format containing a MOF.
mofstructure cif_file
The script will deconstruct the MOF present in the cif file and load the output in a folder called 'MOF_building_units' in the current directory/folder. If you wish to load the output in a specific folder, simply add the path to the folder as follows:
mofstructure cif_file path_to_result_folder
For multiple cif files. Simply run a loop and all the Results will be saved in the
for cifs in ciffiles:
mofstructure cifs path_to_result
Creating a database
If you have a folder containg many cif files for different MOF, you could easily create a database. To create such a database, simply run the following command.
mofstructure_database ciffolder
Here the 'ciffolder' is the folder containing the cif files. The ouput will be saved in the default folder called 'MOFDb' in the current folder. Again you can choose the path to the save folder by simply listing it at the end of the command.
mofstructure_database ciffolder path_to_result
Use as a libray
from mofstructure import mofdeconstructor
from mofstructure import porosity
from mofstructure import buildingunits
from ase.io import read, write
import pandas as pd
# Read cif file using ase
ase_atom = read(cif_file)
# check and remove unbound guest molecules
no_guest_indices = mofdeconstructor.remove_unbound_guest(ase_atom)
no_guest_atom = ase_atom[no_guest_indices]
# compute porosit and write output to csv
pores = porosity.zeo_calculation(ase_ato, probe_radius=1.86, number_of_steps=5000)
df = pd.DataFrame(pores, index=[0])
df.to_csv('pore.csv')
sbus and linkers
Compute sbus and linkers
connected_components, atoms_indices_at_breaking_point, porpyrin_checker, all_regions = MOF_deconstructor.secondary_building_units(ase_atom)
metal_sbus, organic_sbus, building_unit_regions = MOF_deconstructor.find_unique_building_units(
connected_components,
atoms_indices_at_breaking_point,
ase_atom,
porpyrin_checker,
all_regions,
cheminfo=True
)
when cheminfo = True
openbabel is called to compute all chemifomatic information, which are all stored on the ase_atom.info metal_sbus and organic_sbus list that contains all the unique instances of the metal sbus and organic sbus.
extracting cheminfor
For each instance in a building unit the various chemiformatic informations are as follows.
for i, sbu in enumerate(metal_sbu):
smi = sbu.info['smi']
inchi = sbu.info['inchi']
inchikey = sbu.info['inchikey]
# for sbus only
number_of_point_of_extension = sbu.info['point_of_extension']
#for metal sbus only
sbu_type = sbu.info['sbu_type'] # sbu_type :rodlike, irmof, uoi66, paddlewheel e.t.c
# write
sbu.write('metal_sbu_'+str(i)+'.cif')
Summary
All of above codes can be run by a single function found mofstructure.buildingunits
for a single cif file
from mofstructure import buildingunits
buildingunits.work_flow(ciffile, result_folder)
for multiple cif files
from mofstructure import buildingunits
import glob
all_cifs = glob.glob(folder/*cif)
# folder corresponds to the folder containing all the cif files
for cif_files in all_cifs:
buildingunits.work_flow(cif_files, result_folder)
# Note that result_folder can be any path. If the path does not exist, it will create one and populate it with all the data.
Documentation
You can access the full project documentation on docs
Roadmap
In the future the code should be able to:
- Compute rcsr topological code
- Subsitutue building units in a MOF to enable framework functionalisation
- Automatic curation of cifs
- Decontsruction of COFs into their building units
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for mofstructure-0.1.8-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 2d98f88dc34fdcf6e2c709c33221d60d825b53aa0082c88c4f76f6831593c13e |
|
MD5 | e17603b59bf6f6a85eaadc2c2efdb061 |
|
BLAKE2b-256 | 20d92d715f11fd61f83d2f4254ba616bb4557b93f0d8314920a6df551b36a9cf |