refactor: remove dead cross-module includes on reverse-dependency edges

Static scan + per-target compile verification identified 9 unused
includes that create reverse/lateral dependency edges between modules.
Removing them weakens the coupling without any behavior change:

  io   -> md    input_conv.cpp        (md_func.h)
  estate -> lcao  elecstate_energy_terms.cpp, elecstate_print.cpp
                                       (module_deepks/LCAO_deepks.h)
  lcao -> pw    rdmft_tools.cpp (structure_factor.h),
                wavefunc_in_pw.cpp (soc.h)
  lcao -> io    FORCE_gamma.cpp, FORCE_k.cpp (module_hs/write_HS.h)
  pw   -> io    forces_cc.cpp, forces_scc.cpp (module_output/output_log.h)

Verified by building io_basic, elecstate, rdmft, hamilt_lcao and
module_pwdft (ENABLE_LCAO=ON) after removal; all link targets compile.

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>

Author: Critsium-xy

source/source_estate/elecstate_energy_terms.cpp

@@ -3,7 +3,6 @@
 #include "source_estate/module_pot/H_Hartree_pw.h"
 #include "source_estate/module_pot/efield.h"
 #include "source_estate/module_pot/gatefield.h"
-#include "source_lcao/module_deepks/LCAO_deepks.h"
 #include "source_lcao/module_deltaspin/spin_constrain.h"
 #include "source_lcao/module_dftu/dftu.h" // mohan add 2025-11-06
 

source/source_estate/elecstate_print.cpp

@@ -6,7 +6,6 @@
 #include "source_estate/module_pot/efield.h"
 #include "source_estate/module_pot/gatefield.h"
 #include "source_hamilt/module_xc/xc_functional.h"
-#include "source_lcao/module_deepks/LCAO_deepks.h"
 #include "source_io/module_parameter/parameter.h"
 #include "occupy.h"
 namespace elecstate

source/source_io/module_parameter/input_conv.cpp

@@ -42,7 +42,6 @@
 #include "source_estate/module_pot/gatefield.h"
 #include "source_hsolver/hsolver_lcao.h"
 #include "source_hsolver/hsolver_pw.h"
-#include "source_md/md_func.h"
 #include "source_relax/bfgs_basic.h"
 #include "source_relax/ions_move_cg.h"
 

source/source_lcao/FORCE_gamma.cpp

@@ -13,7 +13,6 @@
 #include "source_estate/elecstate_lcao.h"
 #include "source_lcao/LCAO_domain.h"
 #include "source_lcao/pulay_fs.h"
-#include "source_io/module_hs/write_HS.h"
 
 template <>
 void Force_LCAO<double>::allocate(const UnitCell& ucell,

source/source_lcao/FORCE_k.cpp

@@ -10,7 +10,6 @@
 #include "source_estate/module_dm/cal_dm_psi.h"
 #include "source_lcao/LCAO_domain.h"
 #include "source_lcao/pulay_fs.h"
-#include "source_io/module_hs/write_HS.h"
 #include "source_io/module_parameter/parameter.h"
 
 #include <map>

source/source_lcao/module_rdmft/rdmft_tools.cpp

@@ -10,7 +10,6 @@
 #include "source_estate/module_pot/H_Hartree_pw.h"
 #include "source_estate/module_pot/pot_local.h"
 #include "source_estate/module_pot/pot_xc.h"
-#include "source_pw/module_pwdft/structure_factor.h"
 #include "source_lcao/module_gint/gint_interface.h"
 #include "source_io/module_parameter/parameter.h"
 

source/source_lcao/wavefunc_in_pw.cpp

@@ -5,7 +5,6 @@
 #include "source_base/math_sphbes.h"
 #include "source_base/math_polyint.h"
 #include "source_base/math_ylmreal.h"
-#include "source_pw/module_pwdft/soc.h"
 
 void Wavefunc_in_pw::make_table_q(
 	const UnitCell &ucell,

source/source_pw/module_pwdft/forces_cc.cpp

@@ -2,7 +2,6 @@
 #include "stress_func.h"
 #include "source_base/parallel_reduce.h"
 #include "source_io/module_parameter/parameter.h"
-#include "source_io/module_output/output_log.h"
 // new
 #include "source_base/complexmatrix.h"
 #include "source_base/libm/libm.h"

source/source_pw/module_pwdft/forces_scc.cpp

@@ -1,6 +1,5 @@
 #include "forces.h"
 #include "source_base/parallel_reduce.h"
-#include "source_io/module_output/output_log.h"
 #include "stress_func.h"
 // new
 #include "source_base/complexmatrix.h"