The topological charge distribution P (Q) is calculated for lattice CP N−1 models. In order to suppress lattice cut-off effects, we employ a fixed point (FP) action. Through transformation of P (Q), we calculate the free energy F (θ) as a function of the θ parameter. For N=4, scaling behavior is observed for P (Q) and F (θ), as well as the correlation lengths ξ(Q). For N=2, however, scaling behavior is not observed, as expected. For comparison, we also make a calculation for the CP 3 model with a standard action. We furthermore pay special attention to the behavior of P (Q) in order to investigate the dynamics of instantons. For this purpose, we carefully consider the behavior of γ eff , which is an effective power of P (Q) (∼ exp(−CQ γ eff )), and reflects the local behavior of P (Q) as a function of Q. We study γ eff for two cases, the dilute gas approximation based on the Poisson distribution of instantons and the Debye-Hückel approximation of instanton quarks. In both cases, we find behavior similar to that observed in numerical simulations. §1. IntroductionIt is interesting to study the phase structure of asymptotic free theories such as QCD and the CP N −1 model. Non-perturbative studies of the phase structure of such theories are necessary in order to understand why effects of the topological term (θ term) are suppressed in Nature. The θ term affects the dynamics at low energy and is expected to lead to rich phase structures. 1) Actually, in the Z(N) gauge model, it has been shown by use of free energy arguments that oblique confinement phases could emerge and that an interesting phase structure may be realized. 2) In this paper we are concerned with the dynamics of the θ vacuum of CP N −1 models with a topological term, which have several dynamical properties in common with QCD. We believe that study of the two-dimensional model will be useful in acquiring information about realistic physics.From the numerical point of view, the topological term introduces a complex Boltzmann weight in the Euclidean lattice path integral formalism. The complex nature of the weight prevents one from straightforwardly applying the standard algorithm used for Monte Carlo simulations. This problem can be circumvented * )