deepmodeling / abacus-develop

An electronic structure package based on either plane wave basis or numerical atomic orbitals.
http://abacus.ustc.edu.cn
GNU Lesser General Public License v3.0
151 stars 122 forks source link

Request: Transition State Optimization in ABACUS itself #4327

Open QuantumMisaka opened 1 month ago

QuantumMisaka commented 1 month ago

Background

There is no transition state optimization (or, saddle-point optimization) implementation in ABACUS itself. Although we can use ASE-ABACUS interfaces to use the TS method in ASE or coded by ASE, but:

  1. There are some problem in ASE-ABACUS, with all single-end TS method (e.g. Dimer) used by ASE-ABACUS interfaces cannot give right result and effective process.
  2. It is needed for a DFT calculation program to have (at least on) saddle-point search implementation.

Describe the solution you'd like

ABACUS should have (at least one) TS optimization method. which can be:

  1. NEB
  2. Dimer
  3. TS based on P-RFO (Partitional Rational Function Optimization) method, namely, Sella

The first and second one can be implemented easily, but the last one is the most recommended. Since in PBC, internal coordinates cannot be effectively used, we can only use the PBC part of Sella. ref: https://pubs.acs.org/doi/full/10.1021/acs.jctc.9b00869

Task list only for developers

Notice Possible Changes of Behavior (Reminder only for developers)

No response

Notice any changes of core modules (Reminder only for developers)

No response

Notice Possible Changes of Core Modules (Reminder only for developers)

No response

Additional Context

No response

Task list for Issue attackers (only for developers)

kirk0830 commented 1 month ago

I doubt the robustness of DIMER, even in other softwares I know people will first use NEB then DIMER. If possible, I think who is willing to implement TS searching algo. should only consider 1 and/or 3.