Dear CP2k community,<div><br></div><div>  I performed geometry optimization for C2H4 at first singlet excited state (S1) using ROKS, and found that the optimized geometry is exactly same as the ground state (S0). Can anyone help me to figure out the way to optimize geometry at low-spin excited state? I also ran an geometry optimization at first triplet excited state (T1) , and obtained a distinct geometry as expected. I guess I might miss some key options in the input file shown below.</div>
<div><br></div><div>Thanks.</div><div><br></div><div>Hanning Chen</div><div>Department of Chemistry</div><div>Northwestern University</div><div>Evanston, IL 60208</div><div>      </div><div><br></div><div><br></div><div>Input file:</div>
<div><br></div><div><br></div><div>  &FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    ROKS</div><div>    MULTIPLICITY 1</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-8</div><div>    &END QS</div>
<div>    &SCF</div><div>      MAX_SCF 100</div><div>      SCF_GUESS ATOMIC</div><div>      &DIAGONALIZATION</div><div>       ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>    &END SCF</div>
<div>    &XC</div><div>      &XC_FUNCTIONAL PADE</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 6.0 6.0 8.0</div>
<div>    &END CELL</div><div>    &COORD</div><div>    C         -0.2558206925       -0.0439868417        0.6850619693</div><div>    C         -0.2214915625       -0.0386665897       -0.6514432536</div><div>    H         -0.4100752923        0.8714166677        1.2581094627</div>
<div>    H         -0.1401850725       -0.9650616559        1.2579495323</div><div>    H         -0.0795226532       -0.9560364450       -1.2252531634</div><div>    H         -0.3528336094        0.8817708920       -1.2228724740</div>
<div>    &END COORD</div><div>        &KIND H</div><div>      BASIS_SET DZVP-GTH-PADE</div><div>      POTENTIAL GTH-PADE-q1</div><div>    &END KIND</div><div>    &KIND C</div><div>      BASIS_SET DZVP-GTH-PADE</div>
<div>      POTENTIAL GTH-PADE-q4</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&GLOBAL</div><div>  PROJECT C2H4</div><div>  RUN_TYPE GEOMETRY_OPTIMIZATION</div><div>  PRINT_LEVEL MEDIUM</div>
<div>&END GLOBAL</div><div><br></div>