diff --git a/source/module_io/output_mulliken.h b/source/module_io/output_mulliken.h
index a5c0c21bbc..4611c124d1 100644
--- a/source/module_io/output_mulliken.h
+++ b/source/module_io/output_mulliken.h
@@ -138,7 +138,7 @@ void cal_mag(Parallel_Orbitals* pv,
             auto sc_lambda
                 = new hamilt::DeltaSpin<hamilt::OperatorLCAO<TK, double>>(nullptr,
                                                                           kv.kvec_d,
-                                                                          nullptr,
+                                                                          dynamic_cast<hamilt::HamiltLCAO<TK, double>*>(p_ham)->getHR(),
                                                                           ucell,
                                                                           &gd,
                                                                           two_center_bundle.overlap_orb_onsite.get(),
@@ -164,7 +164,7 @@ void cal_mag(Parallel_Orbitals* pv,
             auto sc_lambda = new hamilt::DeltaSpin<hamilt::OperatorLCAO<std::complex<double>, std::complex<double>>>(
                 nullptr,
                 kv.kvec_d,
-                nullptr,
+                dynamic_cast<hamilt::HamiltLCAO<std::complex<double>, std::complex<double>>*>(p_ham)->getHR(),
                 ucell,
                 &gd,
                 two_center_bundle.overlap_orb_onsite.get(),