modified on 11 September 2014 at 15:03 ••• 1,207 views

Release66:QMMM Excited States

From NWChem

Revision as of 15:03, 11 September 2014 by Silvia (Talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Template:Release66:QMMM Single Point Calculations The excited state QM/MM energy calculations can be performed with TCE

task qmmm tce energy

or TDDFT

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