Our research interest is the theoretical study of strongly correlated condensed matter systems such as realized in high-temperature superconductors, frustrated quantum spins, nano-tubes and graphene, and ultracold atoms trapped in optical lattices. In particular, in recent years we have been developing a novel many-body method, the entanglement perturbation theory (EPT), for calculating the ground states as well as thermodynamic properties of some fundamental models for these strongly correlated condensed matter systems, such as the Heisenberg model, the Hubbard model and the Kondo-lattice model. The success so far is the calculation of spontaneous magnetization and the specific heat of the 3D Ising model , the ground state of the 1D Hubbard model , and a general description of the ground states in 2 dimensions . Our current efforts include (1) the calculation of the elementary excitations in solids such as spin-wave spectrum in Heisenberg antiferromagnets and fermionic spectrum in Hubbard models, (2) the Kondo impurities and quantum transport through them, (3) thermodynamics as well as dynamics (i.e. density of states) and (4) developments of robust numerical algorithms therein. To put EPT in a word, it is simple, general and exact, and essentially different from the currently available many-body methods such as exact diagonalization, Monte Carlo, numerical renormalization group and various mean field theories. In addition to the EPT, we also pursue a novel path-integral method, the cluster transfer matrix (CTM) [4, 5], for thermodynamics and dynamics (i.e., transport) of Kondo-like systems such as quantum dots, single electron transistors and composites thereof.