Dear cp2k users,<div><br></div><div>I am intending to do a tddfpt calculation to get the 10 possible first excitations for negatively charged-Nitrogen-vacancy center in diamond. The system has negative charge and as a result has got triplet excited states as well. Really appreciate it if someone could help me to make the input file as it should be.  currently I get the error " XC_DERIV method not implemented for GPW-LSD!"</div><div><br></div><div>Input file</div><div>===================================</div><div>&GLOBAL</div><div>  PROJECT nvne</div><div>  RUN_TYPE ENERGY</div><div>  PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div><br></div><div>   &PROPERTIES</div><div>    &TDDFPT</div><div>!       RES_ETYPE   triplets</div><div>       NSTATES     10            # number of excited states</div><div>       MAX_ITER    200           # maximum number of Davidson iterations</div><div>       CONVERGENCE [eV] 1.0e-5   # convergence on maximum energy change between iterations</div><div>       &MGRID</div><div>          CUTOFF 300 # separate cutoff for TDDFPT calc</div><div>       &END</div><div><br></div><div>       RESTART     .TRUE.</div><div>       WFN_RESTART_FILE_NAME RESTART.tdwfn</div><div><br></div><div>    &END TDDFPT</div><div>  &END PROPERTIES</div><div><br></div><div>  !&PRINT</div><div>  !  &FORCES</div><div>  !  &END</div><div>  !  &TOTAL_NUMBERS  ON</div><div>  !  &END TOTAL_NUMBERS</div><div>  !&END PRINT</div><div><br></div><div>  &DFT</div><div><br></div><div>    UKS .TRUE.</div><div>    MULTIPLICITY 3</div><div>    CHARGE -1</div><div>    &PRINT</div><div>      &PDOS SILENT</div><div>        COMPONENTS F</div><div>        NLUMO  -1</div><div>      &END PDOS</div><div>      &MO_CUBES</div><div>        NLUMO -1</div><div>        NHOMO 50</div><div>        WRITE_CUBE .FALSE.</div><div>        STRIDE 1 1 1</div><div>      &END</div><div>    &END PRINT</div><div><br></div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME POTENTIAL</div><div>    WFN_RESTART_FILE_NAME slg-RESTART.wfn</div><div><br></div><div>    &MGRID</div><div>      CUTOFF 300</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-12</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS RESTART</div><div>      EPS_SCF 1.0E-5</div><div>      MAX_SCF 300</div><div><br></div><div>      !&OT</div><div>      !  PRECONDITIONER FULL_ALL</div><div>      !  MINIMIZER DIIS</div><div>      !&END OT</div><div><br></div><div>      ADDED_MOS 30</div><div>      !CHOLESKY INVERSE</div><div>      !&SMEAR ON</div><div>      !  METHOD FERMI_DIRAC</div><div>      !  ELECTRONIC_TEMPERATURE [K] 300</div><div>      !&END SMEAR</div><div>      &DIAGONALIZATION</div><div>        ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &MIXING</div><div>        METHOD BROYDEN_MIXING</div><div>        ALPHA 0.5</div><div>        BETA  0.5</div><div>        NBROYDEN 8</div><div>      &END MIXING</div><div><br></div><div>    &END SCF</div><div><br></div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>      &VDW_POTENTIAL</div><div>         POTENTIAL_TYPE PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>            PARAMETER_FILE_NAME dftd3.dat</div><div>            TYPE DFTD3</div><div>            REFERENCE_FUNCTIONAL PBE</div><div>            R_CUTOFF [angstrom] 16</div><div>         &END</div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div>!    &XC_DERIV COLLOCATE</div><div>!    &END XC_DERIV</div><div><br></div><div>  &END DFT</div><div>  &SUBSYS</div><div><br></div><div>    &KIND C</div><div>      ELEMENT C</div><div>      BASIS_SET ORB DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND</div><div>    &KIND N</div><div>      ELEMENT N</div><div>      BASIS_SET ORB DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q5</div><div>    &END KIND</div><div><br></div><div>    &CELL</div><div>      PERIODIC XYZ</div><div>      A  7.015032481   0.007230903  -0.007210904</div><div>      B  0.007230725   7.015042835  -0.007225784</div><div>      C -0.008733070  -0.008757378  10.521451073</div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME coord.xyz</div><div>      COORD_FILE_FORMAT XYZ</div><div>    &END TOPOLOGY</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div><br></div><div>============================== </div>