def test_pcm_init(self):
test_molecule = QCInput.from_file(os.path.join(test_dir, "new_qchem_files/pcm.qin")).molecule
test_forceset = ForceSet(molecule=test_molecule, pcm_dielectric=10.0)
self.assertEqual(
test_forceset.rem,
{
"job_type": "force",
"gen_scfman": "true",
"basis": "def2-tzvppd",
"max_scf_cycles": "200",
"method": "wb97xd",
"scf_algorithm": "diis",
"xc_grid": "3",
"solvent_method": "pcm",
"symmetry": "false",
"sym_ignore": "true",
"resp_charges": "true",
},
)
self.assertEqual(
test_forceset.pcm,
{"heavypoints": "194", "hpoints": "194", "radii": "uff", "theory": "cpcm", "vdwscale": "1.1"},
)
self.assertEqual(test_forceset.solvent, {"dielectric": 10.0})
self.assertEqual(test_forceset.molecule, test_molecule)
def test_smd_init(self):
test_molecule = QCInput.from_file(
os.path.join(test_dir, "new_qchem_files/pcm.qin")).molecule
test_forceset = ForceSet(molecule=test_molecule, smd_solvent="water")
self.assertEqual(
test_forceset.rem,
{
"job_type": "force",
"gen_scfman": "true",
"basis": "def2-tzvppd",
"max_scf_cycles": "100",
"method": "wb97xd",
"scf_algorithm": "diis",
"xc_grid": "3",
"thresh": "14",
"s2thresh": "16",
"solvent_method": "smd",
"ideriv": "1",
"symmetry": "false",
"sym_ignore": "true",
"resp_charges": "true",
},
)
self.assertEqual(test_forceset.smx, {"solvent": "water"})
self.assertEqual(test_forceset.molecule, test_molecule)
def test_init(self):
test_molecule = QCInput.from_file(os.path.join(test_dir, "new_qchem_files/pcm.qin")).molecule
test_forceset = ForceSet(molecule=test_molecule)
self.assertEqual(
test_forceset.rem,
{
"job_type": "force",
"gen_scfman": "true",
"basis": "def2-tzvppd",
"max_scf_cycles": "200",
"method": "wb97xd",
"scf_algorithm": "diis",
"xc_grid": "3",
"symmetry": "false",
"sym_ignore": "true",
"resp_charges": "true",
},
)
self.assertEqual(test_forceset.pcm, {})
self.assertEqual(test_forceset.solvent, {})
self.assertEqual(test_forceset.molecule, test_molecule)