Feature: Add supersine electric fields and refactor RT-TDDFT field handling

Author: AsTonyshment

docs/advanced/input_files/input-main.md

@@ -8,6 +8,13 @@
 #include "source_lcao/module_rt/td_moving_gauge.h"
 #include "source_lcao/module_rt/velocity_op.h"
 
+#include <memory>
+
+namespace elecstate
+{
+class TDFieldManager;
+}
+
 namespace ModuleESolver
 {
 
@@ -28,12 +35,7 @@ class ESolver_KS_LCAO_TDDFT : public ESolver_KS_LCAO<std::complex<double>, TR>
 
     void store_h_s_psi(UnitCell& ucell, const int istep, const int iter, const bool conv_esolver);
 
-    void iter_finish(UnitCell& ucell,
-                     const int istep,
-                     const int estep,
-                     const int estep_max,
-                     int& iter,
-                     bool& conv_esolver);
+    void iter_finish(UnitCell& ucell, const int istep, const int estep, const int estep_max, int& iter, bool& conv_esolver);
 
     virtual void after_scf(UnitCell& ucell, const int istep, const bool conv_esolver) override;
 
@@ -72,6 +74,8 @@ class ESolver_KS_LCAO_TDDFT : public ESolver_KS_LCAO<std::complex<double>, TR>
     //! Restart flag
     bool restart_done = false;
 
+    std::shared_ptr<elecstate::TDFieldManager> td_field_manager_;
+
   private:
     void weight_dm_rho(const UnitCell& ucell);
 };

docs/parameters.yaml

@@ -47,6 +47,7 @@ if(ENABLE_LCAO)
   list(APPEND objects
       elecstate_lcao.cpp
       module_pot/H_TDDFT_pw.cpp
+      module_pot/td_field.cpp
       module_dm/init_dm.cpp
       module_dm/density_matrix.cpp
       module_dm/density_matrix_io.cpp