Density matrix of anharmonic oscillators by a thermodynamic variation-perturbation technique, a closed form approach

A general approximation technique for the density matrix of coupled anharmo nic oscillators is developed. Starting with the Moyal-Weyl trace formalism for the Wigner phase space distribution function, a combination of the Gibb s-Bogoliubov inequality for the partition function and an operator perturba tion technique is applied to the Moyal-Weyl trace. Introduction of generati ng functions and of operator techniques allows the closed form evaluation o f the Moyal-Weyl formula for anharmonic oscillators. The Fourier transforms relating the Moyal-Weyl formula and the Wigner function as well as those r elating the Wigner function and the density matrix can be given in terms of Hermite polynomials. The Moyal-Weyl formula as the crucial part of the tec hnique is worked out in the first order thermodynamic variation perturbatio n approximation for the anharmonic quartic oscillator as well as for two co upled anharmonic oscillators as an example for arbitrary coupled oscillator s. (C) 2001 Elsevier Science B.V. All rights reserved.