Using the particle-mesh method with 1283 grids and 643 particles, we have carried out a number of N-body simulations of the large scale structure for the cosmological model proposed by Fukuyama et al.(1996): this model contains the matter(Ω0) as well as a scalar field(Φ) with a finite mass that couples non-minimally with the scalar curvature R through the form of 1/2ηΦ2R, where η(= −80) is the coupling constant.

For simplicity, we have adopted the same values as those employed by Fukuyama et al.(1996) for all the parameters other than Ω0, fo which we have varied from 0.001 to 0.15. In performing our simulation, we further assume that only the density of matter ρ spatially fluctuates. The initial condition is created by perturbing the homogeneous distribution of the particles by means of the random Gaussian Harrison-Zeldovich spectrum.

The two-point correlation function ξ(r) is then compputed for each of the resulting structures to compare with the observational data(Davis and Peebles, 1983). It is interesting to note that the structure formation appears to be achieved rather straightforwardly with the scalar field model. However, in order to yield an agreement between the theoretical and the observational two-point correlation functions, we seem to require the Ω0 value much larger than 0.01 adopted by Fukuyama et al.(1996).