In [ ]:
from rmgpy.molecule import Molecule
from rmgpy.species import Species
from rmgpy.molecule.resonance import *
In [ ]:
import rmgpy
from rmgpy.data.rmg import RMGDatabase
databasePath = rmgpy.settings['database.directory']
database = RMGDatabase()
database.load(
path = databasePath,
thermoLibraries = ['Narayanaswamy'],
reactionLibraries = [],
seedMechanisms = [],
kineticsFamilies = 'none'
)
In [ ]:
mol = Molecule(SMILES='N=C[O]')
In [ ]:
mol
In [ ]:
out = mol.generateResonanceIsomers()
out
In [ ]:
out[0]
In [ ]:
out[1]
In [ ]:
out[0].getFormula()
In [ ]:
out[1].getFormula()
In [ ]:
mol2 = Molecule(SMILES='CCC[CH]c1ccccc1')
In [ ]:
mol2
In [ ]:
out2 = mol2.generateResonanceIsomers()
In [ ]:
spec = Species(molecule=[mol2])
spec.generateResonanceIsomers()
spec.molecule
In [ ]:
thermo = database.thermo.getThermoData(spec)
In [ ]:
thermo
In [ ]:
spec.molecule
In [ ]: