You do not have permission to delete messages in this group
Copy link
Report message
Sign in to report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to cp2k
Dear developers,
I know that is not possible to perform SCF calculations using the OT scheme when extra molecular orbitals are requested.
The problem is that for many not trivial systems the OT scheme is always able to reach a SCF convergence, while the standard DIAGONALIZATION algorithms ( no matter which combinations of parameters are used for the mixing ) fail, even starting from a fully converged electron density reached with the OT scheme. Strict diagonalization algorithms are also required in various subroutines like KPOINT, XAS and GW for example.
Having the possibility to calculate added MOS using the OT scheme would be enormously beneficial in terms of applicability of the CP2K software and a considerable gain in computational resources.
So, I was wondering if there is a specific reason why extra MO cannot be computed with the OT scheme, or more simply it has not been implemented yet.
Thank you in advance for the help and clarifications.
Kind regards, Lorenzo Agosta
Jürg Hutter
unread,
May 26, 2023, 4:50:31 AM5/26/23
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Sign in to report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to cp...@googlegroups.com
Hi
OT is a direct energy optimization method. As the energy doesn't depend on the virtual states, they cannot
be calculated directly within the optimization.
However, you can calculate additional states (LUMO +++) after an OT optimization. Several property options
allow for such calculations (e.g. PRINT/MO_CUBES). Another possibility is to restart an optimized structure
and requesting a single diagonalization step.