The temperature dependence of electric dipolar susceptibility χP is discussed on the basis of the Anderson-Holstein model with the use of a numerical renormalization group (NRG) technique. Note that χP is related with phonon Green's function D. In order to obtain correct temperature dependence of χP at low temperatures, we propose a method to evaluate χP through the Dyson equation from charge susceptibility χc calculated by the NRG, in contrast to the direct NRG calculation of D. We find that the irreducible charge susceptibility estimated from χc agree with the perturbation calculation, suggesting that our method works well. In the research field of condensed matter physics, exotic magnetism in cage structure materials such as filled skutterudites has attracted much attention due to the interests on new electronic properties caused by oscillation of a guest atom in a cage composed of relatively light atoms [1]. Such oscillation with large amplitude is frequently called rattling and it is considered to play crucial roles for the formation of magnetically-robust heavy electron state in SmOs 4 Sb 12 [2]. This peculiar heavyelectron state has been theoretically investigated from various aspects by several groups [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17].Recently, the Kondo effect of a vibrating magnetic ion in a cage has been theoretically discussed on the basis of a two-channel conduction electron system hybridized with a vibrating magnetic ion [18][19][20]. Note here that a vibrating ion inevitably induces electric dipole moment. Then, it has been found that magnetic and non-magnetic Kondo effects alternatively occur due to the screening of spin moment and electric dipole moment of vibrating ion [21]. In particular, electric dipolar two-channel Kondo effect has been found to occur for weak Coulomb interaction. Then, it has been proposed that magnetically robust heavy-electron state appears near the fixed point of electric dipolar two-channel Kondo effect.In this paper, in order to promote our understanding on the Kondo effect concerning electric dipole moment P , we analyze the temperature dependence of electric dipolar susceptibility χ P on the basis of the AndersonHolstein Hamiltonian with the use of a numerical renormalization group method. For the reproduction of correct temperature dependence of χ P at low temperatures, we propose a method to evaluate χ P through the Dyson equation from charge susceptibility χ c . This method is found to provide correct results in the temperature region lower than the Kondo temperature, in sharp contrast to the numerical evaluation of the phonon Green's function which is directly related to χ P . Now we explain the model Hamiltonian. We consider a conduction electron system in which an impurity ion * Electronic address: fuse-takahiro@tmu.ac.jp is embedded. On the impurity site, localized electrons are coupled with ion vibration. The situation is well described by the Anderson-Holstein model, given by [22] (1) where ε k is the dispersion of conduction electron, c ...