/openCOSMO-RS_conformer_pipeline

Primary LanguagePythonGNU General Public License v2.0GPL-2.0

openCOSMO-RS_conformer_pipeline

Important for the execution are following requirements:

  • Working installation of ORCA 6.0+ (please make sure to download and extract all available compressed files for your OS)
  • The folder containing the ORCA installation has to be added to the path environment variable.
  • Install the xtb executable from the official repo releases. To do this copy the main binary from the 'bin' folder inside the compressed file into the ORCA installation directory and rename it to 'otool_xtb' or 'otool_xtb.exe' depending on your OS.
  • In case you want to use balloon to generate conformers instead of rdkit, you can get balloon here.
  • If you want to run parallelized calculations please check the ORCA manual to see how to setup this properly (installing MPI, setting paths, etc.)

The idea is that you can run the conformer generator as follows:

python ConformerGenerator.py --structures_file file.inp --cpcm_radii cpcm_radii.inp --n_cores 2

with file.inp beeing a TAB separated file similar to the following:

name [TAB] SMILES [TAB] optional xyz file [TAB] charge

Examples:

methane[TAB]C[TAB]xyz_file[TAB]0
ethanol[TAB]CCO[TAB][TAB]0

This is to calculate and plot the sigma profile:

from matplotlib import pyplot as plt
from input_parsers import SigmaProfileParser

spp = SigmaProfileParser(r'path.orcacosmo')
sigmas, areas = spp.cluster_and_create_sigma_profile()

plt.plot(sigmas, areas)
plt.xlim([-0.03, 0.03])
plt.savefig('sigma_profile.png')