Abstract. The effect of the electron-electron (e-e) interaction on the Peierls lattice distortions due to the electron-lattice (e-l) interaction is studied in the two-dimensional Peierls-Hubbard model, treating the fluctuation of e-e interaction around the Hartree-Fock solution within the 2nd order perturbation theory. In our previous work, using the Hartree-Fock approximation, we found multimode Peierls lattice distortions with wave vectors, the nesting vector Q and those parallel to it, are not affected by an e-e interaction if it is weak compared with the e-l coupling. The phase transition between the BOW (bond order wave) with multimode lattice distortions and the SDW (spin density wave) with the wave vector Q behaves as the 1st order transition. The property of the multimode BOW is found to change drastically when we consider the fluctuation effect within the 2nd order perturbation theory. The Fourier components of the multimode BOW increase gradually as the e-e interaction parameter is increased. Especially the Fourier component with the wave vector equal to the smallest reciprocal lattice vector in the presence of the multimode BOW is most strongly affected by the fluctuation effect.Keywords: multimode Peierls phase, 2D Peierls-Hubbard model, electron correlation, BOW-SDW transition, perturbation theory PACS: 71.30.+h, 71.45.Lr, 75.30.Fv, 63.20.Kr Competition between the Peierls lattice distortion and the electron-electron (e-e) interaction is widely discussed in the 1D Peierls-Hubbard model with a half-filled electronic band. In this system the lattice distortion which has a wave number Q π plays a very important role in the all range of e-e interaction parameter [1].We are interested in the same problem in the 2D system, since the ground state of the 2D electron-lattice (e-l) system, especially in the 2D SSH (Su-SchriefferHeeger) model with a half-filled electronic band, is an unusual Peierls state [2] in contrast to the Peierls state expected from that of the 1D system [3]. The ground state of 2D SSH model has multimode lattice distortions with the wave vector Q π πµ and those parallel to it. In addition, there are infinite number of degenerate states with non-equivalent distortion patterns corresponding to many different combinations of the Fourier components of the distortion for the wave vectors parallel to Q [4].The effects of the e-e interaction on this multimode Peierls state can be discussed using the 2D PeierlsHubbard (2D PH) model. The Hamiltonian of the 2D PH model is described as follows and e x´yµ the lattice displacements at the site r and the unit lattice vector for respective direction, respectively, and the parameters t 0 , α, K and U are the electron transfer integral between nearest-neighbor sites, the e-l coupling constant, the force constant describing the ionic coupling strength and the strength of the on-site e-e interaction between up and down spins, respectively.In our previous paper [6] we used the Hartree-Fock approximation (HFA) as the easiest approach to discuss the e-e i...