Hidemi NAGAO, Masayoahi NAKANO* (*Osaka Univ.), Satoru YAMADA*, Koji OHTA (Osaka Natl. Res. Inst.) and Kizashi YAMAGUCHI*
[Synthetic Metals 85, 1159 (1997)]
Based on the path integral method formulated in terms of the Thouless parametrization, numerical calculations of the ensemble average of the polarizability and the second polarizability are carried out by means of the Monte Carlo method. An H3+ system is examined to elucidate the performance of the present method. We discuss the dependence of the calculated (hyper)polarizabilities on the number of the division of the partition function and the accuracy of the results obtained by the present method.