QMMM Excited States
The excited state QM/MM energy calculations can be performed with tce
task qmmm tce energy
task qmmm tddft energy
levels of theory. The excited state QM/MM gradient energy calculations can be performed only numerically.
The example below illustrates single point excited state energy calculation at DFT/B3LYP level for quantum water molecule embedded into 20 angstrom box of classical SPCE/E water molecules. The preparation stage that involves definition of the QM region and solvation is performed as part of the calculation. Note that water fragment file wtr.frg is required in this calculation. Prepare run will generate restart file (wtr_ref.rst) and topology file (wtr.top)
In the QM/MM interface block the use of bq_zone value of 3.0 Angstrom is specified.
start wtr permanent_dir ./perm scratch_dir ./data prepare source wtr0.pdb new_top new_seq new_rst modify segment 1 quantum center orient solvate box 3.0 update lists ignore write wtr_ref.rst write wtr_ref.pdb end task prepare md system wtr_ref end basis * library "6-31G" end dft xc b3lyp end tddft nroots 4 notriplet end qmmm bq_zone 3.0 end task qmmm tddft energy