Finally, we note that the plasmon resonance effect of MNPs has not been taken into account in this work. We expect to include the resonance effect by adopting the frequencydependent parameters for metal atoms in future publications, which requires solving complex-valued response equations. The occupied MO coefficients could be expanded to the second-order in orbital rotations Θ vo as Ref. 90 suggested.