Density operator cluster expansion

Obviously we may expect that the simple two- and three-particle collision approximation discussed in the previous sections is not appropriate, because a large number of particles always interact simultaneously. Formally this approximation leads to divergencies. In the previous sections we used in a systematic way cluster expansions for the two- and three-particle density operator in order to include two-particle bound states and their relevant interaction in three- and four-particle clusters. In the framework of that consideration we started with the elementary particles (e, p) and their interactions. The bound states turned out to be special states, and, especially, scattering states were dealt with in a consistent manner. 

Here we assume that it is possible to expand these reduced distribution functions as a power series in the density at time t = 0. If this is not possible, one can still derive equations similar to (197) and (198) by making a cluster expansion of the N-particle streaming operator exp(-i, ). In the modified cluster expansions only the functions Fi(x,-, 0), F2(x,-, Xj, 0),..., appear. tAs a result dt this property, it follows that if at t = 0 all the s particles are more than a few molecular diameters from each other then A. = 1. 

Solution of the Kohn-Sham equations as outlined above are done within the static limit, i.e. use of the Born-Oppenheimer approximation, which implies that the motions of the nuclei and electrons are solved separately. It should however in many cases be of interest to include the dynamics of, for example, the reaction of molecules with clusters or surfaces. A combined ab initio method for solving both the geometric and electronic problem simultaneously is the Car-Parrinello method, which is a DFT dynamics method [52]. This method uses a plane wave expansion for the density, and the inner ions are replaced by pseudo-potentials [53]. Today this method has been extensively used for studies of dynamic problems in solids, clusters, fullerenes etc [54-61]. We have recently in a co-operation project with Andreoni at IBM used this technique for studying the existence of different isomers of transition metal clusters [62,63]. 

---