def test_gauss_optts_calc():
calc = Calculation(name='test_ts_reopt_optts', molecule=test_mol,
method=method, keywords=optts_keywords,
bond_ids_to_add=[(0, 1)])
calc.run()
print(calc.input.added_internals)
assert os.path.exists('test_ts_reopt_optts_g09.com')
bond_added = False
for line in open('test_ts_reopt_optts_g09.com', 'r'):
if 'B' in line and len(line.split()) == 3:
bond_added = True
assert line.split()[0] == 'B'
assert line.split()[1] == '1'
assert line.split()[2] == '2'
assert bond_added
assert calc.get_normal_mode_displacements(mode_number=6) is not None
assert calc.terminated_normally()
assert calc.optimisation_converged()
assert calc.optimisation_nearly_converged() is False
assert len(calc.get_imaginary_freqs()) == 1
assert -40.324 < calc.get_free_energy() < -40.322
assert -40.301 < calc.get_enthalpy() < -40.299
def test_opt_calc():
calc = Calculation(name='opt',
molecule=test_mol,
method=method,
keywords=opt_keywords)
calc.run()
assert os.path.exists('opt_nwchem.nw')
assert os.path.exists('opt_nwchem.out')
final_atoms = calc.get_final_atoms()
assert len(final_atoms) == 5
assert type(final_atoms[0]) is Atom
assert -40.4165 < calc.get_energy() < -40.4164
assert calc.output.exists()
assert calc.output.file_lines is not None
assert calc.get_imaginary_freqs() == []
assert calc.input.filename == 'opt_nwchem.nw'
assert calc.output.filename == 'opt_nwchem.out'
assert calc.terminated_normally()
assert calc.optimisation_converged()
assert calc.optimisation_nearly_converged() is False
charges = calc.get_atomic_charges()
assert len(charges) == 5
assert all(-1.0 < c < 1.0 for c in charges)
# Optimisation should result in small gradients
gradients = calc.get_gradients()
assert len(gradients) == 5
assert all(-0.1 < np.linalg.norm(g) < 0.1 for g in gradients)
def test_psi4_opt_calculation():
methylchloride = Molecule(name='CH3Cl',
smiles='[H]C([H])(Cl)[H]',
solvent_name='water')
calc = Calculation(name='opt',
molecule=methylchloride,
method=method,
keywords=opt_keywords)
calc.run()
assert os.path.exists('opt_psi4.inp') is True
assert os.path.exists('opt_orca.out') is True
assert len(calc.get_final_atoms()) == 5
assert -499.735 < calc.get_energy() < -499.730
assert calc.output.exists()
assert calc.output.file_lines is not None
assert calc.get_imaginary_freqs() == []
assert calc.input.filename == 'opt_psi4.inp'
assert calc.output.filename == 'opt_psi4.out'
assert calc.terminated_normally()
assert calc.optimisation_converged()
calc = Calculation(name='opt',
molecule=methylchloride,
method=method,
keywords=opt_keywords)
# If the calculation is not run with calc.run() then there should be no
# input and the calc should raise that there is no input
with pytest.raises(NoInputError):
execute_calc(calc)
def test_mopac_opt_calculation():
os.chdir(os.path.join(here, 'data'))
calc = Calculation(name='opt',
molecule=methylchloride,
method=method,
keywords=Config.MOPAC.keywords.opt)
calc.run()
assert os.path.exists('opt_mopac.mop') is True
assert os.path.exists('opt_mopac.out') is True
assert len(calc.get_final_atoms()) == 5
# Actual energy in Hartrees
energy = Constants.eV2ha * -430.43191
assert energy - 0.0001 < calc.get_energy() < energy + 0.0001
assert calc.output.exists()
assert calc.output.file_lines is not None
assert calc.input.filename == 'opt_mopac.mop'
assert calc.output.filename == 'opt_mopac.out'
assert calc.terminated_normally()
assert calc.optimisation_converged() is True
with pytest.raises(NotImplementedError):
_ = calc.optimisation_nearly_converged()
with pytest.raises(NotImplementedError):
_ = calc.get_imaginary_freqs()
with pytest.raises(NotImplementedError):
_ = calc.get_normal_mode_displacements(4)
os.remove('opt_mopac.mop')
os.chdir(here)
def test_termination_short():
calc = Calculation(name='test', molecule=methylchloride,
method=method, keywords=Config.MOPAC.keywords.sp)
calc.output.filename = 'test.out'
calc.output.file_lines = ['JOB ENDED NORMALLY', 'another line']
assert calc.terminated_normally()
def test_get_version_no_output():
calc = Calculation(name='test',
molecule=methylchloride,
method=method,
keywords=method.keywords.sp)
calc.output.filename = 'test.out'
calc.output.file_lines = ['some', 'incorrect', 'lines']
assert not calc.terminated_normally()
version = method.get_version(calc)
assert version == '???'
def test_bad_geometry():
# Calculation with the wrong spin state should fail
calc = Calculation(name='h2_overlap_opt',
molecule=Molecule(atoms=[Atom('H'), Atom('H')]),
method=method,
keywords=Config.MOPAC.keywords.opt)
calc.output.filename = 'h2_overlap_opt_mopac.out'
calc.output.file_lines = open(calc.output.filename, 'r').readlines()
assert not calc.terminated_normally()
assert calc.get_energy() is None
assert not calc.optimisation_converged()
def test_fix_angle_error():
os.chdir(os.path.join(here, 'data', 'g09'))
mol = Molecule(smiles='CC/C=C/CO')
mol.name = 'molecule'
calc = Calculation(name='angle_fail', molecule=mol, method=method,
keywords=opt_keywords)
calc.run()
assert os.path.exists('angle_fail_g09_cartesian.com') is True
assert os.path.exists('angle_fail_g09_internal.com') is True
assert calc.output.filename == 'angle_fail_g09_internal.log'
assert calc.terminated_normally()
def test_bad_orca_output():
calc = Calculation(name='no_output', molecule=test_mol, method=method,
keywords=opt_keywords)
calc.output.file_lines = []
calc.output.rev_file_lines = []
assert calc.get_energy() is None
with pytest.raises(AtomsNotFound):
calc.get_final_atoms()
with pytest.raises(NoInputError):
calc.execute_calculation()
calc.output_file_lines = None
assert calc.terminated_normally() is False
def test_gauss_opt_calc():
os.chdir(os.path.join(here, 'data'))
methylchloride = Molecule(name='CH3Cl',
smiles='[H]C([H])(Cl)[H]',
solvent_name='water')
calc = Calculation(name='opt',
molecule=methylchloride,
method=method,
keywords=opt_keywords)
calc.run()
assert os.path.exists('opt_g09.com')
assert os.path.exists('opt_g09.log')
assert len(calc.get_final_atoms()) == 5
assert os.path.exists('opt_g09.xyz')
assert calc.get_energy() == -499.729222331
assert calc.output.exists()
assert calc.output.file_lines is not None
assert calc.get_imaginary_freqs() == []
with pytest.raises(NoNormalModesFound):
calc.get_normal_mode_displacements(mode_number=1)
assert calc.input.filename == 'opt_g09.com'
assert calc.output.filename == 'opt_g09.log'
assert calc.terminated_normally()
assert calc.optimisation_converged()
assert calc.optimisation_nearly_converged() is False
charges = calc.get_atomic_charges()
assert len(charges) == methylchloride.n_atoms
# Should be no very large atomic charges in this molecule
assert all(-1.0 < c < 1.0 for c in charges)
gradients = calc.get_gradients()
assert len(gradients) == methylchloride.n_atoms
assert len(gradients[0]) == 3
# Should be no large forces for an optimised molecule
assert sum(gradients[0]) < 0.1
os.remove('opt_g09.com')
os.chdir(here)
def test_orca_opt_calculation():
os.chdir(os.path.join(here, 'data'))
methylchloride = Molecule(name='CH3Cl',
smiles='[H]C([H])(Cl)[H]',
solvent_name='water')
calc = Calculation(name='opt', molecule=methylchloride, method=method,
keywords=opt_keywords)
calc.run()
assert os.path.exists('opt_orca.inp') is True
assert os.path.exists('opt_orca.out') is True
assert len(calc.get_final_atoms()) == 5
assert -499.735 < calc.get_energy() < -499.730
assert calc.output.exists()
assert calc.output.file_lines is not None
assert calc.get_imaginary_freqs() == []
assert calc.input.filename == 'opt_orca.inp'
assert calc.output.filename == 'opt_orca.out'
assert calc.terminated_normally()
assert calc.optimisation_converged()
assert calc.optimisation_nearly_converged() is False
with pytest.raises(NoNormalModesFound):
calc.get_normal_mode_displacements(mode_number=0)
# Should have a partial atomic charge for every atom
charges = calc.get_atomic_charges()
assert len(charges) == 5
assert type(charges[0]) == float
assert -1.0 < charges[0] < 1.0
calc = Calculation(name='opt', molecule=methylchloride, method=method,
keywords=opt_keywords)
# If the calculation is not run with calc.run() then there should be no
# input and the calc should raise that there is no input
with pytest.raises(NoInputError):
execute_calc(calc)
os.remove('opt_orca.inp')
os.chdir(here)
