def test_spin_multiplicity_from_qchem_log(self):
"""
Uses a QChem log file for npropyl to test that its
molecular degrees of freedom can be properly read.
"""
log = QChemLog(os.path.join(os.path.dirname(__file__), 'data', 'npropyl.out'))
conformer, unscaled_frequencies = log.load_conformer()
self.assertEqual(conformer.spin_multiplicity, 2)
log = QChemLog(os.path.join(os.path.dirname(__file__), 'data', 'co.out'))
conformer, unscaled_frequencies = log.load_conformer()
self.assertEqual(conformer.spin_multiplicity, 1)
def test_load_vibrations_from_qchem_log(self):
"""
Uses a QChem log files to test that
molecular energies can be properly read.
"""
log = QChemLog(os.path.join(os.path.dirname(__file__), 'data', 'npropyl.out'))
conformer, unscaled_frequencies = log.load_conformer()
self.assertEqual(len(conformer.modes[2]._frequencies.value), 24)
self.assertEqual(conformer.modes[2]._frequencies.value[5], 881.79)
log = QChemLog(os.path.join(os.path.dirname(__file__), 'data', 'co.out'))
conformer, unscaled_frequencies = log.load_conformer()
self.assertEqual(len(conformer.modes[2]._frequencies.value), 1)
self.assertEqual(conformer.modes[2]._frequencies.value, 2253.16)
def test_load_npropyl_modes_from_qchem_log(self):
"""
Uses a QChem log file for npropyl to test that its
molecular modes can be properly read.
"""
log = QChemLog(
os.path.join(os.path.dirname(__file__), 'data', 'npropyl.out'))
conformer, unscaled_frequencies = log.load_conformer()
self.assertTrue(
len([
mode for mode in conformer.modes
if isinstance(mode, IdealGasTranslation)
]) == 1)
self.assertTrue(
len([
mode for mode in conformer.modes
if isinstance(mode, NonlinearRotor)
]) == 1)
self.assertTrue(
len([
mode for mode in conformer.modes
if isinstance(mode, HarmonicOscillator)
]) == 1)
self.assertTrue(
len([
mode for mode in conformer.modes
if isinstance(mode, HinderedRotor)
]) == 0)
