Skip to content

Latest commit

 

History

History
99 lines (66 loc) · 2.39 KB

readme.rst

File metadata and controls

99 lines (66 loc) · 2.39 KB
Raw

Psi Phi Package

Do you need rdkit? Do you need quantum chemistry? We got you bro.

Examples

Assume all codesnippets below are using RDKit molecule objs

molecule = Chem.MolFromSmiles("O")
Chem.AddHydrogens(molecule)
AllChem.UFFOptimizeMolecule(molecule)

Optimize using XTB

from ppqm import XtbCalculator
xc = XtbCalculator()
molecule2 = xc.optimize(molecule, return_copy=True)

Example of using GAMESS calculator and using specific options. As you notice, GAMESS needs a lot of settings to work with.

from ppqm import GamessCalculator

# Let's set the complicated GAMESS settings
gamess_options = {
    "scr": "/tmp/node/scr/space/slurm/id",
    "cmd": "/opt/gamess/rungms",
    "gamess_scr": "~/scr",
    "gamess_userscr": "~/userscr",
}
gc = GamessCalculator(**gamess_options)

# Now that we have gamess setup, we can then the GAMESS options we all know
# and love. Knowing exactly what keywords to set in GAMESS, you'll have to
# read the manual
calculation_option = {
    "runtyp": "optimize",
    "statpt": {
        "opttol": 0.005,
        "nstep": 300,
        "projct": False,
    }
}

# We then use the options to get properties for the molecule.
# The return will be a list of dictionaries, per conformer in the molobj.
results = gc.calculate(molecule, calculation_option)

for properties in results:
    print(properties)

Notes on Jupyter usage

Example notebooks uses nglview to visualize the molecules in notebooks

# install nglview
conda install nglview -c conda-forge
# install plugin for jupyter lab
jupyter labextension install nglview-js-widgets

Please note, if you are using Jupyter Lab (not notebook) there are som additional notes to consider

Contributions

Fork, branch and use pre-commit.

Other code bases

"is this the first python wrapper for quantum chemistry?" No, check the others and find the one right for your project. Know one, not on the list? Add it.