Energy decomposition analysis (EDA) is a very important class of methods to investigate interaction between different fragments in a chemical system. By decomposing the total interaction energy into different easily understandable physical components, chemists can readily recognize which factors play a major and minor role in the interaction of interest, this is evidently very important in understanding the nature of the interaction. EDA methods have been widely employed in theoretical chemical studies. There are many EDA methods in the field of quantum chemistry. The present article will not make a comprehensive overview, interested readers are recommended to read relevant reviews. Briefly speaking, interfragment EDA can be mainly divided into two categories. The methods in the first category are dependent of variationally solved orbitals of every fragment and entire system, such as Kitaura-Morokuma, LMOEDA, GKS-EDA , NEDA and ALMO-EDA. The second category includes the methods that only depend on fragment wavefunctions, the most representative one is the symmetry-adapted perturbation theory (SAPT) method, which rigorously derive various interaction components between molecules through intermolecular perturbation theory. There are also many other forms of EDA. For example, Mayer's energy decomposition calculates various interaction terms between atoms based on electron integrals and density matrix of the entire system. The EDA-FF method proposed by us provides interaction information between atoms or fragments based on forcefields and atomic charges. ETS-NOCV method can further decompose orbital interaction into components corresponding to different features and make them visible. The scheme proposed by Liu decomposes total energy into steric, electrostatic, and fermionic quantum effects to explore the physical essence that causes the difference in total energy between different status. Although, as mentioned above, there are many interfragment EDA methods, each of them has its own limitations. For example, Kitaura-Morokuma is defined based on the now obsolete Hartree-Fock (HF) method. The use of LMOEDA and GKS-EDA mainly depends on the modified version of the GAMESS-US program, whose use is relatively complicated, and the source code needs to be compiled by users, which is somewhat difficult for junior researchers in quantum chemistry. ALMO-EDA is only implemented in the commercial Q-Chem program, thus limiting its accessibility. The well-known SAPT is mostly suitable for the study of weak interactions between two molecules but can hardly study the nature of chemical bonds. In addition, the accurate high-order SAPT realizations such as SAPT2+(3) are fairly expensive and implementation for open-shell systems are not available in most codes.